Update deformation gradient logic

[CompRhys]

Initially we focussed on matching the Frechet filter in ASE. That implementation fixes the reference cell as the initial cell. This can lead to blowup of optimization if the cells change too much.

Look at alternatives such as trust regions