Merge pull request #1384 from ChrisRackauckas-Claude/add-irkgausslegendre-benchmarks

Add IRKGaussLegendre.jl methods to Non-stiff ODE benchmarks

Author: ChrisRackauckas

benchmarks/NonStiffODE/LotkaVolterra_wpd.jmd

@@ -9,7 +9,7 @@ The purpose of this problem is to test the performance on easy problems. Since i
 
 ```julia
 using OrdinaryDiffEq, ParameterizedFunctions, ODEInterfaceDiffEq, LSODA,
-      Sundials, DiffEqDevTools, StaticArrays, OrdinaryDiffEqSIMDRK
+      Sundials, DiffEqDevTools, StaticArrays, OrdinaryDiffEqSIMDRK, IRKGaussLegendre
 
 f = @ode_def LotkaVolterra begin
     dx = a*x - b*x*y
@@ -42,6 +42,8 @@ setups = [Dict(:alg=>DP5())
           Dict(:alg=>dopri5())
           Dict(:alg=>Tsit5())
           Dict(:alg=>Vern6())
+          Dict(:alg=>IRKGL16(simd = false))
+          Dict(:alg=>IRKGL16(simd = true))
           Dict(:alg=>Tsit5(), :prob_choice => 2)
           Dict(:alg=>Vern6(), :prob_choice => 2)
           Dict(:alg=>MER5v2(), :prob_choice => 2)
@@ -63,9 +65,12 @@ setups = [Dict(:alg=>DP5())
           #Dict(:alg=>ode45())
           Dict(:alg=>Tsit5())
           Dict(:alg=>Vern6())
+          Dict(:alg=>IRKGL16(simd = false))
+          Dict(:alg=>IRKGL16(simd = true))
           Dict(:alg=>Tsit5(), :prob_choice => 2)
           Dict(:alg=>Vern6(), :prob_choice => 2)]
-wp = WorkPrecisionSet(probs, abstols, reltols, setups; appxsol = test_sol, maxiters = 10000,
+wp = WorkPrecisionSet(
+    probs, abstols, reltols, setups; appxsol = test_sol, maxiters = 10000,
     error_estimate = :L2, dense_errors = true, numruns = 100)
 plot(wp)
 ```
@@ -82,6 +87,8 @@ setups = [Dict(:alg=>DP8())
           Dict(:alg=>Vern7())
           Dict(:alg=>Vern8())
           Dict(:alg=>Vern9())
+          Dict(:alg=>IRKGL16(simd = false))
+          Dict(:alg=>IRKGL16(simd = true))
           Dict(:alg=>Vern6(), :prob_choice => 2)
           Dict(:alg=>Vern7(), :prob_choice => 2)
           Dict(:alg=>Vern8(), :prob_choice => 2)
@@ -98,6 +105,8 @@ setups = [Dict(:alg=>odex())
           Dict(:alg=>Vern7())
           Dict(:alg=>Vern8())
           Dict(:alg=>Vern9())
+          Dict(:alg=>IRKGL16(simd = false))
+          Dict(:alg=>IRKGL16(simd = true))
           Dict(:alg=>Vern6(), :prob_choice => 2)
           Dict(:alg=>Vern7(), :prob_choice => 2)
           Dict(:alg=>Vern8(), :prob_choice => 2)
@@ -119,6 +128,8 @@ setups = [Dict(:alg=>DP8())
           Dict(:alg=>Vern7())
           Dict(:alg=>Vern8())
           Dict(:alg=>Vern9())
+          Dict(:alg=>IRKGL16(simd = false))
+          Dict(:alg=>IRKGL16(simd = true))
           Dict(:alg=>Vern6(), :prob_choice => 2)
           Dict(:alg=>Vern7(), :prob_choice => 2)
           Dict(:alg=>Vern8(), :prob_choice => 2)
@@ -141,14 +152,17 @@ reltols = 1.0 ./ 10.0 .^ (8:13)
 setups = [Dict(:alg=>Tsit5())
           Dict(:alg=>Vern9())
           Dict(:alg=>VCABM())
+          Dict(:alg=>IRKGL16(simd = false))
+          Dict(:alg=>IRKGL16(simd = true))
           Dict(:alg=>Vern9(), :prob_choice => 2)
           Dict(:alg=>VCABM(), :prob_choice => 2)
           Dict(:alg=>AitkenNeville(min_order = 1, max_order = 9, init_order = 4, threading = true))
           Dict(:alg=>ExtrapolationMidpointDeuflhard(
               min_order = 1, max_order = 9, init_order = 4, threading = true))
           Dict(:alg=>ExtrapolationMidpointHairerWanner(
               min_order = 2, max_order = 11, init_order = 4, threading = true))]
-solnames = ["Tsit5", "Vern9", "VCABM", "Vern9 Static", "VCABM Static",
+solnames = ["Tsit5", "Vern9", "VCABM", "IRKGL16", "IRKGL16 SIMD",
+    "Vern9 Static", "VCABM Static",
     "AitkenNeville", "Midpoint Deuflhard", "Midpoint Hairer Wanner"]
 wp = WorkPrecisionSet(
     probs, abstols, reltols, setups; appxsol = test_sol, names = solnames,