Re: [firedrake] Jacobian wrong with -snes_fd_color

On 23/11/15 08:55, Buesing, Henrik wrote:

So, the only way to deal with this is the following?

Use my functions to calculate the values pointwise. Then, provide the correct derivatives via the coefficient_derivatives keyword in the derivative() function.

I think so. If you need to do this update once per nonlinear step, at the moment you will have to write your own nonlinear iteration (e.g. fixed point). Since we do not really have support for hooking arbitrary callbacks into the nonlinear solver pipeline. Lawrence