open-space-toolkit-astrodynamics 2.2.0__py311-none-any.whl → 5.0.0__py311-none-any.whl

ostk/astrodynamics/test/flight/system/test_satellite_system.py

@@ -1,77 +0,0 @@
-# Apache License 2.0
-
-import pytest
-
-import numpy as np
-
-import ostk.mathematics as mathematics
-
-import ostk.physics as physics
-
-import ostk.astrodynamics as astrodynamics
-
-Cuboid = mathematics.geometry.d3.objects.Cuboid
-Composite = mathematics.geometry.d3.objects.Composite
-Point = mathematics.geometry.d3.objects.Point
-
-Mass = physics.units.Mass
-
-SatelliteSystem = astrodynamics.flight.system.SatelliteSystem
-
-
-@pytest.fixture
-def satellite_system_default_inputs():
-    mass = Mass(90.0, Mass.Unit.Kilogram)
-    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))
-    surface_area = 0.8
-    drag_coefficient = 2.2
-
-    return (mass, satellite_geometry, inertia_tensor, surface_area, drag_coefficient)
-
-
-@pytest.fixture
-def satellite_system(satellite_system_default_inputs) -> SatelliteSystem:
-    return SatelliteSystem(*satellite_system_default_inputs)
-
-
-class TestSatelliteSystem:
-    def test_constructors(
-        self,
-        satellite_system: SatelliteSystem,
-    ):
-        assert satellite_system is not None
-        assert isinstance(satellite_system, SatelliteSystem)
-        assert satellite_system.is_defined()
-
-    def test_comparators(
-        self,
-        satellite_system: SatelliteSystem,
-    ):
-        assert (satellite_system == satellite_system) is True
-        assert (satellite_system != satellite_system) is False
-
-    def test_getters(
-        self,
-        satellite_system_default_inputs,
-        satellite_system: SatelliteSystem,
-    ):
-        (
-            mass,
-            satellite_geometry,
-            inertia_tensor,
-            surface_area,
-            drag_coefficient,
-        ) = satellite_system_default_inputs
-
-        assert satellite_system.get_mass() == mass
-        assert satellite_system.get_geometry() == satellite_geometry
-        assert np.array_equal(satellite_system.get_inertia_tensor(), inertia_tensor)
-        assert satellite_system.get_cross_sectional_surface_area() == surface_area
-        assert satellite_system.get_drag_coefficient() == drag_coefficient

ostk/astrodynamics/test/flight/test_profile.py

@@ -1,153 +0,0 @@
-# Apache License 2.0
-
-from datetime import datetime
-
-import pytest
-
-from ostk.mathematics.geometry.d3.transformations.rotations import Quaternion
-
-from ostk.physics import Environment
-from ostk.physics.time import DateTime
-from ostk.physics.time import Time
-from ostk.physics.time import Scale
-from ostk.physics.time import Instant
-from ostk.physics.units import Length
-from ostk.physics.coordinate import Transform
-from ostk.physics.coordinate import Position
-from ostk.physics.coordinate import Velocity
-from ostk.physics.coordinate import Frame
-from ostk.physics.coordinate import Axes
-from ostk.physics.coordinate.frame.providers import Dynamic as DynamicProvider
-
-from ostk.astrodynamics import Trajectory
-from ostk.astrodynamics.trajectory import Orbit
-from ostk.astrodynamics.flight import Profile
-from ostk.astrodynamics.flight.profile import State
-from ostk.astrodynamics.flight.profile.models import Transform as TransformModel
-from ostk.astrodynamics.flight.profile.models import Tabulated as TabulatedModel
-
-
-@pytest.fixture
-def instant() -> Instant:
-    return Instant.date_time(DateTime(2020, 1, 3), Scale.UTC)
-
-
-@pytest.fixture
-def profile() -> Profile:
-    def dynamic_provider_generator(instant: Instant):
-        return Transform.identity(instant)
-
-    return Profile(
-        model=TransformModel(
-            dynamic_provider=DynamicProvider(dynamic_provider_generator),
-            frame=Frame.GCRF(),
-        ),
-    )
-
-
-class TestProfile:
-    def test_constructors(self, profile: Profile):
-        assert profile is not None
-        assert isinstance(profile, Profile)
-
-    def test_get_state_at(self, profile: Profile, instant: Instant):
-        state: State = profile.get_state_at(instant)
-
-        assert state is not None
-        assert isinstance(state, State)
-        state.is_defined()
-
-    def test_get_states_at(self, profile: Profile, instant: Instant):
-        states = profile.get_states_at([instant, instant])
-
-        assert states is not None
-
-    def test_get_axes_at(self, profile: Profile, instant: Instant):
-        axes = profile.get_axes_at(instant)
-
-        assert axes is not None
-        assert isinstance(axes, Axes)
-
-    def test_get_body_frame(self, profile: Profile, instant: Instant):
-        frame = profile.get_body_frame("Name")
-
-        assert frame is not None
-        assert isinstance(frame, Frame)
-
-    def test_undefined(self):
-        profile: Profile = Profile.undefined()
-
-        assert profile is not None
-        assert isinstance(profile, Profile)
-        assert profile.is_defined() is False
-
-    def test_inertial_pointing(self):
-        quaternion: Quaternion = Quaternion([0.0, 0.0, 0.0, 1.0], Quaternion.Format.XYZS)
-
-        trajectory: Trajectory = Trajectory.position(
-            Position.meters((0.0, 0.0, 0.0), Frame.GCRF())
-        )
-
-        profile: Profile = Profile.inertial_pointing(trajectory, quaternion)
-
-        assert profile is not None
-        assert isinstance(profile, Profile)
-        assert profile.is_defined()
-
-    def test_nadir_pointing(self):
-        environment = Environment.default()
-
-        orbit = Orbit.sun_synchronous(
-            epoch=Instant.date_time(datetime(2020, 1, 1, 0, 0, 0), Scale.UTC),
-            altitude=Length.kilometers(500.0),
-            local_time_at_descending_node=Time(14, 0, 0),
-            celestial_object=environment.access_celestial_object_with_name("Earth"),
-        )
-
-        profile: Profile = Profile.nadir_pointing(orbit, Orbit.FrameType.VVLH)
-
-        assert profile is not None
-        assert isinstance(profile, Profile)
-        assert profile.is_defined()
-
-    def test_tabulated(self):
-        profile = Profile(
-            model=TabulatedModel(
-                states=[
-                    State(
-                        instant=Instant.date_time(
-                            datetime(2020, 1, 1, 0, 0, 0), Scale.UTC
-                        ),
-                        position=Position.meters((0.0, 0.0, 0.0), Frame.GCRF()),
-                        velocity=Velocity.meters_per_second(
-                            (0.0, 0.0, 0.0), Frame.GCRF()
-                        ),
-                        attitude=Quaternion.unit(),
-                        angular_velocity=(0.0, 0.0, 0.0),
-                        reference_frame=Frame.GCRF(),
-                    ),
-                    State(
-                        instant=Instant.date_time(
-                            datetime(2020, 1, 1, 0, 1, 0), Scale.UTC
-                        ),
-                        position=Position.meters((0.0, 0.0, 0.0), Frame.GCRF()),
-                        velocity=Velocity.meters_per_second(
-                            (0.0, 0.0, 0.0), Frame.GCRF()
-                        ),
-                        attitude=Quaternion.unit(),
-                        angular_velocity=(0.0, 0.0, 0.0),
-                        reference_frame=Frame.GCRF(),
-                    ),
-                ],
-            ),
-        )
-
-        assert isinstance(profile, Profile)
-        assert profile.is_defined()
-
-        assert (
-            profile.get_state_at(
-                Instant.date_time(datetime(2020, 1, 1, 0, 0, 30), Scale.UTC)
-            )
-            is not None
-        )

ostk/astrodynamics/test/flight/test_system.py

@@ -1,55 +0,0 @@
-# Apache License 2.0
-
-import pytest
-
-import numpy as np
-
-import ostk.mathematics as mathematics
-
-import ostk.physics as physics
-
-import ostk.astrodynamics as astrodynamics
-
-Cuboid = mathematics.geometry.d3.objects.Cuboid
-Composite = mathematics.geometry.d3.objects.Composite
-Point = mathematics.geometry.d3.objects.Point
-
-Mass = physics.units.Mass
-
-System = astrodynamics.flight.System
-
-
-@pytest.fixture
-def system_default_inputs():
-    mass = Mass(90.0, Mass.Unit.Kilogram)
-    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],
-        )
-    )
-
-    return (mass, geometry)
-
-
-@pytest.fixture
-def system(system_default_inputs) -> System:
-    return System(*system_default_inputs)
-
-
-class TestSatelliteSystem:
-    def test_constructors(self, system: System):
-        assert system is not None
-        assert isinstance(system, System)
-        assert system.is_defined()
-
-    def test_comparators(self, system: System):
-        assert (system == system) is True
-        assert (system != system) is False
-
-    def test_getters(self, system_default_inputs, system: System):
-        (mass, geometry) = system_default_inputs
-
-        assert system.get_mass() == mass
-        assert system.get_geometry() == geometry

ostk/astrodynamics/test/solvers/__init__.py

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

ostk/astrodynamics/test/solvers/test_temporal_condition_solver.py

@@ -1,153 +0,0 @@
-# Apache License 2.0
-
-import pytest
-
-from ostk.physics.units import Length
-from ostk.physics.units import Angle
-from ostk.physics.time import Instant
-from ostk.physics.time import Duration
-from ostk.physics.time import Interval
-from ostk.physics.time import DateTime
-from ostk.physics.time import Scale
-from ostk.physics import Environment
-
-from ostk.astrodynamics.trajectory import Orbit
-from ostk.astrodynamics.trajectory.orbit.models import Kepler
-from ostk.astrodynamics.trajectory.orbit.models.kepler import COE
-from ostk.astrodynamics.access import Generator
-from ostk.astrodynamics.solvers import TemporalConditionSolver
-
-
-@pytest.fixture
-def temporal_condition_solver() -> TemporalConditionSolver:
-    return TemporalConditionSolver(
-        time_step=Duration.seconds(30.0),
-        tolerance=Duration.milliseconds(1.0),
-        maximum_iteration_count=1234,
-    )
-
-
-@pytest.fixture
-def interval() -> Interval:
-    return Interval.closed(
-        Instant.date_time(DateTime(2018, 1, 1, 0, 0, 0), Scale.UTC),
-        Instant.date_time(DateTime(2018, 1, 1, 2, 0, 0), Scale.UTC),
-    )
-
-
-class TestTemporalConditionSolver:
-    def test_constructor_success(self):
-        temporal_condition_solver = TemporalConditionSolver(
-            time_step=Duration.seconds(30.0),
-            tolerance=Duration.milliseconds(1.0),
-            maximum_iteration_count=1234,
-        )
-
-        assert isinstance(temporal_condition_solver, TemporalConditionSolver)
-
-    def test_getters_success(
-        self,
-        temporal_condition_solver: TemporalConditionSolver,
-    ):
-        assert temporal_condition_solver.get_time_step() == Duration.seconds(30.0)
-        assert temporal_condition_solver.get_tolerance() == Duration.milliseconds(1.0)
-        assert temporal_condition_solver.get_maximum_iteration_count() == 1234
-
-    def test_solve_success_one_condition_always_true(
-        self,
-        temporal_condition_solver: TemporalConditionSolver,
-        interval: Interval,
-    ):
-        solution: list[Interval] = temporal_condition_solver.solve(
-            condition=lambda _: True,
-            interval=interval,
-        )
-
-        assert isinstance(solution, list)
-        assert solution == [interval]
-
-    def test_solve_success_one_condition_always_false(
-        self,
-        temporal_condition_solver: TemporalConditionSolver,
-        interval: Interval,
-    ):
-        solution: list[Interval] = temporal_condition_solver.solve(
-            condition=lambda _: False,
-            interval=interval,
-        )
-
-        assert isinstance(solution, list)
-        assert solution == []
-
-    def test_solve_success_multiple_conditions_always_true(
-        self,
-        temporal_condition_solver: TemporalConditionSolver,
-        interval: Interval,
-    ):
-        solution: list[Interval] = temporal_condition_solver.solve(
-            conditions=[
-                lambda _: True,
-                lambda _: True,
-            ],
-            interval=interval,
-        )
-
-        assert isinstance(solution, list)
-        assert solution == [interval]
-
-    def test_solve_success_multiple_conditions_one_always_false(
-        self,
-        temporal_condition_solver: TemporalConditionSolver,
-        interval: Interval,
-    ):
-        solution: list[Interval] = temporal_condition_solver.solve(
-            conditions=[
-                lambda _: True,
-                lambda _: False,
-            ],
-            interval=interval,
-        )
-
-        assert isinstance(solution, list)
-        assert solution == []
-
-    def test_solve_success_using_access_generator(
-        self,
-        temporal_condition_solver: TemporalConditionSolver,
-        interval: Interval,
-    ):
-        environment = Environment.default()
-
-        generator = Generator(
-            environment=environment,
-        )
-
-        earth = environment.access_celestial_object_with_name("Earth")
-
-        trajectory = Orbit(
-            model=Kepler(
-                coe=COE(
-                    semi_major_axis=Length.kilometers(7000.0),
-                    eccentricity=0.0,
-                    inclination=Angle.degrees(45.0),
-                    raan=Angle.degrees(0.0),
-                    aop=Angle.degrees(0.0),
-                    true_anomaly=Angle.degrees(0.0),
-                ),
-                epoch=Instant.date_time(DateTime(2018, 1, 1, 0, 0, 0), Scale.UTC),
-                celestial_object=earth,
-                perturbation_type=Kepler.PerturbationType.No,
-            ),
-            celestial_object=earth,
-        )
-
-        solution: list[Interval] = temporal_condition_solver.solve(
-            condition=generator.get_condition_function(
-                from_trajectory=trajectory,
-                to_trajectory=trajectory,
-            ),
-            interval=interval,
-        )
-
-        assert isinstance(solution, list)
-        assert solution == [interval]

ostk/astrodynamics/test/trajectory/__init__.py

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

ostk/astrodynamics/test/trajectory/orbit/__init__.py

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

ostk/astrodynamics/test/trajectory/orbit/messages/__init__.py

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

ostk/astrodynamics/test/trajectory/orbit/messages/spacex/__init__.py

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

ostk/astrodynamics/test/trajectory/orbit/messages/spacex/conftest.py

@@ -1,18 +0,0 @@
-# Apache License 2.0
-
-import pytest
-
-import pathlib
-
-from ostk.core.filesystem import Path
-from ostk.core.filesystem import File
-
-
-@pytest.fixture
-def data_directory_path() -> str:
-    return f"{pathlib.Path(__file__).parent.absolute()}/data"
-
-
-@pytest.fixture
-def opm_file(data_directory_path: str) -> File:
-    return File.path(Path.parse(f"{data_directory_path}/opm_1.yaml"))

ostk/astrodynamics/test/trajectory/orbit/messages/spacex/data/opm_1.yaml

@@ -1,44 +0,0 @@
-# Dummy SpaceX OPM output
-
-# Notes:
-# - ECEF velocity is Earth relative
-# - Apogee/Perigee altitude assumes a spherical Earth, 6378.137 km radius
-# - Orbital elements are computed in an inertial frame realized by inertially
-#   freezing the WGS84 ECEF frame at time of current state
-# - State is post-deployment, so includes separation delta-velocity
-
-
-header:
-  generation_date: 2020-01-01T12:34:56.789Z
-  launch_date: 2020-01-02T12:34:56.789Z
-
-
-deployments:
-
-- name: satellite_a
-  sequence_number: 1
-  mission_time_s: 3600.0
-  date: 2020-01-02T13:34:56.789Z
-  r_ecef_m: [693289.644, 6876578.628, -133035.288]
-  v_ecef_m_per_s: [1305.783, 39.783, 7525.920]
-  mean_perigee_altitude_km: 526.768
-  mean_apogee_altitude_km: 568.430
-  mean_inclination_deg: 97.123
-  mean_argument_of_perigee_deg: -179.513
-  mean_longitude_ascending_node_deg: 85.057
-  mean_mean_anomaly_deg: 179.263
-  ballistic_coef_kg_per_m2: 47.55
-
-- name: satellite_b
-  sequence_number: 2
-  mission_time_s: 7200.0
-  date: 2020-01-02T14:34:56.789Z
-  r_ecef_m: [699863.059, 6875647.517, -123777.595]
-  v_ecef_m_per_s: [1504.658, 6.705, 7538.669]
-  mean_perigee_altitude_km: 536.779
-  mean_apogee_altitude_km: 529.851
-  mean_inclination_deg: 97.124
-  mean_argument_of_perigee_deg: 136.875
-  mean_longitude_ascending_node_deg: 85.032
-  mean_mean_anomaly_deg: -127.164
-  ballistic_coef_kg_per_m2: 44.26

ostk/astrodynamics/test/trajectory/orbit/messages/spacex/test_opm.py

@@ -1,108 +0,0 @@
-# Apache License 2.0
-
-import pytest
-
-from datetime import datetime
-
-from ostk.core.containers import Dictionary
-from ostk.core.filesystem import Path
-from ostk.core.filesystem import File
-
-from ostk.physics.units import Length
-from ostk.physics.units import Angle
-from ostk.physics.time import Instant
-from ostk.physics.time import Scale
-from ostk.physics.time import Duration
-from ostk.physics.coordinate import Position
-from ostk.physics.coordinate import Velocity
-from ostk.physics.coordinate import Frame
-
-from ostk.astrodynamics.trajectory.orbit.messages.spacex import OPM
-
-
-@pytest.fixture
-def opm() -> OPM:
-    return OPM(
-        header=OPM.Header(
-            generation_date=Instant.date_time(datetime(2020, 1, 2, 3, 4, 5), Scale.UTC),
-            launch_date=Instant.date_time(datetime(2020, 1, 3, 3, 4, 5), Scale.UTC),
-        ),
-        deployments=[
-            OPM.Deployment(
-                name="A",
-                sequence_number=1,
-                mission_time=Duration.hours(3600.0),
-                date=Instant.date_time(datetime(2020, 1, 4, 3, 4, 5), Scale.UTC),
-                position=Position.meters((1.0, 2.0, 3.0), Frame.ITRF()),
-                velocity=Velocity.meters_per_second((4.0, 5.0, 6.0), Frame.ITRF()),
-                mean_perigee_altitude=Length.kilometers(500.0),
-                mean_apogee_altitude=Length.kilometers(500.0),
-                mean_inclination=Angle.degrees(1.0),
-                mean_argument_of_perigee=Angle.degrees(2.0),
-                mean_longitude_ascending_node=Angle.degrees(3.0),
-                mean_mean_anomaly=Angle.degrees(4.0),
-                ballistic_coefficient=123.456,
-            ),
-        ],
-    )
-
-
-class TestOPM:
-    def test_constructor(self, opm: OPM):
-        assert opm is not None
-
-    def test_get_header(self, opm: OPM):
-        assert opm.get_header() is not None
-
-    def test_get_deployments(self, opm: OPM):
-        assert len(opm.get_deployments()) == 1
-
-    def test_get_deployment_at(self, opm: OPM):
-        assert opm.get_deployment_at(index=0).name == "A"
-
-    def test_get_deployment_with_name(self, opm: OPM):
-        assert opm.get_deployment_with_name(name="A").name == "A"
-
-    def test_undefined(self):
-        assert OPM.undefined().is_defined() is False
-
-    def test_dictionary(self):
-        dictionary = Dictionary(
-            {
-                "header": {
-                    "generation_date": "2020-01-01T12:34:56.789Z",
-                    "launch_date": "2020-01-02T12:34:56.789Z",
-                },
-                "deployments": [
-                    {
-                        "name": "satellite_a",
-                        "sequence_number": 1,
-                        "mission_time_s": 3600.0,
-                        "date": "2020-01-02T13:34:56.789Z",
-                        "r_ecef_m": [693289.644, 6876578.628, -133035.288],
-                        "v_ecef_m_per_s": [1305.783, 39.783, 7525.920],
-                        "mean_perigee_altitude_km": 526.768,
-                        "mean_apogee_altitude_km": 568.430,
-                        "mean_inclination_deg": 97.123,
-                        "mean_argument_of_perigee_deg": -179.513,
-                        "mean_longitude_ascending_node_deg": 85.057,
-                        "mean_mean_anomaly_deg": 179.263,
-                        "ballistic_coef_kg_per_m2": 47.55,
-                    }
-                ],
-            }
-        )
-
-        assert OPM.dictionary(dictionary) is not None
-
-    def test_parse(self, opm_file: File):
-        with open(str(opm_file.get_path().to_string()), "r") as stream:
-            assert OPM.parse(string=stream.read()) is not None
-
-    def test_load(self, opm_file: File):
-        assert OPM.load(file=opm_file) is not None
-
-
-class TestOPMDeployment:
-    def test_to_state(self, opm: OPM):
-        assert opm.get_deployment_with_name(name="A").to_state() is not None

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

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

ostk/astrodynamics/test/trajectory/orbit/models/kepler/__init__.py

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