petsc-3.8.4 2018-03-24

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.

Examples