Fix SteadyState Adjoint for OOP jacobian

SciML · Mar 28, 2024 · 9202505 · 9202505
1 parent 10b49a4
commit 9202505
Showing 1 changed file with 5 additions and 1 deletion.
diff --git a/src/steadystate_adjoint.jl b/src/steadystate_adjoint.jl
@@ -55,7 +55,11 @@ end
 
     if needs_jac
         if DiffEqBase.has_jac(f)
-            f.jac(diffcache.J, y, p, nothing)
+            if DiffEqBase.isinplace(sol.prob)
+                f.jac(diffcache.J, y, p, nothing)
+            else
+                copyto!(diffcache.J, f.jac(y, p, nothing))
+            end
         else
             if DiffEqBase.isinplace(sol.prob)
                 jacobian!(diffcache.J, diffcache.uf, y, diffcache.f_cache,