petsc-3.4.5 2014-06-29

DMPlexComputeJacobianFEM

Form the local portion of the Jacobian matrix J at the local solution X using pointwise functions specified by the user.

Synopsis

#include "petscdmplex.h"   
PetscErrorCode DMPlexComputeJacobianFEM(DM dm, Vec X, Mat Jac, Mat JacP, MatStructure *str,void *user)

Input Parameters

	dm	- The mesh
	X	- Local input vector
	user	- The user context

Output Parameter

Jac -Jacobian matrix

Note

The second member of the user context must be an FEMContext.

We form the residual one batch of elements at a time. This allows us to offload work onto an accelerator, like a GPU, or vectorize on a multicore machine.

Examples