Hello, I was trying to test NonlinearSolve.jl with automatic differentiation like Enzyme.jl or Zygote.jl. But the following simple code fails f(x, p) = p.a * x^2 - p.b prob = NonlinearProblem(f, 0.0, (a=1.0, b=2.0)) solve(prob, NewtonRaphson(autodiff=AutoEnzyme())) It returns zero, which is not t…

Using NonlinearSolve with autodiff

isaacsas October 15, 2024, 5:46pm 5

Note that f'(x_0) = 0 here so you might want to use another initial guess. Even using the in-place version I get that NewtonRaphson stalls on this problem.

Using NonlinearSolve with autodiff