#include "petscdmplextransform.h" PetscErrorCode DMPlexRefineRegularGetAffineFaceTransforms(DMPlexTransform tr, DMPolytopeType ct, PetscInt *Nf, PetscReal *v0[], PetscReal *J[], PetscReal *invJ[], PetscReal *detJ[])
cr | - The DMPlexCellRefiner object | |
ct | - The cell type |
Nf | - The number of faces for this cell type | |
v0 | - The translation of the first vertex for each face | |
J | - The Jacobian for each face (map from original cell to subcell) | |
invJ | - The inverse Jacobian for each face | |
detJ | - The determinant of the Jacobian for each face |
Note: The Jacobian and inverse Jacboian will be rectangular, and the inverse is really a generalized inverse.
v0 + j x_face = x_cell
invj (x_cell - v0) = x_face