open-space-toolkit-astrodynamics 9.4.1__py38-none-manylinux2014_x86_64.whl

ostk/astrodynamics/test/trajectory/orbit/models/test_kepler.py

@@ -0,0 +1,130 @@
+# Apache License 2.0
+
+import ostk.physics as physics
+
+import ostk.astrodynamics as astrodynamics
+
+Length = physics.unit.Length
+Angle = physics.unit.Angle
+Scale = physics.time.Scale
+Instant = physics.time.Instant
+Interval = physics.time.Interval
+DateTime = physics.time.DateTime
+Position = physics.coordinate.Position
+Velocity = physics.coordinate.Velocity
+Frame = physics.coordinate.Frame
+Environment = physics.Environment
+
+Trajectory = astrodynamics.Trajectory
+Model = astrodynamics.trajectory.Model
+Orbit = astrodynamics.trajectory.Orbit
+Pass = astrodynamics.trajectory.orbit.Pass
+Kepler = astrodynamics.trajectory.orbit.model.Kepler
+COE = astrodynamics.trajectory.orbit.model.kepler.COE
+SGP4 = astrodynamics.trajectory.orbit.model.SGP4
+TLE = astrodynamics.trajectory.orbit.model.sgp4.TLE
+State = astrodynamics.trajectory.State
+Access = astrodynamics.Access
+
+earth = Environment.default().access_celestial_object_with_name("Earth")
+
+
+def construct_kepler():
+    a = Length.kilometers(7000.0)
+    e = 0.1
+    i = Angle.degrees(35.0)
+    raan = Angle.degrees(40.0)
+    aop = Angle.degrees(50.0)
+    nu = Angle.degrees(60.0)
+
+    coe = COE(a, e, i, raan, aop, nu)
+
+    assert coe is not None
+    assert isinstance(coe, COE)
+
+    epoch = Instant.date_time(DateTime(2018, 1, 1, 0, 0, 0), Scale.UTC)
+
+    kepler = Kepler(coe, epoch, earth, Kepler.PerturbationType.No)
+
+    return kepler
+
+
+def test_trajectory_orbit_models_kepler_constructors():
+    a = Length.kilometers(7000.0)
+    e = 0.1
+    i = Angle.degrees(35.0)
+    raan = Angle.degrees(40.0)
+    aop = Angle.degrees(50.0)
+    nu = Angle.degrees(60.0)
+
+    coe = COE(a, e, i, raan, aop, nu)
+
+    assert coe is not None
+    assert isinstance(coe, COE)
+
+    epoch = Instant.date_time(DateTime(2018, 1, 1, 0, 0, 0), Scale.UTC)
+
+    kepler = Kepler(coe, epoch, earth, Kepler.PerturbationType.No)
+
+    assert kepler is not None
+    assert isinstance(kepler, Kepler)
+    assert kepler.is_defined()
+
+
+def test_trajectory_orbit_models_kepler_comparators():
+    kepler: Kepler = construct_kepler()
+
+    assert kepler == kepler
+    assert (kepler != kepler) is False
+
+
+def test_trajectory_orbit_models_kepler_getters():
+    kepler: Kepler = construct_kepler()
+
+    # get_classical_orbital_elements
+
+    assert kepler.get_classical_orbital_elements() is not None
+
+    # get_epoch
+
+    assert kepler.get_epoch() is not None
+
+    # get_revolution_number_at_epoch()
+
+    assert kepler.get_revolution_number_at_epoch() is not None
+
+    # get_gravitational_parameter()
+
+    assert kepler.get_gravitational_parameter() is not None
+
+    # get_equatorial_radius()
+
+    assert kepler.get_equatorial_radius() is not None
+
+    # get_j2()
+
+    assert kepler.get_j2() is not None
+
+    # get_j4()
+
+    assert kepler.get_j4() is not None
+
+    # get_perturbation_type()
+
+    assert kepler.get_perturbation_type() is not None
+
+
+def test_trajectory_orbit_models_kepler_calculate_state_at_epoch():
+    kepler: Kepler = construct_kepler()
+
+    epoch: Instant = Instant.date_time(DateTime(2018, 1, 1, 0, 0, 0), Scale.UTC)
+
+    assert kepler.calculate_state_at(epoch) is not None
+
+
+def test_trajectory_orbit_models_kepler_calculate_rev_number_at_epoch():
+    kepler: Kepler = construct_kepler()
+
+    epoch: Instant = Instant.date_time(DateTime(2018, 1, 1, 0, 0, 0), Scale.UTC)
+
+    assert kepler.calculate_revolution_number_at(epoch) is not None

ostk/astrodynamics/test/trajectory/orbit/models/test_propagated.py

@@ -0,0 +1,234 @@
+# 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.object.celestial import Earth
+
+from ostk.astrodynamics.trajectory.state import NumericalSolver
+from ostk.astrodynamics.trajectory import State
+from ostk.astrodynamics.trajectory import Orbit
+from ostk.astrodynamics.trajectory import Propagator
+from ostk.astrodynamics import Dynamics
+from ostk.astrodynamics.dynamics import CentralBodyGravity
+from ostk.astrodynamics.dynamics import PositionDerivative
+from ostk.astrodynamics.trajectory.orbit.model import Propagated
+
+
+@pytest.fixture
+def state() -> State:
+    frame: Frame = Frame.GCRF()
+    position: Position = Position.meters([7500000.0, 0.0, 0.0], frame)
+    velocity: Velocity = Velocity.meters_per_second(
+        [0.0, 5335.865450622126, 5335.865450622126], frame
+    )
+
+    instant: Instant = Instant.date_time(DateTime(2018, 1, 1, 0, 0, 0), Scale.UTC)
+    return State(instant, position, velocity)
+
+
+@pytest.fixture
+def central_body_dynamics() -> CentralBodyGravity:
+    return CentralBodyGravity(Earth.WGS84(20, 0))
+
+
+@pytest.fixture
+def position_derivative() -> PositionDerivative:
+    return PositionDerivative()
+
+
+@pytest.fixture
+def dynamics(
+    position_derivative: PositionDerivative, central_body_dynamics: CentralBodyGravity
+) -> list:
+    return [position_derivative, central_body_dynamics]
+
+
+@pytest.fixture
+def numerical_solver() -> NumericalSolver:
+    return NumericalSolver(
+        NumericalSolver.LogType.NoLog,
+        NumericalSolver.StepperType.RungeKuttaFehlberg78,
+        5.0,
+        1.0e-15,
+        1.0e-15,
+    )
+
+
+@pytest.fixture
+def propagator(numerical_solver: NumericalSolver, dynamics: list[Dynamics]) -> Propagator:
+    return Propagator(numerical_solver, dynamics)
+
+
+@pytest.fixture
+def revolution_number() -> int:
+    return 5
+
+
+@pytest.fixture
+def propagated(
+    propagator: Propagator, state: State, revolution_number: int
+) -> Propagated:
+    return Propagated(propagator, state, revolution_number)
+
+
+@pytest.fixture
+def earth() -> Earth:
+    return Earth.spherical()
+
+
+@pytest.fixture
+def orbit(propagated: Propagated, earth: Earth) -> Orbit:
+    return Orbit(propagated, earth)
+
+
+class TestPropagated:
+    def test_constructors(
+        self,
+        propagated: Propagated,
+        earth: Earth,
+        dynamics: list,
+        numerical_solver: NumericalSolver,
+        state: State,
+    ):
+        assert propagated is not None
+        assert isinstance(propagated, Propagated)
+        assert propagated.is_defined()
+
+        orbit = Orbit(propagated, earth)
+
+        assert orbit is not None
+        assert isinstance(orbit, Orbit)
+        assert orbit.is_defined()
+
+        state_array = [state, state]
+        propagated_with_state_array = Propagated(
+            Propagator(numerical_solver, dynamics), state_array
+        )
+
+        assert propagated_with_state_array is not None
+        assert isinstance(propagated_with_state_array, Propagated)
+        assert propagated_with_state_array.is_defined()
+
+    def test_comparators(
+        self,
+        propagated: Propagated,
+    ):
+        assert (propagated == propagated) is True
+        assert (propagated != propagated) is False
+
+    def test_getters(
+        self,
+        propagated: Propagated,
+        state: State,
+        revolution_number: int,
+    ):
+        assert propagated.get_epoch() == state.get_instant()
+
+        assert propagated.get_revolution_number_at_epoch() == revolution_number
+
+        assert propagated.access_cached_state_array() == [state]
+
+    def test_calculate_state_at(
+        self,
+        propagated: Propagated,
+        orbit: Orbit,
+    ):
+        instant: Instant = Instant.date_time(DateTime(2018, 1, 1, 0, 10, 0), Scale.UTC)
+
+        propagated_state = propagated.calculate_state_at(instant)
+        propagated_state_orbit = orbit.get_state_at(instant)
+
+        assert propagated_state == propagated_state_orbit
+
+        propagated_state_position_ref = np.array(
+            [6265892.25765909, 3024770.94961259, 3024359.72137468]
+        )
+        propagated_state_velocity_ref = np.array(
+            [-3974.49168221, 4468.16996776, 4466.19232746]
+        )
+
+        propagated_state_position = propagated_state.get_position().get_coordinates()
+        propagated_state_velocity = propagated_state.get_velocity().get_coordinates()
+
+        assert all(
+            [
+                round(propagated_state_position[i], 8)
+                == round(propagated_state_position_ref[i], 8)
+                for i in range(0, len(propagated_state_position_ref))
+            ]
+        )
+        assert all(
+            [
+                round(propagated_state_velocity[i], 8)
+                == round(propagated_state_velocity_ref[i], 8)
+                for i in range(0, len(propagated_state_velocity_ref))
+            ]
+        )
+        assert propagated_state.get_instant() == instant
+
+    def test_calculate_states_at(
+        self,
+        propagated: Propagated,
+        orbit: Orbit,
+    ):
+        instant_array = [
+            Instant.date_time(DateTime(2018, 1, 1, 0, 10, 0), Scale.UTC),
+            Instant.date_time(DateTime(2018, 1, 1, 0, 20, 0), Scale.UTC),
+        ]
+
+        propagated_state_array = propagated.calculate_states_at(instant_array)
+        propagated_state_array_orbit = orbit.get_states_at(instant_array)
+
+        assert propagated_state_array_orbit == propagated_state_array
+
+        assert propagated_state_array_orbit[0].get_instant() == instant_array[0]
+        assert propagated_state_array_orbit[1].get_instant() == instant_array[1]
+
+    def test_calculate_revolution_number_at(
+        self,
+        propagated: Propagated,
+        orbit: Orbit,
+        revolution_number: int,
+    ):
+        instant: Instant = Instant.date_time(DateTime(2018, 1, 1, 0, 40, 0), Scale.UTC)
+
+        assert propagated.calculate_revolution_number_at(instant) == revolution_number + 1
+        assert orbit.get_revolution_number_at(instant) == revolution_number + 1
+
+    def test_access_cached_state_array(
+        self,
+        propagated: Propagated,
+        state: State,
+    ):
+        assert len(propagated.access_cached_state_array()) == 1
+        assert propagated.access_cached_state_array()[0] == state
+
+    def test_access_propagator(
+        self,
+        propagated: Propagated,
+    ):
+        assert propagated.access_propagator() is not None
+        assert isinstance(propagated.access_propagator(), Propagator)
+
+    def test_set_cached_state_array(
+        self,
+        propagated: Propagated,
+        state: State,
+    ):
+        assert len(propagated.access_cached_state_array()) == 1
+
+        propagated.set_cached_state_array([state, state, state])
+
+        assert len(propagated.access_cached_state_array()) == 1
+        assert propagated.access_cached_state_array()[0] == state
+
+        with pytest.raises(Exception) as e:
+            propagated.set_cached_state_array([])

ostk/astrodynamics/test/trajectory/orbit/models/test_sgp4.py

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