Merge pull request #69 from SciML/myb/fix

Fix constructDormandLockyerMcCorriganPrince6

Author: YingboMa

src/ode_tableaus.jl

@@ -2381,7 +2381,7 @@ function constructDormandLockyerMcCorriganPrince6(T::Type = Float64)
   α = map(T,α)
   c = map(T,c)
   αEEst = map(T,αEEst)
-  return(ExplicitRKTableau(A,c,α,6,fsal=true))
+  return(ExplicitRKTableau(A,c,α,6,adaptiveorder=5,αEEst=αEEst,fsal=true))
 end
 
 

test/ode_tableau_convergence_tests.jl

@@ -18,7 +18,7 @@ dts = 1 .//2 .^(8:-1:4)
 testTol = 0.3
 superduperbool = Vector{Bool}(undef, 2)
 
-for tab in [constructRalston4(), constructTsitourasPapakostas6()]
+for tab in [constructRalston4(Rational{BigInt}), constructTsitourasPapakostas6(Rational{BigInt}), constructDormandLockyerMcCorriganPrince6(Rational{BigInt})]
   @test all(i -> residual_order_condition(tab, i, +, abs) < 10eps(1.0), 1:tab.order)
   if tab.adaptiveorder != 0
     @test all(i -> residual_order_condition(tab, i, +, abs; embedded=true) < 10eps(1.0), tab.adaptiveorder)