remove @inline

Author: Datseris

src/continuous_famous_systems.jl

@@ -27,14 +27,14 @@ function's documentation string.
 function lorenz(u0=[0.0, 10.0, 0.0]; σ = 10.0, ρ = 28.0, β = 8/3)
     return CDS(loop, u0, [σ, ρ, β], loop_jac)
 end
-@inline @inbounds function loop(u, p, t)
+@inbounds function loop(u, p, t)
     σ = p[1]; ρ = p[2]; β = p[3]
     du1 = σ*(u[2]-u[1])
     du2 = u[1]*(ρ-u[3]) - u[2]
     du3 = u[1]*u[2] - β*u[3]
     return SVector{3}(du1, du2, du3)
 end
-@inline @inbounds function loop_jac(u, p, t)
+@inbounds function loop_jac(u, p, t)
     σ, ρ, β = p
     J = @SMatrix [-σ  σ  0;
     ρ - u[3]  (-1)  (-u[1]);
@@ -45,14 +45,14 @@ end
 function lorenz_iip(u0=[0.0, 10.0, 0.0]; σ = 10.0, ρ = 28.0, β = 8/3)
     return CDS(liip, u0, [σ, ρ, β], liip_jac)
 end
-@inline @inbounds function liip(du, u, p, t)
+@inbounds function liip(du, u, p, t)
     σ = p[1]; ρ = p[2]; β = p[3]
     du[1] = σ*(u[2]-u[1])
     du[2] = u[1]*(ρ-u[3]) - u[2]
     du[3] = u[1]*u[2] - β*u[3]
     return nothing
 end
-@inline @inbounds function liip_jac(J, u, p, t)
+@inbounds function liip_jac(J, u, p, t)
     σ, ρ, β = p
     J[1,1] = -σ; J[1, 2] = σ; J[1,3] = 0
     J[2,1] = ρ - u[3]; J[2,2] = -1; J[2,3] = -u[1]
@@ -87,14 +87,14 @@ function's documentation string.
 function roessler(u0=rand(3); a = 0.2, b = 0.2, c = 5.7)
     return CDS(roessler_eom, u0, [a, b, c], roessler_jacob)
 end
-@inline @inbounds function roessler_eom(u, p, t)
+@inbounds function roessler_eom(u, p, t)
     a, b, c = p
     du1 = -u[2]-u[3]
     du2 = u[1] + a*u[2]
     du3 = b + u[3]*(u[1] - c)
     return SVector{3, Float64}(du1, du2, du3)
 end
-@inline @inbounds function roessler_jacob(u, p, t)
+@inbounds function roessler_jacob(u, p, t)
     a, b, c = p
     return @SMatrix [0 (-1) (-1);
                      1 a 0;
@@ -168,10 +168,10 @@ function henonheiles(u0=[0, -0.25, 0.42081, 0]#=; conserveE::Bool = true=#)
     o = zero(eltype(u0))
     J = zeros(eltype(u0), 4, 4)
 
-    # @inline Vhh(q1, q2) = 1//2 * (q1^2 + q2^2 + 2q1^2 * q2 - 2//3 * q2^3)
-    # @inline Thh(p1, p2) = 1//2 * (p1^2 + p2^2)
-    # @inline Hhh(q1, q2, p1, p2) = Thh(p1, p2) + Vhh(q1, q2)
-    # @inline Hhh(u::AbstractVector) = Hhh(u...)
+    # Vhh(q1, q2) = 1//2 * (q1^2 + q2^2 + 2q1^2 * q2 - 2//3 * q2^3)
+    # Thh(p1, p2) = 1//2 * (p1^2 + p2^2)
+    # Hhh(q1, q2, p1, p2) = Thh(p1, p2) + Vhh(q1, q2)
+    # Hhh(u::AbstractVector) = Hhh(u...)
     #
     # E = Hhh(u0)
     #

src/discrete.jl

@@ -30,7 +30,7 @@ function reinit!(integ::MDI, u = integ.u, Q0 = nothing)
     return
 end
 
-@inline function (integ::MinimalDiscreteIntegrator)(t::Real)
+function (integ::MinimalDiscreteIntegrator)(t::Real)
     if t == integ.t
         return integ.u
     else

src/discrete_famous_systems.jl

@@ -23,13 +23,13 @@ Default values are the ones used in the original paper.
 function towel(u0=[0.085, -0.121, 0.075])
     return DDS(eom_towel, u0, nothing, jacob_towel)
 end# should result in lyapunovs: [0.432207,0.378834,-3.74638]
-@inline function eom_towel(x, p, n)
+function eom_towel(x, p, n)
     @inbounds x1, x2, x3 = x[1], x[2], x[3]
     SVector( 3.8*x1*(1-x1) - 0.05*(x2+0.35)*(1-2*x3),
     0.1*( (x2+0.35)*(1-2*x3) - 1 )*(1 - 1.9*x1),
     3.78*x3*(1-x3)+0.2*x2 )
 end
-@inline function jacob_towel(x, p, n)
+function jacob_towel(x, p, n)
     @SMatrix [3.8*(1 - 2x[1]) -0.05*(1-2x[3]) 0.1*(x[2] + 0.35);
     -0.19((x[2] + 0.35)*(1-2x[3]) - 1)  0.1*(1-2x[3])*(1-1.9x[1])  -0.2*(x[2] + 0.35)*(1-1.9x[1]);
     0.0  0.2  3.78(1-2x[3]) ]
@@ -96,7 +96,7 @@ Nuclear Physics, Novosibirsk (1969)
 function standardmap(u0=0.001rand(2); k = 0.971635)
     return DDS(standardmap_eom, u0, [k], standardmap_jacob)
 end
-@inline @inbounds function standardmap_eom(x, par, n)
+@inbounds function standardmap_eom(x, par, n)
     theta = x[1]; p = x[2]
     p += par[1]*sin(theta)
     theta += p
@@ -106,7 +106,7 @@ end
     while p < 0; p += twopi; end
     return SVector(theta, p)
 end
-@inline @inbounds standardmap_jacob(x, p, n) =
+@inbounds standardmap_jacob(x, p, n) =
 @SMatrix [1 + p[1]*cos(x[1])    1;
           p[1]*cos(x[1])        1]
 
@@ -221,8 +221,8 @@ function's documentation string.
 function henon(u0=zeros(2); a = 1.4, b = 0.3)
     return DDS(hoop, u0, [a,b], hoop_jac)
 end # should give lyapunov exponents [0.4189, -1.6229]
-@inline hoop(x, p, n) = SVector{2}(1.0 - p[1]*x[1]^2 + x[2], p[2]*x[1])
-@inline hoop_jac(x, p, n) = @SMatrix [-2*p[1]*x[1] 1.0; p[2] 0.0]
+hoop(x, p, n) = SVector{2}(1.0 - p[1]*x[1]^2 + x[2], p[2]*x[1])
+hoop_jac(x, p, n) = @SMatrix [-2*p[1]*x[1] 1.0; p[2] 0.0]
 
 function henon_iip(u0=zeros(2); a = 1.4, b = 0.3)
     return DDS(hiip, u0, [a, b], hiip_jac)
@@ -265,5 +265,5 @@ function's documentation string.
 function logistic(x0=rand(); r = 4.0)
     return DDS(logistic_eom, x0, [r], logistic_jacob)
 end
-@inline logistic_eom(x, p, n) = p[1]*x*(1-x)
-@inline logistic_jacob(x, p, n) = p[1]*(1-2x)
+logistic_eom(x, p, n) = p[1]*x*(1-x)
+logistic_jacob(x, p, n) = p[1]*(1-2x)

src/dynamicalsystem.jl

@@ -140,10 +140,10 @@ const DDS = DiscreteDynamicalSystem
 systemtype(::DDS) = "discrete"
 
 
-@inline isinplace(::DS{IIP}) where {IIP} = IIP
-@inline statetype(::DS{IIP, S}) where {IIP, S} = S
-@inline stateeltype(::DS{IIP, S}) where {IIP, S} = eltype(S)
-@inline isautodiff(::DS{IIP, S, D, F, P, JAC, JM, IAD}) where
+isinplace(::DS{IIP}) where {IIP} = IIP
+statetype(::DS{IIP, S}) where {IIP, S} = S
+stateeltype(::DS{IIP, S}) where {IIP, S} = eltype(S)
+isautodiff(::DS{IIP, S, D, F, P, JAC, JM, IAD}) where
 {IIP, S, D, F, P, JAC, JM, IAD} = IAD
 
 get_state(ds::DS) = ds.u0
@@ -288,7 +288,7 @@ end
 #######################################################################################
 #                                    Jacobians                                        #
 #######################################################################################
-@inline function create_jacobian(
+function create_jacobian(
     f::F, ::Val{IIP}, s::S, p::P, t::T, ::Val{D}) where {F, IIP, S, P, T, D}
     if IIP
         dum = deepcopy(s)
@@ -311,7 +311,7 @@ end
     end
 end
 
-@inline function get_J(jacob, u0, p, t0, iip) where {JAC}
+function get_J(jacob, u0, p, t0, iip) where {JAC}
     D = length(u0)
     if iip
         J = similar(u0, (D,D))
@@ -340,7 +340,7 @@ ds.jacobian(u, ds.p, t)
 #                                 Tanget Dynamics                                     #
 #######################################################################################
 # IIP Tangent Space dynamics
-@inline function create_tangent(f::F, jacobian::JAC, J::JM,
+function create_tangent(f::F, jacobian::JAC, J::JM,
     ::Val{true}, ::Val{k}) where {F, JAC, JM, k}
     J = deepcopy(J)
     tangentf = (du, u, p, t) -> begin
@@ -376,7 +376,7 @@ end
 
 
 # OOP Tangent Space dynamics
-@inline function create_tangent(f::F, jacobian::JAC, J::JM,
+function create_tangent(f::F, jacobian::JAC, J::JM,
     ::Val{false}, ::Val{k}) where {F, JAC, JM, k}
 
     ws_index = SVector{k, Int}(2:(k+1)...)
@@ -390,8 +390,8 @@ struct TangentOOP{F, JAC, k} <: Function
 end
 function (tan::TangentOOP)(u, p, t)
     @inbounds s = u[:, 1]
-    @inbounds du = tan.f(s, p, t)
-    @inbounds J = tan.jacobian(s, p, t)
+    du = tan.f(s, p, t)
+    J = tan.jacobian(s, p, t)
     @inbounds dW = J*u[:, tan.ws]
     return hcat(du, dW)
 end