Update ode_firk_tests.jl

Shreyas-Ekanathan · Aug 27, 2024 · 064c77d · 064c77d
1 parent aaeed16
commit 064c77d
Showing 1 changed file with 7 additions and 3 deletions.
diff --git a/lib/OrdinaryDiffEqFIRK/test/ode_firk_tests.jl b/lib/OrdinaryDiffEqFIRK/test/ode_firk_tests.jl
@@ -8,9 +8,8 @@ for prob in [prob_ode_linear, prob_ode_2Dlinear]
     @test sim21.𝒪est[:final]≈5 atol=testTol
 end
 
-
 sol = solve(prob_ode_2Dlinear, AdaptiveRadau(), adaptive = false, dt = 1e-2)
-sol = solve(prob_ode_2Dlinear, RadauIIA9(), adaptive = false, dt = 1e-2)
+sol = solve(prob_ode_linear, RadauIIA9(), adaptive = false, dt = 1e-2)
 sol = solve(prob_ode_2Dlinear, RadauIIA5(), adaptive = false, dt = 1e-2)
 
 sim21 = test_convergence(1 ./ 2 .^ (2.5:-1:0.5), prob_ode_linear, RadauIIA9())
@@ -19,9 +18,14 @@ sim21 = test_convergence(1 ./ 2 .^ (2.5:-1:0.5), prob_ode_linear, RadauIIA9())
 sim21 = test_convergence(1 ./ 2 .^ (2.5:-1:0.5), prob_ode_2Dlinear, RadauIIA9())
 @test sim21.𝒪est[:final]≈9 atol=testTol
 
-sim21 = test_convergence(1 ./ 2 .^ (2.25:-1:0.25), prob_ode_linear, AdaptiveRadau())
+prob_ode_linear_big = remake(prob_ode_linear, u0 = big.(prob_ode_linear.u0), tspan = big.(prob_ode_linear.tspan))
+prob_ode_2Dlinear_big = remake(prob_ode_2Dlinear, u0 = big.(prob_ode_2Dlinear.u0), tspan = big.(prob_ode_2Dlinear.tspan))
+
+sol = solve(prob_ode_linear_big, AdaptiveRadau(), adaptive=false, dt = 1e-2)
+sim21 = test_convergence(1 ./ 2 .^ (2.25:-1:0.25), prob_ode_linear_big, AdaptiveRadau())
 @test sim21.𝒪est[:final]≈9 atol=testTol
 
+sol = solve(prob_ode_2Dlinear_big, AdaptiveRadau(), adaptive=false, dt = 1e-2)
 sim21 = test_convergence(1 ./ 2 .^(2.25:-1:0.25), prod_ode_2Dlinear, AdaptiveRadau())
 @test sim21.𝒪est[:final]≈9 atol=testTol