Bump sparspak dependency to 0.3.4, add factorize=false in sparspaklu call.

Author: j-fu

Project.toml

@@ -39,7 +39,7 @@ Reexport = "1"
 SciMLBase = "1.68"
 Setfield = "0.7, 0.8, 1"
 SnoopPrecompile = "1"
-Sparspak = "0.3"
+Sparspak = "0.3.4"
 UnPack = "1"
 julia = "1.6"
 

src/factorization.jl

@@ -511,7 +511,7 @@ function init_cacheval(::SparspakFactorization, A, b, u, Pl, Pr, maxiters::Int,
                        reltol,
                        verbose::Bool, assumptions::OperatorAssumptions)
     A = convert(AbstractMatrix, A)
-    return sparspaklu(A)
+    return sparspaklu(A, factorize=false)
 end
 
 function SciMLBase.solve(cache::LinearCache, alg::SparspakFactorization; kwargs...)