[Libmesh-users] Working with 2 meshes

[Harshad S. <hsa...@pu...>]

Hi,

I need help with transferring data between 2 meshes. I have one mesh for an
equation system and another mesh (which is on a sub-domain of the first
mesh) for a second equation system. Both the meshes are parallelized over
all the available processors. Now, I need to get the solution from second
mesh to first mesh on the correct processors. What is the best way to do
this?

MeshFunction can calculate the solution of the second system on the mesh
points of the first system, but it can only do this for local elements.

If anyone wants to know the details of what I'm trying to do, read further.
I'm working on a Schrodinger-Poisson self-consistent simulation. Basically,
I solve Poisson equation on a large (and coarse) mesh. Then I feed the
potential (solution of Poisson system) to the Schrodinger system which
resides on a smaller mesh contained within the region occupied by the
Poisson system. Then I solve the Schrodinger system and feed the charge
density (solution of Schrodinger system) back to the Poisson system. This
loop runs till the potential and charge converge.

This runs fine serially, but when I try to parallelize it the data required
on a certain processor resides on some other processor.

Any help would be appreciated!

Thanks!
Harshad