From 6bc1bd0839ed6ca7d5755ab71f5fb12de7627c0e Mon Sep 17 00:00:00 2001
From: oscarddssmith <oscar.smith@juliacomputing.com>
Date: Mon, 28 Nov 2022 11:15:29 -0500
Subject: [PATCH] use undefmatrix to init jacobian

---
 Project.toml   | 2 +-
 src/raphson.jl | 2 +-
 2 files changed, 2 insertions(+), 2 deletions(-)

diff --git a/Project.toml b/Project.toml
index 7639b7e0a..8228a0d69 100644
--- 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
index d61878683..f5f507dcf 100644
--- 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)