open-space-toolkit-astrodynamics 5.1.5__py38-none-any.whl → 5.2.0__py38-none-any.whl

ostk/astrodynamics/test/conjunction/messages/ccsds/test_cdm.py

@@ -0,0 +1,416 @@
+# Apache License 2.0
+
+import pytest
+
+from datetime import datetime
+
+from ostk.core.containers import Dictionary
+from ostk.core.filesystem import File
+
+from ostk.physics.units import Length
+from ostk.physics.units import Mass
+from ostk.physics.time import Instant
+from ostk.physics.time import Scale
+from ostk.physics.time import Duration
+from ostk.physics.coordinate import Frame
+from ostk.physics.coordinate import Position
+from ostk.physics.coordinate import Velocity
+
+from ostk.astrodynamics.trajectory import State
+from ostk.astrodynamics.conjunction.messages.ccsds import CDM
+
+
+@pytest.fixture
+def cdm() -> CDM:
+    return CDM(
+        header=CDM.Header(
+            ccsds_cdm_version="1.0",
+            creation_date=Instant.date_time(datetime(2010, 3, 12, 22, 31, 12), Scale.UTC),
+            originator="JSPOC",
+            message_for="SATELLITE A",
+            message_id="201113719185",
+        ),
+        relative_metadata=CDM.RelativeMetadata(
+            comment="",
+            time_of_closest_approach=Instant.date_time(
+                datetime(2010, 3, 13, 22, 37, 52, 618), Scale.UTC
+            ),
+            miss_distance=Length.meters(715.0),
+            start_screen_period=Instant.date_time(
+                datetime(2010, 3, 12, 18, 29, 32, 212), Scale.UTC
+            ),
+            end_screen_period=Instant.date_time(
+                datetime(2010, 3, 15, 18, 29, 32, 212), Scale.UTC
+            ),
+            screen_volume_frame="RTN",
+            screen_volume_shape="ELLIPSOID",
+            screen_volume_x=200.0,
+            screen_volume_y=1000.0,
+            screen_volume_z=1000.0,
+            screen_entry_time=Instant.date_time(
+                datetime(2010, 3, 13, 22, 37, 52, 222), Scale.UTC
+            ),
+            screen_exit_time=Instant.date_time(
+                datetime(2010, 3, 13, 22, 37, 52, 824), Scale.UTC
+            ),
+            collision_probability=0.0000435,
+            collision_probability_method="FOSTER-1992",
+        ),
+        objects_metadata_array=[
+            CDM.Metadata(
+                comment="",
+                object="OBJECT1",
+                object_designator=12345,
+                catalog_name="SATCAT",
+                object_name="SATELLITE A",
+                international_designator="1997−030E",
+                object_type=CDM.ObjectType.Payload,
+                operator_contact_position="OSA",
+                operator_organization="EUTMETSAT",
+                operator_phone="+49615130312",
+                operator_email="JOHN.DOE@SOMEWHERE.NET",
+                ephemeris_name="EPHEMERIS SATELLITE A",
+                covariance_method="CALCULATED",
+                maneuverable="YES",
+                orbit_center="",
+                reference_frame="EME2000",
+                gravity_model="EGM-96: 36D 36O",
+                atmospheric_model="JACCHIA 70 DCA",
+                n_body_perturbations="MOON, SUN",
+                solar_radiation_pressure=False,
+                earth_tides=False,
+                in_track_thrust=False,
+            ),
+            CDM.Metadata(
+                object="OBJECT2",
+                object_designator=12346,
+                catalog_name="SATCAT",
+                object_name="SATELLITE B",
+                international_designator="1997−030F",
+                object_type=CDM.ObjectType.Payload,
+                ephemeris_name="EPHEMERIS SATELLITE B",
+                covariance_method="CALCULATED",
+                maneuverable="YES",
+                reference_frame="EME2000",
+            ),
+        ],
+        objects_data_array=[
+            CDM.Data(
+                time_last_observation_start=Instant.date_time(
+                    datetime(2010, 3, 12, 2, 14, 12, 746), Scale.UTC
+                ),
+                time_last_observation_end=Instant.date_time(
+                    datetime(2010, 3, 12, 2, 14, 12, 746), Scale.UTC
+                ),
+                recommended_od_span=Duration.days(7.88),
+                actual_od_span=Duration.days(5.50),
+                observations_available=592,
+                observations_used=579,
+                tracks_available=123,
+                tracks_used=119,
+                residuals_accepted=97.8,
+                weighted_rms=0.864,
+                area_pc=5.2,
+                area_drag=0.0,
+                area_srp=0.0,
+                mass=Mass.kilograms(251.6),
+                cd_area_over_mass=0.0,
+                cr_area_over_mass=0.0,
+                thrust_acceleration=0.0,
+                sedr=0.0,
+            ),
+            CDM.Data(
+                time_last_observation_start=Instant.date_time(
+                    datetime(2010, 3, 12, 2, 14, 12, 746), Scale.UTC
+                ),
+                time_last_observation_end=Instant.date_time(
+                    datetime(2010, 3, 12, 2, 14, 12, 746), Scale.UTC
+                ),
+                recommended_od_span=Duration.days(7.88),
+                actual_od_span=Duration.days(5.50),
+                observations_available=592,
+                observations_used=579,
+                tracks_available=123,
+                tracks_used=119,
+                residuals_accepted=97.8,
+                weighted_rms=0.864,
+                area_pc=5.2,
+                area_drag=0.0,
+                area_srp=0.0,
+                mass=Mass.kilograms(251.6),
+                cd_area_over_mass=0.0,
+                cr_area_over_mass=0.0,
+                thrust_acceleration=0.0,
+                sedr=0.0,
+            ),
+        ],
+    )
+
+
+class TestCDM:
+    def test_constructor(self, cdm: CDM):
+        assert cdm is not None
+
+    def test_get_header(self, cdm: CDM):
+        cdm_header: CDM.Header = cdm.get_header()
+
+        assert cdm_header is not None
+        assert cdm_header.ccsds_cdm_version == "1.0"
+        assert cdm_header.creation_date == Instant.date_time(
+            datetime(2010, 3, 12, 22, 31, 12), Scale.UTC
+        )
+        assert cdm_header.originator == "JSPOC"
+        assert cdm_header.message_for == "SATELLITE A"
+        assert cdm_header.message_id == "201113719185"
+
+    def test_get_relative_metadata(self, cdm: CDM):
+        cdm_relative_metadata: CDM.RelativeMetadata = cdm.get_relative_metadata()
+
+        assert cdm_relative_metadata is not None
+        assert cdm_relative_metadata.comment == ""
+        assert cdm_relative_metadata.time_of_closest_approach == Instant.date_time(
+            datetime(2010, 3, 13, 22, 37, 52, 618), Scale.UTC
+        )
+        assert cdm_relative_metadata.miss_distance == Length.meters(715.0)
+        assert cdm_relative_metadata.start_screen_period == Instant.date_time(
+            datetime(2010, 3, 12, 18, 29, 32, 212), Scale.UTC
+        )
+        assert cdm_relative_metadata.end_screen_period == Instant.date_time(
+            datetime(2010, 3, 15, 18, 29, 32, 212), Scale.UTC
+        )
+        assert cdm_relative_metadata.screen_volume_frame == "RTN"
+        assert cdm_relative_metadata.screen_volume_shape == "ELLIPSOID"
+        assert cdm_relative_metadata.screen_volume_x == 200.0
+        assert cdm_relative_metadata.screen_volume_y == 1000.0
+        assert cdm_relative_metadata.screen_volume_z == 1000.0
+        assert cdm_relative_metadata.screen_entry_time == Instant.date_time(
+            datetime(2010, 3, 13, 22, 37, 52, 222), Scale.UTC
+        )
+        assert cdm_relative_metadata.screen_exit_time == Instant.date_time(
+            datetime(2010, 3, 13, 22, 37, 52, 824), Scale.UTC
+        )
+        assert cdm_relative_metadata.collision_probability == 0.0000435
+        assert cdm_relative_metadata.collision_probability_method == "FOSTER-1992"
+
+    def test_get_metadata_array(self, cdm: CDM):
+        assert cdm.get_metadata_array() is not None
+        assert len(cdm.get_metadata_array()) == 2
+
+    def test_get_data_array(self, cdm: CDM):
+        assert cdm.get_data_array() is not None
+        assert len(cdm.get_data_array()) == 2
+
+    def test_get_object_metadata_at(self, cdm: CDM):
+        assert cdm.get_object_metadata_at(0) is not None
+        assert cdm.get_object_metadata_at(1) is not None
+
+        cdm_object_1_metadata: CDM.Metadata = cdm.get_object_metadata_at(0)
+
+        assert cdm_object_1_metadata.comment == ""
+        assert cdm_object_1_metadata.object == "OBJECT1"
+        assert cdm_object_1_metadata.object_designator == 12345
+        assert cdm_object_1_metadata.catalog_name == "SATCAT"
+        assert cdm_object_1_metadata.object_name == "SATELLITE A"
+        assert cdm_object_1_metadata.international_designator == "1997−030E"
+        assert cdm_object_1_metadata.object_type == CDM.ObjectType.Payload
+        assert cdm_object_1_metadata.operator_contact_position == "OSA"
+        assert cdm_object_1_metadata.operator_organization == "EUTMETSAT"
+        assert cdm_object_1_metadata.operator_phone == "+49615130312"
+        assert cdm_object_1_metadata.operator_email == "JOHN.DOE@SOMEWHERE.NET"
+        assert cdm_object_1_metadata.ephemeris_name == "EPHEMERIS SATELLITE A"
+        assert cdm_object_1_metadata.covariance_method == "CALCULATED"
+        assert cdm_object_1_metadata.maneuverable == "YES"
+        assert cdm_object_1_metadata.orbit_center == ""
+        assert cdm_object_1_metadata.reference_frame == "EME2000"
+        assert cdm_object_1_metadata.gravity_model == "EGM-96: 36D 36O"
+        assert cdm_object_1_metadata.atmospheric_model == "JACCHIA 70 DCA"
+        assert cdm_object_1_metadata.n_body_perturbations == "MOON, SUN"
+        assert cdm_object_1_metadata.solar_radiation_pressure == False
+        assert cdm_object_1_metadata.earth_tides == False
+        assert cdm_object_1_metadata.in_track_thrust == False
+
+        cdm_object_2_metadata: CDM.Metadata = cdm.get_object_metadata_at(1)
+
+        assert cdm_object_2_metadata.comment == ""
+        assert cdm_object_2_metadata.object == "OBJECT2"
+        assert cdm_object_2_metadata.object_designator == 12346
+        assert cdm_object_2_metadata.catalog_name == "SATCAT"
+        assert cdm_object_2_metadata.object_name == "SATELLITE B"
+        assert cdm_object_2_metadata.international_designator == "1997−030F"
+        assert cdm_object_2_metadata.object_type == CDM.ObjectType.Payload
+        assert cdm_object_2_metadata.operator_contact_position == ""
+        assert cdm_object_2_metadata.operator_organization == ""
+        assert cdm_object_2_metadata.operator_phone == ""
+        assert cdm_object_2_metadata.operator_email == ""
+        assert cdm_object_2_metadata.ephemeris_name == "EPHEMERIS SATELLITE B"
+        assert cdm_object_2_metadata.covariance_method == "CALCULATED"
+        assert cdm_object_2_metadata.maneuverable == "YES"
+        assert cdm_object_2_metadata.orbit_center == ""
+        assert cdm_object_2_metadata.reference_frame == "EME2000"
+        assert cdm_object_2_metadata.gravity_model == ""
+        assert cdm_object_2_metadata.atmospheric_model == ""
+        assert cdm_object_2_metadata.n_body_perturbations == ""
+        assert cdm_object_2_metadata.solar_radiation_pressure == False
+        assert cdm_object_2_metadata.earth_tides == False
+        assert cdm_object_2_metadata.in_track_thrust == False
+
+        with pytest.raises(Exception) as e:
+            cdm.get_object_metadata_at(2)
+
+    def test_get_object_data_at(self, cdm: CDM):
+        assert cdm.get_object_data_at(0) is not None
+        assert cdm.get_object_data_at(1) is not None
+
+        cdm_object_1_data: CDM.Data = cdm.get_object_data_at(0)
+
+        assert cdm_object_1_data.time_last_observation_start == Instant.date_time(
+            datetime(2010, 3, 12, 2, 14, 12, 746), Scale.UTC
+        )
+        assert cdm_object_1_data.time_last_observation_end == Instant.date_time(
+            datetime(2010, 3, 12, 2, 14, 12, 746), Scale.UTC
+        )
+        assert cdm_object_1_data.recommended_od_span == Duration.days(7.88)
+        assert cdm_object_1_data.actual_od_span == Duration.days(5.50)
+        assert cdm_object_1_data.observations_available == 592
+        assert cdm_object_1_data.observations_used == 579
+        assert cdm_object_1_data.tracks_available == 123
+        assert cdm_object_1_data.tracks_used == 119
+        assert cdm_object_1_data.residuals_accepted == 97.8
+        assert cdm_object_1_data.weighted_rms == 0.864
+        assert cdm_object_1_data.area_pc == 5.2
+        assert cdm_object_1_data.area_drag == 0.0
+        assert cdm_object_1_data.area_srp == 0.0
+        assert cdm_object_1_data.mass == Mass.kilograms(251.6)
+        assert cdm_object_1_data.cd_area_over_mass == 0.0
+        assert cdm_object_1_data.cr_area_over_mass == 0.0
+        assert cdm_object_1_data.thrust_acceleration == 0.0
+        assert cdm_object_1_data.sedr == 0.0
+
+        cdm_object_2_data: CDM.Data = cdm.get_object_data_at(1)
+
+        assert cdm_object_2_data.time_last_observation_start == Instant.date_time(
+            datetime(2010, 3, 12, 2, 14, 12, 746), Scale.UTC
+        )
+        assert cdm_object_2_data.time_last_observation_end == Instant.date_time(
+            datetime(2010, 3, 12, 2, 14, 12, 746), Scale.UTC
+        )
+        assert cdm_object_2_data.recommended_od_span == Duration.days(7.88)
+        assert cdm_object_2_data.actual_od_span == Duration.days(5.50)
+        assert cdm_object_2_data.observations_available == 592
+        assert cdm_object_2_data.observations_used == 579
+        assert cdm_object_2_data.tracks_available == 123
+        assert cdm_object_2_data.tracks_used == 119
+        assert cdm_object_2_data.residuals_accepted == 97.8
+        assert cdm_object_2_data.weighted_rms == 0.864
+        assert cdm_object_2_data.area_pc == 5.2
+        assert cdm_object_2_data.area_drag == 0.0
+        assert cdm_object_2_data.area_srp == 0.0
+        assert cdm_object_2_data.mass == Mass.kilograms(251.6)
+        assert cdm_object_2_data.cd_area_over_mass == 0.0
+        assert cdm_object_2_data.cr_area_over_mass == 0.0
+        assert cdm_object_2_data.thrust_acceleration == 0.0
+        assert cdm_object_2_data.sedr == 0.0
+
+        with pytest.raises(Exception) as e:
+            cdm.get_object_data_at(2)
+
+    def test_get_ccsds_cdm_version(self, cdm: CDM):
+        assert cdm.get_ccsds_cdm_version() == "1.0"
+
+    def test_get_creation_instant(self, cdm: CDM):
+        assert cdm.get_creation_instant() == Instant.date_time(
+            datetime(2010, 3, 12, 22, 31, 12), Scale.UTC
+        )
+
+    def test_get_originator(self, cdm: CDM):
+        assert cdm.get_originator() == "JSPOC"
+
+    def test_get_message_for(self, cdm: CDM):
+        assert cdm.get_message_for() == "SATELLITE A"
+
+    def test_get_message_id(self, cdm: CDM):
+        assert cdm.get_message_id() == "201113719185"
+
+    def test_get_time_of_closest_approach(self, cdm: CDM):
+        assert cdm.get_time_of_closest_approach() == Instant.date_time(
+            datetime(2010, 3, 13, 22, 37, 52, 618), Scale.UTC
+        )
+
+    def test_get_miss_distance(self, cdm: CDM):
+        assert cdm.get_miss_distance() == Length.meters(715.0)
+
+    def test_get_relative_position(self, cdm: CDM):
+        assert cdm.get_relative_position().is_defined() is False
+
+    def test_get_relative_velocity(self, cdm: CDM):
+        assert cdm.get_relative_velocity().is_defined() is False
+
+    def test_get_collision_probability(self, cdm: CDM):
+        assert cdm.get_collision_probability() == 0.0000435
+
+    def test_get_collision_probability_method(self, cdm: CDM):
+        assert cdm.get_collision_probability_method() == "FOSTER-1992"
+
+    def test_undefined(self):
+        assert CDM.undefined().is_defined() is False
+
+    def test_dictionary(self, cdm_spacetrack_dictionary: dict):
+        dictionary = Dictionary(cdm_spacetrack_dictionary)
+
+        cdm: CDM = CDM.dictionary(dictionary)
+
+        assert cdm is not None
+        assert cdm.get_header() is not None
+        assert cdm.get_ccsds_cdm_version() == "1.0"
+        assert cdm.get_creation_instant() == Instant.date_time(
+            datetime(2022, 12, 25, 0, 33, 16), Scale.UTC
+        )
+        assert cdm.get_originator() == "CSpOC"
+        assert cdm.get_message_for() == "YAM-2"
+        assert cdm.get_message_id() == "406320986"
+        assert cdm.get_time_of_closest_approach() == Instant.date_time(
+            datetime(2022, 12, 27, 13, 28, 59, 516000), Scale.UTC
+        )
+        assert cdm.get_miss_distance() == Length.meters(974.0)
+        assert cdm.get_relative_position().is_defined() is False
+        assert cdm.get_relative_velocity().is_defined() is False
+        assert cdm.get_collision_probability() == 5.162516e-16
+        assert cdm.get_collision_probability_method() == "FOSTER-1992"
+
+        object_1_metadata: CDM.Metadata = cdm.get_object_metadata_at(0)
+
+        assert object_1_metadata.object == "OBJECT1"
+        assert object_1_metadata.object_designator == 48911
+        assert object_1_metadata.catalog_name == "SATCAT"
+        assert object_1_metadata.object_name == "YAM-2"
+        assert object_1_metadata.international_designator == "2021-059AJ"
+        assert object_1_metadata.object_type == CDM.ObjectType.Payload
+        assert object_1_metadata.ephemeris_name == "NONE"
+        assert object_1_metadata.covariance_method == "CALCULATED"
+        assert object_1_metadata.maneuverable == "N/A"
+        assert object_1_metadata.reference_frame == "ITRF"
+
+        # object_1_data: CDM.Data = cdm.get_object_data_at(0)
+
+        # assert object_1_data.state == State(
+        #     Instant.date_time(datetime(2022, 12, 27, 13, 28, 59, 516000), Scale.UTC),
+        #     Position.meters([-4988150.231999999844, -1691825.955000000075, -4469421.481999999844], Frame.ITRF()),
+        #     Velocity.meters_per_second([-5122.248844000000, 54.300816000000, 5699.434412000000], Frame.ITRF())
+        # )
+
+    def test_parse(self, cdm_file: File):
+        with open(str(cdm_file.get_path().to_string()), "r") as stream:
+            assert CDM.parse(string=stream.read()) is not None
+
+    def test_load(self, cdm_file: File):
+        assert CDM.load(file=cdm_file) is not None
+
+    def test_object_type_from_string(self):
+        assert CDM.object_type_from_string("PAYLOAD") == CDM.ObjectType.Payload
+        assert CDM.object_type_from_string("ROCKET BODY") == CDM.ObjectType.RocketBody
+        assert CDM.object_type_from_string("DEBRIS") == CDM.ObjectType.Debris
+        assert CDM.object_type_from_string("UNKNOWN") == CDM.ObjectType.Unknown
+        assert CDM.object_type_from_string("OTHER") == CDM.ObjectType.Other
+
+        with pytest.raises(Exception) as e:
+            CDM.object_type_from_string("ANOTHERSUPPORTEDTYPE")

ostk/astrodynamics/test/dynamics/__init__.py

@@ -0,0 +1 @@
+# Apache License 2.0

ostk/astrodynamics/test/dynamics/test_atmospheric_drag.py

@@ -0,0 +1,128 @@
+# Apache License 2.0
+
+import pytest
+
+import numpy as np
+
+from ostk.mathematics.geometry.d3.objects import Cuboid
+from ostk.mathematics.geometry.d3.objects import Composite
+from ostk.mathematics.geometry.d3.objects import Point
+
+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.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.astrodynamics.trajectory.state import CoordinatesSubset
+from ostk.astrodynamics.trajectory.state.coordinates_subset import CartesianPosition
+from ostk.astrodynamics.trajectory.state.coordinates_subset import CartesianVelocity
+from ostk.astrodynamics.trajectory.state import CoordinatesBroker
+
+from ostk.astrodynamics.trajectory import State
+from ostk.astrodynamics.trajectory.state import CoordinatesSubset, CoordinatesBroker
+from ostk.astrodynamics.trajectory.state.coordinates_subset import (
+    CartesianPosition,
+    CartesianVelocity,
+)
+from ostk.astrodynamics import Dynamics
+from ostk.astrodynamics.dynamics import AtmosphericDrag
+
+
+@pytest.fixture
+def earth() -> Earth:
+    return Earth.from_models(
+        EarthGravitationalModel(EarthGravitationalModel.Type.Undefined),
+        EarthMagneticModel(EarthMagneticModel.Type.Undefined),
+        EarthAtmosphericModel(EarthAtmosphericModel.Type.Exponential),
+    )
+
+
+@pytest.fixture
+def dry_mass() -> float:
+    return 100.0
+
+
+@pytest.fixture
+def surface_area() -> float:
+    return 500.0
+
+
+@pytest.fixture
+def drag_coefficient() -> float:
+    return 2.1
+
+
+@pytest.fixture
+def dynamics(earth: Earth) -> AtmosphericDrag:
+    return AtmosphericDrag(earth)
+
+
+@pytest.fixture
+def coordinates_broker() -> CoordinatesBroker:
+    return CoordinatesBroker(
+        [
+            CartesianPosition.default(),
+            CartesianVelocity.default(),
+            CoordinatesSubset.mass(),
+            CoordinatesSubset.surface_area(),
+            CoordinatesSubset.drag_coefficient(),
+        ]
+    )
+
+
+@pytest.fixture
+def instant() -> Instant:
+    return Instant.date_time(DateTime(2021, 3, 20, 12, 0, 0), Scale.UTC)
+
+
+@pytest.fixture
+def position_coordinates() -> list:
+    return [7000000.0, 0.0, 0.0]
+
+
+@pytest.fixture
+def velocity_coordinates() -> list:
+    return [0.0, 7546.05329, 0.0]
+
+
+@pytest.fixture
+def state(
+    instant: Instant,
+    position_coordinates: list,
+    velocity_coordinates: list,
+    dry_mass: float,
+    surface_area: float,
+    drag_coefficient: float,
+    coordinates_broker: CoordinatesBroker,
+) -> State:
+    coordinates = (
+        position_coordinates
+        + velocity_coordinates
+        + [dry_mass, surface_area, drag_coefficient]
+    )
+    return State(instant, coordinates, Frame.GCRF(), coordinates_broker)
+
+
+class TestAtmosphericDrag:
+    def test_constructors(self, dynamics: AtmosphericDrag):
+        assert dynamics is not None
+        assert isinstance(dynamics, AtmosphericDrag)
+        assert isinstance(dynamics, Dynamics)
+        assert dynamics.is_defined()
+
+    def test_getters(self, dynamics: AtmosphericDrag, earth: Earth):
+        assert dynamics.get_celestial() == earth
+
+    def test_compute_contribution(self, dynamics: AtmosphericDrag, state: State):
+        contribution = dynamics.compute_contribution(
+            state.get_instant(), state.get_coordinates(), state.get_frame()
+        )
+
+        assert len(contribution) == 3

ostk/astrodynamics/test/dynamics/test_central_body_gravity.py

@@ -0,0 +1,58 @@
+# Apache License 2.0
+
+import pytest
+
+import numpy as np
+
+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.physics.environment.objects.celestial_bodies import Earth
+
+from ostk.astrodynamics.trajectory import State
+from ostk.astrodynamics import Dynamics
+from ostk.astrodynamics.dynamics import CentralBodyGravity
+
+
+@pytest.fixture
+def earth() -> Earth:
+    return Earth.spherical()
+
+
+@pytest.fixture
+def dynamics(earth: Earth) -> CentralBodyGravity:
+    return CentralBodyGravity(earth)
+
+
+@pytest.fixture
+def state() -> State:
+    frame: Frame = Frame.GCRF()
+    position: Position = Position.meters([7000000.0, 0.0, 0.0], frame)
+    velocity: Velocity = Velocity.meters_per_second([0.0, 0.0, 0.0], frame)
+    instant = Instant.date_time(DateTime(2021, 3, 20, 12, 0, 0), Scale.UTC)
+    return State(instant, position, velocity)
+
+
+class TestCentralBodyGravity:
+    def test_constructors(
+        self,
+        dynamics: CentralBodyGravity,
+    ):
+        assert dynamics is not None
+        assert isinstance(dynamics, CentralBodyGravity)
+        assert isinstance(dynamics, Dynamics)
+        assert dynamics.is_defined()
+
+    def test_getters(self, dynamics: CentralBodyGravity, earth: Earth):
+        assert dynamics.get_celestial() == earth
+
+    def test_compute_contribution(self, dynamics: CentralBodyGravity, state: State):
+        contribution = dynamics.compute_contribution(
+            state.get_instant(), state.get_coordinates(), state.get_frame()
+        )
+
+        assert len(contribution) == 3
+        assert contribution == pytest.approx([-8.134702887755102, 0.0, 0.0])

ostk/astrodynamics/test/dynamics/test_dynamics.py

@@ -0,0 +1,50 @@
+# Apache License 2.0
+
+import pytest
+
+from ostk.physics import Environment
+
+from ostk.astrodynamics import Dynamics
+from ostk.astrodynamics.trajectory.state.coordinates_subset import CartesianPosition
+
+
+@pytest.fixture
+def name() -> str:
+    return "MyDynamics"
+
+
+@pytest.fixture
+def environment() -> Environment:
+    return Environment.default()
+
+
+@pytest.fixture
+def dynamics(name: str) -> Dynamics:
+    class MyDynamics(Dynamics):
+        def is_defined():
+            return True
+
+        def get_read_coordinates_subsets(self):
+            return [CartesianPosition.default()]
+
+        def get_write_coordinates_subsets(self):
+            return [CartesianPosition.default()]
+
+        def compute_contribution(self, instant, state_vector, frame):
+            return state_vector
+
+    return MyDynamics(name)
+
+
+class TestDynamics:
+    def test_subclass(self, dynamics: Dynamics):
+        assert dynamics is not None
+
+    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)