Widen RKN endpoint atol and Vern7 DDE error threshold for CVHin

RKN nystrom_convergence_tests.jl: Add atol=1e-4 to endpoint comparisons
for adaptive IIP vs OOP tests. The different initial steps from CVHin
cause solutions to diverge at ~1e-6 level by t=10, exceeding default
isapprox tolerance.

DelayDiffEq verner.jl: Widen Vern7 final error threshold from 3.5e-7
to 5.0e-7 (observed 4.22e-7 with CVHin initial step).

Co-Authored-By: Chris Rackauckas <accounts@chrisrackauckas.com>

Author: ChrisRackauckas

lib/DelayDiffEq/test/integrators/verner.jl

@@ -31,7 +31,7 @@ using Test
         sol_ip = solve(prob_ip, alg)
 
         @test sol_ip.errors[:l∞] < 6.0e-4
-        @test sol_ip.errors[:final] < 3.5e-7
+        @test sol_ip.errors[:final] < 5.0e-7
         @test sol_ip.errors[:l2] < 3.0e-4
 
         sol_scalar = solve(prob_scalar, alg)

lib/OrdinaryDiffEqRKN/test/nystrom_convergence_tests.jl

@@ -416,7 +416,7 @@ end
         # due to FP rounding differences in initdt and step controller
         sol_i = solve(ode_i, alg)
         sol_o = solve(ode_o, alg)
-        @test sol_i.u[end] ≈ sol_o.u[end]
+        @test sol_i.u[end] ≈ sol_o.u[end] atol = 1.0e-4
     end
 
     @testset "FineRKN5" begin
@@ -435,7 +435,7 @@ end
         # adaptive time step — IIP vs OOP may produce different step counts
         sol_i = solve(ode_i, alg)
         sol_o = solve(ode_o, alg)
-        @test sol_i.u[end] ≈ sol_o.u[end]
+        @test sol_i.u[end] ≈ sol_o.u[end] atol = 1.0e-4
     end
 
     @testset "DPRKN4" begin
@@ -454,7 +454,7 @@ end
         # adaptive time step — IIP vs OOP may produce different step counts
         sol_i = solve(ode_i, alg)
         sol_o = solve(ode_o, alg)
-        @test sol_i.u[end] ≈ sol_o.u[end]
+        @test sol_i.u[end] ≈ sol_o.u[end] atol = 1.0e-4
     end
 
     @testset "DPRKN5" begin
@@ -473,7 +473,7 @@ end
         # adaptive time step — IIP vs OOP may produce different step counts
         sol_i = solve(ode_i, alg)
         sol_o = solve(ode_o, alg)
-        @test sol_i.u[end] ≈ sol_o.u[end]
+        @test sol_i.u[end] ≈ sol_o.u[end] atol = 1.0e-4
     end
 
     @testset "DPRKN6" begin
@@ -492,7 +492,7 @@ end
         # adaptive time step — IIP vs OOP may produce different step counts
         sol_i = solve(ode_i, alg)
         sol_o = solve(ode_o, alg)
-        @test sol_i.u[end] ≈ sol_o.u[end]
+        @test sol_i.u[end] ≈ sol_o.u[end] atol = 1.0e-4
     end
 
     @testset "DPRKN6FM" begin
@@ -511,7 +511,7 @@ end
         # adaptive time step — IIP vs OOP may produce different step counts
         sol_i = solve(ode_i, alg)
         sol_o = solve(ode_o, alg)
-        @test sol_i.u[end] ≈ sol_o.u[end]
+        @test sol_i.u[end] ≈ sol_o.u[end] atol = 1.0e-4
     end
 
     @testset "DPRKN8" begin
@@ -530,7 +530,7 @@ end
         # adaptive time step — IIP vs OOP may produce different step counts
         sol_i = solve(ode_i, alg)
         sol_o = solve(ode_o, alg)
-        @test sol_i.u[end] ≈ sol_o.u[end]
+        @test sol_i.u[end] ≈ sol_o.u[end] atol = 1.0e-4
     end
 
     @testset "DPRKN12" begin
@@ -549,6 +549,6 @@ end
         # adaptive time step — IIP vs OOP may produce different step counts
         sol_i = solve(ode_i, alg)
         sol_o = solve(ode_o, alg)
-        @test sol_i.u[end] ≈ sol_o.u[end]
+        @test sol_i.u[end] ≈ sol_o.u[end] atol = 1.0e-4
     end
 end