@@ -88,55 +88,32 @@ plt = heatmap(xs, tslices, ys', xlabel="x", ylabel="t")
8888
8989#### Implicit-Explicit Methods
9090
91- Dense/banded linear solvers.
92-
93- ```julia
94- abstols = 0.1 .^ (7:10)
95- reltols = 0.1 .^ (4:7)
96- setups = [
97- # KenCarp methods take forever with adaptive timestepping for some reason
98- Dict(:alg => KenCarp3()),
99- Dict(:alg => KenCarp4()),
100- Dict(:alg => KenCarp5()),
101- # Dict(:alg => ARKODE(Sundials.Implicit(), order=3, linear_solver=:Band, jac_upper=1, jac_lower=1)),
102- # Dict(:alg => ARKODE(Sundials.Implicit(), order=4, linear_solver=:Band, jac_upper=1, jac_lower=1)),
103- # Dict(:alg => ARKODE(Sundials.Implicit(), order=5, linear_solver=:Band, jac_upper=1, jac_lower=1))
104- ]
105- labels = hcat(
106- "KenCarp3",
107- "KenCarp4",
108- "KenCarp5",
109- # "ARKODE3",
110- # "ARKODE4",
111- # "ARKODE5",
112- )
113- @time wp = WorkPrecisionSet(prob, abstols, reltols, setups;
114- print_names=true, names=labels, numruns=5, error_estimate=:l2,
115- save_everystep=false, appxsol=test_sol, maxiters=Int(2e5));
116-
117- plot(wp, label=labels, markershape=:auto, title="IMEX Methods, Band Linsolve, Low Tolerances")
118- ```
119-
12091Krylov linear solvers.
12192```julia
12293abstols = 0.1 .^ (7:10)
12394reltols = 0.1 .^ (4:7)
12495setups = [
12596 # KenCarp methods take forever with adaptive timestepping for some reason
126- Dict(:alg => KenCarp3(linsolve=KrylovJL_GMRES())),
127- Dict(:alg => KenCarp4(linsolve=KrylovJL_GMRES())),
128- Dict(:alg => KenCarp5(linsolve=KrylovJL_GMRES())),
97+ # Dict(:alg => KenCarp3()),
98+ # Dict(:alg => KenCarp4()),
99+ # Dict(:alg => KenCarp5()),
100+ # Dict(:alg => KenCarp3(linsolve=KrylovJL_GMRES())),
101+ # Dict(:alg => KenCarp4(linsolve=KrylovJL_GMRES())),
102+ # Dict(:alg => KenCarp5(linsolve=KrylovJL_GMRES())),
129103 Dict(:alg => ARKODE(Sundials.Implicit(), order=3, linear_solver=:GMRES)),
130104 Dict(:alg => ARKODE(Sundials.Implicit(), order=4, linear_solver=:GMRES)),
131105 Dict(:alg => ARKODE(Sundials.Implicit(), order=5, linear_solver=:GMRES)),
132106]
133107labels = hcat(
134- "KenCarp3",
135- "KenCarp4",
136- "KenCarp5",
137- "ARKODE3",
138- "ARKODE4",
139- "ARKODE5",
108+ # "KenCarp3",
109+ # "KenCarp4",
110+ # "KenCarp5",
111+ # "KenCarp3 (Krylov)",
112+ # "KenCarp4 (Krylov)",
113+ # "KenCarp5 (Krylov)",
114+ "ARKODE3 (Krylov)",
115+ "ARKODE4 (Krylov)",
116+ "ARKODE5 (Krylov)",
140117)
141118@time wp = WorkPrecisionSet(prob, abstols, reltols, setups;
142119 print_names=true, names=labels, numruns=5, error_estimate=:l2,
0 commit comments