Disable computation of second-order vector field for Lagrangian systems.

JuliaGNI · Apr 16, 2024 · ef810d8 · ef810d8
1 parent 04856f8
commit ef810d8
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 3 deletions.
diff --git a/src/lagrangian.jl b/src/lagrangian.jl
@@ -70,7 +70,7 @@ struct LagrangianSystem
 
         equs_subs = substitute_lagrangian_variables(equs, x, ẋ, v)
         equs_subs = merge(equs_subs, (
-            a = inv(equs_subs.M) * (equs_subs.f - equs_subs.N * V),
+            # a = inv(equs_subs.M) * (equs_subs.f - equs_subs.N * V),
             ϕ = P .- equs_subs.ϑ,
             ψ = F .- equs_subs.g,
             # σ = _simplify(inv(equs_subs.ω), dosimplify),
@@ -86,7 +86,7 @@ struct LagrangianSystem
         code = (
             L  = substitute_parameters(build_function(equs_subs.L,  t, X, V, params...; nanmath = false), params),
             EL = substitute_parameters(build_function(equs_subs.EL, t, X, V, params...; nanmath = false)[2], params),
-            a  = substitute_parameters(build_function(equs_subs.a,  t, X, V, params...; nanmath = false)[2], params),
+            # a  = substitute_parameters(build_function(equs_subs.a,  t, X, V, params...; nanmath = false)[2], params),
             f  = substitute_parameters(build_function(equs_subs.f,  t, X, V, params...; nanmath = false)[2], params),
             g  = substitute_parameters(build_function(equs_subs.g,  t, X, Λ, V, params...; nanmath = false)[2], params),
             p  = substitute_parameters(build_function(equs_subs.ϑ,  t, X, V, params...; nanmath = false)[1], params),

diff --git a/test/lagrangian_particle.jl b/test/lagrangian_particle.jl
@@ -37,7 +37,7 @@ for k in (:L, :EL, :f, :g, :ϑ, :θ, :ω, :Ω, :ϕ, :ψ, :M, :N)
     @test k ∈ keys(EulerLagrange.equations(lag_sys))
 end
 
-for k in (:L, :EL, :a, :f, :g, :p, :ϑ, :θ, :ω, :Ω, :ϕ, :ψ, :M)#, :P
+for k in (:L, :EL, :f, :g, :p, :ϑ, :θ, :ω, :Ω, :ϕ, :ψ, :M)# :a, :P
     @test k ∈ keys(EulerLagrange.functions(lag_sys))
 end