From 14a0f1199859c3438bb206a0b45b695fef3750af Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?J=C3=BCrgen=20Fuhrmann?= Date: Mon, 31 Jan 2022 14:57:20 +0100 Subject: [PATCH] fix default linear solver in DifferentialEquations interface --- Project.toml | 2 +- src/vfvm_diffeq_interface.jl | 4 ++-- 2 files changed, 3 insertions(+), 3 deletions(-) diff --git a/Project.toml b/Project.toml index 44714a404..a4850d2ba 100644 --- a/Project.toml +++ b/Project.toml @@ -1,7 +1,7 @@ name = "VoronoiFVM" uuid = "82b139dc-5afc-11e9-35da-9b9bdfd336f3" authors = ["Juergen Fuhrmann "] -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 index a16c512e8..2c7b9337e 100644 --- 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)