add default

Author: vyudu

test/extensions/lattice_simulation_plotting.jl

@@ -1,7 +1,7 @@
 ### Preparations ###
 
 # Fetch packages.
-using Catalyst, CairoMakie, GraphMakie, Graphs, JumpProcesses, OrdinaryDiffEq, Test
+using Catalyst, CairoMakie, GraphMakie, Graphs, JumpProcesses, OrdinaryDiffEqDefault, OrdinaryDiffEq, Test
 
 
 ### Checks Basic Plot Cases ###
@@ -91,7 +91,7 @@ let
     lattice = CartesianGrid((2,2,2))
     lrs = LatticeReactionSystem(rs, [diffusion_rx], lattice)
     oprob = ODEProblem(lrs, [:X => 1.0], 1.0, [:d => 1.0, :D => 0.2])
-    osol = solve(oprob, Tsit5())
+    osol = solve(oprob)
 
     @test_throws ArgumentError lattice_plot(osol, :X, lrs)
     @test_throws ArgumentError lattice_animation(osol, :X, lrs, "animation_tmp.mp4")

test/network_analysis/conservation_laws.jl

@@ -1,7 +1,7 @@
 ### Prepares Tests ###
 
 # Fetch packages.
-using Catalyst, JumpProcesses, LinearAlgebra, NonlinearSolve, OrdinaryDiffEqTsit5, SteadyStateDiffEq, StochasticDiffEq, Test
+using Catalyst, JumpProcesses, LinearAlgebra, NonlinearSolve, OrdinaryDiffEqTsit5, OrdinaryDiffEqDefault, OrdinaryDiffEqVerner, SteadyStateDiffEq, StochasticDiffEq, Test
 
 # Sets stable rng number.
 using StableRNGs
@@ -192,7 +192,7 @@ let
     oprob1 = ODEProblem(rn, u0_1, 10.0, ps; remove_conserved = true, remove_conserved_warn)
     oprob2 = ODEProblem(rn, u0_2, 10.0, ps; remove_conserved = true, remove_conserved_warn)
     oprob3 = ODEProblem(rn, u0_3, 10.0, ps; remove_conserved = true, remove_conserved_warn)
-    @test solve(oprob1, Tsit5())[sps] ≈ solve(oprob2, Tsit5())[sps] ≈ solve(oprob3, Tsit5())[sps]
+    @test solve(oprob1)[sps] ≈ solve(oprob2)[sps] ≈ solve(oprob3)[sps]
 end
 
 # Tests conservation laws in SDE simulation.
@@ -224,8 +224,8 @@ let
     oprob2 = ODEProblem(osys, u0, 10.0, ps_2)
 
     # Checks that the solutions generates the correct conserved quantities.
-    sol1 = solve(oprob1, Tsit5(); saveat = 1.0)
-    sol2 = solve(oprob2, Tsit5(); saveat = 1.0)
+    sol1 = solve(oprob1; saveat = 1.0)
+    sol2 = solve(oprob2; saveat = 1.0)
     @test all(sol1[osys.X1 + osys.X2] .== 2.0)
     @test all(sol2[osys.X1 + osys.X2] .== 4.0)
 end