Merge pull request #91 from SciML/prob_choice

fix prob_choice filtering and add a better test

Author: ChrisRackauckas

Project.toml

@@ -17,7 +17,7 @@ Statistics = "10745b16-79ce-11e8-11f9-7d13ad32a3b2"
 
 [compat]
 DelayDiffEq = "5.20"
-DiffEqBase = "6.40"
+DiffEqBase = "6.94.4"
 DiffEqNoiseProcess = "5.0"
 DiffEqProblemLibrary = "4.5"
 NLsolve = "4.2"

test/benchmark_tests.jl

@@ -115,6 +115,39 @@ wp = WorkPrecisionSet(probs, abstols, reltols, setups; appxsol = [test_sol, noth
 wp = WorkPrecisionSet(probs, abstols, reltols, setups; appxsol = [test_sol, test_sol1],
                       save_everystep = false, numruns = 20, maxiters = 10000)
 
+# Dual DAE Problems
+function rober(du, u, p, t)
+    y₁, y₂, y₃ = u
+    k₁, k₂, k₃ = p
+    du[1] = -k₁ * y₁ + k₃ * y₂ * y₃
+    du[2] = k₁ * y₁ - k₂ * y₂^2 - k₃ * y₂ * y₃
+    du[3] = k₂ * y₂^2
+    nothing
+end
+prob1 = ODEProblem(rober, [1.0, 0.0, 0.0], (0.0, 1e5), [0.04, 3e7, 1e4])
+
+ode_ref_sol = solve(prob1, Rodas5(), abstol = 1 / 10^14, reltol = 1 / 10^14);
+
+function dae_rober(out, du, u, p, t)
+    out[1] = -0.04u[1] + 1e4 * u[2] * u[3] - du[1]
+    out[2] = +0.04u[1] - 3e7 * u[2]^2 - 1e4 * u[2] * u[3] - du[2]
+    out[3] = u[1] + u[2] + u[3] - 1.0
+end
+u₀ = [1.0, 0, 0]
+du₀ = [-0.04, 0.04, 0.0]
+tspan = (0.0, 100000.0)
+
+differential_vars = [true, true, false]
+prob2 = DAEProblem(dae_rober, du₀, u₀, tspan, differential_vars = differential_vars)
+
+dae_ref_sol = solve(prob2, DFBDF(), abstol = 1 / 10^14, reltol = 1 / 10^14);
+
+probs = [prob1, prob2]
+setups = [Dict(:alg => Rodas5())
+          Dict(:alg => DFBDF(), :prob_choice => 2)]
+wp = WorkPrecisionSet(probs, abstols, reltols, setups; appxsol = [ode_ref_sol, dae_ref_sol],
+                      save_everystep = false, numruns = 20, maxiters = 10000)
+
 # DDE problem
 prob = prob_dde_constant_1delay_ip