fix default linear solver in DifferentialEquations interface

j-fu · Jan 31, 2022 · 14a0f11 · 14a0f11 · j-fu · Jan 31, 2022
1 parent 2af10f2
commit 14a0f11
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 3 deletions.
diff --git a/Project.toml b/Project.toml
@@ -1,7 +1,7 @@
 name = "VoronoiFVM"
 uuid = "82b139dc-5afc-11e9-35da-9b9bdfd336f3"
 authors = ["Juergen Fuhrmann <juergen.fuhrmann@wias-berlin.de>"]
-version = "0.15.2-dev"
+version = "0.15.2"
 
 [deps]
 DiffResults = "163ba53b-c6d8-5494-b064-1a9d43ac40c5"

diff --git a/src/vfvm_diffeq_interface.jl b/src/vfvm_diffeq_interface.jl
@@ -195,7 +195,7 @@ Alias for [`solve(system::VoronoiFVM.AbstractSystem)`](@ref) with the correspond
                tspan; solver=nothing, kwargs...)
 
     if isnothing(solver)
-        solver=DiffEq.Rosenbrock23()
+        solver=DiffEq.Rosenbrock23(linsolve=DiffEq.KLUFactorization())
     end
 
     sys.history=DiffEqHistory()
@@ -205,7 +205,7 @@ Alias for [`solve(system::VoronoiFVM.AbstractSystem)`](@ref) with the correspond
     flush!(sys.matrix)
     f = DiffEq.ODEFunction(eval_rhs!,
                            jac=eval_jacobian!,
-                           jac_prototype=sys.matrix,
+                           jac_prototype=sys.matrix.cscmatrix,
                            mass_matrix=mass_matrix(sys))
 
     prob = DiffEq.ODEProblem(f,vec(inival),tspan,sys)