QNDF Benchmarks for ROBER

using Pkg

pkg"update"

pkg"add BenchmarkTools DiffEqBase DiffEqDevTools DiffEqFlux DiffEqOperators DiffEqProblemLibrary DifferentialEquations FillArrays Flux LSODA ODE ODEInterfaceDiffEq ParameterizedFunctions Test Plots RecursiveArrayTools StaticArrays Sundials LinearAlgebra Random"

pkg"add OrdinaryDiffEq"

using OrdinaryDiffEq, DiffEqDevTools, Sundials, ParameterizedFunctions, Plots, ODE, ODEInterfaceDiffEq, LSODA

gr() # gr(fmt=:png)

using LinearAlgebra

LinearAlgebra.BLAS.set_num_threads(1)

rober = @ode_def begin

  dy₁ = -k₁*y₁+k₃*y₂*y₃

  dy₂ =  k₁*y₁-k₂*y₂^2-k₃*y₂*y₃

  dy₃ =  k₂*y₂^2

end k₁ k₂ k₃

prob = ODEProblem(rober,[1.0,0.0,0.0],(0.0,1e5),(0.04,3e7,1e4))

sol = solve(prob,CVODE_BDF(),abstol=1/10^14,reltol=1/10^14)

test_sol = TestSolution(sol)

plot(sol,dpi=200)

abstols = 1.0 ./ 10.0 .^ (5:8)

reltols = 1.0 ./ 10.0 .^ (1:4);

setups = [Dict(:alg=>Rosenbrock23()),

          Dict(:alg=>TRBDF2()),

          Dict(:alg=>ImplicitEulerExtrapolation()),

          #Dict(:alg=>ImplicitDeuflhardExtrapolation()), # Diverges

          #Dict(:alg=>ImplicitHairerWannerExtrapolation()), # Diverges

          #Dict(:alg=>ABDF2()), # Maxiters

          Dict(:alg=>QNDF()),

          Dict(:alg=>Exprb43()),

          Dict(:alg=>Exprb32()),

]

names= ["Rosenbrock23" "TRBDF2" "ImplicitEulerExtrapolation" "QNDF"  "Exprb43" "Exprb32"]

wp = WorkPrecisionSet(prob,abstols,reltols,setups; 

 save_everystep=false,appxsol=test_sol,maxiters=Int(1e5),numruns=10)

plot(wp)

versioninfo()