densify sparse W for inv

SciML · Nov 27, 2024 · be5fc08 · be5fc08
1 parent 07e9ffd
commit be5fc08
Showing 1 changed file with 2 additions and 1 deletion.
diff --git a/lib/OrdinaryDiffEqDifferentiation/src/derivative_utils.jl b/lib/OrdinaryDiffEqDifferentiation/src/derivative_utils.jl
@@ -17,7 +17,8 @@ struct StaticWOperator{isinv, T, F} <: AbstractSciMLOperator{T}
         # when constructing W for the first time for the type
         # inv(W) can be singular
         _W = if isinv && callinv
-            inv(W)
+            # W may be sparse, needs to be dense to inv
+            inv(Matrix(W))
         else
             W
         end