On Tue, Sep 22, 2015 at 2:42 PM, Harshad Sahasrabudhe <hsaha...@purdue.edu> wrote:
> Hi, > > I'm using ParallelMesh to discretize the Poisson equation. I want to > integrate a certain quantity over each element. That quantity is not > available for all quadrature points of some elements on one process (such > elements are ghosted), so I sum the parts which are available in respective > processes for that particular element and MPI reduce the sum to get the > complete integral. > > Now to be able to do this, I have to identify the elements for which the > integral has to be MPI reduced. Can anyone help me with how the element can > be identified? I tried using elem->id() but I get an exception in > DofObject::id(), since the _id == invalid_id. > Not sure what this means: "quantity is not available for all quadrature points of some elements on one process" You shouldn't be computing anything on ghosted elements... they will be the responsibility of the processor that owns them. -- John ------------------------------------------------------------------------------ _______________________________________________ Libmesh-users mailing list Libmesh-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-users
