erhart1

One could consider the extension of the tensorial models beyond rank 2. This could be used, e.g., to predict the piezoelectric coupling.

It would be useful to include a header in the `thermo.out` file that specifies the meaning of the columns. If these lines are prepended by a widely used comment symbols...