From 47c1c34212dbacd40586a6f842354b179882c212 Mon Sep 17 00:00:00 2001 From: Nathanael Bosch Date: Sat, 4 Nov 2023 16:07:19 +0100 Subject: [PATCH] Change the LotkaVolterra benchmark on the initialization to see what's going on --- docs/src/benchmarks/lotkavolterra.md | 29 ++++++++++++++++------------ 1 file changed, 17 insertions(+), 12 deletions(-) diff --git a/docs/src/benchmarks/lotkavolterra.md b/docs/src/benchmarks/lotkavolterra.md index 9375a816c..ccd13c899 100644 --- a/docs/src/benchmarks/lotkavolterra.md +++ b/docs/src/benchmarks/lotkavolterra.md @@ -205,15 +205,22 @@ DENSE = false; SAVE_EVERYSTEP = false; _setups = [ - "EK1(2) TaylorInit" => Dict(:alg => EK1(order=2, smooth=DENSE, initialization=TaylorModeInit())) - "EK1(3) TaylorInit" => Dict(:alg => EK1(order=3, smooth=DENSE, initialization=TaylorModeInit())) - "EK1(5) TaylorInit" => Dict(:alg => EK1(order=5, smooth=DENSE, initialization=TaylorModeInit())) - "EK1(2) Tsit5Init" => Dict(:alg => EK1(order=2, smooth=DENSE, initialization=ClassicSolverInit())) - "EK1(3) Tsit5Init" => Dict(:alg => EK1(order=3, smooth=DENSE, initialization=ClassicSolverInit())) - "EK1(5) Tsit5Init" => Dict(:alg => EK1(order=5, smooth=DENSE, initialization=ClassicSolverInit())) - "EK1(2) Tsit5Init+ddu" => Dict(:alg => EK1(order=2, smooth=DENSE, initialization=ClassicSolverInit(init_on_ddu=true))) - "EK1(3) Tsit5Init+ddu" => Dict(:alg => EK1(order=3, smooth=DENSE, initialization=ClassicSolverInit(init_on_ddu=true))) - "EK1(5) Tsit5Init+ddu" => Dict(:alg => EK1(order=5, smooth=DENSE, initialization=ClassicSolverInit(init_on_ddu=true))) + "EK1(2) TaylorInit" => Dict(:alg => EK1(order=2, smooth=DENSE, initialization=TaylorModeInit(2))) + "EK1(3) TaylorInit" => Dict(:alg => EK1(order=3, smooth=DENSE, initialization=TaylorModeInit(3))) + "EK1(5) TaylorInit" => Dict(:alg => EK1(order=5, smooth=DENSE, initialization=TaylorModeInit(5))) + "EK1(8) TaylorInit" => Dict(:alg => EK1(order=8, smooth=DENSE, initialization=TaylorModeInit(8))) + "EK1(2) TaylorInit" => Dict(:alg => EK1(order=2, smooth=DENSE, initialization=ForwardDiffInit(2))) + "EK1(3) TaylorInit" => Dict(:alg => EK1(order=3, smooth=DENSE, initialization=ForwardDiffInit(3))) + "EK1(5) TaylorInit" => Dict(:alg => EK1(order=5, smooth=DENSE, initialization=ForwardDiffInit(5))) + "EK1(8) TaylorInit" => Dict(:alg => EK1(order=8, smooth=DENSE, initialization=ForwardDiffInit(8))) + "EK1(2) SimpleInit" => Dict(:alg => EK1(order=2, smooth=DENSE, initialization=SimpleInit())) + "EK1(3) SimpleInit" => Dict(:alg => EK1(order=3, smooth=DENSE, initialization=SimpleInit())) + "EK1(5) SimpleInit" => Dict(:alg => EK1(order=5, smooth=DENSE, initialization=SimpleInit())) + "EK1(8) SimpleInit" => Dict(:alg => EK1(order=8, smooth=DENSE, initialization=SimpleInit())) + "EK1(2) ClassicSolverInit" => Dict(:alg => EK1(order=2, smooth=DENSE, initialization=ClassicSolverInit())) + "EK1(3) ClassicSolverInit" => Dict(:alg => EK1(order=3, smooth=DENSE, initialization=ClassicSolverInit())) + "EK1(5) ClassicSolverInit" => Dict(:alg => EK1(order=5, smooth=DENSE, initialization=ClassicSolverInit())) + "EK1(8) ClassicSolverInit" => Dict(:alg => EK1(order=8, smooth=DENSE, initialization=ClassicSolverInit())) ] labels = first.(_setups) @@ -235,7 +242,7 @@ wp = WorkPrecisionSet( verbose = false, ) -plot(wp, color=[2 2 2 4 4 4 5 5 5], xticks = 10.0 .^ (-16:1:5)) +plot(wp, color=[2 2 2 2 4 4 4 4 5 5 5 5 6 6 6 6], xticks = 10.0 .^ (-16:1:5)) ``` ![](figures/lotkavolterra_7_1.svg) @@ -719,5 +726,3 @@ Info Packages marked with ⌃ and ⌅ have new versions available, but those wi th ⌅ are restricted by compatibility constraints from upgrading. To see why use `status --outdated -m` ``` - -