Skip to content

Commit

Permalink
Fix the ROBER benchmark plot
Browse files Browse the repository at this point in the history
  • Loading branch information
nathanaelbosch committed Oct 12, 2023
1 parent 7e57fba commit ea45a30
Show file tree
Hide file tree
Showing 4 changed files with 471 additions and 312 deletions.
13 changes: 4 additions & 9 deletions benchmarks/rober.jmd
Original file line number Diff line number Diff line change
Expand Up @@ -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;
Expand All @@ -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))
```


Expand Down
232 changes: 116 additions & 116 deletions docs/src/benchmarks/figures/rober_2_1.svg
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
529 changes: 346 additions & 183 deletions docs/src/benchmarks/figures/rober_3_1.svg
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
9 changes: 5 additions & 4 deletions docs/src/benchmarks/rober.md
Original file line number Diff line number Diff line change
Expand Up @@ -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)
Expand All @@ -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;
Expand All @@ -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)
Expand Down

0 comments on commit ea45a30

Please sign in to comment.