some import fixes

Author: ChrisRackauckas

lib/OrdinaryDiffEqExponentialRK/src/OrdinaryDiffEqExponentialRK.jl

@@ -10,7 +10,7 @@ using RecursiveArrayTools
 using MuladdMacro, FastBroadcast
 using LinearAlgebra: axpy!, mul!
 using DiffEqBase, SciMLBase
-using ExponentialUtilities: _phiv_timestep_caches
+using ExponentialUtilities
 
 include("algorithms.jl")
 include("alg_utils.jl")

lib/OrdinaryDiffEqLinear/src/OrdinaryDiffEqLinear.jl

@@ -9,7 +9,6 @@ using LinearAlgebra: mul!
 using DiffEqBase
 using SciMLOperators: AbstractSciMLOperator
 using ExponentialUtilities
-using ExponentialUtilities: _phiv_timestep_caches
 
 include("algorithms.jl")
 include("alg_utils.jl")

lib/OrdinaryDiffEqLinear/src/linear_caches.jl

@@ -562,6 +562,17 @@ end
     exp_cache::expType
 end
 
+function _phiv_timestep_caches(u_prototype, maxiter::Int, p::Int)
+    n = length(u_prototype)
+    T = eltype(u_prototype)
+    u = zero(u_prototype)                         # stores the current state
+    W = Matrix{T}(undef, n, p + 1)                  # stores the w vectors
+    P = Matrix{T}(undef, n, p + 2)                  # stores output from phiv!
+    Ks = KrylovSubspace{T}(n, maxiter)            # stores output from arnoldi!
+    phiv_cache = PhivCache(u_prototype, maxiter, p + 1) # cache used by phiv! (need +1 for error estimation)
+    return u, W, P, Ks, phiv_cache
+end
+
 function alg_cache(alg::LinearExponential, u, rate_prototype, ::Type{uEltypeNoUnits},
         ::Type{uBottomEltypeNoUnits},
         ::Type{tTypeNoUnits}, uprev, uprev2, f, t, dt, reltol, p, calck,
@@ -585,4 +596,4 @@ function alg_cache(alg::LinearExponential, u, rate_prototype, ::Type{uEltypeNoUn
     end
     exp_cache = ExponentialUtilities.alloc_mem(f, ExpMethodGeneric())
     LinearExponentialCache(u, uprev, tmp, rtmp, KsCache, exp_cache)
-end
+end