open-space-toolkit-astrodynamics 2.4.0__py39-none-any.whl → 3.0.0__py39-none-any.whl

{open_space_toolkit_astrodynamics-2.4.0.dist-info → open_space_toolkit_astrodynamics-3.0.0.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: open-space-toolkit-astrodynamics
-Version: 2.4.0
+Version: 3.0.0
 Summary: Orbit, attitude, access.
 Author: Open Space Collective
 Author-email: contact@open-space-collective.org

{open_space_toolkit_astrodynamics-2.4.0.dist-info → open_space_toolkit_astrodynamics-3.0.0.dist-info}/RECORD

@@ -1,9 +1,9 @@
 ostk/__init__.py,sha256=epnVn2PwdQkUDZ1msqBRO5nEZIOUBIq-IfK3IlNPijE,21
-ostk/astrodynamics/OpenSpaceToolkitAstrodynamicsPy.cpython-39-x86_64-linux-gnu.so,sha256=p8Q2kmoC_eidYT-XofoL3k0H5SjWzAHrPZZFxyvJp7A,1414584
+ostk/astrodynamics/OpenSpaceToolkitAstrodynamicsPy.cpython-39-x86_64-linux-gnu.so,sha256=YXAEb6MOEWwW1KOz7sCniXUmonQS-SHeg8LryEQRLCo,1447384
 ostk/astrodynamics/__init__.py,sha256=BE7yyrOnPqFCTPT3GhYEI4IWR5nzs96VSjKphhQ09ew,97
 ostk/astrodynamics/converters.py,sha256=sBZ3_JMjowWbQ9gBPA2cSEkbBfAwXmt1IBbLks76-Vk,3502
 ostk/astrodynamics/display.py,sha256=y9FnoQbPFGM6LzkUdgXgeqtuVGhv57GuTKbeDdcFPgw,6306
-ostk/astrodynamics/libopen-space-toolkit-astrodynamics.so.2,sha256=3Zkok2yopaE5WjopTGYXl_oIh65GrnI7W8unc0jUkCk,69730504
+ostk/astrodynamics/libopen-space-toolkit-astrodynamics.so.3,sha256=6Q6q--9IL9_hOCb-ooRA_zUBP_GsZe33716f0OfkZzU,72418128
 ostk/astrodynamics/utilities.py,sha256=NbeKN_CtcCRHIinhoNWTEosBMR7xc924S-ZFD97mduw,4478
 ostk/astrodynamics/viewer.py,sha256=_4uucNtZrIADtKtEYHbiu91RVVvKRG3vPULUAshPLJU,8912
 ostk/astrodynamics/test/__init__.py,sha256=epnVn2PwdQkUDZ1msqBRO5nEZIOUBIq-IfK3IlNPijE,21
@@ -40,20 +40,21 @@ ostk/astrodynamics/test/flight/system/test_satellite_system.py,sha256=-5KOGVTWRI
 ostk/astrodynamics/test/flight/system/dynamics/__init__.py,sha256=epnVn2PwdQkUDZ1msqBRO5nEZIOUBIq-IfK3IlNPijE,21
 ostk/astrodynamics/test/flight/system/dynamics/test_atmospheric_drag.py,sha256=bTFDvI1riLFlgKzs3TRt8dafFtb--L_AMEDuuRWEgGc,3793
 ostk/astrodynamics/test/flight/system/dynamics/test_central_body_gravity.py,sha256=r3JIIX7cD63BhrgvLtjFaK7ew1BqHdKP5_VxwGmuv1Q,1818
-ostk/astrodynamics/test/flight/system/dynamics/test_dynamics.py,sha256=arNkAOj2y4CDmc9CXxItbFtC4QYjaByuJefdi5S6alw,914
+ostk/astrodynamics/test/flight/system/dynamics/test_dynamics.py,sha256=eD5lpF8C8b-AMrLOXlQSo-j00-5AHTaQOybbSJJmvGc,1254
 ostk/astrodynamics/test/flight/system/dynamics/test_position_derivative.py,sha256=ZVVewpqZwqiBLDqoNQDHDgDodtj6_M0J6ZmgHe_1gz8,1646
 ostk/astrodynamics/test/flight/system/dynamics/test_third_body_gravity.py,sha256=Z8OkM-EFafFpWziu3bfatlSlZA4Tic2fsVLU35rDp5E,1915
-ostk/astrodynamics/test/flight/system/dynamics/thruster/test_constant_thrust.py,sha256=MP6ocA0umkASJfidTni62APX9tmSxxXFQ4dxvMTKrH4,4480
+ostk/astrodynamics/test/flight/system/dynamics/thruster/test_constant_thrust.py,sha256=EdEI0tMXXJVtaj7XogLuCkfK24iX9J1sxx6rr7_y-1o,4464
 ostk/astrodynamics/test/solvers/__init__.py,sha256=epnVn2PwdQkUDZ1msqBRO5nEZIOUBIq-IfK3IlNPijE,21
 ostk/astrodynamics/test/solvers/test_temporal_condition_solver.py,sha256=s7r7VelXx1_b-6wO_GyOlmxlsfA4I_zpxrWojriQls8,4867
 ostk/astrodynamics/test/trajectory/__init__.py,sha256=epnVn2PwdQkUDZ1msqBRO5nEZIOUBIq-IfK3IlNPijE,21
-ostk/astrodynamics/test/trajectory/test_local_orbital_frame_direction.py,sha256=DJltOvRcxaEJbZsH_PQe7oPdN9KRE6ufXmJ1kqUw6O4,2518
+ostk/astrodynamics/test/trajectory/test_local_orbital_frame_direction.py,sha256=d72J50UGG9m0Atei0UQ8sITU1WdJmws5xgUiacLZMbw,2515
 ostk/astrodynamics/test/trajectory/test_local_orbital_frame_factory.py,sha256=FziYh9XxD2FEuj6WhnqOe5mF9twLTQtvcem8KzQZKXg,1864
 ostk/astrodynamics/test/trajectory/test_model.py,sha256=epnVn2PwdQkUDZ1msqBRO5nEZIOUBIq-IfK3IlNPijE,21
 ostk/astrodynamics/test/trajectory/test_orbit.py,sha256=SlCBZ_broPaJ6L8khMuvbdPem9-zmqG6IZwxr9cECB4,3155
-ostk/astrodynamics/test/trajectory/test_propagator.py,sha256=e_fmvvaoMhYs0ApqCgbzfBYrYPaGOIuzouo-6x2jbK4,12504
-ostk/astrodynamics/test/trajectory/test_state.py,sha256=Ohjc_cQ28XoOtEY4mP-EuyrCKirRf0-RxJNzFVFqiJ4,3699
-ostk/astrodynamics/test/trajectory/test_trajectory_segment.py,sha256=eOCI9GLcAnMrwXkdaHbd0thnyNDCyBgAemhuMO6B7xY,3863
+ostk/astrodynamics/test/trajectory/test_propagator.py,sha256=ct4Sg8l3R-PQus0NYBviu5ExUXQC5f8TcPGqJiYUA_0,12208
+ostk/astrodynamics/test/trajectory/test_segment.py,sha256=T4OJ25z4-ogAV8WN9W7CEM_EMQtvVmGD8-Iql4NjJY8,5060
+ostk/astrodynamics/test/trajectory/test_sequence.py,sha256=n-B2_h7EYNd--Slw0680U4IRPnhbNYzcXOzNNBALyiU,10350
+ostk/astrodynamics/test/trajectory/test_state.py,sha256=Kr4IbGBc1TNGdEOsXFLGSQGC9SwCER-lJ6E-xumqgN4,3829
 ostk/astrodynamics/test/trajectory/orbit/__init__.py,sha256=epnVn2PwdQkUDZ1msqBRO5nEZIOUBIq-IfK3IlNPijE,21
 ostk/astrodynamics/test/trajectory/orbit/test_model.py,sha256=epnVn2PwdQkUDZ1msqBRO5nEZIOUBIq-IfK3IlNPijE,21
 ostk/astrodynamics/test/trajectory/orbit/test_pass.py,sha256=MtPRKhVz4hCzfFuE0p1x89Ps68M9s7miT7VtPIDVeOo,1898
@@ -76,8 +77,8 @@ ostk/astrodynamics/test/trajectory/state/test_coordinates_subset.py,sha256=eMmif
 ostk/astrodynamics/test/trajectory/state/test_numerical_solver.py,sha256=ZTWRSVFOJnjzQP7JDsHZxvZ4ldyDHmXTptG2csC7_ic,7480
 ostk/astrodynamics/test/trajectory/state/coordinates_subset/test_cartesian_position.py,sha256=V5d7UnV6KIJjJvh5R0i3imRfyN04JCDvCUsKdvfKmhE,2708
 ostk/astrodynamics/test/trajectory/state/coordinates_subset/test_cartesian_velocity.py,sha256=5a7yeOVt7gY4oQdINSElb_U2jznJmuWPrT_YJrI6KaM,3070
-open_space_toolkit_astrodynamics-2.4.0.dist-info/METADATA,sha256=R8ipuXPGkD7i4U6XXBZxvb8nJrsbXWQIw03RKlBdSyE,1779
-open_space_toolkit_astrodynamics-2.4.0.dist-info/WHEEL,sha256=x0YL6OSV_WitIuMP0HsTpwBgQTKZeHzZKOAyhSZVN4o,93
-open_space_toolkit_astrodynamics-2.4.0.dist-info/top_level.txt,sha256=zOR18699uDYnafgarhL8WU_LmTZY_5NVqutv-flp_x4,5
-open_space_toolkit_astrodynamics-2.4.0.dist-info/zip-safe,sha256=AbpHGcgLb-kRsJGnwFEktk7uzpZOCcBY74-YBdrKVGs,1
-open_space_toolkit_astrodynamics-2.4.0.dist-info/RECORD,,
+open_space_toolkit_astrodynamics-3.0.0.dist-info/METADATA,sha256=tIw96qI_UMBkuf9PFsOjObKhca4j9kWzWQoAOFi2BM4,1779
+open_space_toolkit_astrodynamics-3.0.0.dist-info/WHEEL,sha256=x0YL6OSV_WitIuMP0HsTpwBgQTKZeHzZKOAyhSZVN4o,93
+open_space_toolkit_astrodynamics-3.0.0.dist-info/top_level.txt,sha256=zOR18699uDYnafgarhL8WU_LmTZY_5NVqutv-flp_x4,5
+open_space_toolkit_astrodynamics-3.0.0.dist-info/zip-safe,sha256=AbpHGcgLb-kRsJGnwFEktk7uzpZOCcBY74-YBdrKVGs,1
+open_space_toolkit_astrodynamics-3.0.0.dist-info/RECORD,,

ostk/astrodynamics/test/flight/system/dynamics/test_dynamics.py

@@ -2,6 +2,8 @@
 
 import pytest
 
+from ostk.physics import Environment
+
 from ostk.astrodynamics.flight.system import Dynamics
 from ostk.astrodynamics.trajectory.state.coordinates_subset import CartesianPosition
 
@@ -11,6 +13,11 @@ def name() -> str:
     return "MyDynamics"
 
 
+@pytest.fixture
+def environment() -> Environment:
+    return Environment.default()
+
+
 @pytest.fixture
 def dynamics(name: str) -> Dynamics:
     class MyDynamics(Dynamics):
@@ -35,3 +42,9 @@ class TestDynamics:
 
     def test_get_name(self, dynamics: Dynamics, name: str):
         assert dynamics.get_name() == name
+
+    def test_from_environment(self, environment: Environment):
+        dynamics: list[Dynamics] = Dynamics.from_environment(environment)
+
+        assert dynamics is not None
+        assert isinstance(dynamics, list)

ostk/astrodynamics/test/flight/system/dynamics/thruster/test_constant_thrust.py

@@ -12,8 +12,6 @@ from ostk.physics.units import Mass
 from ostk.physics.time import Instant
 from ostk.physics.time import DateTime
 from ostk.physics.time import Scale
-from ostk.physics.coordinate import Position
-from ostk.physics.coordinate import Velocity
 from ostk.physics.coordinate import Frame
 
 from ostk.astrodynamics.trajectory import LocalOrbitalFrameFactory
@@ -34,8 +32,8 @@ from ostk.astrodynamics.flight.system.dynamics.thruster import ConstantThrust
 @pytest.fixture
 def local_orbital_frame_direction() -> LocalOrbitalFrameDirection:
     return LocalOrbitalFrameDirection(
-        [1.0, 0.0, 0.0],
-        LocalOrbitalFrameFactory.VNC(Frame.GCRF()),
+        vector=[1.0, 0.0, 0.0],
+        local_orbital_frame_factory=LocalOrbitalFrameFactory.VNC(Frame.GCRF()),
     )
 
 
@@ -77,9 +75,9 @@ def dynamics(
     local_orbital_frame_direction: LocalOrbitalFrameDirection,
 ) -> ConstantThrust:
     return ConstantThrust(
-        satellite_system,
-        local_orbital_frame_direction,
-        "Constant Thrust Thruster Dynamics",
+        satellite_system=satellite_system,
+        thrust_direction=local_orbital_frame_direction,
+        name="Constant Thrust Thruster Dynamics",
     )
 
 

ostk/astrodynamics/test/trajectory/test_local_orbital_frame_direction.py

@@ -2,7 +2,6 @@
 
 import pytest
 
-from ostk.physics.time import Instant
 from ostk.physics.coordinate import Frame
 
 from ostk.astrodynamics.trajectory import LocalOrbitalFrameFactory
@@ -25,8 +24,8 @@ def local_orbital_frame_direction(
     local_orbital_frame_factory: LocalOrbitalFrameFactory,
 ) -> LocalOrbitalFrameDirection:
     return LocalOrbitalFrameDirection(
-        direction_vector,
-        local_orbital_frame_factory,
+        vector=direction_vector,
+        local_orbital_frame_factory=local_orbital_frame_factory,
     )
 
 

ostk/astrodynamics/test/trajectory/test_propagator.py

@@ -372,16 +372,9 @@ class TestPropagator:
         self,
         numerical_solver: NumericalSolver,
         environment: Environment,
-        satellite_system: SatelliteSystem,
     ):
         assert Propagator.from_environment(numerical_solver, environment) is not None
 
-        assert (
-            Propagator.from_environment(numerical_solver, environment, satellite_system)
-            is not None
-        )
-
-    def test_default(self, environment: Environment, satellite_system: SatelliteSystem):
+    def test_default(self, environment: Environment):
         assert Propagator.default()
         assert Propagator.default(environment) is not None
-        assert Propagator.default(environment, satellite_system) is not None

ostk/astrodynamics/test/trajectory/{test_trajectory_segment.py → test_segment.py}

@@ -13,10 +13,14 @@ from ostk.physics.coordinate import Frame
 from ostk.physics.environment.objects.celestial_bodies import Earth
 
 from ostk.astrodynamics.trajectory.state import NumericalSolver
+from ostk.astrodynamics.flight.system import SatelliteSystem
 from ostk.astrodynamics.flight.system.dynamics import CentralBodyGravity
 from ostk.astrodynamics.flight.system.dynamics import PositionDerivative
+from ostk.astrodynamics.flight.system.dynamics.thruster import ConstantThrust
+from ostk.astrodynamics.trajectory import LocalOrbitalFrameDirection
+from ostk.astrodynamics.trajectory import LocalOrbitalFrameFactory
 from ostk.astrodynamics.trajectory import State
-from ostk.astrodynamics.trajectory import TrajectorySegment
+from ostk.astrodynamics.trajectory import Segment
 from ostk.astrodynamics.event_condition import InstantCondition
 
 
@@ -84,45 +88,82 @@ def coast_duration_segment(
     dynamics: list,
     numerical_solver: NumericalSolver,
 ):
-    return TrajectorySegment.coast(name, instant_condition, dynamics, numerical_solver)
+    return Segment.coast(name, instant_condition, dynamics, numerical_solver)
 
 
-class TestTrajectorySegment:
-    def test_get_name(self, coast_duration_segment: TrajectorySegment, name: str):
+@pytest.fixture
+def thruster_dynamics() -> ConstantThrust:
+    return ConstantThrust(
+        satellite_system=SatelliteSystem.default(),
+        thrust_direction=LocalOrbitalFrameDirection(
+            vector=[1.0, 0.0, 0.0],
+            local_orbital_frame_factory=LocalOrbitalFrameFactory.VNC(Frame.GCRF()),
+        ),
+    )
+
+
+class TestSegment:
+    def test_get_name(self, coast_duration_segment: Segment, name: str):
         assert coast_duration_segment.get_name() == name
 
     def test_get_event_condition(
         self,
         instant_condition: InstantCondition,
-        coast_duration_segment: TrajectorySegment,
+        coast_duration_segment: Segment,
     ):
         assert coast_duration_segment.get_event_condition() == instant_condition
 
     def test_get_dynamics(
         self,
         dynamics: list,
-        coast_duration_segment: TrajectorySegment,
+        coast_duration_segment: Segment,
     ):
         assert len(coast_duration_segment.get_dynamics()) == len(dynamics)
 
     def test_get_numerical_solver(
         self,
         numerical_solver: NumericalSolver,
-        coast_duration_segment: TrajectorySegment,
+        coast_duration_segment: Segment,
     ):
         assert coast_duration_segment.get_numerical_solver() == numerical_solver
 
     def test_get_type(
         self,
-        coast_duration_segment: TrajectorySegment,
+        coast_duration_segment: Segment,
     ):
-        assert coast_duration_segment.get_type() == TrajectorySegment.Type.Coast
+        assert coast_duration_segment.get_type() == Segment.Type.Coast
+
+    def test_coast(
+        self,
+        name: str,
+        instant_condition: InstantCondition,
+        dynamics: list,
+        numerical_solver: NumericalSolver,
+    ):
+        assert (
+            Segment.coast(name, instant_condition, dynamics, numerical_solver) is not None
+        )
+
+    def test_maneuver(
+        self,
+        name: str,
+        instant_condition: InstantCondition,
+        thruster_dynamics: ConstantThrust,
+        dynamics: list,
+        numerical_solver: NumericalSolver,
+    ):
+        assert (
+            Segment.maneuver(
+                name, instant_condition, thruster_dynamics, dynamics, numerical_solver
+            )
+            is not None
+        )
 
     def test_solve(
         self,
         state: State,
         end_instant: Instant,
-        coast_duration_segment: TrajectorySegment,
+        coast_duration_segment: Segment,
     ):
         solution = coast_duration_segment.solve(state)
 

ostk/astrodynamics/test/trajectory/test_sequence.py

@@ -0,0 +1,379 @@
+# Apache License 2.0
+
+import pytest
+
+import numpy as np
+
+from ostk.mathematics.geometry.d3.objects import Composite
+from ostk.mathematics.geometry.d3.objects import Cuboid
+from ostk.mathematics.geometry.d3.objects import Point
+
+from ostk.physics import Environment
+from ostk.physics.coordinate import Frame
+from ostk.physics.environment.atmospheric import Earth as EarthAtmosphericModel
+from ostk.physics.environment.gravitational import Earth as EarthGravitationalModel
+from ostk.physics.environment.magnetic import Earth as EarthMagneticModel
+from ostk.physics.environment.objects.celestial_bodies import Earth
+from ostk.physics.time import DateTime
+from ostk.physics.time import Duration
+from ostk.physics.time import Instant
+from ostk.physics.time import Scale
+from ostk.physics.units import Derived
+from ostk.physics.units import Length
+from ostk.physics.units import Mass
+
+from ostk.astrodynamics.event_condition import COECondition
+from ostk.astrodynamics.event_condition import InstantCondition
+from ostk.astrodynamics.event_condition import RealCondition
+from ostk.astrodynamics.flight.system import Dynamics
+from ostk.astrodynamics.flight.system import PropulsionSystem
+from ostk.astrodynamics.flight.system import SatelliteSystem
+from ostk.astrodynamics.flight.system.dynamics.thruster import ConstantThrust
+from ostk.astrodynamics.trajectory import LocalOrbitalFrameDirection
+from ostk.astrodynamics.trajectory import LocalOrbitalFrameFactory
+from ostk.astrodynamics.trajectory import Segment
+from ostk.astrodynamics.trajectory import Sequence
+from ostk.astrodynamics.trajectory import State
+from ostk.astrodynamics.trajectory.state import CoordinatesBroker
+from ostk.astrodynamics.trajectory.state import CoordinatesSubset
+from ostk.astrodynamics.trajectory.state import NumericalSolver
+from ostk.astrodynamics.trajectory.state.coordinates_subset import CartesianPosition
+from ostk.astrodynamics.trajectory.state.coordinates_subset import CartesianVelocity
+
+
+@pytest.fixture
+def propulsion_system() -> PropulsionSystem:
+    return PropulsionSystem(
+        1.0,
+        1500.0,
+    )
+
+
+@pytest.fixture
+def dry_mass() -> Mass:
+    return Mass.kilograms(100.0)
+
+
+@pytest.fixture
+def wet_mass() -> Mass:
+    return Mass.kilograms(10.0)
+
+
+@pytest.fixture
+def cross_sectional_surface_area() -> float:
+    return 1.0
+
+
+@pytest.fixture
+def drag_coefficient() -> float:
+    return 2.2
+
+
+@pytest.fixture
+def satellite_system(
+    dry_mass: Mass,
+    cross_sectional_surface_area: float,
+    drag_coefficient: float,
+    propulsion_system: PropulsionSystem,
+) -> SatelliteSystem:
+    satellite_geometry = Composite(
+        Cuboid(
+            Point(0.0, 0.0, 0.0),
+            [[1.0, 0.0, 0.0], [0.0, 1.0, 0.0], [0.0, 0.0, 1.0]],
+            [1.0, 0.0, 0.0],
+        )
+    )
+    inertia_tensor = np.ndarray(shape=(3, 3))
+
+    return SatelliteSystem(
+        dry_mass,
+        satellite_geometry,
+        inertia_tensor,
+        cross_sectional_surface_area,
+        drag_coefficient,
+        propulsion_system,
+    )
+
+
+@pytest.fixture
+def environment() -> Environment:
+    return Environment(
+        Instant.J2000(),
+        [
+            Earth.from_models(
+                EarthGravitationalModel(EarthGravitationalModel.Type.Spherical),
+                EarthMagneticModel(EarthMagneticModel.Type.Undefined),
+                EarthAtmosphericModel(EarthAtmosphericModel.Type.Exponential),
+            )
+        ],
+    )
+
+
+@pytest.fixture
+def coordinates_broker() -> CoordinatesBroker:
+    return CoordinatesBroker(
+        [
+            CartesianPosition.default(),
+            CartesianVelocity.default(),
+            CoordinatesSubset.mass(),
+            CoordinatesSubset.surface_area(),
+            CoordinatesSubset.drag_coefficient(),
+        ]
+    )
+
+
+@pytest.fixture
+def cartesian_coordinates() -> list[float]:
+    return [
+        717094.039086306,
+        -6872433.2241124,
+        46175.9696673281,
+        -970.650826004612,
+        -45.4598114773158,
+        7529.82424886455,
+    ]
+
+
+@pytest.fixture
+def state(
+    cartesian_coordinates: list[float],
+    dry_mass: Mass,
+    wet_mass: Mass,
+    cross_sectional_surface_area: float,
+    drag_coefficient: float,
+    coordinates_broker: CoordinatesBroker,
+) -> State:
+    frame: Frame = Frame.GCRF()
+    instant: Instant = Instant.date_time(DateTime(2018, 1, 1, 0, 0, 0), Scale.UTC)
+    coordinates = [
+        *cartesian_coordinates,
+        dry_mass.in_kilograms() + wet_mass.in_kilograms(),
+        cross_sectional_surface_area,
+        drag_coefficient,
+    ]
+
+    return State(
+        instant,
+        [
+            717094.039086306,
+            -6872433.2241124,
+            46175.9696673281,
+            -970.650826004612,
+            -45.4598114773158,
+            7529.82424886455,
+            dry_mass.in_kilograms() + wet_mass.in_kilograms(),
+            cross_sectional_surface_area,
+            drag_coefficient,
+        ],
+        frame,
+        coordinates_broker,
+    )
+
+
+@pytest.fixture
+def dynamics(environment: Environment) -> list:
+    return Dynamics.from_environment(environment)
+
+
+@pytest.fixture
+def numerical_solver() -> NumericalSolver:
+    return NumericalSolver.default_conditional()
+
+
+@pytest.fixture
+def duration() -> Duration:
+    return Duration.minutes(5.0)
+
+
+@pytest.fixture
+def instant_condition(state: State, duration: Duration) -> InstantCondition:
+    return InstantCondition(
+        RealCondition.Criterion.AnyCrossing, state.get_instant() + duration
+    )
+
+
+@pytest.fixture
+def sma() -> Length:
+    return Length.kilometers(6907.000)
+
+
+@pytest.fixture
+def gravitational_parameter() -> Derived:
+    return EarthGravitationalModel.spherical.gravitational_parameter
+
+
+@pytest.fixture
+def constant_thrust(satellite_system: SatelliteSystem) -> ConstantThrust:
+    return ConstantThrust(
+        satellite_system=satellite_system,
+        thrust_direction=LocalOrbitalFrameDirection(
+            vector=[1.0, 0.0, 0.0],
+            local_orbital_frame_factory=LocalOrbitalFrameFactory.VNC(Frame.GCRF()),
+        ),
+    )
+
+
+@pytest.fixture
+def sma_condition(sma: Length, gravitational_parameter: Derived) -> COECondition:
+    return COECondition.semi_major_axis(
+        criterion=RealCondition.Criterion.AnyCrossing,
+        frame=Frame.GCRF(),
+        semi_major_axis=sma,
+        gravitational_parameter=gravitational_parameter,
+    )
+
+
+@pytest.fixture
+def coast_duration_segment(
+    instant_condition: InstantCondition,
+    dynamics: list[Dynamics],
+    numerical_solver: NumericalSolver,
+):
+    return Segment.coast(
+        name="duration coast",
+        event_condition=instant_condition,
+        dynamics=dynamics,
+        numerical_solver=numerical_solver,
+    )
+
+
+@pytest.fixture
+def coast_sma_segment(
+    sma_condition: COECondition,
+    dynamics: list[Dynamics],
+    numerical_solver: NumericalSolver,
+):
+    return Segment.coast(
+        name="sma coast",
+        event_condition=sma_condition,
+        dynamics=dynamics,
+        numerical_solver=numerical_solver,
+    )
+
+
+@pytest.fixture
+def thrust_segment(
+    sma_condition: RealCondition,
+    constant_thrust: ConstantThrust,
+    dynamics: list[Dynamics],
+    numerical_solver: NumericalSolver,
+):
+    return Segment.maneuver(
+        name="duration thrust",
+        event_condition=sma_condition,
+        thruster_dynamics=constant_thrust,
+        dynamics=dynamics,
+        numerical_solver=numerical_solver,
+    )
+
+
+@pytest.fixture
+def segments(
+    coast_duration_segment: Segment,
+    thrust_segment: Segment,
+) -> list[Segment]:
+    return [coast_duration_segment, thrust_segment]
+
+
+@pytest.fixture
+def maximum_propagation_duration() -> Duration:
+    return Duration.days(2.0)
+
+
+@pytest.fixture
+def repetition_count() -> int:
+    return 1
+
+
+@pytest.fixture
+def sequence(
+    segments: list[Segment],
+    repetition_count: int,
+    numerical_solver: NumericalSolver,
+    dynamics: list,
+    maximum_propagation_duration: Duration,
+):
+    return Sequence(
+        segments=segments,
+        repetition_count=repetition_count,
+        dynamics=dynamics,
+        numerical_solver=numerical_solver,
+        maximum_propagation_duration=maximum_propagation_duration,
+    )
+
+
+class TestSequence:
+    def test_get_segments(
+        self,
+        sequence: Sequence,
+        segments: list[Segment],
+    ):
+        assert len(sequence.get_segments()) == len(segments)
+
+    def test_get_numerical_solver(
+        self,
+        sequence: Sequence,
+        numerical_solver: NumericalSolver,
+    ):
+        assert sequence.get_numerical_solver() == numerical_solver
+
+    def test_get_dynamics(
+        self,
+        sequence: Sequence,
+        dynamics: list,
+    ):
+        assert len(sequence.get_dynamics()) == len(dynamics)
+
+    def test_get_maximum_propagation_duration(
+        self,
+        sequence: Sequence,
+        maximum_propagation_duration: Duration,
+    ):
+        assert sequence.get_maximum_propagation_duration() == maximum_propagation_duration
+
+    def test_add_segment(
+        self,
+        sequence: Sequence,
+        coast_duration_segment: Segment,
+    ):
+        segments_count: int = len(sequence.get_segments())
+
+        sequence.add_segment(coast_duration_segment)
+
+        assert len(sequence.get_segments()) == segments_count + 1
+
+        segments_count = len(sequence.get_segments())
+
+        sequence.add_segment([coast_duration_segment, coast_duration_segment])
+
+        assert len(sequence.get_segments()) == segments_count + 2
+
+    def test_add_coast_segment(
+        self,
+        sequence: Sequence,
+        instant_condition: InstantCondition,
+    ):
+        segments_count: int = len(sequence.get_segments())
+
+        sequence.add_coast_segment(instant_condition)
+
+        assert len(sequence.get_segments()) == segments_count + 1
+
+    def test_add_maneuver_segment(
+        self,
+        sequence: Sequence,
+        instant_condition: InstantCondition,
+        constant_thrust: ConstantThrust,
+    ):
+        segments_count: int = len(sequence.get_segments())
+
+        sequence.add_maneuver_segment(instant_condition, constant_thrust)
+
+        assert len(sequence.get_segments()) == segments_count + 1
+
+    def test_solve(self, state: State, sequence: Sequence, segments: list[Segment]):
+        solution = sequence.solve(state)
+
+        assert len(solution.segment_solutions) == len(segments)
+
+        assert solution.get_states() is not None

ostk/astrodynamics/test/trajectory/test_state.py

@@ -102,6 +102,7 @@ class TestState:
         position: Position,
         velocity: Velocity,
         frame: Frame,
+        coordinates_broker: CoordinatesBroker,
     ):
         assert state.get_instant() == instant
         assert state.get_position() == position
@@ -111,6 +112,7 @@ class TestState:
             state.get_coordinates()
             == np.append(position.get_coordinates(), velocity.get_coordinates())
         ).all()
+        assert state.get_coordinates_subsets() == coordinates_broker.get_subsets()
 
     def test_in_frame(
         self,