update DDEProblems to have no parameters

Author: ChrisRackauckas

src/dde_premade_problems.jl

@@ -47,7 +47,7 @@ function f_1delay(::Type{Val{:analytic}}, u₀, p, t)
 end
 
 build_prob_dde_1delay(u₀, ::T=u₀) where {T} =
-    DDEProblem(f_1delay, t->[zero(u₀)], [u₀], (zero(T), T(10)), [oneunit(T)])
+    DDEProblem(f_1delay, t->[zero(u₀)], [u₀], (zero(T), T(10)), nothing, [oneunit(T)])
 
 """
     prob_dde_1delay
@@ -81,12 +81,12 @@ function f_1delay_notinplace(u, h, p, t)
     - h(t - oneunit(t)) ./ oneunit(t)
 end
 
-f_1delay_notinplace(::Type{Val{:analytic}}, t, u0) = f_1delay(Val{:analytic}, t, u0)
+f_1delay_notinplace(::Type{Val{:analytic}}, u0, p, t) = f_1delay(Val{:analytic}, u0, p, t)
 
 #### Vectorized history function
 
 build_prob_dde_1delay_notinplace(u₀, ::T=u₀) where {T} =
-    DDEProblem(f_1delay_notinplace, t->[zero(u₀)], [u₀], (zero(T), T(10)), [oneunit(T)])
+    DDEProblem(f_1delay_notinplace, t->[zero(u₀)], [u₀], (zero(T), T(10)), nothing, [oneunit(T)])
 
 """
     prob_dde_1delay_notinplace
@@ -99,7 +99,7 @@ prob_dde_1delay_notinplace = build_prob_dde_1delay_notinplace(1.0)
 #### Scalar history function
 
 build_prob_dde_1delay_scalar_notinplace(u₀, ::T=u₀) where {T} =
-    DDEProblem(f_1delay_notinplace, t -> zero(u₀), u₀, (zero(T), T(10)), [oneunit(T)])
+    DDEProblem(f_1delay_notinplace, t -> zero(u₀), u₀, (zero(T), T(10)), nothing, [oneunit(T)])
 
 """
     prob_dde_1delay_scalar_notinplace
@@ -113,7 +113,7 @@ prob_dde_1delay_scalar_notinplace = build_prob_dde_1delay_scalar_notinplace(1.0)
 
 ### In-place function
 
-function f_2delays(du::T, u, h, p, t) where T
+function f_2delays(du, u, h, p, t::T) where T
     du[1] = (- h(t - T(1//3))[1] - h(t - T(1//5))[1]) / oneunit(t)
 end
 
@@ -154,7 +154,7 @@ function f_2delays(::Type{Val{:analytic}}, u₀, p, t)
 end
 
 build_prob_dde_2delays(u₀, ::T=u₀) where {T} =
-    DDEProblem(f_2delays, t -> [zero(u₀)], [u₀], (zero(T), oneunit(T)), [T(1//3), T(1//5)])
+    DDEProblem(f_2delays, t -> [zero(u₀)], [u₀], (zero(T), oneunit(T)), nothing, [T(1//3), T(1//5)])
 
 """
     prob_dde_2delays
@@ -184,18 +184,18 @@ prob_dde_2delays = build_prob_dde_2delays(1.0)
 
 ### Not in-place function
 
-function f_2delays_notinplace(u::T, h, p, t) where T
+function f_2delays_notinplace(u, h, p, t::T) where T
     (- h(t - T(1//3)) .- h(t - T(1//5))) ./ oneunit(t)
 end
 
-f_2delays_notinplace(::Type{Val{:analytic}}, t, u0) =
-    f_2delays(Val{:analytic}, t, u0)
+f_2delays_notinplace(::Type{Val{:analytic}}, u0, p, t) =
+    f_2delays(Val{:analytic}, u0, p, t)
 
 #### Vectorized history function
 
 build_prob_dde_2delays_notinplace(u₀, ::T=u₀) where {T} =
     DDEProblem(f_2delays_notinplace, t -> [zero(u₀)], [u₀], (zero(T), oneunit(T)),
-               [T(1//3), T(1//5)])
+               nothing, [T(1//3), T(1//5)])
 
 """
     prob_dde_2delays_notinplace
@@ -209,7 +209,7 @@ prob_dde_2delays_notinplace = build_prob_dde_2delays_notinplace(1.0)
 
 build_prob_dde_2delays_scalar_notinplace(u₀, ::T=u₀) where {T} =
     DDEProblem(f_2delays_notinplace, t -> zero(u₀), u₀, (zero(T), oneunit(T)),
-               [T(1//3), T(1//5)])
+               nothing, [T(1//3), T(1//5)])
 
 """
     prob_dde_2delays_scalar_notinplace
@@ -225,12 +225,12 @@ prob_dde_2delays_scalar_notinplace = build_prob_dde_2delays_scalar_notinplace(1.
 
 ### In-place function
 
-function f_1delay_long(du::T, u, h, p, t) where T
+function f_1delay_long(du, u, h, p, t::T) where T
     du[1] = (- h(t - T(1//5))[1] + u[1]) / oneunit(t)
 end
 
 build_prob_dde_1delay_long(u₀, ::T=u₀) where {T} =
-    DDEProblem(f_1delay_long, t -> [zero(u₀)], [u₀], (zero(T), T(100)), [T(1//5)])
+    DDEProblem(f_1delay_long, t -> [zero(u₀)], [u₀], (zero(T), T(100)), nothing, [T(1//5)])
 
 """
     prob_dde_1delay_long
@@ -258,13 +258,13 @@ prob_dde_1delay_long = build_prob_dde_1delay_long(1.0)
 
 ### Not in-place function
 
-function f_1delay_long_notinplace(u::T, h, p, t) where T
+function f_1delay_long_notinplace(u, h, p, t::T) where T
     (- h(t - T(1//5)) .+ u ) ./ oneunit(t)
 end
 
 build_prob_dde_1delay_long_notinplace(u₀, ::T=u₀) where {T} =
     DDEProblem(f_1delay_long_notinplace, t -> [zero(u₀)], [u₀], (zero(T), T(100)),
-               [T(1//5)])
+               nothing, [T(1//5)])
 
 """
     prob_dde_1delay_long_notinplace
@@ -275,7 +275,7 @@ in-place function.
 prob_dde_1delay_long_notinplace = build_prob_dde_1delay_long_notinplace(1.0)
 
 build_prob_dde_1delay_long_scalar_notinplace(u₀, ::T=u₀) where {T} =
-    DDEProblem(f_1delay_long_notinplace, t -> zero(u₀), u₀, (zero(T), T(100)), [T(1//5)])
+    DDEProblem(f_1delay_long_notinplace, t -> zero(u₀), u₀, (zero(T), T(100)), nothing, [T(1//5)])
 
 """
     prob_dde_1delay_long_scalar_notinplace
@@ -289,12 +289,12 @@ prob_dde_1delay_long_scalar_notinplace = build_prob_dde_1delay_long_scalar_notin
 
 ### In-place function
 
-function f_2delays_long(du::T, u, h, p, t) where T
+function f_2delays_long(du, u, h, p, t::T) where T
     du[1] = (- h(t - T(1//3))[1] - h(t - T(1//5))[1]) / oneunit(t)
 end
 
 build_prob_dde_2delays_long(u₀, ::T=u₀) where {T} =
-    DDEProblem(f_2delays_long, t -> [zero(u₀)], [u₀], (zero(T), T(100)), [T(1//3), T(1//5)])
+    DDEProblem(f_2delays_long, t -> [zero(u₀)], [u₀], (zero(T), T(100)), nothing, [T(1//3), T(1//5)])
 
 """
     prob_dde_2delays_long
@@ -322,15 +322,15 @@ prob_dde_2delays_long = build_prob_dde_2delays_long(1.0)
 
 ### Not in-place function
 
-function f_2delays_long_notinplace(u::T, h, p, t) where T
+function f_2delays_long_notinplace(u, h, p, t::T) where T
     (- h(t - T(1//3)) .- h(t - T(1//5))) ./ oneunit(t)
 end
 
 #### Vectorized history function
 
 build_prob_dde_2delays_long_notinplace(u₀, ::T=u₀) where {T} =
      DDEProblem(f_2delays_long_notinplace, t -> [zero(u₀)], [u₀], (zero(T), T(100)),
-                [T(1//3), T(1//5)])
+                nothing, [T(1//3), T(1//5)])
 
 """
     prob_dde_2delays_long_notinplace
@@ -344,7 +344,7 @@ prob_dde_2delays_long_notinplace = build_prob_dde_2delays_long_notinplace(1.0)
 
 build_prob_dde_2delays_long_scalar_notinplace(u₀, ::T=u₀) where {T} =
      DDEProblem(f_2delays_long_notinplace, t -> zero(u₀), u₀, (zero(T), T(100)),
-                [T(1//3), T(1//5)])
+                nothing, [T(1//3), T(1//5)])
 
 """
     prob_dde_2delays_long_scalar_notinplace
@@ -370,7 +370,7 @@ end
 Model of blood production with constant delay (M. C. Mackey and L. Glass, Oscillation and
 chaos in physiological control systems, 1977).
 """
-prob_dde_mackey = DDEProblem(f_dde_mackey, t -> [0.5], [0.5], (0.0, 500.0), [14])
+prob_dde_mackey = DDEProblem(f_dde_mackey, t -> [0.5], [0.5], (0.0, 500.0), nothing, [14])
 
 function f_dde_wheldon(du, u, h, p, t)
     du[1] = 1.1/(1 + sqrt(10)*(h(t-20)[1])^(5/4)) - 10*u[1]/(1 + 40*u[2])
@@ -402,7 +402,7 @@ DDE with vanishing time dependent delay at ``t = 1`` (K. W. Neves, Automatic int
 of functional differential equations: An approach, 1975).
 """
 prob_dde_neves_1 = DDEProblem(f_dde_neves1, t -> [log(t)], [log(0.1)], (0.1, 10.), [],
-                              [(t, u) -> t - exp(1 - 1/t)])
+                              nothing, [(t, u) -> t - exp(1 - 1/t)])
 
 function f_dde_neves_thompson(du, u, h, p, t)
     if h(t/2)[1] < 0
@@ -447,7 +447,7 @@ u(t) = 1
 for ``t \\leq 0``.
 """
 prob_dde_neves_thompson = DDEProblem(f_dde_neves_thompson, t -> [1.], [1.], (0., 2*log(66)),
-                                     [], [(t, u) -> t/2])
+                                     nothing, [], [(t, u) -> t/2])
 
 function f_dde_paul1(du, u, h, p, t)
     du[1] = - 2*h(t - 1 - abs(u[1]))[1]*(1 - u[1]^2)
@@ -471,7 +471,7 @@ u(t) = 1/2
 
 for ``t \\leq 0``.
 """
-prob_dde_paul1 = DDEProblem(f_dde_paul1, t -> [0.5], [0.5], (0., 30.), [],
+prob_dde_paul1 = DDEProblem(f_dde_paul1, t -> [0.5], [0.5], (0., 30.), nothing, [],
                              [(t, u) -> 1 + abs(u[1])])
 
 function f_dde_paul2(du, u, h, p, t)
@@ -485,7 +485,7 @@ end
 
 DDE with state dependent delay (C. A. H. Paul, A test set of functional differential equations, 1994).
 """
-prob_dde_paul2 = DDEProblem(f_dde_paul2, t -> [1; 0.5], [1; 0.5], (0., 30.), [],
+prob_dde_paul2 = DDEProblem(f_dde_paul2, t -> [1; 0.5], [1; 0.5], (0., 30.), nothing, [],
                              [(t, u) -> u[2]])
 
 function build_prob_dde_mahaffy(tspan, h, σ₀, T₁, γ, Q, k, a, K, r)
@@ -495,7 +495,7 @@ function build_prob_dde_mahaffy(tspan, h, σ₀, T₁, γ, Q, k, a, K, r)
         du[3] = 1 - Q*exp(γ*u[3])/(σ₀*h(t-T₁-u[3])[2])
     end
 
-    DDEProblem(f, h, h(0), tspan, [T₁], [(t, u) -> T₁ + u[3]])
+    DDEProblem(f, h, h(0), tspan, nothing, [T₁], [(t, u) -> T₁ + u[3]])
 end
 
 function h_mahaffy1(t)
@@ -535,7 +535,7 @@ DDE with vanishing state dependent delay at ``t = 1`` (K. W. Neves, Automatic in
 of functional differential equations: An approach, 1975).
 """
 prob_dde_neves2 = DDEProblem(f_dde_neves2, t -> [log(t); 1/t], [log(0.1); 10], (0.1, 5.),
-                              [], [(t, u) -> t - exp(1 - u[2])])
+                             nothing, [], [(t, u) -> t - exp(1 - u[2])])
 
 function build_f_dde_gatica(r₁, r₂, α, δ)
     function f_dde_gatica(du, u, h, p, t)
@@ -559,7 +559,7 @@ Model of antigen antibody dynamics with fading memory, with vanishing state depe
 at ``t = 0`` (J. Gatica and P. Waltman, A threshold model of antigen antibody dynamics with fading memory, 1982).
 """
 prob_dde_gatica = DDEProblem(build_f_dde_gatica(0.02, 0.005, 3, 0.01), t -> [5; 0.1; 0; 0],
-                             [5; 0.1; 0; 0], (0., 40.), [], [(t, u) -> u[4]])
+                             [5; 0.1; 0; 0], (0., 40.), nothing, [], [(t, u) -> u[4]])
 
 #=
 Quorum Sensing model
@@ -611,7 +611,7 @@ function build_prob_dde_qs(u₀, tspan, τ, D, γₛ, Kₘ, nₛ, a, αₐ, β
         end
     end
 
-    DDEProblem(f_dde_qs, t -> u₀, u₀, tspan, [τ])
+    DDEProblem(f_dde_qs, t -> u₀, u₀, tspan, nothing, [τ])
 end
 
 """