[petsc-users] snes set function and Jacobian concurrently

[Hoang Giang Bui]

Hello

Is it possible to combine the SNESSetFunction and SNESSetJacobian into one
call, i.e. SNESSetFunctionAndJacobian which takes unique function pointer
to compute rhs and Jacobian concurrently? The reason for it is that the
finite element code I'm working on contribute to global rhs and stiffness
at one call, so to cater for SNES would break the existing interface.

Another supporting point for that feature is that for irreversible problem,
the rhs (equilirium) and Jacobian shall not be decoupled, since they are
not able to return to the past state after proceed.

Giang
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20160418/1b16f8a5/attachment-0001.html>