Description

Improving the cell dynamics geometry optimization

The relaxation is now done with the strain and the stress directly. Because of the symmetry, only 6 components are optimized and not 9. Moreover, this makes the code more robust and efficient.

The internal variable is the strain tensor, which has the advantage to start from identity and to be properly normalized, thus avoid small numbers and instabilities in the optimization procedure.