[Libmesh-users] L2 Projection between meshes with split dofs

[Renato P. <re...@gm...>]

Hi all,

I have a doubt on L2 projection between meshes with split dofs.

I am currently using the following code to project values from one mesh to
another more refined one (I perform refinement myself):

      MeshFunction * mesh_func = new MeshFunction( . . . );
      mesh_func->init();
      mesh_func->enable_out_of_mesh_mode( 0 );
      sys.project_solution( ... )

I can see it interpolates the projected values based on the coordinate of
the dofs. Is that so?

It happens that in some coordinates I have duplicated dofs representing
fractures. Displacements are different on both sides of the fracture even
though their coordinates are the same.
Can you see any smart-hack to avoid the interpolation _only_ on fracture
nodes?
Or should I rewrite the whole 'project_solution' with my own stuff?

Thanks,
Renato