From ea45a30e4ee5a0e455521928f4fa9118560f9cfe Mon Sep 17 00:00:00 2001 From: Nathanael Bosch Date: Thu, 12 Oct 2023 16:20:23 +0200 Subject: [PATCH] Fix the ROBER benchmark plot --- benchmarks/rober.jmd | 13 +- docs/src/benchmarks/figures/rober_2_1.svg | 232 +++++----- docs/src/benchmarks/figures/rober_3_1.svg | 529 ++++++++++++++-------- docs/src/benchmarks/rober.md | 9 +- 4 files changed, 471 insertions(+), 312 deletions(-) diff --git a/benchmarks/rober.jmd b/benchmarks/rober.jmd index 0ef839f77..5c5025557 100644 --- a/benchmarks/rober.jmd +++ b/benchmarks/rober.jmd @@ -46,18 +46,14 @@ SAVE_EVERYSTEP = false; _setups = [ "EK1($order)" => Dict(:alg => EK1(order=order, smooth=DENSE)) - for order in 2:6 + for order in 1:4 ] labels = first.(_setups) setups = last.(_setups) -# works: -# abstols = 1.0 ./ 10.0 .^ (4:10) -# reltols = 1.0 ./ 10.0 .^ (1:7) -# test: -abstols = 1.0 ./ 10.0 .^ (4:9) -reltols = 1.0 ./ 10.0 .^ (1:6) +abstols = 1.0 ./ 10.0 .^ (3:9) +reltols = 1.0 ./ 10.0 .^ (1:7) wp = WorkPrecisionSet( mmprob, abstols, reltols, setups; @@ -72,8 +68,7 @@ wp = WorkPrecisionSet( verbose = false, ) -plot(wp, palette=Plots.palette([:blue, :red], length(_setups)), xticks = 10.0 .^ (-16:1:5), - xlims = (2e-15, 3e-7), ylims = (1e-2, 6e-1)) +plot(wp, palette=Plots.palette([:blue, :red], length(_setups)), xticks = 10.0 .^ (-16:1:5)) ``` diff --git a/docs/src/benchmarks/figures/rober_2_1.svg b/docs/src/benchmarks/figures/rober_2_1.svg index 6cf14dc73..19c5a15f4 100644 --- a/docs/src/benchmarks/figures/rober_2_1.svg +++ b/docs/src/benchmarks/figures/rober_2_1.svg @@ -1,131 +1,131 @@ - + - + - + - + - + - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + diff --git a/docs/src/benchmarks/figures/rober_3_1.svg b/docs/src/benchmarks/figures/rober_3_1.svg index ffbff4443..46b2ca1b4 100644 --- a/docs/src/benchmarks/figures/rober_3_1.svg +++ b/docs/src/benchmarks/figures/rober_3_1.svg @@ -1,197 +1,360 @@ - + - + - + - + - + - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + diff --git a/docs/src/benchmarks/rober.md b/docs/src/benchmarks/rober.md index 7fc790941..454735b4a 100644 --- a/docs/src/benchmarks/rober.md +++ b/docs/src/benchmarks/rober.md @@ -53,7 +53,7 @@ SAVE_EVERYSTEP = false; _setups = [ "EK1($order)" => Dict(:alg => EK1(order=order, smooth=DENSE)) - for order in 2:6 + for order in 1:4 ] labels = first.(_setups) @@ -63,8 +63,8 @@ setups = last.(_setups) # abstols = 1.0 ./ 10.0 .^ (4:10) # reltols = 1.0 ./ 10.0 .^ (1:7) # test: -abstols = 1.0 ./ 10.0 .^ (4:9) -reltols = 1.0 ./ 10.0 .^ (1:6) +abstols = 1.0 ./ 10.0 .^ (3:9) +reltols = 1.0 ./ 10.0 .^ (1:7) wp = WorkPrecisionSet( mmprob, abstols, reltols, setups; @@ -80,7 +80,8 @@ wp = WorkPrecisionSet( ) plot(wp, palette=Plots.palette([:blue, :red], length(_setups)), xticks = 10.0 .^ (-16:1:5), - xlims = (2e-15, 3e-7), ylims = (1e-2, 6e-1)) + #xlims = (2e-15, 3e-7), ylims = (1e-2, 6e-1) + ) ``` ![](figures/rober_3_1.svg)