[chombier]

nice work! you could probably speed the computation quite a lot by using a small deformation basis (e.g. using modal analysis) and optimize in this smaller space instead.