Implemented dense_output in python engines.

Added quadratic interpolators to `BaseTrajData` and `TrajectoryData`.
Moved all `TrajectoryData` instantiators to trajectory_data.py.
Rewrote `TrajectoryDataFilter` to interpolate for all requested points, and to generate multiple `TrajectoryData` points in a single call when required.

Author: dbookstaber

examples/Examples.ipynb

@@ -115,7 +115,7 @@ def _integrate(self, shot_info: Shot, maximum_range: float, record_step: float,
             if filter_flags:  # require check before call to improve performance
 
                 # Record TrajectoryData row
-                if (data := data_filter.should_record(range_vector, velocity_vector, mach, time)) is not None:
+                if (data := data_filter.record(range_vector, velocity_vector, mach, time)) is not None:
                     ranges.append(create_trajectory_row(data.time, data.position, data.velocity,
                                                         data.velocity.magnitude(), data.mach,
                                                         self.spin_drift(data.time), self.look_angle_rad,

examples/ExtremeExamples.ipynb

@@ -143,7 +143,7 @@ def _integrate(self, shot_info: Shot, maximum_range: float, record_step: float,
             # region Check whether to record TrajectoryData row at current point
             if filter_flags:  # require check before call to improve performance
                 # Record TrajectoryData row
-                if (data := data_filter.should_record(range_vector, velocity_vector, mach, time)) is not None:
+                if (data := data_filter.record(range_vector, velocity_vector, mach, time)) is not None:
                     ranges.append(create_trajectory_row(data.time, data.position, data.velocity,
                                                         data.velocity.magnitude(), data.mach,
                                                         self.spin_drift(data.time), self.look_angle_rad,

examples/integrators/euler_inline_op.py

@@ -151,7 +151,7 @@ def _integrate(self, shot_info: Shot, maximum_range: float, record_step: float,
             # region Check whether to record TrajectoryData row at current point
             if filter_flags:  # require check before call to improve performance
                 # Record TrajectoryData row
-                if (data := data_filter.should_record(range_vector, velocity_vector, mach, time)) is not None:
+                if (data := data_filter.record(range_vector, velocity_vector, mach, time)) is not None:
                     ranges.append(create_trajectory_row(data.time, data.position, data.velocity,
                                                         data.velocity.magnitude(), data.mach,
                                                         self.spin_drift(data.time), self.look_angle_rad,

examples/integrators/leapfrog.py

@@ -150,7 +150,7 @@ def _integrate(self, shot_info: Shot, maximum_range: float, record_step: float,
             # region Check whether to record TrajectoryData row at current point
             if filter_flags:  # require check before call to improve performance
                 # Record TrajectoryData row
-                if (data := data_filter.should_record(range_vector, velocity_vector, mach, time)) is not None:
+                if (data := data_filter.record(range_vector, velocity_vector, mach, time)) is not None:
                     ranges.append(create_trajectory_row(data.time, data.position, data.velocity,
                                                         data.velocity.magnitude(), data.mach,
                                                         self.spin_drift(data.time), self.look_angle_rad,

examples/integrators/leapfrog2.py

@@ -111,7 +111,7 @@ def _integrate(self, shot_info: Shot, maximum_range: float, record_step: float,
             # region Check whether to record TrajectoryData row at current point
             if filter_flags:  # require check before call to improve performance
                 # Record TrajectoryData row
-                if (data := data_filter.should_record(range_vector, velocity_vector, mach, time)) is not None:
+                if (data := data_filter.record(range_vector, velocity_vector, mach, time)) is not None:
                     ranges.append(create_trajectory_row(data.time, data.position, data.velocity,
                                                         data.velocity.magnitude(), data.mach,
                                                         self.spin_drift(data.time), self.look_angle_rad,

examples/integrators/leapfrog3_time_step.py

@@ -127,7 +127,7 @@ def _integrate(self, shot_info: Shot, maximum_range: float, record_step: float,
             # region Check whether to record TrajectoryData row at current point
             if filter_flags:  # require check before call to improve performance
                 # Record TrajectoryData row
-                if (data := data_filter.should_record(range_vector, velocity_vector, mach, time)) is not None:
+                if (data := data_filter.record(range_vector, velocity_vector, mach, time)) is not None:
                     ranges.append(create_trajectory_row(data.time, data.position, data.velocity,
                                                         data.velocity.magnitude(), data.mach,
                                                         self.spin_drift(data.time), self.look_angle_rad,

examples/integrators/rk4_inline_op.py

@@ -267,7 +267,7 @@ cdef class CythonizedBaseIntegrationEngine:
         self._free_trajectory()
         props = ShotProps.from_shot(shot_info)
         props.filter_flags = filter_flags
-        return HitResult(props, object[0], filter_flags > 0, object[1])
+        return HitResult(props, object[0], None, filter_flags > 0, object[1])
 
     cdef void _free_trajectory(CythonizedBaseIntegrationEngine self):
         if self._wind_sock is not NULL:

examples/integrators/verlet3.py

@@ -178,7 +178,6 @@ def _load_mixed_units():
     'calculate_ogw',
     'get_correction',
     'create_trajectory_row',
-    'make_trajectory_row',
     'calculate_curve',
     'CurvePoint',
     'exceptions',