Merge pull request #2538 from mcarmesin/master

Adapt KrylovSubspace Vtype to u
SciML · Nov 19, 2024 · 9ab927b · 9ab927b
2 parents 52bcf16 + 28313cd
commit 9ab927b
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 2 deletions.
diff --git a/lib/OrdinaryDiffEqLinear/Project.toml b/lib/OrdinaryDiffEqLinear/Project.toml
@@ -18,7 +18,7 @@ SciMLOperators = "c0aeaf25-5076-4817-a8d5-81caf7dfa961"
 [compat]
 DiffEqBase = "6.152.2"
 DiffEqDevTools = "2.44.4"
-ExponentialUtilities = "1.26.1"
+ExponentialUtilities = "1.27"
 LinearAlgebra = "<0.0.1, 1"
 OrdinaryDiffEqCore = "1.1"
 OrdinaryDiffEqRosenbrock = "<0.0.1, 1"

diff --git a/lib/OrdinaryDiffEqLinear/src/linear_caches.jl b/lib/OrdinaryDiffEqLinear/src/linear_caches.jl
@@ -591,7 +591,7 @@ function alg_cache(alg::LinearExponential, u, rate_prototype, ::Type{uEltypeNoUn
     if alg.krylov == :off
         KsCache = nothing
     elseif alg.krylov == :simple
-        Ks = KrylovSubspace{T}(n, m)
+        Ks = KrylovSubspace{T,T,typeof(similar(u,size(u,1),2))}(n, m)
         expv_cache = ExpvCache{T}(m)
         KsCache = (Ks, expv_cache)
     elseif alg.krylov == :adaptive