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