I don't understand this. An element (and therefore all quadrature
points on that element) is either entirely owned by a processor or is
entirely not. Do you have a coarser grid that you're integrating on?
> Each element is entirely on one process, but the quantity that I'm
integrating isn't. This is why I ghosted the element and I'm
performing part of the integration one multiple processes.
How are you finding the elements which contain each of your quadrature
points? I.e. when you call elem->id(), where did "elem" come from?
I initially use the Poisson mesh point locator on quadrature points of the
> second mesh to find the elements of Poisson mesh and ghost them.
> Then I run an integral over the quadrature points of these Poisson
> elements, while using interpolated quantity from second mesh.
I forgot to mention that I call allgather() on the Poisson mesh and
reinitialize the point locator before calling it.
------------------------------------------------------------------------------
_______________________________________________
Libmesh-users mailing list
Libmesh-users@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/libmesh-users
