add stiff to nelson

Author: ChrisRackauckas

benchmarks/AstroChem/nelson.jmd

@@ -157,17 +157,29 @@ plot(sol; yscale=:log10, idxs = (0,5))
 abstols = 1.0 ./ 10.0 .^ (8:10)
 reltols = 1.0 ./ 10.0 .^ (8:10)
 
-sol = solve(prob, DP8(), abstol=1e-8, reltol=1e-8)
+sol = solve(prob, CVODE_BDF(), abstol=1e-14, reltol=1e-14)
+
+setups = [
+          Dict(:alg=>FBDF()),
+          Dict(:alg=>QNDF()),
+          #Dict(:alg=>Rodas4P()),
+          Dict(:alg=>CVODE_BDF()),
+          #Dict(:alg=>ddebdf()),
+          #Dict(:alg=>Rodas4()),
+          Dict(:alg=>Rodas5P()),
+          Dict(:alg=>KenCarp4()),
+          Dict(:alg=>KenCarp47()),
+          Dict(:alg=>RadauIIA9()),
+		  Dict(:alg=>lsoda()),
+          #Dict(:alg=>rodas()),
+          #Dict(:alg=>radau()),
+          #Dict(:alg=>lsoda()),
+          #Dict(:alg=>ImplicitEulerExtrapolation(min_order = 5, init_order = 3,threading = OrdinaryDiffEqCore.PolyesterThreads())),
+          #Dict(:alg=>ImplicitEulerExtrapolation(min_order = 5, init_order = 3,threading = false)),
+          #Dict(:alg=>ImplicitEulerBarycentricExtrapolation(min_order = 5, threading = OrdinaryDiffEqCore.PolyesterThreads())),
+          #Dict(:alg=>ImplicitEulerBarycentricExtrapolation(min_order = 5, threading = false)),
+          ]
 
-
-setups = [Dict(:alg=>VCABM())
-          Dict(:alg=>Vern7())
-          Dict(:alg=>CVODE_Adams())
-          Dict(:alg=>CVODE_BDF())
-          Dict(:alg=>lsoda())
-          Dict(:alg=>odex())
-          Dict(:alg=>ddeabm())
-]
 wp = WorkPrecisionSet(prob,abstols,reltols,setups;appxsol=refsol,save_everystep=false, print_names = true)
 plot(wp)
 ```