Make package compatible with OrdinaryDiffEq v6.6 and v6.7 (#130)

Author: nathanaelbosch

Project.toml

@@ -26,7 +26,7 @@ DiffEqBase = "6"
 DiffEqDevTools = "2"
 ForwardDiff = "0.10"
 GaussianDistributions = "0.5"
-OrdinaryDiffEq = "5.71 - 6.5"
+OrdinaryDiffEq = "6.2 - 6.7"
 RecipesBase = "1"
 RecursiveArrayTools = "2"
 Reexport = "0.2, 1.0"

src/alg_utils.jl

@@ -8,6 +8,10 @@ OrdinaryDiffEq.alg_autodiff(alg::EK1{CS,AD}) where {CS,AD} = AD
 OrdinaryDiffEq.alg_difftype(alg::EK1{CS,AD,DiffType}) where {CS,AD,DiffType} = DiffType
 OrdinaryDiffEq.standardtag(alg::AbstractEK) = false
 OrdinaryDiffEq.standardtag(alg::EK1{CS,AD,DiffType,ST}) where {CS,AD,DiffType,ST} = ST
+OrdinaryDiffEq.concrete_jac(alg::AbstractEK) = nothing
+OrdinaryDiffEq.concrete_jac(alg::EK1{CS,AD,DiffType,ST,CJ}) where {CS,AD,DiffType,ST,CJ} =
+    CJ
+
 @inline DiffEqBase.get_tmp_cache(integ, alg::EK1, cache) = (cache.tmp, cache.atmp)
 OrdinaryDiffEq.get_chunksize(alg::AbstractEK) = Val(0)
 OrdinaryDiffEq.isfsal(alg::AbstractEK) = false

src/algorithms.jl

@@ -48,7 +48,7 @@ ForwardDiff.jl to compute Jacobians.
 
 # [References](@ref references)
 """
-struct EK1{CS,AD,DiffType,ST,DT,IT} <: AbstractEK
+struct EK1{CS,AD,DiffType,ST,CJ,DT,IT} <: AbstractEK
     order::Int
     diffusionmodel::DT
     smooth::Bool
@@ -63,11 +63,13 @@ EK1(;
     autodiff=true,
     diff_type=Val{:forward},
     standardtag=true,
+    concrete_jac=nothing,
 ) = EK1{
     chunk_size,
     autodiff,
     diff_type,
     standardtag,
+    concrete_jac,
     typeof(diffusionmodel),
     typeof(initialization),
 }(