Merge pull request #49 from dpad/fix/ModelingToolkit-breaks-model-arrays

[Bugfix] Temporary fix for dynamical models not working under ModelingToolkit v5.21.0.
dpad · Jun 28, 2021 · 84261f8 · 84261f8 · dpad · Jun 28, 2021
2 parents 486e3a3 + e47c3a6
commit 84261f8
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 2 deletions.
diff --git a/Project.toml b/Project.toml
@@ -1,7 +1,7 @@
 name = "OrbitalTrajectories"
 uuid = "2b613a20-8d2a-5290-b19f-e06f4bcc2e7d"
 authors = ["Dan Padilha"]
-version = "0.1.13"
+version = "0.1.14"
 
 [deps]
 DataStructures = "864edb3b-99cc-5e75-8d2d-829cb0a9cfe8"

diff --git a/src/dynamics/dynamical_models.jl b/src/dynamics/dynamical_models.jl
@@ -57,7 +57,7 @@ has_jacobian(X::Type{<:DiffEqBase.ODEFunction}) = !isnothing(fieldtype(X, :jac))
     # The State Transition Matrix (STM) ODE function is defined as follows, including the N^2 Jacobian equations +
     # the N first-order equations of motion. [Koon 2011]
     # NOTE: the Differential is defined element-wise and flattened to a list.
-    stm_eqs = simplify.(D.(ϕ) .~ A * ϕ)
+    stm_eqs = collect(simplify.(D.(ϕ) .~ A * ϕ))
 
     # Create the ODE system and generate its functions
     stm_ode = ODESystem(