<div class="gmail_quote">On Sat, Mar 3, 2012 at 18:51, Barry Smith <span dir="ltr"><<a href="mailto:bsmith@mcs.anl.gov">bsmith@mcs.anl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
 If no compute Jacobian routine was provided it just always uses the original matrix that was set into the SNESSetJacobian.  I can live with removing this feature, I agree it is a little dangerous.</blockquote></div><div>
<br></div><div>What about making it explicit?</div><br><div>SNESSetJacobian(...,SNESDefaultComputeJacobianConstant,PETSC_NULL);</div><div><br></div><div>This routine would explicitly leave the Jacobian unchanged.</div>