From 064c77d6d01fa0c45fe88012c9d9453c27c235db Mon Sep 17 00:00:00 2001 From: Shreyas Ekanathan Date: Tue, 27 Aug 2024 17:20:25 -0400 Subject: [PATCH] Update ode_firk_tests.jl --- lib/OrdinaryDiffEqFIRK/test/ode_firk_tests.jl | 10 +++++++--- 1 file changed, 7 insertions(+), 3 deletions(-) diff --git a/lib/OrdinaryDiffEqFIRK/test/ode_firk_tests.jl b/lib/OrdinaryDiffEqFIRK/test/ode_firk_tests.jl index bead50a768..f8eeac4831 100644 --- 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