<div dir="ltr">Is there a way to tell SNES to simultaneously compute both the residual and the Jacobian in one callback?<div><br></div><div>My code can compute both simultaneously and it will be more efficient (think FE where you can reuse the shape-functions, variables, material properties, etc. for both residual and Jacobian computation).  In addition, I also have automatic differentiation as an option which _definitely_ computes both efficiently (and actually computes residuals, by themselves, much slower).</div><div><br></div><div>I was thinking that I may just save off the Jacobian whenever the initial residual computation is asked for by SNES... and then just return that Jacobian when SNES asks for it.  This may be a bit dicey though as SNES can ask for residual computations at many different points during the solve.</div><div><br></div><div>Thanks for any help!</div><div><br></div><div>Derek</div></div>