petsc-3.14.6 2021-03-30

DMPlexGetTransitiveClosure

Return the points on the transitive closure of the in-edges or out-edges for this point in the DAG

#include "petscdmplex.h"   
PetscErrorCode DMPlexGetTransitiveClosure(DM dm, PetscInt p, PetscBool useCone, PetscInt *numPoints, PetscInt *points[])

Not collective

	mesh	- The DMPlex
	p	- The point, which must lie in the chart set with DMPlexSetChart()
	useCone	- PETSC_TRUE for in-edges, otherwise use out-edges
	points	- If points is NULL on input, internal storage will be returned, otherwise the provided array is used

	numPoints	- The number of points in the closure, so points[] is of size 2*numPoints
	points	- The points and point orientations, interleaved as pairs [p0, o0, p1, o1, ...]

If using internal storage (points is NULL on input), each call overwrites the last output.

Since it returns an array, this routine is only available in Fortran 90, and you must include petsc.h90 in your code.

The numPoints argument is not present in the Fortran 90 binding since it is internal to the array.

beginner