On Sat, Mar 3, 2012 at 18:51, Barry Smith <bsmith at mcs.anl.gov> wrote:
> 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. What about making it explicit? SNESSetJacobian(...,SNESDefaultComputeJacobianConstant,PETSC_NULL); This routine would explicitly leave the Jacobian unchanged. -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20120303/160bab37/attachment.html>
