Disable AD for newton / continuation?

[dchang24]

Is it possible to disable automatic differentiation and use finite difference instead when finding equilibrium (i.e. find u such that F(u, par)=0) using newton or continuation?

[rveltz]

Yes, you pass the jacobian to the problem. Example:

prob = BifurcationProblem(F_chan, sol0, par, (@lens _.α); 
      J = (x,p) -> BK.finiteDifferences(z -> F_chan(z,p),x), )

but you can pass your own.

More help with ?BifurcationProblem