On Tue, 22 Sep 2015, Harshad Sahasrabudhe wrote: > I forgot to mention that I call allgather() on the Poisson mesh and > reinitialize the point locator before calling it.
At that point there's little reason to bother using a ParallelMesh at all; after an allgather it's become a less-efficient-than-SerialMesh serial mesh. You can use a point locator on the distributed mesh without serializing first; you just have to account for the fact that it will return NULL when you don't have any semilocal elements containing the point you requested. This is sounding like you might be triggering a libMesh bug, though. An allgather() should leave every newly gathered element with the correct id. Any chance you can boil the problem down to a test case we can replicate? --- Roy ------------------------------------------------------------------------------ _______________________________________________ Libmesh-users mailing list Libmesh-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-users
