open-space-toolkit-mathematics 4.6.0__py39-none-manylinux2014_aarch64.whl

ostk/mathematics/test/object/test_interval.py

@@ -0,0 +1,515 @@
+# Apache License 2.0
+
+import pytest
+
+from ostk.core.type import Real
+from ostk.core.type import String
+
+from ostk.mathematics.object import RealInterval
+
+
+class TestInterval:
+    def test_type(self):
+        enum_members = RealInterval.Type.__members__
+
+        assert list(enum_members.keys()) == [
+            "Undefined",
+            "Closed",
+            "Open",
+            "HalfOpenLeft",
+            "HalfOpenRight",
+        ]
+        assert list(enum_members.values()) == [
+            RealInterval.Type.Undefined,
+            RealInterval.Type.Closed,
+            RealInterval.Type.Open,
+            RealInterval.Type.HalfOpenLeft,
+            RealInterval.Type.HalfOpenRight,
+        ]
+
+    def test_default_constructor(self):
+        # Input types for RealInterval
+        interval_1 = RealInterval(-4.31, 1.0, RealInterval.Type.Open)
+        interval_2 = RealInterval(-2.0, -1.0, RealInterval.Type.Closed)
+        interval_3 = RealInterval(3.5, 4567.35566, RealInterval.Type.HalfOpenRight)
+        interval_4 = RealInterval(1.45, 1.45, RealInterval.Type.Closed)
+
+        assert isinstance(interval_1, RealInterval)
+        assert isinstance(interval_2, RealInterval)
+        assert isinstance(interval_3, RealInterval)
+        assert isinstance(interval_4, RealInterval)
+        assert interval_1 is not None
+        assert interval_2 is not None
+        assert interval_3 is not None
+        assert interval_4 is not None
+
+        with pytest.raises(TypeError):
+            interval = RealInterval(3.0, 1, RealInterval.Type.Closed)
+
+        interval_5 = RealInterval(Real(-4.31), Real(1.0), RealInterval.Type.Open)
+        interval_6 = RealInterval(Real(-2.0), Real(-1.0), RealInterval.Type.Closed)
+        interval_7 = RealInterval(
+            Real(3.5), Real(4567.35566), RealInterval.Type.HalfOpenRight
+        )
+        interval_8 = RealInterval(Real(1.45), Real(1.45), RealInterval.Type.Closed)
+
+        assert isinstance(interval_5, RealInterval)
+        assert isinstance(interval_6, RealInterval)
+        assert isinstance(interval_7, RealInterval)
+        assert isinstance(interval_8, RealInterval)
+        assert interval_5 is not None
+        assert interval_6 is not None
+        assert interval_7 is not None
+        assert interval_8 is not None
+
+        # Interval Bounds
+        a = -4.31
+        b = 3.0
+
+        # Types of RealInterval
+        interval_undefined = RealInterval(a, b, RealInterval.Type.Undefined)
+        interval_closed = RealInterval(a, b, RealInterval.Type.Closed)
+        interval_open = RealInterval(a, b, RealInterval.Type.Open)
+        interval_halfopenleft = RealInterval(a, b, RealInterval.Type.HalfOpenLeft)
+        interval_halfopenright = RealInterval(a, b, RealInterval.Type.HalfOpenRight)
+
+        assert isinstance(interval_undefined, RealInterval)
+        assert isinstance(interval_closed, RealInterval)
+        assert isinstance(interval_open, RealInterval)
+        assert isinstance(interval_halfopenleft, RealInterval)
+        assert isinstance(interval_halfopenright, RealInterval)
+
+        # Unvalid interval definition
+
+        with pytest.raises(RuntimeError):
+            invalid_interval_1 = RealInterval(4.8, 3.5, RealInterval.Type.Open)
+
+        with pytest.raises(RuntimeError):
+            invalid_interval_2 = RealInterval(4.8, 3.5, RealInterval.Type.Closed)
+
+        with pytest.raises(RuntimeError):
+            invalid_interval_3 = RealInterval(4.8, 3.5, RealInterval.Type.HalfOpenLeft)
+
+        with pytest.raises(RuntimeError):
+            invalid_interval_4 = RealInterval(4.8, 3.5, RealInterval.Type.HalfOpenRight)
+
+    def test_undefined_constructor(self):
+        undefined_interval = RealInterval.undefined()
+
+        assert isinstance(undefined_interval, RealInterval)
+        assert undefined_interval is not None
+        assert undefined_interval.is_defined() is False
+
+    @pytest.mark.parametrize(
+        "static_func",
+        [
+            (RealInterval.closed),
+            (RealInterval.open),
+            (RealInterval.half_open_left),
+            (RealInterval.half_open_right),
+        ],
+    )
+    def test_static_constructors(self, static_func):
+        a = -3.1
+        b = 45.6
+
+        closed_interval = static_func(a, b)
+
+        assert isinstance(closed_interval, RealInterval)
+        assert closed_interval is not None
+        assert closed_interval.is_defined()
+        assert closed_interval.get_lower_bound() == a
+        assert closed_interval.get_upper_bound() == b
+
+    def test_is_defined(self):
+        # Interval Bounds
+        a = -4.31
+        b = 3.0
+
+        # Types of RealInterval
+        interval_undefined = RealInterval(a, b, RealInterval.Type.Undefined)
+        interval_closed = RealInterval(a, b, RealInterval.Type.Closed)
+        interval_open = RealInterval(a, b, RealInterval.Type.Open)
+        interval_halfopenleft = RealInterval(a, b, RealInterval.Type.HalfOpenLeft)
+        interval_halfopenright = RealInterval(a, b, RealInterval.Type.HalfOpenRight)
+
+        assert interval_undefined.is_defined() is False
+        assert interval_closed.is_defined() is True
+        assert interval_open.is_defined() is True
+        assert interval_halfopenleft.is_defined() is True
+        assert interval_halfopenright.is_defined() is True
+        assert RealInterval.undefined().is_defined() is False
+
+    # def test_is_degenerate ():
+
+    def test_intersects(self):
+        # Interval Bounds
+        a = -4.31
+        b = 3.0
+
+        # Types of RealInterval
+        interval_undefined = RealInterval(a, b, RealInterval.Type.Undefined)
+        interval_closed = RealInterval(a, b, RealInterval.Type.Closed)
+        interval_open = RealInterval(a, b, RealInterval.Type.Open)
+        interval_halfopenleft = RealInterval(a, b, RealInterval.Type.HalfOpenLeft)
+        interval_halfopenright = RealInterval(a, b, RealInterval.Type.HalfOpenRight)
+
+        # Define Test Intervals
+        interval_left = RealInterval(-5.0, -4.5, RealInterval.Type.Closed)
+        interval_intersects_left = RealInterval(-5.0, -4.26, RealInterval.Type.Closed)
+        interval_right = RealInterval(4.56, 4.67, RealInterval.Type.Closed)
+        interval_intersects_right = RealInterval(2.78, 46.09, RealInterval.Type.Closed)
+        interval_between = RealInterval(-3.4, 2.45, RealInterval.Type.Closed)
+        interval_bigger = RealInterval(-45.0, 34.12, RealInterval.Type.Closed)
+
+        # Add test cases with contained intervals not RealInterval.Type.Closed...
+        # Add test cases with open intervals and half open intervals...
+
+        # Test intersects on undefined
+
+        # with pytest.raises(RuntimeError):
+
+        #     interval_undefined.intersects
+
+        # Test intersects on closed
+
+        assert interval_closed.intersects(interval_left) is False
+        assert interval_closed.intersects(interval_intersects_left) is True
+        assert interval_closed.intersects(interval_right) is False
+        assert interval_closed.intersects(interval_intersects_right) is True
+        assert interval_closed.intersects(interval_between) is True
+        assert interval_closed.intersects(interval_bigger) is True
+
+        # Test intersects on open
+
+        assert interval_open.intersects(interval_left) is False
+        assert interval_open.intersects(interval_intersects_left) is True
+        assert interval_open.intersects(interval_right) is False
+        assert interval_open.intersects(interval_intersects_right) is True
+        assert interval_open.intersects(interval_between) is True
+        assert interval_open.intersects(interval_bigger) is True
+
+        # Test intersects on halfopenleft
+
+        assert interval_halfopenleft.intersects(interval_left) is False
+        assert interval_halfopenleft.intersects(interval_intersects_left) is True
+        assert interval_halfopenleft.intersects(interval_right) is False
+        assert interval_halfopenleft.intersects(interval_intersects_right) is True
+        assert interval_halfopenleft.intersects(interval_between) is True
+        assert interval_halfopenleft.intersects(interval_bigger) is True
+
+        # Test intersects on halfopenright
+
+        assert interval_halfopenright.intersects(interval_left) is False
+        assert interval_halfopenright.intersects(interval_intersects_left) is True
+        assert interval_halfopenright.intersects(interval_right) is False
+        assert interval_halfopenright.intersects(interval_intersects_right) is True
+        assert interval_halfopenright.intersects(interval_between) is True
+        assert interval_halfopenright.intersects(interval_bigger) is True
+
+    def test_contains_real(self):
+        # Interval Bounds
+        a = -4.31
+        b = 3.0
+
+        # Types of RealInterval
+        interval_undefined = RealInterval(a, b, RealInterval.Type.Undefined)
+        interval_closed = RealInterval(a, b, RealInterval.Type.Closed)
+        interval_open = RealInterval(a, b, RealInterval.Type.Open)
+        interval_halfopenleft = RealInterval(a, b, RealInterval.Type.HalfOpenLeft)
+        interval_halfopenright = RealInterval(a, b, RealInterval.Type.HalfOpenRight)
+
+        # Define Reals
+        real_left = -5.43
+        real_right = 3.1
+        real_leftbound = a
+        real_rightbound = b
+        real_between_1 = 2.3
+        real_between_2 = 2.999999999
+
+        # Test contains on undefined
+
+        with pytest.raises(RuntimeError):
+            interval_undefined.contains(real_left)
+
+        # Test contains on closed
+
+        assert interval_closed.contains(real_left) is False
+        assert interval_closed.contains(real_right) is False
+        assert interval_closed.contains(real_leftbound) is True
+        assert interval_closed.contains(real_rightbound) is True
+        assert interval_closed.contains(real_between_1) is True
+        assert interval_closed.contains(real_between_2) is True
+
+        # Test contains on open
+
+        assert interval_open.contains(real_left) is False
+        assert interval_open.contains(real_right) is False
+        assert interval_open.contains(real_leftbound) is False
+        assert interval_open.contains(real_rightbound) is False
+        assert interval_open.contains(real_between_1) is True
+        assert interval_open.contains(real_between_2) is True
+
+        # Test contains on halfopenleft
+
+        assert interval_halfopenleft.contains(real_left) is False
+        assert interval_halfopenleft.contains(real_right) is False
+        assert interval_halfopenleft.contains(real_leftbound) is False
+        assert interval_halfopenleft.contains(real_rightbound) is True
+        assert interval_halfopenleft.contains(real_between_1) is True
+        assert interval_halfopenleft.contains(real_between_2) is True
+
+        # Test contains on halfopenright
+
+        assert interval_halfopenright.contains(real_left) is False
+        assert interval_halfopenright.contains(real_right) is False
+        assert interval_halfopenright.contains(real_leftbound) is True
+        assert interval_halfopenright.contains(real_rightbound) is False
+        assert interval_halfopenright.contains(real_between_1) is True
+        assert interval_halfopenright.contains(real_between_2) is True
+
+    def test_contains_interval(self):
+        # Interval Bounds
+        a = -4.31
+        b = 3.0
+
+        # Types of RealInterval
+        interval_undefined = RealInterval(a, b, RealInterval.Type.Undefined)
+        interval_closed = RealInterval(a, b, RealInterval.Type.Closed)
+        interval_open = RealInterval(a, b, RealInterval.Type.Open)
+        interval_halfopenleft = RealInterval(a, b, RealInterval.Type.HalfOpenLeft)
+        interval_halfopenright = RealInterval(a, b, RealInterval.Type.HalfOpenRight)
+
+        # Define Test Intervals
+        interval_left = RealInterval(-5.0, -4.5, RealInterval.Type.Closed)
+        interval_intersects_left = RealInterval(-5.0, -4.56, RealInterval.Type.Closed)
+        interval_right = RealInterval(4.56, 4.67, RealInterval.Type.Closed)
+        interval_intersects_right = RealInterval(2.78, 46.09, RealInterval.Type.Closed)
+        interval_between = RealInterval(-3.4, 2.45, RealInterval.Type.Closed)
+        interval_bigger = RealInterval(-45.0, 34.12, RealInterval.Type.Closed)
+
+        # Add test cases with contained intervals not RealInterval.Type.Closed...
+
+        # Test contains on undefined
+
+        with pytest.raises(RuntimeError):
+            interval_undefined.contains(interval_left)
+
+        # Test contains on closed
+
+        assert interval_closed.contains(interval_left) is False
+        assert interval_closed.contains(interval_intersects_left) is False
+        assert interval_closed.contains(interval_right) is False
+        assert interval_closed.contains(interval_intersects_right) is False
+        assert interval_closed.contains(interval_between) is True
+        assert interval_closed.contains(interval_bigger) is False
+
+        # Test contains on open
+
+        assert interval_open.contains(interval_left) is False
+        assert interval_open.contains(interval_intersects_left) is False
+        assert interval_open.contains(interval_right) is False
+        assert interval_open.contains(interval_intersects_right) is False
+        assert interval_open.contains(interval_between) is True
+        assert interval_open.contains(interval_bigger) is False
+
+        # Test contains on halfopenleft
+
+        assert interval_halfopenleft.contains(interval_left) is False
+        assert interval_halfopenleft.contains(interval_intersects_left) is False
+        assert interval_halfopenleft.contains(interval_right) is False
+        assert interval_halfopenleft.contains(interval_intersects_right) is False
+        assert interval_halfopenleft.contains(interval_between) is True
+        assert interval_halfopenleft.contains(interval_bigger) is False
+
+        # Test contains on halfopenright
+
+        assert interval_halfopenright.contains(interval_left) is False
+        assert interval_halfopenright.contains(interval_intersects_left) is False
+        assert interval_halfopenright.contains(interval_right) is False
+        assert interval_halfopenright.contains(interval_intersects_right) is False
+        assert interval_halfopenright.contains(interval_between) is True
+        assert interval_halfopenright.contains(interval_bigger) is False
+
+    def test_get_bounds(self):
+        # Interval Bounds
+        a = -4.31
+        b = 3.0
+
+        # Types of RealInterval
+        interval_undefined = RealInterval(a, b, RealInterval.Type.Undefined)
+        interval_closed = RealInterval(a, b, RealInterval.Type.Closed)
+        interval_open = RealInterval(a, b, RealInterval.Type.Open)
+        interval_halfopenleft = RealInterval(a, b, RealInterval.Type.HalfOpenLeft)
+        interval_halfopenright = RealInterval(a, b, RealInterval.Type.HalfOpenRight)
+
+        # get_lower_bound
+        with pytest.raises(RuntimeError):
+            interval_undefined.get_lower_bound()
+
+        assert interval_closed.get_lower_bound() == a
+        assert interval_open.get_lower_bound() == a
+        assert interval_halfopenleft.get_lower_bound() == a
+        assert interval_halfopenright.get_lower_bound() == a
+
+        # get_upper_bound
+        with pytest.raises(RuntimeError):
+            interval_undefined.get_upper_bound()
+
+        assert interval_closed.get_upper_bound() == b
+        assert interval_open.get_upper_bound() == b
+        assert interval_halfopenleft.get_upper_bound() == b
+        assert interval_halfopenright.get_upper_bound() == b
+
+    def test_to_string(self):
+        # Interval Bounds
+        a = -4.31
+        b = 3.0
+
+        # Types of RealInterval
+        interval_undefined = RealInterval(a, b, RealInterval.Type.Undefined)
+        interval_closed = RealInterval(a, b, RealInterval.Type.Closed)
+        interval_open = RealInterval(a, b, RealInterval.Type.Open)
+        interval_halfopenleft = RealInterval(a, b, RealInterval.Type.HalfOpenLeft)
+        interval_halfopenright = RealInterval(a, b, RealInterval.Type.HalfOpenRight)
+
+        with pytest.raises(RuntimeError):
+            interval_undefined.to_string()
+
+        assert isinstance(interval_closed.to_string(), String)
+        assert interval_closed.to_string() == "[-4.3099999999999996, 3.0]"
+        assert interval_open.to_string() == "]-4.3099999999999996, 3.0["
+        assert interval_halfopenleft.to_string() == "]-4.3099999999999996, 3.0]"
+        assert interval_halfopenright.to_string() == "[-4.3099999999999996, 3.0["
+
+    def test_get_intersection_with(self):
+        test_cases = [
+            (
+                RealInterval(0.0, 10.0, RealInterval.Type.HalfOpenLeft),
+                RealInterval(5.0, 7.0, RealInterval.Type.HalfOpenLeft),
+                RealInterval(5.0, 7.0, RealInterval.Type.HalfOpenLeft),
+            ),
+            (
+                RealInterval(0.0, 10.0, RealInterval.Type.HalfOpenRight),
+                RealInterval(-15.0, 25.0, RealInterval.Type.HalfOpenRight),
+                RealInterval(0.0, 10.0, RealInterval.Type.HalfOpenRight),
+            ),
+            (
+                RealInterval(0.0, 10.0, RealInterval.Type.Open),
+                RealInterval(-5.0, 7.0, RealInterval.Type.Open),
+                RealInterval(0.0, 7.0, RealInterval.Type.Open),
+            ),
+            (
+                RealInterval(0.0, 10.0, RealInterval.Type.Closed),
+                RealInterval(5.0, 15.0, RealInterval.Type.Closed),
+                RealInterval(5.0, 10.0, RealInterval.Type.Closed),
+            ),
+        ]
+
+        for test_case in test_cases:
+            first_interval = test_case[0]
+            second_interval = test_case[1]
+            expected_interval = test_case[2]
+
+            assert (
+                first_interval.get_intersection_with(second_interval) == expected_interval
+            )
+
+    def test_get_union_with(self):
+        test_cases = [
+            (
+                RealInterval(0.0, 10.0, RealInterval.Type.HalfOpenLeft),
+                RealInterval(5.0, 7.0, RealInterval.Type.HalfOpenLeft),
+                RealInterval(0.0, 10.0, RealInterval.Type.HalfOpenLeft),
+            ),
+            (
+                RealInterval(0.0, 10.0, RealInterval.Type.HalfOpenRight),
+                RealInterval(-15.0, 25.0, RealInterval.Type.HalfOpenRight),
+                RealInterval(-15.0, 25.0, RealInterval.Type.HalfOpenRight),
+            ),
+            (
+                RealInterval(0.0, 10.0, RealInterval.Type.Open),
+                RealInterval(-5.0, 7.0, RealInterval.Type.Open),
+                RealInterval(-5.0, 10.0, RealInterval.Type.Open),
+            ),
+            (
+                RealInterval(0.0, 10.0, RealInterval.Type.Closed),
+                RealInterval(5.0, 15.0, RealInterval.Type.Closed),
+                RealInterval(0.0, 15.0, RealInterval.Type.Closed),
+            ),
+        ]
+
+        for test_case in test_cases:
+            first_interval = test_case[0]
+            second_interval = test_case[1]
+            expected_interval = test_case[2]
+
+            assert first_interval.get_union_with(second_interval) == expected_interval
+
+    def test_clip(self):
+        intervals: list[RealInterval] = [
+            RealInterval(0.0, 10.0, RealInterval.Type.HalfOpenLeft),
+            RealInterval(6.5, 8.0, RealInterval.Type.HalfOpenLeft),
+        ]
+        clipping_interval: RealInterval = RealInterval.closed(6.0, 7.0)
+
+        assert RealInterval.clip(intervals, clipping_interval) is not None
+
+    def test_sort(self):
+        intervals: list[RealInterval] = [
+            RealInterval(0.0, 10.0, RealInterval.Type.HalfOpenLeft),
+            RealInterval(10.0, 20.0, RealInterval.Type.HalfOpenLeft),
+        ]
+        assert RealInterval.sort(intervals) is not None
+        assert RealInterval.sort(intervals, False) is not None
+        assert RealInterval.sort(intervals, False, False) is not None
+
+    def test_merge(self):
+        intervals: list[RealInterval] = [
+            RealInterval(0.0, 10.0, RealInterval.Type.HalfOpenLeft),
+            RealInterval(5.0, 7.0, RealInterval.Type.HalfOpenLeft),
+            RealInterval(0.0, 10.0, RealInterval.Type.HalfOpenLeft),
+        ]
+
+        assert RealInterval.merge(intervals) is not None
+
+    def test_get_gaps(self):
+        intervals: list[RealInterval] = [
+            RealInterval(0.0, 10.0, RealInterval.Type.HalfOpenLeft),
+            RealInterval(5.0, 7.0, RealInterval.Type.HalfOpenLeft),
+            RealInterval(5.0, 15.0, RealInterval.Type.HalfOpenLeft),
+        ]
+
+        assert RealInterval.get_gaps(intervals) is not None
+
+        assert (
+            RealInterval.get_gaps(intervals, RealInterval.closed(-5.0, 25.0)) is not None
+        )
+
+    def test_logical_or(self):
+        intervals_1: list[RealInterval] = [
+            RealInterval(8.0, 9.0, RealInterval.Type.Closed),
+            RealInterval(0.0, 1.0, RealInterval.Type.HalfOpenLeft),
+            RealInterval(2.0, 3.0, RealInterval.Type.Open),
+        ]
+
+        intervals_2: list[RealInterval] = [
+            RealInterval(0.5, 3.5, RealInterval.Type.Open),
+            RealInterval(5.0, 7.0, RealInterval.Type.HalfOpenLeft),
+        ]
+
+        assert RealInterval.logical_or(intervals_1, intervals_2) is not None
+
+    def test_logical_and(self):
+        intervals_1: list[RealInterval] = [
+            RealInterval(8.0, 9.0, RealInterval.Type.Closed),
+            RealInterval(0.0, 1.0, RealInterval.Type.HalfOpenLeft),
+            RealInterval(2.0, 3.0, RealInterval.Type.Open),
+        ]
+
+        intervals_2: list[RealInterval] = [
+            RealInterval(0.5, 3.5, RealInterval.Type.Open),
+            RealInterval(5.0, 7.0, RealInterval.Type.HalfOpenLeft),
+        ]
+
+        assert RealInterval.logical_and(intervals_1, intervals_2) is not None

ostk/mathematics/test/object/test_vector.py

@@ -0,0 +1,5 @@
+# Apache License 2.0
+
+import pytest
+
+import ostk.mathematics as mathematics

ostk/mathematics/test/solver/test_numerical_solver.py

@@ -0,0 +1,176 @@
+# Apache License 2.0
+
+import pytest
+
+import numpy as np
+import math
+
+from ostk.mathematics.solver import NumericalSolver
+
+
+def oscillator(x, dxdt, _):
+    dxdt[0] = x[1]
+    dxdt[1] = -x[0]
+    return dxdt
+
+
+def get_state_vec(time: float) -> np.ndarray:
+    return np.array([math.sin(time), math.cos(time)])
+
+
+@pytest.fixture
+def initial_state_vec() -> np.ndarray:
+    return get_state_vec(0.0)
+
+
+@pytest.fixture
+def numerical_solver_default_inputs() -> (
+    tuple[NumericalSolver.LogType, NumericalSolver.StepperType, float, float, float]
+):
+    log_type = NumericalSolver.LogType.NoLog
+    stepper_type = NumericalSolver.StepperType.RungeKuttaCashKarp54
+    initial_time_step = 5.0
+    relative_tolerance = 1.0e-12
+    absolute_tolerance = 1.0e-12
+
+    return (
+        log_type,
+        stepper_type,
+        initial_time_step,
+        relative_tolerance,
+        absolute_tolerance,
+    )
+
+
+@pytest.fixture
+def numerical_solver(numerical_solver_default_inputs) -> NumericalSolver:
+    return NumericalSolver(*numerical_solver_default_inputs)
+
+
+@pytest.fixture
+def numerical_solver_conditional() -> NumericalSolver:
+    return NumericalSolver(
+        NumericalSolver.LogType.NoLog,
+        NumericalSolver.StepperType.RungeKuttaDopri5,
+        5.0,
+        1.0e-12,
+        1.0e-12,
+    )
+
+
+class TestNumericalSolver:
+    def test_constructor_success(self, numerical_solver: NumericalSolver):
+        assert numerical_solver is not None
+        assert isinstance(numerical_solver, NumericalSolver)
+        assert numerical_solver.is_defined()
+
+    def test_comparators(self, numerical_solver: NumericalSolver):
+        assert numerical_solver == numerical_solver
+        assert (numerical_solver != numerical_solver) is False
+
+    def test_get_types(
+        self,
+        numerical_solver_default_inputs: tuple[
+            NumericalSolver.LogType, NumericalSolver.StepperType, float, float, float
+        ],
+        numerical_solver: NumericalSolver,
+    ):
+        (
+            log_type,
+            stepper_type,
+            initial_time_step,
+            relative_tolerance,
+            absolute_tolerance,
+        ) = numerical_solver_default_inputs
+
+        assert numerical_solver.get_log_type() == log_type
+        assert numerical_solver.get_stepper_type() == stepper_type
+        assert numerical_solver.get_time_step() == initial_time_step
+        assert numerical_solver.get_relative_tolerance() == relative_tolerance
+        assert numerical_solver.get_absolute_tolerance() == absolute_tolerance
+        assert numerical_solver.get_observed_state_vectors() is not None
+
+    def test_get_string_from_types(self):
+        assert (
+            NumericalSolver.string_from_stepper_type(
+                NumericalSolver.StepperType.RungeKutta4
+            )
+            == "RungeKutta4"
+        )
+        assert (
+            NumericalSolver.string_from_stepper_type(
+                NumericalSolver.StepperType.RungeKuttaCashKarp54
+            )
+            == "RungeKuttaCashKarp54"
+        )
+        assert (
+            NumericalSolver.string_from_stepper_type(
+                NumericalSolver.StepperType.RungeKuttaFehlberg78
+            )
+            == "RungeKuttaFehlberg78"
+        )
+        assert (
+            NumericalSolver.string_from_log_type(NumericalSolver.LogType.NoLog) == "NoLog"
+        )
+        assert (
+            NumericalSolver.string_from_log_type(NumericalSolver.LogType.LogConstant)
+            == "LogConstant"
+        )
+        assert (
+            NumericalSolver.string_from_log_type(NumericalSolver.LogType.LogAdaptive)
+            == "LogAdaptive"
+        )
+
+    def test_integrate_duration(
+        self, numerical_solver: NumericalSolver, initial_state_vec: np.ndarray
+    ):
+        integration_duration: float = 100.0
+
+        state_vector, _ = numerical_solver.integrate_duration(
+            initial_state_vec, integration_duration, oscillator
+        )
+
+        assert 5e-9 >= abs(state_vector[0] - math.sin(integration_duration))
+        assert 5e-9 >= abs(state_vector[1] - math.cos(integration_duration))
+
+        integration_durations = np.arange(100.0, 200.0, 50.0).astype(float)
+        solutions = numerical_solver.integrate_duration(
+            initial_state_vec, integration_durations, oscillator
+        )
+
+        for solution, integration_duration in zip(solutions, integration_durations):
+            state_vector, _ = solution
+
+            assert 5e-9 >= abs(state_vector[0] - math.sin(integration_duration))
+            assert 5e-9 >= abs(state_vector[1] - math.cos(integration_duration))
+
+    def test_integrate_time(self, numerical_solver: NumericalSolver):
+        start_time: float = 500.0
+        end_time: float = start_time + 100.0
+
+        initial_state_vec = get_state_vec(start_time)
+
+        state_vector, _ = numerical_solver.integrate_time(
+            initial_state_vec, start_time, end_time, oscillator
+        )
+
+        assert 5e-9 >= abs(state_vector[0] - math.sin(end_time))
+        assert 5e-9 >= abs(state_vector[1] - math.cos(end_time))
+
+        end_times = np.arange(600.0, 700.0, 50.0).astype(float)
+        solutions = numerical_solver.integrate_time(
+            initial_state_vec, start_time, end_times, oscillator
+        )
+
+        for solution, end_time in zip(solutions, end_times):
+            state_vector, _ = solution
+
+            assert 5e-9 >= abs(state_vector[0] - math.sin(end_time))
+            assert 5e-9 >= abs(state_vector[1] - math.cos(end_time))
+
+    def test_default(self):
+        assert NumericalSolver.default() is not None
+
+    def test_undefined(self):
+        assert NumericalSolver.undefined() is not None
+        assert NumericalSolver.undefined().is_defined() is False