Specialize UMFPack initialization

Author: ChrisRackauckas

Project.toml

@@ -1,7 +1,7 @@
 name = "LinearSolve"
 uuid = "7ed4a6bd-45f5-4d41-b270-4a48e9bafcae"
 authors = ["SciML"]
-version = "0.3.0"
+version = "0.3.1"
 
 [deps]
 ArrayInterface = "4fba245c-0d91-5ea0-9b3e-6abc04ee57a9"

src/factorization.jl

@@ -45,6 +45,16 @@ Base.@kwdef struct UMFPACKFactorization <: AbstractFactorization
     reuse_symbolic::Bool = true
 end
 
+function init_cacheval(alg::UMFPACKFactorization, A, b, u, Pl, Pr, maxiters, abstol, reltol, verbose)
+    zerobased = SparseArrays.getcolptr(A)[1] == 0
+    res = SuiteSparse.UMFPACK.UmfpackLU(C_NULL, C_NULL, size(A, 1), size(A, 2),
+                    zerobased ? copy(SparseArrays.getcolptr(A)) : SuiteSparse.decrement(SparseArrays.getcolptr(A)),
+                    zerobased ? copy(rowvals(A)) : SuiteSparse.decrement(rowvals(A)),
+                    copy(nonzeros(A)), 0)
+    finalizer(SuiteSparse.UMFPACK.umfpack_free_symbolic, res)
+    res
+end
+
 function do_factorization(::UMFPACKFactorization, A, b, u)
     if A isa AbstractDiffEqOperator
         A = A.A