Hi, If I have a pointwise Jacobian function f, I know that I can call DMGetDS() and pass f to PetscDSSetJacobian(), and that f will be used by PetscFEIntegrate() and thus by DMPlexSNESComputeJacobianFEM(). It looks like PetscFEIntegrate() is only used by plex and not da. Is there anyway that I can use f to construct a Jacobian with a da object?
Toby
signature.asc
Description: Digital signature
