Adds support for MTK v9.3

Author: joeycarpinelli

Project.toml

@@ -24,7 +24,7 @@ SPICEBodiesExt = "SPICEBodies"
 AstrodynamicalCalculations = "^0.5.1"
 DocStringExtensions = "0.9"
 Memoize = "0.4"
-ModelingToolkit = "8"
+ModelingToolkit = "~9.3"
 SciMLBase = "2"
 StaticArrays = "1"
 Symbolics = "5"

src/Attitude.jl

@@ -155,8 +155,9 @@ end
 @memoize function AttitudeFunction(; stm=false, name=:Attitude, kwargs...)
     defaults = (; jac=true)
     options = merge(defaults, kwargs)
+    sys = complete(AttitudeSystem(; stm=stm, name=name); split=false)
     return ODEFunction{true,SciMLBase.FullSpecialize}(
-        complete(AttitudeSystem(; stm=stm, name=name));
+        sys, ModelingToolkit.unknowns(sys), ModelingToolkit.parameters(sys);
         options...
     )
 end

src/CR3BP.jl

@@ -107,8 +107,9 @@ end
 @memoize function CR3BFunction(; stm=false, name=:CR3B, kwargs...)
     defaults = (; jac=true)
     options = merge(defaults, kwargs)
+    sys = complete(CR3BSystem(; stm=stm, name=name); split=false)
     return ODEFunction{true,SciMLBase.FullSpecialize}(
-        complete(CR3BSystem(; stm=stm, name=name));
+        sys, ModelingToolkit.unknowns(sys), ModelingToolkit.parameters(sys);
         options...
     )
 end

src/Entry.jl

@@ -119,8 +119,9 @@ end
 @memoize function PlanarEntryFunction(; name=:PlanarEntry, kwargs...)
     defaults = (; jac=true)
     options = merge(defaults, kwargs)
+    sys = complete(PlanarEntrySystem(; name=name); split=false)
     return ODEFunction{true,SciMLBase.FullSpecialize}(
-        complete(PlanarEntrySystem(; name=name));
+        sys, ModelingToolkit.unknowns(sys), ModelingToolkit.parameters(sys);
         options...
     )
 end

src/NBP.jl

@@ -155,8 +155,9 @@ end
         calculations! Consider setting `jac=false`, `stm=false`, or both.
         """
     end
+    sys = complete(NBSystem(N; stm=stm, name=name); split=false)
     return ODEFunction{true,SciMLBase.FullSpecialize}(
-        complete(NBSystem(N; stm=stm, name=name));
+        sys, ModelingToolkit.unknowns(sys), ModelingToolkit.parameters(sys);
         options...
     )
 end

src/R2BP.jl

@@ -107,8 +107,9 @@ end
 @memoize function R2BFunction(; stm=false, name=:R2B, kwargs...)
     defaults = (; jac=true)
     options = merge(defaults, kwargs)
+    sys = complete(R2BSystem(; stm=stm, name=name); split=false)
     return ODEFunction{true,SciMLBase.FullSpecialize}(
-        complete(R2BSystem(; stm=stm, name=name));
+        sys, ModelingToolkit.unknowns(sys), ModelingToolkit.parameters(sys);
         options...
     )
 end