Generalize the input for the stress-free transformation strains
The Fortran code assumes that the (unrotated) SFTS tensor is diagonal. However in general it isn't (and can't be diagonalized), as in the beta1 precipitates in Mg-Nd. The input for the SFTS should be a 3x3 matrix.