use undefmatrix to init jacobian

SciML · Nov 28, 2022 · 4115605 · 4115605
1 parent d721e0a
commit 4115605
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 2 deletions.
diff --git a/Project.toml b/Project.toml
@@ -19,7 +19,7 @@ StaticArraysCore = "1e83bf80-4336-4d27-bf5d-d5a4f845583c"
 UnPack = "3a884ed6-31ef-47d7-9d2a-63182c4928ed"
 
 [compat]
-ArrayInterfaceCore = "0.1.1"
+ArrayInterfaceCore = "0.1.22"
 DiffEqBase = "6"
 FiniteDiff = "2"
 ForwardDiff = "0.10.3"

diff --git a/src/raphson.jl b/src/raphson.jl
@@ -54,7 +54,7 @@ end
 
 function jacobian_caches(alg::NewtonRaphson, f, u, p, ::Val{true})
     uf = JacobianWrapper(f, p)
-    J = false .* u .* u'
+    J = ArrayInterfaceCore.undefmatrix(u)
 
     linprob = LinearProblem(J, _vec(zero(u)); u0 = _vec(zero(u)))
     weight = similar(u)