Interpolating to globally refined distributed mesh

[Konrad Simon]

Dear deal.ii community,

I am trying to interpolate a distributed solution vector onto a globally refined distributed mesh. Actually I need to globally refine more than once but I believe I can iterate this procedure (is there a more efficient way? The finite element can be vector or scalar valued, and the solution vector a block vector)

However, despite looking at tutorials and the documentation I get an error (segmentation fault) without further explanation. I did not manage to track down the error - probably I am using the interface not correctly. Anyone has experience with this?

A minimal example is attached. (I used deal.ii 9.2.0 but it should work with 9.1.1, too.)

Best regards,
Konrad

[Konrad Simon]

Add-on: The problem is in the function RefineInterpolate<dim>::refine_and_interpolate_on_distributed_mesh()

[Konrad Simon]

Dear deal.ii community,

I solved it, sorry for bugging you with it but simple mistakes can bug you for long... 

I simply forgot to re-distribute the dofs for the finite element after refining the mesh. Ooofff :-/

Best,
Konrad