I'm fine with that. On Mar 3, 2012, at 8:06 PM, Jed Brown wrote:
> 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.
