petsc-3.7.3 2016-08-01

DMPlexGetTransitiveClosure_Internal

Return the points on the transitive closure of the in-edges or out-edges for this point in the Sieve DAG with a specified initial orientation

Synopsis

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

Not collective

Input Parameters

	mesh	- The DMPlex
	p	- The Sieve point, which must lie in the chart set with DMPlexSetChart()
	orientation	- The orientation of the point
	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

Output Parameters

	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, ...]

Note

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

Fortran Notes

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