simplify ODE test dependency

Author: Datseris

src/core_systems/continuous_time_ode.jl

@@ -141,7 +141,7 @@ set_state!(ds::CoupledODEs, u::AbstractArray) = (set_state!(ds.integ, u); ds)
 # so that `ds` is printed
 SciMLBase.step!(ds::CoupledODEs, args...) = (step!(ds.integ, args...); ds)
 
-function SciMLBase.reinit!(ds::ContinuousTimeDynamicalSystem, u = initial_state(ds);
+function SciMLBase.reinit!(ds::ContinuousTimeDynamicalSystem, u::AbstractArray = initial_state(ds);
         p = current_parameters(ds), t0 = initial_time(ds), kw...
     )
     set_parameters!(ds, p)

test/Project.toml

@@ -3,7 +3,8 @@ DiffEqNoiseProcess = "77a26b50-5914-5dd7-bc55-306e6241c503"
 DynamicalSystemsBase = "6e36e845-645a-534a-86f2-f5d4aa5a06b4"
 LinearAlgebra = "37e2e46d-f89d-539d-b4ee-838fcccc9c8e"
 ModelingToolkit = "961ee093-0014-501f-94e3-6117800e7a78"
-OrdinaryDiffEq = "1dea7af3-3e70-54e6-95c3-0bf5283fa5ed"
+OrdinaryDiffEqRosenbrock = "43230ef6-c299-4910-a778-202eb28ce4ce"
+OrdinaryDiffEqTsit5 = "b1df2697-797e-41e3-8120-5422d3b24e4a"
 Statistics = "10745b16-79ce-11e8-11f9-7d13ad32a3b2"
 StochasticDiffEq = "789caeaf-c7a9-5a7d-9973-96adeb23e2a0"
 Test = "8dfed614-e22c-5e08-85e1-65c5234f0b40"

test/continuous.jl

@@ -1,6 +1,7 @@
 using DynamicalSystemsBase, Test
 
-using OrdinaryDiffEq: Vern9, ODEProblem, Rodas5, Tsit5
+using OrdinaryDiffEqTsit5: ODEProblem, Tsit5
+using OrdinaryDiffEqRosenbrock: Rodas5
 
 include("test_system_function.jl")
 
@@ -25,14 +26,9 @@ p0 = [10, 28, 8/3]
 
 lorenz_oop = CoupledODEs(lorenz_rule, u0, p0)
 lorenz_iip = CoupledODEs(ODEProblem(lorenz_rule_iip, copy(u0), (0.0, Inf), p0))
-lorenz_vern = CoupledODEs(lorenz_rule, u0, p0;
-    diffeq = (alg = Vern9(), abstol = 1e-9, reltol = 1e-9)
-)
 
-for (ds, iip) in zip((lorenz_oop, lorenz_iip, lorenz_vern), (false, true, false))
-
-    name = (ds === lorenz_vern) ? "lorvern" : "lorenz"
-    @testset "$(name) IIP=$(iip)" begin
+for (ds, iip) in zip((lorenz_oop, lorenz_iip), (false, true))
+    @testset "IIP=$(iip)" begin
         @test dynamic_rule(ds) == (iip ? lorenz_rule_iip : lorenz_rule)
         test_dynamical_system(ds, u0, p0; idt = false, iip=iip)
     end
@@ -42,12 +38,6 @@ end
     lorenz_oop = CoupledODEs(lorenz_rule, u0, p0)
     @test lorenz_oop.integ.alg isa Tsit5
 
-    lorenz_vern = CoupledODEs(lorenz_rule, u0, p0;
-        diffeq = (alg = Vern9(), verbose = false, abstol = 1e-9, reltol = 1e-9)
-    )
-    @test lorenz_vern.integ.alg isa Vern9
-    @test lorenz_vern.integ.opts.verbose == false
-
     # also test ODEproblem creation
     prob = lorenz_vern.integ.sol.prob