fixes

Author: oscardssmith

Project.toml

@@ -121,8 +121,8 @@ OrdinaryDiffEqQPRK = "1"
 OrdinaryDiffEqRKN = "1"
 OrdinaryDiffEqRosenbrock = "1"
 OrdinaryDiffEqSDIRK = "1"
-OrdinaryDiffEqStabilizedIRK = "1"
 OrdinaryDiffEqSSPRK = "1"
+OrdinaryDiffEqStabilizedIRK = "1"
 OrdinaryDiffEqStabilizedRK = "1"
 OrdinaryDiffEqSymplecticRK = "1"
 OrdinaryDiffEqTsit5 = "1"

lib/OrdinaryDiffEqNonlinearSolve/src/newton.jl

@@ -220,15 +220,9 @@ end
         reltol = eps(eltype(dz))
     end
 
-    if is_always_new(nlsolver) || (iter == 1 && new_W)
-        linres = dolinsolve(integrator, linsolve; A = W, b = _vec(b), linu = _vec(dz),
-            reltol = reltol)
-    else
-        linres = dolinsolve(
-            integrator, linsolve; A = nothing, b = _vec(b), linu = _vec(dz),
-            reltol = reltol)
-    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),
+        reltol)
     cache.linsolve = linres.cache
 
     if DiffEqBase.has_stats(integrator)