On Mon, Apr 18, 2016 at 9: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. > The Jacobian and Residual are not always formed together, so you would have to stash your computation in a user context and then check that the same X had been passed to both functions. For example, the Residual can be called many times in a line search routine. Matt > 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 > > -- What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead. -- Norbert Wiener
