VdP _again_

nathanaelbosch · Nov 14, 2023 · 0b4536d · 0b4536d
1 parent a127df1
commit 0b4536d
Show file tree

Hide file tree

Showing 15 changed files with 4,384 additions and 4,495 deletions.
diff --git a/docs/src/benchmarks/figures/vanderpol_10_1.svg b/docs/src/benchmarks/figures/vanderpol_10_1.svg
diff --git a/docs/src/benchmarks/figures/vanderpol_11_1.svg b/docs/src/benchmarks/figures/vanderpol_11_1.svg
diff --git a/docs/src/benchmarks/figures/vanderpol_12_1.svg b/docs/src/benchmarks/figures/vanderpol_12_1.svg
diff --git a/docs/src/benchmarks/figures/vanderpol_13_1.svg b/docs/src/benchmarks/figures/vanderpol_13_1.svg
diff --git a/docs/src/benchmarks/figures/vanderpol_14_1.svg b/docs/src/benchmarks/figures/vanderpol_14_1.svg
diff --git a/docs/src/benchmarks/figures/vanderpol_15_1.svg b/docs/src/benchmarks/figures/vanderpol_15_1.svg
diff --git a/docs/src/benchmarks/figures/vanderpol_2_1.svg b/docs/src/benchmarks/figures/vanderpol_2_1.svg
diff --git a/docs/src/benchmarks/figures/vanderpol_3_1.svg b/docs/src/benchmarks/figures/vanderpol_3_1.svg
diff --git a/docs/src/benchmarks/figures/vanderpol_4_1.svg b/docs/src/benchmarks/figures/vanderpol_4_1.svg
diff --git a/docs/src/benchmarks/figures/vanderpol_5_1.svg b/docs/src/benchmarks/figures/vanderpol_5_1.svg
diff --git a/docs/src/benchmarks/figures/vanderpol_6_1.svg b/docs/src/benchmarks/figures/vanderpol_6_1.svg
diff --git a/docs/src/benchmarks/figures/vanderpol_7_1.svg b/docs/src/benchmarks/figures/vanderpol_7_1.svg
diff --git a/docs/src/benchmarks/figures/vanderpol_8_1.svg b/docs/src/benchmarks/figures/vanderpol_8_1.svg
diff --git a/docs/src/benchmarks/figures/vanderpol_9_1.svg b/docs/src/benchmarks/figures/vanderpol_9_1.svg
diff --git a/docs/src/benchmarks/vanderpol.md b/docs/src/benchmarks/vanderpol.md
@@ -38,7 +38,7 @@ u0 = [2.0, 0.0]
 prob = ODEProblem(vanderpol!, u0, tspan, p)
 
 test_sol = solve(prob, RadauIIA5(), abstol=1/10^14, reltol=1/10^14)
-plot(test_sol, title="Van der Pol Solution", legend=false, ylims=(-2.5, 2.5), xticks=:auto)
+plot(test_sol, title="Van der Pol Solution", legend=false, ylims=(-5, 5), xticks=:auto)
 ```
 
 ![](figures/vanderpol_2_1.svg)
@@ -259,7 +259,7 @@ du0 = [0.0]
 prob2 = SecondOrderODEProblem(vanderpol2!, du0, u0, tspan, p)
 
 test_sol2 = solve(prob2, RadauIIA5(), abstol=1/10^14, reltol=1/10^14)
-plot(test_sol2, title="Van der Pol Solution (2nd order)", legend=false, ylims=(-2.5, 2.5), xticks=:auto)
+plot(test_sol2, title="Van der Pol Solution (2nd order)", legend=false, ylims=(-5, 5), xticks=:auto)
 ```
 
 ![](figures/vanderpol_11_1.svg)
@@ -282,8 +282,9 @@ _setups = [
 labels = first.(_setups)
 setups = last.(_setups)
 
-abstols = 1.0 ./ 10.0 .^ (6:11)
-reltols = 1.0 ./ 10.0 .^ (3:8)
+abstols = 1.0 ./ 10.0 .^ (5:8)
+reltols = 1.0 ./ 10.0 .^ (2:5)
+
 
 wp = WorkPrecisionSet(
     [prob, prob2], abstols, reltols, setups;