Read in deformation function on each process  and adapt riemannian trustregion solver - parallel assembly works again!

See merge request !44