http://www.mcs.anl.gov/petsc/documentation/faq.html#functionjacobian
> On Apr 18, 2016, at 8:34 AM, Hoang Giang Bui <[email protected]> wrote: > > 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 >
