diff --git a/lib/OrdinaryDiffEqNonlinearSolve/src/newton.jl b/lib/OrdinaryDiffEqNonlinearSolve/src/newton.jl index 20c92ae5aa..b3cb37a5da 100644 --- a/lib/OrdinaryDiffEqNonlinearSolve/src/newton.jl +++ b/lib/OrdinaryDiffEqNonlinearSolve/src/newton.jl @@ -222,7 +222,7 @@ end make_new_W = is_always_new(nlsolver) || (iter == 1 && new_W) linres = dolinsolve(integrator, linsolve; A = make_new_W ? W : nothing, b = _vec(b), - linu = dz, reltol) + linu = _vec(dz), reltol) cache.linsolve = linres.cache if DiffEqBase.has_stats(integrator)