petsc-3.8.4 2018-03-24
Report Typos and Errors

DMPlexComputeCellGeometryFEM

Compute the Jacobian, inverse Jacobian, and Jacobian determinant at each quadrature point in the given cell

Synopsis

#include "petscdmplex.h"   
#include "petscfe.h"       
PetscErrorCode DMPlexComputeCellGeometryFEM(DM dm, PetscInt cell, PetscQuadrature quad, PetscReal *v, PetscReal *J, PetscReal *invJ, PetscReal *detJ)
Collective on DM

Input Arguments

dm - the DM
cell - the cell
quad - the quadrature containing the points in the reference element where the geometry will be evaluated. If quad == NULL, geometry will be evaluated at the first vertex of the reference element

Output Arguments

v - the image of the transformed quadrature points, otherwise the image of the first vertex in the closure of the reference element
J - the Jacobian of the transform from the reference element at each quadrature point
invJ - the inverse of the Jacobian at each quadrature point
detJ - the Jacobian determinant at each quadrature point

Fortran Notes

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

See Also

DMGetCoordinateSection(), DMGetCoordinateVec()

Level:advanced
Location:
src/dm/impls/plex/plexgeometry.c
Index of all DMPLEX routines
Table of Contents for all manual pages
Index of all manual pages

Examples

src/dm/impls/plex/examples/tutorials/ex2.c.html