<div dir="ltr">Oh man!  Sorry Barry!  I swear I looked around before I sent the email.  I should have checked the FAQ a little more closely!<div><br></div><div>I can understand the reasoning in the FAQ... but I still wonder if it might not be useful to provide all three options (Function, Jacobian, FunctionJacobian).  In my case I could fill in each one to do the right thing.  That way PETSc could call the "FunctionJacobian" one when it knew it needed both (by default that could just farm out to the individual calls).  But you guys have definitely thought a lot more about this than I have.</div><div><br></div><div>So, do you still recommend what's suggested in the FAQ?  Save off the Jacobian computation during the residual computation and then use that when SNES asks for a Jacobian?  In the case of automatic differentiation this could make a pretty huge difference in time...</div><div><br></div><div>Derek</div><div><div><br><div class="gmail_quote"><div dir="ltr">On Fri, Dec 9, 2016 at 1:49 PM Barry Smith <<a href="mailto:bsmith@mcs.anl.gov">bsmith@mcs.anl.gov</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><br class="gmail_msg">
  Sorry the title in the FAQ is a bit tongue-in-check.<br class="gmail_msg">
<br class="gmail_msg">
  <a href="http://www.mcs.anl.gov/petsc/documentation/faq.html#functionjacobian" rel="noreferrer" class="gmail_msg" target="_blank">http://www.mcs.anl.gov/petsc/documentation/faq.html#functionjacobian</a><br class="gmail_msg">
<br class="gmail_msg">
<br class="gmail_msg">
> On Dec 9, 2016, at 12:45 PM, Derek Gaston <<a href="mailto:friedmud@gmail.com" class="gmail_msg" target="_blank">friedmud@gmail.com</a>> wrote:<br class="gmail_msg">
><br class="gmail_msg">
> Is there a way to tell SNES to simultaneously compute both the residual and the Jacobian in one callback?<br class="gmail_msg">
><br class="gmail_msg">
> 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).<br class="gmail_msg">
><br class="gmail_msg">
> 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.<br class="gmail_msg">
><br class="gmail_msg">
> Thanks for any help!<br class="gmail_msg">
><br class="gmail_msg">
> Derek<br class="gmail_msg">
<br class="gmail_msg">
</blockquote></div></div></div></div>