Skip to content

Commit

Permalink
Clean up the lotkavolterra benchmark
Browse files Browse the repository at this point in the history
  • Loading branch information
nathanaelbosch committed Nov 4, 2023
1 parent 00332f5 commit a726d7b
Showing 1 changed file with 1 addition and 47 deletions.
48 changes: 1 addition & 47 deletions benchmarks/lotkavolterra.jmd
Original file line number Diff line number Diff line change
Expand Up @@ -175,53 +175,7 @@ wp = WorkPrecisionSet(
plot(wp, color=[2 2 2 3 3 3], xticks = 10.0 .^ (-16:1:5))
```

## TaylorModeInit vs ClassicSolverInit

```julia
DENSE = false;
SAVE_EVERYSTEP = false;

labels = first.(_setups)
setups = last.(_setups)

abstols = 1.0 ./ 10.0 .^ (4:14)
reltols = 1.0 ./ 10.0 .^ (1:11)

_setups = [
"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) ForwardDiffInit" => Dict(:alg => EK1(order=2, smooth=DENSE, initialization=ForwardDiffInit(2)))
"EK1(3) ForwardDiffInit" => Dict(:alg => EK1(order=3, smooth=DENSE, initialization=ForwardDiffInit(3)))
"EK1(5) ForwardDiffInit" => Dict(:alg => EK1(order=5, smooth=DENSE, initialization=ForwardDiffInit(5)))
"EK1(8) ForwardDiffInit" => 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()))
]

# wp = WorkPrecisionSet(
# prob, abstols, reltols, setups;
# names = labels,
# #print_names = true,
# appxsol = test_sol,
# dense = DENSE,
# save_everystep = SAVE_EVERYSTEP,
# numruns = 10,
# maxiters = Int(1e7),
# timeseries_errors = false,
# verbose = false,
# )

# 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))
```

## Comparison of the different initialization schemes

```julia
DENSE = false;
Expand Down

0 comments on commit a726d7b

Please sign in to comment.