Change the LotkaVolterra benchmark on the initialization to see what'…

…s going on
nathanaelbosch · Nov 4, 2023 · 47c1c34 · 47c1c34
1 parent 32c4c2c
commit 47c1c34
Showing 1 changed file with 17 additions and 12 deletions.
diff --git 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`
 ```
-
-