#include "petscdmplex.h" PetscErrorCode DMPlexGetClosureIndices(DM dm, PetscSection section, PetscSection idxSection, PetscInt point, PetscInt *numIndices, PetscInt **indices, PetscInt *outOffsets)Not collective
dm | - The DM | |
section | - The section describing the points (a local section) | |
idxSection | - The section on which to obtain indices (may be local or global) | |
point | - The mesh point |
numIndices | - The number of indices | |
indices | - The indices | |
outOffsets | - Field offset if not NULL |
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.