diff --git a/src/chaotic_attractors.jl b/src/chaotic_attractors.jl index 674a82b..9c20fa5 100644 --- a/src/chaotic_attractors.jl +++ b/src/chaotic_attractors.jl @@ -476,49 +476,17 @@ function Duffing() return prob end -function MackeyGlass end -originalcode(::typeof(MackeyGlass)) = """ -class MackeyGlass(DynSysDelay): - @staticjit - def _rhs(x, xt, t, beta, gamma, n, tau): - xdot = beta * (xt / (1 + xt ** n)) - gamma * x - return xdot -""" -@doc make_docstring(MackeyGlass) MackeyGlass -function MackeyGlass() - function rhs(du, u, p, t) - beta, gamma, n = p - du[1] = beta * (u[2] / (1 + u[2]^n)) - gamma * u[1] - end - u0 = Float64.(ATTRACTOR_DATA["MackeyGlass"]["initial_conditions"]) - p = dict_to_componentarray(ATTRACTOR_DATA["MackeyGlass"]["parameters"]) - tspan = (0.0, 1.0) - f = ODEFunction(rhs) - prob = DDEProblem(f, u0, tspan, p, constant_lags = [1.0]) - return prob -end +# class MackeyGlass(DynSysDelay): +# @staticjit +# def _rhs(x, xt, t, beta, gamma, n, tau): +# xdot = beta * (xt / (1 + xt ** n)) - gamma * x +# return xdot -function IkedaDelay end -originalcode(::typeof(IkedaDelay)) = """ -class IkedaDelay(DynSysDelay): - @staticjit - def _rhs(x, xt, t, c, mu, tau, x0): - xdot = mu * np.sin(xt - x0) - c * x - return xdot -""" -@doc make_docstring(IkedaDelay) IkedaDelay -function IkedaDelay() - function rhs(du, u, p, t) - c, mu, x0 = p - du[1] = mu * sin(u[2] - x0) - c * u[1] - end - u0 = Float64.(ATTRACTOR_DATA["IkedaDelay"]["initial_conditions"]) - p = dict_to_componentarray(ATTRACTOR_DATA["IkedaDelay"]["parameters"]) - tspan = (0.0, 1.0) - f = ODEFunction(rhs) - prob = DDEProblem(f, u0, tspan, p, constant_lags = [1.0]) - return prob -end +# class IkedaDelay(DynSysDelay): +# @staticjit +# def _rhs(x, xt, t, c, mu, tau, x0): +# xdot = mu * np.sin(xt - x0) - c * x +# return xdot # class SprottDelay(IkedaDelay): # pass