Skip to content

Commit 0f296e1

Browse files
ChrisRackauckasChrisRackauckas
committed
remove broken and unnecessary plot limit calculation
1 parent 52a3401 commit 0f296e1

File tree

5 files changed

+24
-113
lines changed

5 files changed

+24
-113
lines changed

‎benchmarks/Bio/BCR.jmd‎

Lines changed: 5 additions & 23 deletions
Original file line numberDiff line numberDiff line change
@@ -99,18 +99,6 @@ test_sol = TestSolution(sol);
9999
default(legendfontsize=7,framestyle=:box,gridalpha=0.3,gridlinewidth=2.5)
100100
```
101101

102-
#### Declares a plotting helper function
103-
104-
```julia
105-
function plot_settings(wp)
106-
times = vcat(map(wp -> wp.times, wp.wps)...)
107-
errors = vcat(map(wp -> wp.errors, wp.wps)...)
108-
xlimit = 10 .^ (floor(log10(minimum(errors))), ceil(log10(maximum(errors))))
109-
ylimit = 10 .^ (floor(log10(minimum(times))), ceil(log10(maximum(times))))
110-
return xlimit,ylimit
111-
end
112-
```
113-
114102
#### Declare pre-conditioners
115103
```julia
116104
using IncompleteLU, LinearAlgebra
@@ -189,8 +177,7 @@ wp = WorkPrecisionSet([oprob,oprob_sparse,sparsejacprob],abstols,reltols,setups;
189177
saveat=tf/10000.,appxsol=[test_sol,test_sol,test_sol],maxiters=Int(1e6),numruns=1)
190178

191179
names = ["lsoda" "CVODE_BDF" "CVODE_BDF (LapackDense)" "CVODE_BDF (GMRES)" "CVODE_BDF (GMRES, iLU)" "CVODE_BDF (KLU, sparse jac)"]
192-
xlimit,ylimit = plot_settings(wp)
193-
plot(wp;label=names,xlimit=xlimit,ylimit=ylimit)
180+
plot(wp;label=names)
194181
```
195182

196183
## Work-Precision Diagrams (various Julia solvers)
@@ -213,8 +200,7 @@ wp = WorkPrecisionSet(oprob,abstols,reltols,setups;error_estimate=:l2,
213200
saveat=tf/10000.,appxsol=test_sol,maxiters=Int(1e6),numruns=1)
214201

215202
names = ["TRBDF2" "QNDF" "FBDF" "KenCarp4"]
216-
xlimit,ylimit = plot_settings(wp)
217-
plot(wp;label=names,xlimit=xlimit,ylimit=ylimit)
203+
plot(wp;label=names)
218204
```
219205

220206
#### Declare solvers (using GMRES linear solver).
@@ -235,8 +221,7 @@ wp = WorkPrecisionSet(oprob,abstols,reltols,setups;error_estimate=:l2,
235221
saveat=tf/10000.,appxsol=test_sol,maxiters=Int(1e6),numruns=1)
236222

237223
names = ["TRBDF2 (GMRES)" "QNDF (GMRES)" "FBDF (GMRES)" "KenCarp4 (GMRES)"]
238-
xlimit,ylimit = plot_settings(wp)
239-
plot(wp;label=names,xlimit=xlimit,ylimit=ylimit)
224+
plot(wp;label=names)
240225
```
241226

242227
#### Declare solvers (using GMRES linear solver, with pre-conditioner).
@@ -257,8 +242,7 @@ wp = WorkPrecisionSet(sparsejacprob,abstols,reltols,setups;error_estimate=:l2,
257242
saveat=tf/10000.,appxsol=test_sol,maxiters=Int(1e6),numruns=1)
258243

259244
names = ["TRBDF2 (GMRES, iLU)" "QNDF (GMRES, iLU)" "FBDF (GMRES, iLU)" "KenCarp4 (GMRES, iLU)"]
260-
xlimit,ylimit = plot_settings(wp)
261-
plot(wp;label=names,xlimit=xlimit,ylimit=ylimit)
245+
plot(wp;label=names)
262246
```
263247

264248
#### Declare solvers (using sparse jacobian)
@@ -282,8 +266,7 @@ wp = WorkPrecisionSet(sparsejacprob,abstols,reltols,setups;error_estimate=:l2,
282266
saveat=tf/10000.,appxsol=test_sol,maxiters=Int(1e6),numruns=1)
283267

284268
names = ["TRBDF2 (KLU, sparse jac)" "QNDF (KLU, sparse jac)" "FBDF (KLU, sparse jac)" "KenCarp4 (KLU, sparse jac)"]
285-
xlimit,ylimit = plot_settings(wp)
286-
plot(wp;label=names,xlimit=xlimit,ylimit=ylimit)
269+
plot(wp;label=names)
287270
```
288271

289272
## Summary of results
@@ -312,7 +295,6 @@ wp = WorkPrecisionSet([oprob,oprob_sparse,sparsejacprob],abstols,reltols,setups;
312295
names = ["CVODE_BDF (GMRES, iLU)" "QNDF (GMRES, iLU)" "FBDF (GMRES, iLU)" "QNDF (KLU, sparse jac)" "FBDF (KLU, sparse jac)" "KenCarp4 (KLU, sparse jac)"]
313296
colors = [:green :deepskyblue1 :dodgerblue2 :royalblue2 :slateblue3 :lightskyblue]
314297
markershapes = [:octagon :hexagon :rtriangle :pentagon :ltriangle :star5]
315-
xlimit,ylimit = plot_settings(wp)
316298
xlimit = xlimit .* [0.95,1/0.95]; ylimit = ylimit .* [0.95,1/0.95];
317299
plot(wp;label=names,left_margin=10Plots.mm,right_margin=10Plots.mm,xlimit=xlimit,ylimit=ylimit,xticks=[1e-3,1e-2,1e-1,1e0,1e1,1e2,1e3],yticks=[1e0,1e1,1e2,1e3],color=colors,markershape=markershapes,legendfontsize=15,tickfontsize=15,guidefontsize=15, legend=:topright, lw=20, la=0.8, markersize=20,markerstrokealpha=1.0, markerstrokewidth=1.5, gridalpha=0.3, gridlinewidth=7.5,size=(1100,1000))
318300
```

‎benchmarks/Bio/egfr_net.jmd‎

Lines changed: 5 additions & 23 deletions
Original file line numberDiff line numberDiff line change
@@ -93,18 +93,6 @@ test_sol = TestSolution(sol);
9393
default(legendfontsize=7,framestyle=:box,gridalpha=0.3,gridlinewidth=2.5)
9494
```
9595

96-
#### Declares a plotting helper function
97-
98-
```julia
99-
function plot_settings(wp)
100-
times = vcat(map(wp -> wp.times, wp.wps)...)
101-
errors = vcat(map(wp -> wp.errors, wp.wps)...)
102-
xlimit = 10 .^ (floor(log10(minimum(errors))), ceil(log10(maximum(errors))))
103-
ylimit = 10 .^ (floor(log10(minimum(times))), ceil(log10(maximum(times))))
104-
return xlimit,ylimit
105-
end
106-
```
107-
10896
#### Sets tolerances
10997

11098
```julia
@@ -144,8 +132,7 @@ wp = WorkPrecisionSet(oprob,abstols,reltols,setups;error_estimate=:l2,
144132
saveat=tf/10000.,appxsol=test_sol,maxiters=Int(1e9),numruns=100)
145133

146134
names = ["lsoda" "CVODE_BDF" "CVODE_BDF (Lapack Dense)" "CVODE_Adams" "TRBDF2" "QNDF" "FBDF" "KenCarp4" "Rosenbrock23" "Rodas4" "Rodas5P"]
147-
xlimit,ylimit = plot_settings(wp)
148-
plot(wp;label=names,xlimit=xlimit,ylimit=ylimit)
135+
plot(wp;label=names)
149136
```
150137

151138
#### Declare solvers (using GMRES linear solver)
@@ -174,8 +161,7 @@ wp = WorkPrecisionSet(oprob,abstols,reltols,setups;error_estimate=:l2,
174161
saveat=tf/10000.,appxsol=test_sol,maxiters=Int(1e9),numruns=100)
175162

176163
names = ["lsoda" "CVODE_BDF (GMRES)" "TRBDF2 (GMRES)" "QNDF (GMRES)" "FBDF (GMRES)" "KenCarp4 (GMRES)" "Rosenbrock23 (GMRES)" "Rodas4 (GMRES)" "Rodas5P (GMRES)"]
177-
xlimit,ylimit = plot_settings(wp)
178-
plot(wp;label=names,xlimit=xlimit,ylimit=ylimit)
164+
plot(wp;label=names)
179165
```
180166

181167
#### Declare solvers (using sparse jacobian)
@@ -203,8 +189,7 @@ wp = WorkPrecisionSet(sparsejacprob,abstols,reltols,setups;error_estimate=:l2,
203189
saveat=tf/10000.,appxsol=test_sol,maxiters=Int(1e9),numruns=100)
204190

205191
names = ["CVODE_BDF (KLU, sparse jac)" "TRBDF2 (KLU, sparse jac)" "QNDF (KLU, sparse jac)" "FBDF (KLU, sparse jac)" "KenCarp4 (KLU, sparse jac)" "Rosenbrock23 (KLU, sparse jac)" "Rodas4 (KLU, sparse jac)" "Rodas5P (KLU, sparse jac)"]
206-
xlimit,ylimit = plot_settings(wp)
207-
plot(wp;label=names,xlimit=xlimit,ylimit=ylimit)
192+
plot(wp;label=names)
208193
```
209194

210195
## Explicit Work-Precision Diagram
@@ -236,8 +221,7 @@ wp = WorkPrecisionSet(oprob,abstols,reltols,setups;error_estimate=:l2,
236221
saveat=tf/10000.,appxsol=test_sol,maxiters=Int(1e9),numruns=200)
237222

238223
names = ["lsoda" "CVODE_Adams" "Tsit5" "BS5" "VCABM" "Vern6" "Vern7" "Vern8" "Vern9" "ROCK4"]
239-
xlimit,ylimit = plot_settings(wp)
240-
plot(wp;label=names,xlimit=xlimit,ylimit=ylimit)
224+
plot(wp;label=names)
241225
```
242226

243227
#### Additional explicit solvers
@@ -247,8 +231,7 @@ setups = [Dict(:alg=>ROCK2())];
247231
wp = WorkPrecisionSet(oprob,abstols,reltols,setups;error_estimate=:l2,
248232
saveat=tf/10000.,appxsol=test_sol,maxiters=Int(1e9),numruns=200)
249233
names = ["ROCK2"]
250-
xlimit,ylimit = plot_settings(wp)
251-
plot(wp;label=names,xlimit=xlimit,ylimit=ylimit)
234+
plot(wp;label=names)
252235
```
253236

254237
## Summary of results
@@ -278,7 +261,6 @@ wp = WorkPrecisionSet([oprob,sparsejacprob],abstols,reltols,setups;error_estimat
278261
names = ["lsoda" "CVODE_BDF (GMRES)" "QNDF (GMRES)" "QNDF (KLU)" "BS5" "Vern6" "ROCK4"]
279262
colors = [:seagreen1 :darkgreen :deepskyblue1 :deepskyblue4 :thistle2 :lightslateblue :purple4]
280263
markershapes = [:star4 :rect :hexagon :octagon :star8 :rtriangle :square]
281-
xlimit,ylimit = plot_settings(wp)
282264
xlimit = xlimit .* [0.95,1/0.95]; ylimit = ylimit .* [0.95,1/0.95];
283265
plot(wp;label=names,left_margin=10Plots.mm,right_margin=10Plots.mm,xlimit=xlimit,ylimit=ylimit,xticks=[1e-9,1e-8,1e-7,1e-6,1e-5,1e-4,1e-3,1e-2],yticks=[1e-2,1e-1],color=colors,markershape=markershapes,legendfontsize=15,tickfontsize=15,guidefontsize=15, legend=:topright, lw=20, la=0.8, markersize=20,markerstrokealpha=1.0, markerstrokewidth=1.5, gridalpha=0.3, gridlinewidth=7.5,size=(1100,1000))
284266
```

‎benchmarks/Bio/fceri_gamma2.jmd‎

Lines changed: 6 additions & 25 deletions
Original file line numberDiff line numberDiff line change
@@ -89,18 +89,6 @@ test_sol = TestSolution(sol);
8989
default(legendfontsize=7,framestyle=:box,gridalpha=0.3,gridlinewidth=2.5)
9090
```
9191

92-
#### Declares a plotting helper function
93-
94-
```julia
95-
function plot_settings(wp)
96-
times = vcat(map(wp -> wp.times, wp.wps)...)
97-
errors = vcat(map(wp -> wp.errors, wp.wps)...)
98-
xlimit = 10 .^ (floor(log10(minimum(errors))), ceil(log10(maximum(errors))))
99-
ylimit = 10 .^ (floor(log10(minimum(times))), ceil(log10(maximum(times))))
100-
return xlimit,ylimit
101-
end
102-
```
103-
10492
#### Declare pre-conditioners
10593
```julia
10694
using IncompleteLU, LinearAlgebra
@@ -169,8 +157,7 @@ wp = WorkPrecisionSet([oprob,oprob_sparse,sparsejacprob],abstols,reltols,setups;
169157
saveat=tf/10000.,appxsol=[test_sol,test_sol,test_sol],maxiters=Int(1e9),numruns=10)
170158

171159
names = ["lsoda" "CVODE_BDF" "CVODE_BDF (LapackDense)" "CVODE_BDF (GMRES)" "CVODE_BDF (GMRES, iLU)" "CVODE_BDF (KLU, sparse jac)"]
172-
xlimit,ylimit = plot_settings(wp)
173-
plot(wp;label=names,xlimit=xlimit,ylimit=ylimit)
160+
plot(wp;label=names)
174161
```
175162

176163
## Work-Precision Diagrams (various Julia solvers)
@@ -193,8 +180,7 @@ wp = WorkPrecisionSet(oprob,abstols,reltols,setups;error_estimate=:l2,
193180
saveat=tf/10000.,appxsol=test_sol,maxiters=Int(1e12),dtmin=1e-18,numruns=10)
194181

195182
names = ["TRBDF2" "QNDF" "FBDF" "KenCarp4"]
196-
xlimit,ylimit = plot_settings(wp)
197-
plot(wp;label=names,xlimit=xlimit,ylimit=ylimit)
183+
plot(wp;label=names)
198184
```
199185

200186
#### Declare solvers (using GMRES linear solver).
@@ -215,8 +201,7 @@ wp = WorkPrecisionSet(oprob,abstols,reltols,setups;error_estimate=:l2,
215201
saveat=tf/10000.,appxsol=test_sol,maxiters=Int(1e12),dtmin=1e-18,numruns=10)
216202

217203
names = ["TRBDF2 (GMRES)" "QNDF (GMRES)" "FBDF (GMRES)" "KenCarp4 (GMRES)"]
218-
xlimit,ylimit = plot_settings(wp)
219-
plot(wp;label=names,xlimit=xlimit,ylimit=ylimit)
204+
plot(wp;label=names)
220205
```
221206

222207
#### Declare solvers (using GMRES linear solver, with pre-conditioner).
@@ -237,8 +222,7 @@ wp = WorkPrecisionSet(sparsejacprob,abstols,reltols,setups;error_estimate=:l2,
237222
saveat=tf/10000.,appxsol=test_sol,maxiters=Int(1e12),dtmin=1e-18,numruns=10)
238223

239224
names = ["TRBDF2 (GMRES, iLU)" "QNDF (GMRES, iLU)" "FBDF (GMRES, iLU)" "KenCarp4 (GMRES, iLU)"]
240-
xlimit,ylimit = plot_settings(wp)
241-
plot(wp;label=names,xlimit=xlimit,ylimit=ylimit)
225+
plot(wp;label=names)
242226
```
243227

244228
#### Declare solvers (using sparse jacobian)
@@ -262,8 +246,7 @@ wp = WorkPrecisionSet(sparsejacprob,abstols,reltols,setups;error_estimate=:l2,
262246
saveat=tf/10000.,appxsol=test_sol,maxiters=Int(1e12),dtmin=1e-18,numruns=10)
263247

264248
names = ["TRBDF2 (KLU, sparse jac)" "QNDF (KLU, sparse jac)" "FBDF (KLU, sparse jac)" "KenCarp4 (KLU, sparse jac)"]
265-
xlimit,ylimit = plot_settings(wp)
266-
plot(wp;label=names,xlimit=xlimit,ylimit=ylimit)
249+
plot(wp;label=names)
267250
```
268251

269252
## Explicit Work-Precision Diagram
@@ -293,8 +276,7 @@ wp = WorkPrecisionSet(oprob,abstols,reltols,setups;error_estimate=:l2,
293276
saveat=tf/10000.,appxsol=test_sol,maxiters=Int(1e9),numruns=10)
294277

295278
names = ["CVODE_Adams" "Tsit5" "BS5" "Vern6" "Vern7" "Vern8" "Vern9"]
296-
xlimit,ylimit = plot_settings(wp)
297-
plot(wp;label=names,xlimit=xlimit,ylimit=ylimit)
279+
plot(wp;label=names)
298280
```
299281

300282
## Summary of results
@@ -322,7 +304,6 @@ wp = WorkPrecisionSet([oprob,oprob_sparse,sparsejacprob],abstols,reltols,setups;
322304
names = ["CVODE_BDF (GMRES)" "CVODE_BDF (GMRES, iLU)" "QNDF (GMRES, iLU)" "FBDF (GMRES, iLU)" "Tsit5"]
323305
colors = [:darkgreen :green :deepskyblue1 :dodgerblue2 :orchid2]
324306
markershapes = [:rect :octagon :hexagon :rtriangle :ltriangle]
325-
xlimit,ylimit = plot_settings(wp)
326307
xlimit = xlimit .* [0.95,1/0.95]; ylimit = ylimit .* [0.95,1/0.95];
327308
plot(wp;label=names,left_margin=10Plots.mm,right_margin=10Plots.mm,xlimit=xlimit,ylimit=ylimit,xticks=[1e-9,1e-8,1e-7,1e-6,1e-5,1e-4,1e-3,1e-2,1e-1],yticks=[1e-1,1e0,1e1,1e2],color=colors,markershape=markershapes,legendfontsize=15,tickfontsize=15,guidefontsize=15, legend=:topright, lw=20, la=0.8, markersize=20,markerstrokealpha=1.0, markerstrokewidth=1.5, gridalpha=0.3, gridlinewidth=7.5,size=(1100,1000))
328309
```

‎benchmarks/Bio/multisite2.jmd‎

Lines changed: 4 additions & 21 deletions
Original file line numberDiff line numberDiff line change
@@ -92,18 +92,6 @@ test_sol = TestSolution(sol);
9292
default(legendfontsize=7,framestyle=:box,gridalpha=0.3,gridlinewidth=2.5)
9393
```
9494

95-
#### Declares a plotting helper function
96-
97-
```julia
98-
function plot_settings(wp)
99-
times = vcat(map(wp -> wp.times, wp.wps)...)
100-
errors = vcat(map(wp -> wp.errors, wp.wps)...)
101-
xlimit = 10 .^ (floor(log10(minimum(errors))), ceil(log10(maximum(errors))))
102-
ylimit = 10 .^ (floor(log10(minimum(times))), ceil(log10(maximum(times))))
103-
return xlimit,ylimit
104-
end
105-
```
106-
10795
#### Sets tolerances
10896

10997
```julia
@@ -134,8 +122,7 @@ wp = WorkPrecisionSet(oprob,abstols,reltols,setups;error_estimate=:l2,
134122
saveat=tf/10000.,appxsol=test_sol,maxiters=Int(1e9),numruns=200)
135123

136124
names = ["lsoda" "CVODE_BDF" "CVODE_BDF (LapackDense)" "CVODE_BDF (GMRES)"]
137-
xlimit,ylimit = plot_settings(wp)
138-
plot(wp;label=names,xlimit=xlimit,ylimit=ylimit)
125+
plot(wp;label=names)
139126
```
140127

141128
## Implicit Work-Precision Diagram
@@ -166,8 +153,7 @@ wp = WorkPrecisionSet(oprob,abstols,reltols,setups;error_estimate=:l2,
166153
saveat=tf/10000.,appxsol=test_sol,maxiters=Int(1e12),dtmin=1e-18,numruns=200)
167154

168155
names = ["TRBDF2" "QNDF" "FBDF" "KenCarp4" "Rosenbrock23" "Rodas4" "Rodas5P"]
169-
xlimit,ylimit = plot_settings(wp)
170-
plot(wp;label=names,xlimit=xlimit,ylimit=ylimit)
156+
plot(wp;label=names)
171157
```
172158

173159
Implicit methods doing poorly suggests it's non-stiff.
@@ -201,8 +187,7 @@ wp = WorkPrecisionSet(oprob,abstols,reltols,setups;error_estimate=:l2,
201187
saveat=tf/10000.,appxsol=test_sol,maxiters=Int(1e9),numruns=200)
202188

203189
names = ["lsoda" "CVODE_Adams" "Tsit5" "BS5" "VCABM" "Vern6" "Vern7" "Vern8" "Vern9" "ROCK4"]
204-
xlimit,ylimit = plot_settings(wp)
205-
plot(wp;label=names,xlimit=xlimit,ylimit=ylimit)
190+
plot(wp;label=names)
206191
```
207192

208193
#### Additional explicit solvers
@@ -212,8 +197,7 @@ setups = [Dict(:alg=>ROCK2())];
212197
wp = WorkPrecisionSet(oprob,abstols,reltols,setups;error_estimate=:l2,
213198
saveat=tf/10000.,appxsol=test_sol,maxiters=Int(1e9),numruns=200)
214199
names = ["ROCK2"]
215-
xlimit,ylimit = plot_settings(wp)
216-
plot(wp;label=names,xlimit=xlimit,ylimit=ylimit)
200+
plot(wp;label=names)
217201
```
218202

219203
## Summary of results
@@ -245,7 +229,6 @@ wp = WorkPrecisionSet(oprob,abstols,reltols,setups;error_estimate=:l2,
245229
names = ["lsoda" "CVODE_BDF (GMRES)" "QNDF" "FBDF" "Rodas5P" "BS5" "VCABM" "Vern6" "ROCK4"]
246230
colors = [:seagreen1 :darkgreen :deepskyblue1 :dodgerblue2 :blue :thistle2 :lightsteelblue2 :lightslateblue :purple4]
247231
markershapes = [:star4 :rect :hexagon :rtriangle :heptagon :star8 :heptagon :rtriangle :square]
248-
xlimit,ylimit = plot_settings(wp)
249232
xlimit = xlimit .* [0.95,1/0.95]; ylimit = ylimit .* [0.95,1/0.95];
250233
plot(wp;label=names,left_margin=10Plots.mm,right_margin=10Plots.mm,xlimit=xlimit,ylimit=ylimit,xticks=[1e-10,1e-9,1e-8,1e-7,1e-6,1e-5,1e-4,1e-3],yticks=[1e-3,1e-2,1e-1],color=colors,markershape=markershapes,legendfontsize=15,tickfontsize=15,guidefontsize=15, legend=:topright, lw=20, la=0.8, markersize=20,markerstrokealpha=1.0, markerstrokewidth=1.5, gridalpha=0.3, gridlinewidth=7.5,size=(1100,1000))
251234
```

‎benchmarks/Bio/multistate.jmd‎

Lines changed: 4 additions & 21 deletions
Original file line numberDiff line numberDiff line change
@@ -91,18 +91,6 @@ test_sol = TestSolution(sol);
9191
default(legendfontsize=7,framestyle=:box,gridalpha=0.3,gridlinewidth=2.5)
9292
```
9393

94-
#### Declares a plotting helper function
95-
96-
```julia
97-
function plot_settings(wp)
98-
times = vcat(map(wp -> wp.times, wp.wps)...)
99-
errors = vcat(map(wp -> wp.errors, wp.wps)...)
100-
xlimit = 10 .^ (floor(log10(minimum(errors))), ceil(log10(maximum(errors))))
101-
ylimit = 10 .^ (floor(log10(minimum(times))), ceil(log10(maximum(times))))
102-
return xlimit,ylimit
103-
end
104-
```
105-
10694
#### Sets tolerances
10795

10896
```julia
@@ -133,8 +121,7 @@ wp = WorkPrecisionSet(oprob,abstols,reltols,setups;error_estimate=:l2,
133121
saveat=tf/10000.,appxsol=test_sol,maxiters=Int(1e9),numruns=200)
134122

135123
names = ["lsoda" "CVODE_BDF" "CVODE_BDF (LapackDense)" "CVODE_BDF (GMRES)"]
136-
xlimit,ylimit = plot_settings(wp)
137-
plot(wp;label=names,xlimit=xlimit,ylimit=ylimit)
124+
plot(wp;label=names)
138125
```
139126

140127
## Implicit Work-Precision Diagram
@@ -165,8 +152,7 @@ wp = WorkPrecisionSet(oprob,abstols,reltols,setups;error_estimate=:l2,
165152
saveat=tf/10000.,appxsol=test_sol,maxiters=Int(1e12),dtmin=1e-18,numruns=200)
166153

167154
names = ["TRBDF2" "QNDF" "FBDF" "KenCarp4" "Rosenbrock23" "Rodas4" "Rodas5P"]
168-
xlimit,ylimit = plot_settings(wp)
169-
plot(wp;label=names,xlimit=xlimit,ylimit=ylimit)
155+
plot(wp;label=names)
170156
```
171157

172158
Implicit methods doing poorly suggests it's non-stiff.
@@ -199,8 +185,7 @@ wp = WorkPrecisionSet(oprob,abstols,reltols,setups;error_estimate=:l2,
199185
saveat=tf/10000.,appxsol=test_sol,maxiters=Int(1e9),numruns=200)
200186

201187
names = ["lsoda" "Tsit5" "BS5" "VCABM" "Vern6" "Vern7" "Vern8" "Vern9" "ROCK4"]
202-
xlimit,ylimit = plot_settings(wp)
203-
plot(wp;label=names,xlimit=xlimit,ylimit=ylimit)
188+
plot(wp;label=names)
204189
```
205190

206191
#### Additional explicit solvers
@@ -210,8 +195,7 @@ setups = [Dict(:alg=>CVODE_Adams()), Dict(:alg=>ROCK2())];
210195
wp = WorkPrecisionSet(oprob,abstols,reltols,setups;error_estimate=:l2,
211196
saveat=tf/10000.,appxsol=test_sol,maxiters=Int(1e9),numruns=200)
212197
names = ["CVODE_Adams" "ROCK2"]
213-
xlimit,ylimit = plot_settings(wp)
214-
plot(wp;label=names,xlimit=xlimit,ylimit=ylimit)
198+
plot(wp;label=names)
215199
```
216200

217201
## Summary of results
@@ -243,7 +227,6 @@ wp = WorkPrecisionSet(oprob,abstols,reltols,setups;error_estimate=:l2,
243227
names = ["lsoda" "CVODE_BDF" "QNDF" "KenCarp4" "Rodas5P" "Tsit5" "BS5" "VCABM" "Vern7"]
244228
colors = [:seagreen1 :chartreuse1 :deepskyblue1 :lightskyblue :blue :orchid2 :thistle2 :lightsteelblue2 :mediumpurple1]
245229
markershapes = [:star4 :circle :hexagon :star5 :heptagon :ltriangle :star8 :heptagon :star6]
246-
xlimit,ylimit = plot_settings(wp)
247230
xlimit = xlimit .* [0.95,1/0.95]; ylimit = ylimit .* [0.95,1/0.95];
248231
plot(wp;label=names,left_margin=10Plots.mm,right_margin=10Plots.mm,xlimit=xlimit,ylimit=ylimit,xticks=[1e-12,1e-11,1e-10,1e-9,1e-8,1e-7,1e-6,1e-5,1e-4,1e-3,1e-2],yticks=[1e-3,1e-2],color=colors,markershape=markershapes,legendfontsize=15,tickfontsize=15,guidefontsize=15, legend=:topright, lw=20, la=0.8, markersize=20,markerstrokealpha=1.0, markerstrokewidth=1.5, gridalpha=0.3, gridlinewidth=7.5,size=(1100,1000))
249232
```

0 commit comments

Comments
 (0)