[Libmesh-users] Data projection on mesh refinement

[<an...@tn...>]

Hi all

Either I am again missing something obvious (in that case: sorry for
the noise), or I feel something is missing:

Upon mesh refinement, the mesh.data is not projected on the new mesh.
To my understanding this could be accomplished in a similar way as the
projection of solution vectors. Indeed the only workaround I could
come up with was to add a trivial system, solve it for the nodal data
on the RHS, and subsequently use the solution of this system to
obtain the data values on a refined mesh.

I think that such a routine would be either explicitly called as
mesh.data.reinit(), or would automatically be called by       
mesh_refinement.refine_and_coarsen_elements();

Does this make sense?

Best, Martin