> > 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. So the reason why that quantity isn't available on the process which owns the element is because the quantity is interpolated from another mesh. This other mesh is smaller spatially than the Poisson mesh and thus has a different distribution. On Tue, Sep 22, 2015 at 4:55 PM, John Peterson <jwpeter...@gmail.com> wrote: > > > 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
