Newton Raphson ignores jac-prototype

avik-pal · avik-pal · commit b0aa0afb46f9 · 2023-08-18T12:37:44.000-04:00
diff --git a/.gitignore b/.gitignore
@@ -23,3 +23,5 @@ docs/site/
 # environment.
 Manifest.toml
 docs/src/assets/Project.toml
+
+.vscode
diff --git a/Project.toml b/Project.toml
@@ -1,7 +1,7 @@
 name = "NonlinearSolve"
 uuid = "8913a72c-1f9b-4ce2-8d82-65094dcecaec"
 authors = ["SciML"]
-version = "1.9.0"
+version = "1.10.0"
 
 [deps]
 ArrayInterface = "4fba245c-0d91-5ea0-9b3e-6abc04ee57a9"
diff --git a/src/raphson.jl b/src/raphson.jl
@@ -106,7 +106,11 @@ end
 
 function jacobian_caches(alg::NewtonRaphson, f, u, p, ::Val{true})
     uf = JacobianWrapper(f, p)
-    J = ArrayInterface.undefmatrix(u)
+    J = if f.jac_prototype === nothing
+        ArrayInterface.undefmatrix(u)
+    else
+        f.jac_prototype
+    end
 
     linprob = LinearProblem(J, _vec(zero(u)); u0 = _vec(zero(u)))
     weight = similar(u)
