#include "petscdmplextransform.h" PetscErrorCode DMPlexTransformCellTransform(DMPlexTransform tr, DMPolytopeType source, PetscInt p, PetscInt *rt, PetscInt *Nt, DMPolytopeType *target[], PetscInt *size[], PetscInt *cone[], PetscInt *ornt[])
tr | - The DMPlexTransform object | |
source | - The source cell type | |
p | - The source point, which can also determine the refine type |
rt | - The refine type for this point | |
Nt | - The number of types produced by this point | |
target | - An array of length Nt giving the types produced | |
size | - An array of length Nt giving the number of cells of each type produced | |
cone | - An array of length Nt*size[t]*coneSize[t] giving the cell type for each point in the cone of each produced point | |
ornt | - An array of length Nt*size[t]*coneSize[t] giving the orientation for each point in the cone of each produced point |
Note: The cone array gives the cone of each subcell listed by the first three outputs. For each cone point, we need the cell type, point identifier, and orientation within the subcell. The orientation is with respect to the canonical division (described in these outputs) of the cell in the original mesh. The point identifier is given by
the number of cones to be taken, or 0 for the current cell
the cell cone point number at each level from which it is subdivided
the replica number r of the subdivision.The orientation is with respect to the canonical cone orientation. For example, the prescription for edge division is
Nt = 2
target = {DM_POLYTOPE_POINT, DM_POLYTOPE_SEGMENT}
size = {1, 2}
cone = {DM_POLYTOPE_POINT, 1, 0, 0, DM_POLYTOPE_POINT, 0, 0, DM_POLYTOPE_POINT, 0, 0, DM_POLYTOPE_POINT, 1, 1, 0}
ornt = { 0, 0, 0, 0}