I’m struggling to locate the solver stats in NonlinearSolve.jl. After running the basic example from the docs,

using NonlinearSolve function f(u, p) u .* u .- p end u0 = [1.0, 1.0] p = 2.0 prob = NonlinearProblem(f, u0, p) sol = solve(prob, NewtonRaphson())

I’m finding that sol.stats is empty. Is there a keyword I need to set?

What version of NonlinearSolve and SciMLBase?

SciMLBase v1.97.0 NonlinearSolve v1.6.0

It was added in NonlinearSolve.jl v3.

I see, thank you.

it’s also probably worth noting that prior to version 3, Nonlinear solve had a number of bugs that made it fail to converge or converge much slower than it should. As such, I strongly recommend upgrading.