Overlapping mesh - integral calculation

You cannot use direct interpolation: this is a sequential kernel, you’d need to have the same subdomain meshes, which is no the case given that you partition both global meshes without any constraints.

That being said, you are using the low-level API, which I explicitly recommended not to do in a previous post.

Instead, you can use higher-level macros, more user friendly, like buildDmesh() or createMat() . It will allow you to do much more powerful operations like interpolation in parallel […].

You problem is in fact addressed in this other post. You can also have a look at this other example. Let me know if something is not clear to you.
What you do could work but you are not using the partition of unity so everything is over-summed on the overlap. Plus, you have to go from local to global to local again, which is rather inefficient.

1 Like