> Do not pass -snes_mf or -snes_mf_operator with the -snes_type test option Sorry, I am lost. I provide an analytical Jacobian to a polynomial function. Sure I can make mistakes even in such trivial derivations, but how would I get the results I expect from a book without the option? Dominik