> > 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. Yes, this is the reason why I allgather() first, because there will definitely be points which I want to locate outside the scope of the process. But I plan to delete_remote_elements() after I have identified the required elements. The mesh doesn't change afterwards, which is why I think this is feasible. > 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? Yes, I will try to get to the bottom of this issue. If I find it's a LibMesh bug, then I'll let you know. On Tue, Sep 22, 2015 at 5:34 PM, Roy Stogner <royst...@ices.utexas.edu> wrote: > > 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
