DMPlexGetClosureIndices#

Gets the global dof indices associated with the closure of the given point within the provided sections.

Synopsis#

#include "petscdmplex.h"   
PetscErrorCode DMPlexGetClosureIndices(DM dm, PetscSection section, PetscSection idxSection, PetscInt point, PetscBool useClPerm, PetscInt *numIndices, PetscInt *indices[], PetscInt outOffsets[], PetscScalar *values[])

Not collective

Input Parameters#

dm - The DM
section - The PetscSection describing the points (a local section)
idxSection - The PetscSection from which to obtain indices (may be local or global)
point - The point defining the closure
useClPerm - Use the closure point permutation if available

Output Parameters#

numIndices - The number of dof indices in the closure of point with the input sections
indices - The dof indices
outOffsets - Array to write the field offsets into, or NULL
values - The input values, which may be modified if sign flips are induced by the point symmetries, or NULL

Notes#

Must call DMPlexRestoreClosureIndices() to free allocated memory

If idxSection is global, any constrained dofs (see DMAddBoundary(), for example) will get negative indices. The value of those indices is not significant. If idxSection is local, the constrained dofs will yield the involution -(idx+1) of their index in a local vector. A caller who does not wish to distinguish those points may recover the nonnegative indices via involution, -(-(idx+1)+1)==idx. Local indices are provided when idxSection == section, otherwise global indices (with the above semantics) are implied.