[petsc-users] -snes_mf and -snes_fd

Barry Smith bsmith at mcs.anl.gov
Fri May 6 08:58:03 CDT 2016


> On May 5, 2016, at 11:03 PM, Jed Brown <jed at jedbrown.org> wrote:
> 
> Matthew Knepley <knepley at gmail.com> writes:
>>> without -snes_mf AT ALL but with -snes_fd 0 I get RHS called 3 times [1]
>>> 
>> 
>> I need to see the output of -snes_view for these cases. Its possible that we
>> have a bug related to the argument.
> 
> This case looks correct to me.  It's the same as passing no options, so
> the analytic Jacobian gets used.  That said, our options scheme is lame
> to use separate options when only one can possibly be used.  We should
> have -snes_jacobian <analytic,mf,mf_operator,fd,fd_color> or similar.

   Perhaps it should have two parts; one for the mat and one for the pmat. So -snes_jacobian X[:Y] where X can be analytic, mf,fd,fd_color and Y can be analytic,fd,fd_color,none

   Make it an issue.

Barry




More information about the petsc-users mailing list