brahe 1.4.0

/*!
 * Pre-made event detectors
 *
 * Domain-specific event detectors for altitude, orbital elements, and eclipse detection.
 */

use super::common::{DBinaryEvent, DValueEvent, SBinaryEvent, SValueEvent};
use super::traits::{
    DEventCallback, DEventDetector, EdgeType, EventAction, EventDirection, SEventCallback,
    SEventDetector,
};
use crate::access::PolygonLocation;
use crate::constants::AngleFormat;
use crate::coordinates::{point_in_polygon, position_ecef_to_geodetic, state_eci_to_koe};
use crate::frames::position_eci_to_ecef;
use crate::orbit_dynamics::{eclipse_conical, sun_position};
use crate::orbits::{anomaly_mean_to_eccentric, anomaly_mean_to_true};
use crate::propagators::EphemerisSource;
use crate::spice::{SPKKernel, sun_position_de};
use crate::time::Epoch;
use nalgebra::{DVector, SVector, Vector3, Vector6};
use std::sync::Arc;

/// Macro to implement SEventDetector trait delegation for wrapper types (static-sized)
macro_rules! impl_sevent_detector_delegate {
    ($struct:ty, $inner:ident, $s:ident, $p:ident) => {
        impl<const $s: usize, const $p: usize> SEventDetector<$s, $p> for $struct {
            fn evaluate(
                &self,
                t: Epoch,
                state: &SVector<f64, $s>,
                params: Option<&SVector<f64, $p>>,
            ) -> f64 {
                self.$inner.evaluate(t, state, params)
            }

            fn target_value(&self) -> f64 {
                self.$inner.target_value()
            }

            fn name(&self) -> &str {
                self.$inner.name()
            }

            fn direction(&self) -> EventDirection {
                self.$inner.direction()
            }

            fn time_tolerance(&self) -> f64 {
                self.$inner.time_tolerance()
            }

            fn value_tolerance(&self) -> f64 {
                self.$inner.value_tolerance()
            }

            fn step_reduction_factor(&self) -> f64 {
                self.$inner.step_reduction_factor()
            }

            fn callback(&self) -> Option<&SEventCallback<$s, $p>> {
                self.$inner.callback()
            }

            fn action(&self) -> EventAction {
                self.$inner.action()
            }
        }
    };
}

/// Macro to implement DEventDetector trait delegation for wrapper types (dynamic-sized)
macro_rules! impl_devent_detector_delegate {
    ($struct:ty, $inner:ident) => {
        impl DEventDetector for $struct {
            fn evaluate(
                &self,
                t: Epoch,
                state: &DVector<f64>,
                params: Option<&DVector<f64>>,
            ) -> f64 {
                self.$inner.evaluate(t, state, params)
            }

            fn target_value(&self) -> f64 {
                self.$inner.target_value()
            }

            fn name(&self) -> &str {
                self.$inner.name()
            }

            fn direction(&self) -> EventDirection {
                self.$inner.direction()
            }

            fn time_tolerance(&self) -> f64 {
                self.$inner.time_tolerance()
            }

            fn value_tolerance(&self) -> f64 {
                self.$inner.value_tolerance()
            }

            fn step_reduction_factor(&self) -> f64 {
                self.$inner.step_reduction_factor()
            }

            fn callback(&self) -> Option<&DEventCallback> {
                self.$inner.callback()
            }

            fn action(&self) -> EventAction {
                self.$inner.action()
            }
        }
    };
}

/// Altitude-based event detector (for orbital states, static-sized)
///
/// Convenience wrapper around [`SValueEvent`] for detecting geodetic altitude
/// crossings. Assumes the first 3 elements of the state vector are ECI position
/// in meters. Transforms to geodetic coordinates to compute altitude above the
/// WGS84 ellipsoid.
///
/// # Examples
/// ```
/// use brahe::events::SAltitudeEvent;
/// use brahe::events::EventDirection;
///
/// // Detect when altitude drops below 500 km (6D state)
/// let event = SAltitudeEvent::<6, 0>::new(
///     500e3,
///     "Low Altitude Warning",
///     EventDirection::Decreasing
/// );
///
/// // Works with extended states (e.g., 7D with mass)
/// let event = SAltitudeEvent::<7, 4>::new(
///     200e3,
///     "Reentry",
///     EventDirection::Decreasing
/// );
/// ```
pub struct SAltitudeEvent<const S: usize, const P: usize> {
    inner: SValueEvent<S, P>,
}

impl<const S: usize, const P: usize> SAltitudeEvent<S, P> {
    /// Create new altitude event
    ///
    /// # Arguments
    /// * `value_altitude` - Geodetic altitude value in meters above WGS84 ellipsoid
    /// * `name` - Event name
    /// * `direction` - Detection direction
    pub fn new(value_altitude: f64, name: impl Into<String>, direction: EventDirection) -> Self {
        let altitude_fn = |t: Epoch, state: &SVector<f64, S>, _params: Option<&SVector<f64, P>>| {
            // Extract ECI position (first 3 elements)
            let r_eci = state.fixed_rows::<3>(0).into_owned();

            // Transform ECI -> ECEF -> Geodetic
            let r_ecef = position_eci_to_ecef(t, r_eci);
            let geodetic = position_ecef_to_geodetic(r_ecef, AngleFormat::Radians);

            // Return geodetic altitude (index 2: [lon, lat, alt])
            geodetic[2]
        };

        Self {
            inner: SValueEvent::new(name, altitude_fn, value_altitude, direction),
        }
    }

    /// Set instance number for display name
    pub fn with_instance(mut self, instance: usize) -> Self {
        self.inner = self.inner.with_instance(instance);
        self
    }

    /// Set custom tolerances for event detection
    pub fn with_tolerances(mut self, time_tol: f64, value_tol: f64) -> Self {
        self.inner = self.inner.with_tolerances(time_tol, value_tol);
        self
    }

    /// Set custom step reduction factor for bisection search
    ///
    /// When a zero-crossing is detected, the new step size is set to this
    /// factor times the bracket width. Default: 0.2
    pub fn with_step_reduction_factor(mut self, factor: f64) -> Self {
        self.inner = self.inner.with_step_reduction_factor(factor);
        self
    }

    /// Set event callback
    pub fn with_callback(mut self, callback: SEventCallback<S, P>) -> Self {
        self.inner = self.inner.with_callback(callback);
        self
    }

    /// Mark as terminal event (stops propagation)
    pub fn set_terminal(mut self) -> Self {
        self.inner = self.inner.set_terminal();
        self
    }
}

// Delegate EventDetector trait implementation to inner SValueEvent
impl_sevent_detector_delegate!(SAltitudeEvent<S, P>, inner, S, P);

/// Dynamic-sized altitude event detector
///
/// See [`SAltitudeEvent`] for details. This version works with dynamic-sized
/// state vectors. Assumes first 3 elements are ECI position in meters.
pub struct DAltitudeEvent {
    inner: DValueEvent,
}

impl DAltitudeEvent {
    /// Create new altitude event
    ///
    /// # Arguments
    /// * `value_altitude` - Geodetic altitude value in meters above WGS84 ellipsoid
    /// * `name` - Event name
    /// * `direction` - Detection direction
    pub fn new(value_altitude: f64, name: impl Into<String>, direction: EventDirection) -> Self {
        let altitude_fn = |t: Epoch, state: &DVector<f64>, _params: Option<&DVector<f64>>| {
            // Extract ECI position (first 3 elements) and convert to Vector3
            use nalgebra::Vector3;
            let r_eci = Vector3::new(state[0], state[1], state[2]);

            // Transform ECI -> ECEF -> Geodetic
            let r_ecef = position_eci_to_ecef(t, r_eci);
            let geodetic = position_ecef_to_geodetic(r_ecef, AngleFormat::Radians);

            // Return geodetic altitude (index 2: [lon, lat, alt])
            geodetic[2]
        };

        Self {
            inner: DValueEvent::new(name, altitude_fn, value_altitude, direction),
        }
    }

    /// Set instance number for display name
    pub fn with_instance(mut self, instance: usize) -> Self {
        self.inner = self.inner.with_instance(instance);
        self
    }

    /// Set custom tolerances for event detection
    pub fn with_tolerances(mut self, time_tol: f64, value_tol: f64) -> Self {
        self.inner = self.inner.with_tolerances(time_tol, value_tol);
        self
    }

    /// Set custom step reduction factor for bisection search
    ///
    /// When a zero-crossing is detected, the new step size is set to this
    /// factor times the bracket width. Default: 0.2
    pub fn with_step_reduction_factor(mut self, factor: f64) -> Self {
        self.inner = self.inner.with_step_reduction_factor(factor);
        self
    }

    /// Set event callback
    pub fn with_callback(mut self, callback: DEventCallback) -> Self {
        self.inner = self.inner.with_callback(callback);
        self
    }

    /// Mark as terminal event
    pub fn set_terminal(mut self) -> Self {
        self.inner = self.inner.set_terminal();
        self
    }
}

// Delegate EventDetectorD trait implementation to inner DValueEvent
impl_devent_detector_delegate!(DAltitudeEvent, inner);

// =============================================================================
// Orbital Element Events
// =============================================================================

/// Semi-major axis event detector (static-sized)
///
/// Detects when semi-major axis crosses a value value. Computes orbital
/// elements from the Cartesian state (assumes first 6 elements are ECI state).
///
/// # Examples
/// ```
/// use brahe::events::{SSemiMajorAxisEvent, EventDirection};
/// use brahe::constants::R_EARTH;
///
/// // Detect when semi-major axis drops below 6800 km
/// let event = SSemiMajorAxisEvent::<6, 0>::new(
///     R_EARTH + 400e3,
///     "Low SMA",
///     EventDirection::Decreasing
/// );
/// ```
pub struct SSemiMajorAxisEvent<const S: usize, const P: usize> {
    inner: SValueEvent<S, P>,
}

impl<const S: usize, const P: usize> SSemiMajorAxisEvent<S, P> {
    /// Create new semi-major axis event
    ///
    /// # Arguments
    /// * `value` - Semi-major axis value in meters
    /// * `name` - Event name
    /// * `direction` - Detection direction
    pub fn new(value: f64, name: impl Into<String>, direction: EventDirection) -> Self {
        let value_fn = |_t: Epoch, state: &SVector<f64, S>, _params: Option<&SVector<f64, P>>| {
            let state6 = Vector6::new(state[0], state[1], state[2], state[3], state[4], state[5]);
            let koe = state_eci_to_koe(state6, AngleFormat::Radians);
            koe[0] // Semi-major axis
        };

        Self {
            inner: SValueEvent::new(name, value_fn, value, direction),
        }
    }

    /// Set instance number for display name
    pub fn with_instance(mut self, instance: usize) -> Self {
        self.inner = self.inner.with_instance(instance);
        self
    }

    /// Set custom tolerances for event detection
    pub fn with_tolerances(mut self, time_tol: f64, value_tol: f64) -> Self {
        self.inner = self.inner.with_tolerances(time_tol, value_tol);
        self
    }

    /// Set custom step reduction factor for bisection search
    pub fn with_step_reduction_factor(mut self, factor: f64) -> Self {
        self.inner = self.inner.with_step_reduction_factor(factor);
        self
    }

    /// Set event callback
    pub fn with_callback(mut self, callback: SEventCallback<S, P>) -> Self {
        self.inner = self.inner.with_callback(callback);
        self
    }

    /// Mark as terminal event (stops propagation)
    pub fn set_terminal(mut self) -> Self {
        self.inner = self.inner.set_terminal();
        self
    }
}

impl_sevent_detector_delegate!(SSemiMajorAxisEvent<S, P>, inner, S, P);

/// Dynamic-sized semi-major axis event detector
pub struct DSemiMajorAxisEvent {
    inner: DValueEvent,
}

impl DSemiMajorAxisEvent {
    /// Create new semi-major axis event
    pub fn new(value: f64, name: impl Into<String>, direction: EventDirection) -> Self {
        let value_fn = |_t: Epoch, state: &DVector<f64>, _params: Option<&DVector<f64>>| {
            let state6 = Vector6::new(state[0], state[1], state[2], state[3], state[4], state[5]);
            let koe = state_eci_to_koe(state6, AngleFormat::Radians);
            koe[0]
        };

        Self {
            inner: DValueEvent::new(name, value_fn, value, direction),
        }
    }

    /// Set instance number for display name
    pub fn with_instance(mut self, instance: usize) -> Self {
        self.inner = self.inner.with_instance(instance);
        self
    }

    /// Set custom tolerances for event detection
    pub fn with_tolerances(mut self, time_tol: f64, value_tol: f64) -> Self {
        self.inner = self.inner.with_tolerances(time_tol, value_tol);
        self
    }

    /// Set custom step reduction factor for bisection search
    pub fn with_step_reduction_factor(mut self, factor: f64) -> Self {
        self.inner = self.inner.with_step_reduction_factor(factor);
        self
    }

    /// Set event callback
    pub fn with_callback(mut self, callback: DEventCallback) -> Self {
        self.inner = self.inner.with_callback(callback);
        self
    }

    /// Mark as terminal event
    pub fn set_terminal(mut self) -> Self {
        self.inner = self.inner.set_terminal();
        self
    }
}

impl_devent_detector_delegate!(DSemiMajorAxisEvent, inner);

/// Eccentricity event detector (static-sized)
///
/// Detects when eccentricity crosses a value value.
///
/// # Examples
/// ```
/// use brahe::events::{SEccentricityEvent, EventDirection};
///
/// // Detect when eccentricity exceeds 0.1
/// let event = SEccentricityEvent::<6, 0>::new(
///     0.1,
///     "High Eccentricity",
///     EventDirection::Increasing
/// );
/// ```
pub struct SEccentricityEvent<const S: usize, const P: usize> {
    inner: SValueEvent<S, P>,
}

impl<const S: usize, const P: usize> SEccentricityEvent<S, P> {
    /// Create new eccentricity event
    pub fn new(value: f64, name: impl Into<String>, direction: EventDirection) -> Self {
        let value_fn = |_t: Epoch, state: &SVector<f64, S>, _params: Option<&SVector<f64, P>>| {
            let state6 = Vector6::new(state[0], state[1], state[2], state[3], state[4], state[5]);
            let koe = state_eci_to_koe(state6, AngleFormat::Radians);
            koe[1] // Eccentricity
        };

        Self {
            inner: SValueEvent::new(name, value_fn, value, direction),
        }
    }

    /// Set instance number for display name
    pub fn with_instance(mut self, instance: usize) -> Self {
        self.inner = self.inner.with_instance(instance);
        self
    }

    /// Set custom tolerances for event detection
    pub fn with_tolerances(mut self, time_tol: f64, value_tol: f64) -> Self {
        self.inner = self.inner.with_tolerances(time_tol, value_tol);
        self
    }

    /// Set custom step reduction factor for bisection search
    pub fn with_step_reduction_factor(mut self, factor: f64) -> Self {
        self.inner = self.inner.with_step_reduction_factor(factor);
        self
    }

    /// Set event callback
    pub fn with_callback(mut self, callback: SEventCallback<S, P>) -> Self {
        self.inner = self.inner.with_callback(callback);
        self
    }

    /// Mark as terminal event (stops propagation)
    pub fn set_terminal(mut self) -> Self {
        self.inner = self.inner.set_terminal();
        self
    }
}

impl_sevent_detector_delegate!(SEccentricityEvent<S, P>, inner, S, P);

/// Dynamic-sized eccentricity event detector
pub struct DEccentricityEvent {
    inner: DValueEvent,
}

impl DEccentricityEvent {
    /// Create new eccentricity event
    pub fn new(value: f64, name: impl Into<String>, direction: EventDirection) -> Self {
        let value_fn = |_t: Epoch, state: &DVector<f64>, _params: Option<&DVector<f64>>| {
            let state6 = Vector6::new(state[0], state[1], state[2], state[3], state[4], state[5]);
            let koe = state_eci_to_koe(state6, AngleFormat::Radians);
            koe[1]
        };

        Self {
            inner: DValueEvent::new(name, value_fn, value, direction),
        }
    }

    /// Set instance number for display name
    pub fn with_instance(mut self, instance: usize) -> Self {
        self.inner = self.inner.with_instance(instance);
        self
    }

    /// Set custom tolerances for event detection
    pub fn with_tolerances(mut self, time_tol: f64, value_tol: f64) -> Self {
        self.inner = self.inner.with_tolerances(time_tol, value_tol);
        self
    }

    /// Set custom step reduction factor for bisection search
    pub fn with_step_reduction_factor(mut self, factor: f64) -> Self {
        self.inner = self.inner.with_step_reduction_factor(factor);
        self
    }

    /// Set event callback
    pub fn with_callback(mut self, callback: DEventCallback) -> Self {
        self.inner = self.inner.with_callback(callback);
        self
    }

    /// Mark as terminal event
    pub fn set_terminal(mut self) -> Self {
        self.inner = self.inner.set_terminal();
        self
    }
}

impl_devent_detector_delegate!(DEccentricityEvent, inner);

/// Inclination event detector (static-sized)
///
/// Detects when inclination crosses a value value. Returns inclination in radians.
///
/// # Examples
/// ```
/// use brahe::events::{SInclinationEvent, EventDirection};
/// use brahe::AngleFormat;
///
/// // Detect when inclination exceeds 90 degrees
/// let event = SInclinationEvent::<6, 0>::new(
///     90.0,
///     "Retrograde",
///     EventDirection::Increasing,
///     AngleFormat::Degrees
/// );
/// ```
pub struct SInclinationEvent<const S: usize, const P: usize> {
    inner: SValueEvent<S, P>,
}

impl<const S: usize, const P: usize> SInclinationEvent<S, P> {
    /// Create new inclination event
    ///
    /// # Arguments
    /// * `value` - Inclination value
    /// * `name` - Event name
    /// * `direction` - Detection direction
    /// * `angle_format` - Whether value is in degrees or radians
    pub fn new(
        value: f64,
        name: impl Into<String>,
        direction: EventDirection,
        angle_format: AngleFormat,
    ) -> Self {
        // Convert value to radians if needed
        let value_rad = match angle_format {
            AngleFormat::Degrees => value.to_radians(),
            AngleFormat::Radians => value,
        };

        let value_fn = |_t: Epoch, state: &SVector<f64, S>, _params: Option<&SVector<f64, P>>| {
            let state6 = Vector6::new(state[0], state[1], state[2], state[3], state[4], state[5]);
            let koe = state_eci_to_koe(state6, AngleFormat::Radians);
            koe[2] // Inclination
        };

        Self {
            inner: SValueEvent::new(name, value_fn, value_rad, direction),
        }
    }

    /// Set instance number for display name
    pub fn with_instance(mut self, instance: usize) -> Self {
        self.inner = self.inner.with_instance(instance);
        self
    }

    /// Set custom tolerances for event detection
    pub fn with_tolerances(mut self, time_tol: f64, value_tol: f64) -> Self {
        self.inner = self.inner.with_tolerances(time_tol, value_tol);
        self
    }

    /// Set custom step reduction factor for bisection search
    pub fn with_step_reduction_factor(mut self, factor: f64) -> Self {
        self.inner = self.inner.with_step_reduction_factor(factor);
        self
    }

    /// Set event callback
    pub fn with_callback(mut self, callback: SEventCallback<S, P>) -> Self {
        self.inner = self.inner.with_callback(callback);
        self
    }

    /// Mark as terminal event (stops propagation)
    pub fn set_terminal(mut self) -> Self {
        self.inner = self.inner.set_terminal();
        self
    }
}

impl_sevent_detector_delegate!(SInclinationEvent<S, P>, inner, S, P);

/// Dynamic-sized inclination event detector
pub struct DInclinationEvent {
    inner: DValueEvent,
}

impl DInclinationEvent {
    /// Create new inclination event
    ///
    /// # Arguments
    /// * `value` - Inclination value
    /// * `name` - Event name
    /// * `direction` - Detection direction
    /// * `angle_format` - Whether value is in degrees or radians
    pub fn new(
        value: f64,
        name: impl Into<String>,
        direction: EventDirection,
        angle_format: AngleFormat,
    ) -> Self {
        // Convert value to radians if needed
        let value_rad = match angle_format {
            AngleFormat::Degrees => value.to_radians(),
            AngleFormat::Radians => value,
        };

        let value_fn = |_t: Epoch, state: &DVector<f64>, _params: Option<&DVector<f64>>| {
            let state6 = Vector6::new(state[0], state[1], state[2], state[3], state[4], state[5]);
            let koe = state_eci_to_koe(state6, AngleFormat::Radians);
            koe[2]
        };

        Self {
            inner: DValueEvent::new(name, value_fn, value_rad, direction),
        }
    }

    /// Set instance number for display name
    pub fn with_instance(mut self, instance: usize) -> Self {
        self.inner = self.inner.with_instance(instance);
        self
    }

    /// Set custom tolerances for event detection
    pub fn with_tolerances(mut self, time_tol: f64, value_tol: f64) -> Self {
        self.inner = self.inner.with_tolerances(time_tol, value_tol);
        self
    }

    /// Set custom step reduction factor for bisection search
    pub fn with_step_reduction_factor(mut self, factor: f64) -> Self {
        self.inner = self.inner.with_step_reduction_factor(factor);
        self
    }

    /// Set event callback
    pub fn with_callback(mut self, callback: DEventCallback) -> Self {
        self.inner = self.inner.with_callback(callback);
        self
    }

    /// Mark as terminal event
    pub fn set_terminal(mut self) -> Self {
        self.inner = self.inner.set_terminal();
        self
    }
}

impl_devent_detector_delegate!(DInclinationEvent, inner);

/// Argument of perigee event detector (static-sized)
///
/// Detects when argument of perigee crosses a value value. Returns in radians.
pub struct SArgumentOfPerigeeEvent<const S: usize, const P: usize> {
    inner: SValueEvent<S, P>,
}

impl<const S: usize, const P: usize> SArgumentOfPerigeeEvent<S, P> {
    /// Create new argument of perigee event
    ///
    /// # Arguments
    /// * `value` - Argument of perigee value
    /// * `name` - Event name
    /// * `direction` - Detection direction
    /// * `angle_format` - Whether value is in degrees or radians
    pub fn new(
        value: f64,
        name: impl Into<String>,
        direction: EventDirection,
        angle_format: AngleFormat,
    ) -> Self {
        let value_rad = match angle_format {
            AngleFormat::Degrees => value.to_radians(),
            AngleFormat::Radians => value,
        };

        let value_fn = |_t: Epoch, state: &SVector<f64, S>, _params: Option<&SVector<f64, P>>| {
            let state6 = Vector6::new(state[0], state[1], state[2], state[3], state[4], state[5]);
            let koe = state_eci_to_koe(state6, AngleFormat::Radians);
            koe[4] // Argument of perigee
        };

        Self {
            inner: SValueEvent::new(name, value_fn, value_rad, direction),
        }
    }

    /// Set instance number for display name
    pub fn with_instance(mut self, instance: usize) -> Self {
        self.inner = self.inner.with_instance(instance);
        self
    }

    /// Set custom tolerances for event detection
    pub fn with_tolerances(mut self, time_tol: f64, value_tol: f64) -> Self {
        self.inner = self.inner.with_tolerances(time_tol, value_tol);
        self
    }

    /// Set custom step reduction factor for bisection search
    pub fn with_step_reduction_factor(mut self, factor: f64) -> Self {
        self.inner = self.inner.with_step_reduction_factor(factor);
        self
    }

    /// Set event callback
    pub fn with_callback(mut self, callback: SEventCallback<S, P>) -> Self {
        self.inner = self.inner.with_callback(callback);
        self
    }

    /// Mark as terminal event (stops propagation)
    pub fn set_terminal(mut self) -> Self {
        self.inner = self.inner.set_terminal();
        self
    }
}

impl_sevent_detector_delegate!(SArgumentOfPerigeeEvent<S, P>, inner, S, P);

/// Dynamic-sized argument of perigee event detector
pub struct DArgumentOfPerigeeEvent {
    inner: DValueEvent,
}

impl DArgumentOfPerigeeEvent {
    /// Create new argument of perigee event
    ///
    /// # Arguments
    /// * `value` - Argument of perigee value
    /// * `name` - Event name
    /// * `direction` - Detection direction
    /// * `angle_format` - Whether value is in degrees or radians
    pub fn new(
        value: f64,
        name: impl Into<String>,
        direction: EventDirection,
        angle_format: AngleFormat,
    ) -> Self {
        let value_rad = match angle_format {
            AngleFormat::Degrees => value.to_radians(),
            AngleFormat::Radians => value,
        };

        let value_fn = |_t: Epoch, state: &DVector<f64>, _params: Option<&DVector<f64>>| {
            let state6 = Vector6::new(state[0], state[1], state[2], state[3], state[4], state[5]);
            let koe = state_eci_to_koe(state6, AngleFormat::Radians);
            koe[4]
        };

        Self {
            inner: DValueEvent::new(name, value_fn, value_rad, direction),
        }
    }

    /// Set instance number for display name
    pub fn with_instance(mut self, instance: usize) -> Self {
        self.inner = self.inner.with_instance(instance);
        self
    }

    /// Set custom tolerances for event detection
    pub fn with_tolerances(mut self, time_tol: f64, value_tol: f64) -> Self {
        self.inner = self.inner.with_tolerances(time_tol, value_tol);
        self
    }

    /// Set custom step reduction factor for bisection search
    pub fn with_step_reduction_factor(mut self, factor: f64) -> Self {
        self.inner = self.inner.with_step_reduction_factor(factor);
        self
    }

    /// Set event callback
    pub fn with_callback(mut self, callback: DEventCallback) -> Self {
        self.inner = self.inner.with_callback(callback);
        self
    }

    /// Mark as terminal event
    pub fn set_terminal(mut self) -> Self {
        self.inner = self.inner.set_terminal();
        self
    }
}

impl_devent_detector_delegate!(DArgumentOfPerigeeEvent, inner);

/// Mean anomaly event detector (static-sized)
///
/// Detects when mean anomaly crosses a value value. Returns in radians.
pub struct SMeanAnomalyEvent<const S: usize, const P: usize> {
    inner: SValueEvent<S, P>,
}

impl<const S: usize, const P: usize> SMeanAnomalyEvent<S, P> {
    /// Create new mean anomaly event
    ///
    /// # Arguments
    /// * `value` - Mean anomaly value
    /// * `name` - Event name
    /// * `direction` - Detection direction
    /// * `angle_format` - Whether value is in degrees or radians
    pub fn new(
        value: f64,
        name: impl Into<String>,
        direction: EventDirection,
        angle_format: AngleFormat,
    ) -> Self {
        let value_rad = match angle_format {
            AngleFormat::Degrees => value.to_radians(),
            AngleFormat::Radians => value,
        };

        let value_fn = |_t: Epoch, state: &SVector<f64, S>, _params: Option<&SVector<f64, P>>| {
            let state6 = Vector6::new(state[0], state[1], state[2], state[3], state[4], state[5]);
            let koe = state_eci_to_koe(state6, AngleFormat::Radians);
            koe[5] // Mean anomaly
        };

        Self {
            inner: SValueEvent::new(name, value_fn, value_rad, direction),
        }
    }

    /// Set instance number for display name
    pub fn with_instance(mut self, instance: usize) -> Self {
        self.inner = self.inner.with_instance(instance);
        self
    }

    /// Set custom tolerances for event detection
    pub fn with_tolerances(mut self, time_tol: f64, value_tol: f64) -> Self {
        self.inner = self.inner.with_tolerances(time_tol, value_tol);
        self
    }

    /// Set custom step reduction factor for bisection search
    pub fn with_step_reduction_factor(mut self, factor: f64) -> Self {
        self.inner = self.inner.with_step_reduction_factor(factor);
        self
    }

    /// Set event callback
    pub fn with_callback(mut self, callback: SEventCallback<S, P>) -> Self {
        self.inner = self.inner.with_callback(callback);
        self
    }

    /// Mark as terminal event (stops propagation)
    pub fn set_terminal(mut self) -> Self {
        self.inner = self.inner.set_terminal();
        self
    }
}

impl_sevent_detector_delegate!(SMeanAnomalyEvent<S, P>, inner, S, P);

/// Dynamic-sized mean anomaly event detector
pub struct DMeanAnomalyEvent {
    inner: DValueEvent,
}

impl DMeanAnomalyEvent {
    /// Create new mean anomaly event
    ///
    /// # Arguments
    /// * `value` - Mean anomaly value
    /// * `name` - Event name
    /// * `direction` - Detection direction
    /// * `angle_format` - Whether value is in degrees or radians
    pub fn new(
        value: f64,
        name: impl Into<String>,
        direction: EventDirection,
        angle_format: AngleFormat,
    ) -> Self {
        let value_rad = match angle_format {
            AngleFormat::Degrees => value.to_radians(),
            AngleFormat::Radians => value,
        };

        let value_fn = |_t: Epoch, state: &DVector<f64>, _params: Option<&DVector<f64>>| {
            let state6 = Vector6::new(state[0], state[1], state[2], state[3], state[4], state[5]);
            let koe = state_eci_to_koe(state6, AngleFormat::Radians);
            koe[5]
        };

        Self {
            inner: DValueEvent::new(name, value_fn, value_rad, direction),
        }
    }

    /// Set instance number for display name
    pub fn with_instance(mut self, instance: usize) -> Self {
        self.inner = self.inner.with_instance(instance);
        self
    }

    /// Set custom tolerances for event detection
    pub fn with_tolerances(mut self, time_tol: f64, value_tol: f64) -> Self {
        self.inner = self.inner.with_tolerances(time_tol, value_tol);
        self
    }

    /// Set custom step reduction factor for bisection search
    pub fn with_step_reduction_factor(mut self, factor: f64) -> Self {
        self.inner = self.inner.with_step_reduction_factor(factor);
        self
    }

    /// Set event callback
    pub fn with_callback(mut self, callback: DEventCallback) -> Self {
        self.inner = self.inner.with_callback(callback);
        self
    }

    /// Mark as terminal event
    pub fn set_terminal(mut self) -> Self {
        self.inner = self.inner.set_terminal();
        self
    }
}

impl_devent_detector_delegate!(DMeanAnomalyEvent, inner);

/// Eccentric anomaly event detector (static-sized)
///
/// Detects when eccentric anomaly crosses a value value. Returns in radians.
pub struct SEccentricAnomalyEvent<const S: usize, const P: usize> {
    inner: SValueEvent<S, P>,
}

impl<const S: usize, const P: usize> SEccentricAnomalyEvent<S, P> {
    /// Create new eccentric anomaly event
    ///
    /// # Arguments
    /// * `value` - Eccentric anomaly value
    /// * `name` - Event name
    /// * `direction` - Detection direction
    /// * `angle_format` - Whether value is in degrees or radians
    pub fn new(
        value: f64,
        name: impl Into<String>,
        direction: EventDirection,
        angle_format: AngleFormat,
    ) -> Self {
        let value_rad = match angle_format {
            AngleFormat::Degrees => value.to_radians(),
            AngleFormat::Radians => value,
        };

        let value_fn = |_t: Epoch, state: &SVector<f64, S>, _params: Option<&SVector<f64, P>>| {
            let state6 = Vector6::new(state[0], state[1], state[2], state[3], state[4], state[5]);
            let koe = state_eci_to_koe(state6, AngleFormat::Radians);
            let e = koe[1];
            let m = koe[5];
            // Convert mean anomaly to eccentric anomaly
            anomaly_mean_to_eccentric(m, e, AngleFormat::Radians).unwrap_or(m)
        };

        Self {
            inner: SValueEvent::new(name, value_fn, value_rad, direction),
        }
    }

    /// Set instance number for display name
    pub fn with_instance(mut self, instance: usize) -> Self {
        self.inner = self.inner.with_instance(instance);
        self
    }

    /// Set custom tolerances for event detection
    pub fn with_tolerances(mut self, time_tol: f64, value_tol: f64) -> Self {
        self.inner = self.inner.with_tolerances(time_tol, value_tol);
        self
    }

    /// Set custom step reduction factor for bisection search
    pub fn with_step_reduction_factor(mut self, factor: f64) -> Self {
        self.inner = self.inner.with_step_reduction_factor(factor);
        self
    }

    /// Set event callback
    pub fn with_callback(mut self, callback: SEventCallback<S, P>) -> Self {
        self.inner = self.inner.with_callback(callback);
        self
    }

    /// Mark as terminal event (stops propagation)
    pub fn set_terminal(mut self) -> Self {
        self.inner = self.inner.set_terminal();
        self
    }
}

impl_sevent_detector_delegate!(SEccentricAnomalyEvent<S, P>, inner, S, P);

/// Dynamic-sized eccentric anomaly event detector
pub struct DEccentricAnomalyEvent {
    inner: DValueEvent,
}

impl DEccentricAnomalyEvent {
    /// Create new eccentric anomaly event
    ///
    /// # Arguments
    /// * `value` - Eccentric anomaly value
    /// * `name` - Event name
    /// * `direction` - Detection direction
    /// * `angle_format` - Whether value is in degrees or radians
    pub fn new(
        value: f64,
        name: impl Into<String>,
        direction: EventDirection,
        angle_format: AngleFormat,
    ) -> Self {
        let value_rad = match angle_format {
            AngleFormat::Degrees => value.to_radians(),
            AngleFormat::Radians => value,
        };

        let value_fn = |_t: Epoch, state: &DVector<f64>, _params: Option<&DVector<f64>>| {
            let state6 = Vector6::new(state[0], state[1], state[2], state[3], state[4], state[5]);
            let koe = state_eci_to_koe(state6, AngleFormat::Radians);
            let e = koe[1];
            let m = koe[5];
            anomaly_mean_to_eccentric(m, e, AngleFormat::Radians).unwrap_or(m)
        };

        Self {
            inner: DValueEvent::new(name, value_fn, value_rad, direction),
        }
    }

    /// Set instance number for display name
    pub fn with_instance(mut self, instance: usize) -> Self {
        self.inner = self.inner.with_instance(instance);
        self
    }

    /// Set custom tolerances for event detection
    pub fn with_tolerances(mut self, time_tol: f64, value_tol: f64) -> Self {
        self.inner = self.inner.with_tolerances(time_tol, value_tol);
        self
    }

    /// Set custom step reduction factor for bisection search
    pub fn with_step_reduction_factor(mut self, factor: f64) -> Self {
        self.inner = self.inner.with_step_reduction_factor(factor);
        self
    }

    /// Set event callback
    pub fn with_callback(mut self, callback: DEventCallback) -> Self {
        self.inner = self.inner.with_callback(callback);
        self
    }

    /// Mark as terminal event
    pub fn set_terminal(mut self) -> Self {
        self.inner = self.inner.set_terminal();
        self
    }
}

impl_devent_detector_delegate!(DEccentricAnomalyEvent, inner);

/// True anomaly event detector (static-sized)
///
/// Detects when true anomaly crosses a value value. Returns in radians.
pub struct STrueAnomalyEvent<const S: usize, const P: usize> {
    inner: SValueEvent<S, P>,
}

impl<const S: usize, const P: usize> STrueAnomalyEvent<S, P> {
    /// Create new true anomaly event
    ///
    /// # Arguments
    /// * `value` - True anomaly value
    /// * `name` - Event name
    /// * `direction` - Detection direction
    /// * `angle_format` - Whether value is in degrees or radians
    pub fn new(
        value: f64,
        name: impl Into<String>,
        direction: EventDirection,
        angle_format: AngleFormat,
    ) -> Self {
        let value_rad = match angle_format {
            AngleFormat::Degrees => value.to_radians(),
            AngleFormat::Radians => value,
        };

        let value_fn = |_t: Epoch, state: &SVector<f64, S>, _params: Option<&SVector<f64, P>>| {
            let state6 = Vector6::new(state[0], state[1], state[2], state[3], state[4], state[5]);
            let koe = state_eci_to_koe(state6, AngleFormat::Radians);
            let e = koe[1];
            let m = koe[5];
            // Convert mean anomaly to true anomaly
            anomaly_mean_to_true(m, e, AngleFormat::Radians).unwrap_or(m)
        };

        Self {
            inner: SValueEvent::new(name, value_fn, value_rad, direction),
        }
    }

    /// Set instance number for display name
    pub fn with_instance(mut self, instance: usize) -> Self {
        self.inner = self.inner.with_instance(instance);
        self
    }

    /// Set custom tolerances for event detection
    pub fn with_tolerances(mut self, time_tol: f64, value_tol: f64) -> Self {
        self.inner = self.inner.with_tolerances(time_tol, value_tol);
        self
    }

    /// Set custom step reduction factor for bisection search
    pub fn with_step_reduction_factor(mut self, factor: f64) -> Self {
        self.inner = self.inner.with_step_reduction_factor(factor);
        self
    }

    /// Set event callback
    pub fn with_callback(mut self, callback: SEventCallback<S, P>) -> Self {
        self.inner = self.inner.with_callback(callback);
        self
    }

    /// Mark as terminal event (stops propagation)
    pub fn set_terminal(mut self) -> Self {
        self.inner = self.inner.set_terminal();
        self
    }
}

impl_sevent_detector_delegate!(STrueAnomalyEvent<S, P>, inner, S, P);

/// Dynamic-sized true anomaly event detector
pub struct DTrueAnomalyEvent {
    inner: DValueEvent,
}

impl DTrueAnomalyEvent {
    /// Create new true anomaly event
    ///
    /// # Arguments
    /// * `value` - True anomaly value
    /// * `name` - Event name
    /// * `direction` - Detection direction
    /// * `angle_format` - Whether value is in degrees or radians
    pub fn new(
        value: f64,
        name: impl Into<String>,
        direction: EventDirection,
        angle_format: AngleFormat,
    ) -> Self {
        let value_rad = match angle_format {
            AngleFormat::Degrees => value.to_radians(),
            AngleFormat::Radians => value,
        };

        let value_fn = |_t: Epoch, state: &DVector<f64>, _params: Option<&DVector<f64>>| {
            let state6 = Vector6::new(state[0], state[1], state[2], state[3], state[4], state[5]);
            let koe = state_eci_to_koe(state6, AngleFormat::Radians);
            let e = koe[1];
            let m = koe[5];
            anomaly_mean_to_true(m, e, AngleFormat::Radians).unwrap_or(m)
        };

        Self {
            inner: DValueEvent::new(name, value_fn, value_rad, direction),
        }
    }

    /// Set instance number for display name
    pub fn with_instance(mut self, instance: usize) -> Self {
        self.inner = self.inner.with_instance(instance);
        self
    }

    /// Set custom tolerances for event detection
    pub fn with_tolerances(mut self, time_tol: f64, value_tol: f64) -> Self {
        self.inner = self.inner.with_tolerances(time_tol, value_tol);
        self
    }

    /// Set custom step reduction factor for bisection search
    pub fn with_step_reduction_factor(mut self, factor: f64) -> Self {
        self.inner = self.inner.with_step_reduction_factor(factor);
        self
    }

    /// Set event callback
    pub fn with_callback(mut self, callback: DEventCallback) -> Self {
        self.inner = self.inner.with_callback(callback);
        self
    }

    /// Mark as terminal event
    pub fn set_terminal(mut self) -> Self {
        self.inner = self.inner.set_terminal();
        self
    }
}

impl_devent_detector_delegate!(DTrueAnomalyEvent, inner);

/// Argument of latitude event detector (static-sized)
///
/// Detects when argument of latitude (omega + true_anomaly) crosses a value.
/// Returns in radians.
pub struct SArgumentOfLatitudeEvent<const S: usize, const P: usize> {
    inner: SValueEvent<S, P>,
}

impl<const S: usize, const P: usize> SArgumentOfLatitudeEvent<S, P> {
    /// Create new argument of latitude event
    ///
    /// # Arguments
    /// * `value` - Argument of latitude value
    /// * `name` - Event name
    /// * `direction` - Detection direction
    /// * `angle_format` - Whether value is in degrees or radians
    pub fn new(
        value: f64,
        name: impl Into<String>,
        direction: EventDirection,
        angle_format: AngleFormat,
    ) -> Self {
        let value_rad = match angle_format {
            AngleFormat::Degrees => value.to_radians(),
            AngleFormat::Radians => value,
        };

        let value_fn = |_t: Epoch, state: &SVector<f64, S>, _params: Option<&SVector<f64, P>>| {
            let state6 = Vector6::new(state[0], state[1], state[2], state[3], state[4], state[5]);
            let koe = state_eci_to_koe(state6, AngleFormat::Radians);
            let e = koe[1];
            let omega = koe[4]; // Argument of perigee
            let m = koe[5];
            let nu = anomaly_mean_to_true(m, e, AngleFormat::Radians).unwrap_or(m);
            // Argument of latitude = omega + true anomaly
            let u = omega + nu;
            // Normalize to [0, 2*PI)
            u.rem_euclid(2.0 * std::f64::consts::PI)
        };

        Self {
            inner: SValueEvent::new(name, value_fn, value_rad, direction),
        }
    }

    /// Set instance number for display name
    pub fn with_instance(mut self, instance: usize) -> Self {
        self.inner = self.inner.with_instance(instance);
        self
    }

    /// Set custom tolerances for event detection
    pub fn with_tolerances(mut self, time_tol: f64, value_tol: f64) -> Self {
        self.inner = self.inner.with_tolerances(time_tol, value_tol);
        self
    }

    /// Set custom step reduction factor for bisection search
    pub fn with_step_reduction_factor(mut self, factor: f64) -> Self {
        self.inner = self.inner.with_step_reduction_factor(factor);
        self
    }

    /// Set event callback
    pub fn with_callback(mut self, callback: SEventCallback<S, P>) -> Self {
        self.inner = self.inner.with_callback(callback);
        self
    }

    /// Mark as terminal event (stops propagation)
    pub fn set_terminal(mut self) -> Self {
        self.inner = self.inner.set_terminal();
        self
    }
}

impl_sevent_detector_delegate!(SArgumentOfLatitudeEvent<S, P>, inner, S, P);

/// Dynamic-sized argument of latitude event detector
pub struct DArgumentOfLatitudeEvent {
    inner: DValueEvent,
}

impl DArgumentOfLatitudeEvent {
    /// Create new argument of latitude event
    ///
    /// # Arguments
    /// * `value` - Argument of latitude value
    /// * `name` - Event name
    /// * `direction` - Detection direction
    /// * `angle_format` - Whether value is in degrees or radians
    pub fn new(
        value: f64,
        name: impl Into<String>,
        direction: EventDirection,
        angle_format: AngleFormat,
    ) -> Self {
        let value_rad = match angle_format {
            AngleFormat::Degrees => value.to_radians(),
            AngleFormat::Radians => value,
        };

        let value_fn = |_t: Epoch, state: &DVector<f64>, _params: Option<&DVector<f64>>| {
            let state6 = Vector6::new(state[0], state[1], state[2], state[3], state[4], state[5]);
            let koe = state_eci_to_koe(state6, AngleFormat::Radians);
            let e = koe[1];
            let omega = koe[4];
            let m = koe[5];
            let nu = anomaly_mean_to_true(m, e, AngleFormat::Radians).unwrap_or(m);
            let u = omega + nu;
            u.rem_euclid(2.0 * std::f64::consts::PI)
        };

        Self {
            inner: DValueEvent::new(name, value_fn, value_rad, direction),
        }
    }

    /// Set instance number for display name
    pub fn with_instance(mut self, instance: usize) -> Self {
        self.inner = self.inner.with_instance(instance);
        self
    }

    /// Set custom tolerances for event detection
    pub fn with_tolerances(mut self, time_tol: f64, value_tol: f64) -> Self {
        self.inner = self.inner.with_tolerances(time_tol, value_tol);
        self
    }

    /// Set custom step reduction factor for bisection search
    pub fn with_step_reduction_factor(mut self, factor: f64) -> Self {
        self.inner = self.inner.with_step_reduction_factor(factor);
        self
    }

    /// Set event callback
    pub fn with_callback(mut self, callback: DEventCallback) -> Self {
        self.inner = self.inner.with_callback(callback);
        self
    }

    /// Mark as terminal event
    pub fn set_terminal(mut self) -> Self {
        self.inner = self.inner.set_terminal();
        self
    }
}

impl_devent_detector_delegate!(DArgumentOfLatitudeEvent, inner);

/// Ascending node event detector (static-sized)
///
/// Detects when spacecraft crosses the ascending node (equator from south to north).
/// The ascending node occurs when the true argument of latitude crosses 0 with increasing values.
///
/// # Examples
/// ```
/// use brahe::events::SAscendingNodeEvent;
///
/// // Detect ascending node crossing
/// let event = SAscendingNodeEvent::<6, 0>::new("Ascending Node");
/// ```
pub struct SAscendingNodeEvent<const S: usize, const P: usize> {
    inner: SValueEvent<S, P>,
}

impl<const S: usize, const P: usize> SAscendingNodeEvent<S, P> {
    /// Create new ascending node event
    ///
    /// # Arguments
    /// * `name` - Event name
    pub fn new(name: impl Into<String>) -> Self {
        let value_fn = |_t: Epoch, state: &SVector<f64, S>, _params: Option<&SVector<f64, P>>| {
            let state6 = Vector6::new(state[0], state[1], state[2], state[3], state[4], state[5]);
            let koe = state_eci_to_koe(state6, AngleFormat::Radians);
            let e = koe[1];
            let omega = koe[4]; // Argument of perigee
            let m = koe[5]; // Mean anomaly

            // Convert mean anomaly to true anomaly
            let nu = anomaly_mean_to_true(m, e, AngleFormat::Radians).unwrap_or(m);

            // True argument of latitude = omega + nu
            // Normalize to [-pi, pi) for proper zero crossing detection
            // Use rem_euclid to handle negative values correctly
            let u = omega + nu;
            (u + std::f64::consts::PI).rem_euclid(2.0 * std::f64::consts::PI) - std::f64::consts::PI
        };

        Self {
            inner: SValueEvent::new(name, value_fn, 0.0, EventDirection::Increasing),
        }
    }

    /// Set instance number for display name
    pub fn with_instance(mut self, instance: usize) -> Self {
        self.inner = self.inner.with_instance(instance);
        self
    }

    /// Set custom tolerances for event detection
    pub fn with_tolerances(mut self, time_tol: f64, value_tol: f64) -> Self {
        self.inner = self.inner.with_tolerances(time_tol, value_tol);
        self
    }

    /// Set custom step reduction factor for bisection search
    pub fn with_step_reduction_factor(mut self, factor: f64) -> Self {
        self.inner = self.inner.with_step_reduction_factor(factor);
        self
    }

    /// Set event callback
    pub fn with_callback(mut self, callback: SEventCallback<S, P>) -> Self {
        self.inner = self.inner.with_callback(callback);
        self
    }

    /// Mark as terminal event (stops propagation)
    pub fn set_terminal(mut self) -> Self {
        self.inner = self.inner.set_terminal();
        self
    }
}

impl_sevent_detector_delegate!(SAscendingNodeEvent<S, P>, inner, S, P);

/// Dynamic-sized ascending node event detector
pub struct DAscendingNodeEvent {
    inner: DValueEvent,
}

impl DAscendingNodeEvent {
    /// Create new ascending node event
    pub fn new(name: impl Into<String>) -> Self {
        let value_fn = |_t: Epoch, state: &DVector<f64>, _params: Option<&DVector<f64>>| {
            let state6 = Vector6::new(state[0], state[1], state[2], state[3], state[4], state[5]);
            let koe = state_eci_to_koe(state6, AngleFormat::Radians);
            let e = koe[1];
            let omega = koe[4];
            let m = koe[5];

            let nu = anomaly_mean_to_true(m, e, AngleFormat::Radians).unwrap_or(m);
            let u = omega + nu;
            (u + std::f64::consts::PI).rem_euclid(2.0 * std::f64::consts::PI) - std::f64::consts::PI
        };

        Self {
            inner: DValueEvent::new(name, value_fn, 0.0, EventDirection::Increasing),
        }
    }

    /// Set instance number for display name
    pub fn with_instance(mut self, instance: usize) -> Self {
        self.inner = self.inner.with_instance(instance);
        self
    }

    /// Set custom tolerances for event detection
    pub fn with_tolerances(mut self, time_tol: f64, value_tol: f64) -> Self {
        self.inner = self.inner.with_tolerances(time_tol, value_tol);
        self
    }

    /// Set custom step reduction factor for bisection search
    pub fn with_step_reduction_factor(mut self, factor: f64) -> Self {
        self.inner = self.inner.with_step_reduction_factor(factor);
        self
    }

    /// Set event callback
    pub fn with_callback(mut self, callback: DEventCallback) -> Self {
        self.inner = self.inner.with_callback(callback);
        self
    }

    /// Mark as terminal event
    pub fn set_terminal(mut self) -> Self {
        self.inner = self.inner.set_terminal();
        self
    }
}

impl_devent_detector_delegate!(DAscendingNodeEvent, inner);

/// Descending node event detector (static-sized)
///
/// Detects when spacecraft crosses the descending node (equator from north to south).
/// The descending node occurs when the true argument of latitude crosses π with increasing values.
///
/// # Examples
/// ```
/// use brahe::events::SDescendingNodeEvent;
///
/// // Detect descending node crossing
/// let event = SDescendingNodeEvent::<6, 0>::new("Descending Node");
/// ```
pub struct SDescendingNodeEvent<const S: usize, const P: usize> {
    inner: SValueEvent<S, P>,
}

impl<const S: usize, const P: usize> SDescendingNodeEvent<S, P> {
    /// Create new descending node event
    ///
    /// # Arguments
    /// * `name` - Event name
    pub fn new(name: impl Into<String>) -> Self {
        let value_fn = |_t: Epoch, state: &SVector<f64, S>, _params: Option<&SVector<f64, P>>| {
            let state6 = Vector6::new(state[0], state[1], state[2], state[3], state[4], state[5]);
            let koe = state_eci_to_koe(state6, AngleFormat::Radians);
            let e = koe[1];
            let omega = koe[4]; // Argument of perigee
            let m = koe[5]; // Mean anomaly

            // Convert mean anomaly to true anomaly
            let nu = anomaly_mean_to_true(m, e, AngleFormat::Radians).unwrap_or(m);

            // True argument of latitude = omega + nu
            let u = omega + nu;
            // Shift so that π becomes our zero crossing point
            // u - π gives us zero at descending node
            let u_shifted = u - std::f64::consts::PI;
            // Normalize to [-pi, pi) using rem_euclid for correct handling of negative values
            (u_shifted + std::f64::consts::PI).rem_euclid(2.0 * std::f64::consts::PI)
                - std::f64::consts::PI
        };

        Self {
            inner: SValueEvent::new(name, value_fn, 0.0, EventDirection::Increasing),
        }
    }

    /// Set instance number for display name
    pub fn with_instance(mut self, instance: usize) -> Self {
        self.inner = self.inner.with_instance(instance);
        self
    }

    /// Set custom tolerances for event detection
    pub fn with_tolerances(mut self, time_tol: f64, value_tol: f64) -> Self {
        self.inner = self.inner.with_tolerances(time_tol, value_tol);
        self
    }

    /// Set custom step reduction factor for bisection search
    pub fn with_step_reduction_factor(mut self, factor: f64) -> Self {
        self.inner = self.inner.with_step_reduction_factor(factor);
        self
    }

    /// Set event callback
    pub fn with_callback(mut self, callback: SEventCallback<S, P>) -> Self {
        self.inner = self.inner.with_callback(callback);
        self
    }

    /// Mark as terminal event (stops propagation)
    pub fn set_terminal(mut self) -> Self {
        self.inner = self.inner.set_terminal();
        self
    }
}

impl_sevent_detector_delegate!(SDescendingNodeEvent<S, P>, inner, S, P);

/// Dynamic-sized descending node event detector
pub struct DDescendingNodeEvent {
    inner: DValueEvent,
}

impl DDescendingNodeEvent {
    /// Create new descending node event
    pub fn new(name: impl Into<String>) -> Self {
        let value_fn = |_t: Epoch, state: &DVector<f64>, _params: Option<&DVector<f64>>| {
            let state6 = Vector6::new(state[0], state[1], state[2], state[3], state[4], state[5]);
            let koe = state_eci_to_koe(state6, AngleFormat::Radians);
            let e = koe[1];
            let omega = koe[4];
            let m = koe[5];

            let nu = anomaly_mean_to_true(m, e, AngleFormat::Radians).unwrap_or(m);
            let u = omega + nu;
            let u_shifted = u - std::f64::consts::PI;
            (u_shifted + std::f64::consts::PI).rem_euclid(2.0 * std::f64::consts::PI)
                - std::f64::consts::PI
        };

        Self {
            inner: DValueEvent::new(name, value_fn, 0.0, EventDirection::Increasing),
        }
    }

    /// Set instance number for display name
    pub fn with_instance(mut self, instance: usize) -> Self {
        self.inner = self.inner.with_instance(instance);
        self
    }

    /// Set custom tolerances for event detection
    pub fn with_tolerances(mut self, time_tol: f64, value_tol: f64) -> Self {
        self.inner = self.inner.with_tolerances(time_tol, value_tol);
        self
    }

    /// Set custom step reduction factor for bisection search
    pub fn with_step_reduction_factor(mut self, factor: f64) -> Self {
        self.inner = self.inner.with_step_reduction_factor(factor);
        self
    }

    /// Set event callback
    pub fn with_callback(mut self, callback: DEventCallback) -> Self {
        self.inner = self.inner.with_callback(callback);
        self
    }

    /// Mark as terminal event
    pub fn set_terminal(mut self) -> Self {
        self.inner = self.inner.set_terminal();
        self
    }
}

impl_devent_detector_delegate!(DDescendingNodeEvent, inner);

// =============================================================================
// State-Derived Events
// =============================================================================

/// Speed event detector (static-sized)
///
/// Detects when velocity magnitude crosses a value value.
///
/// # Examples
/// ```
/// use brahe::events::{SSpeedEvent, EventDirection};
///
/// // Detect when speed drops below 7 km/s
/// let event = SSpeedEvent::<6, 0>::new(
///     7000.0,
///     "Low Speed",
///     EventDirection::Decreasing
/// );
/// ```
pub struct SSpeedEvent<const S: usize, const P: usize> {
    inner: SValueEvent<S, P>,
}

impl<const S: usize, const P: usize> SSpeedEvent<S, P> {
    /// Create new speed event
    ///
    /// # Arguments
    /// * `value` - Speed value in m/s
    /// * `name` - Event name
    /// * `direction` - Detection direction
    pub fn new(value: f64, name: impl Into<String>, direction: EventDirection) -> Self {
        let value_fn = |_t: Epoch, state: &SVector<f64, S>, _params: Option<&SVector<f64, P>>| {
            // Compute velocity magnitude from state[3:6]
            let v = Vector3::new(state[3], state[4], state[5]);
            v.norm()
        };

        Self {
            inner: SValueEvent::new(name, value_fn, value, direction),
        }
    }

    /// Set instance number for display name
    pub fn with_instance(mut self, instance: usize) -> Self {
        self.inner = self.inner.with_instance(instance);
        self
    }

    /// Set custom tolerances for event detection
    pub fn with_tolerances(mut self, time_tol: f64, value_tol: f64) -> Self {
        self.inner = self.inner.with_tolerances(time_tol, value_tol);
        self
    }

    /// Set custom step reduction factor for bisection search
    pub fn with_step_reduction_factor(mut self, factor: f64) -> Self {
        self.inner = self.inner.with_step_reduction_factor(factor);
        self
    }

    /// Set event callback
    pub fn with_callback(mut self, callback: SEventCallback<S, P>) -> Self {
        self.inner = self.inner.with_callback(callback);
        self
    }

    /// Mark as terminal event (stops propagation)
    pub fn set_terminal(mut self) -> Self {
        self.inner = self.inner.set_terminal();
        self
    }
}

impl_sevent_detector_delegate!(SSpeedEvent<S, P>, inner, S, P);

/// Dynamic-sized speed event detector
pub struct DSpeedEvent {
    inner: DValueEvent,
}

impl DSpeedEvent {
    /// Create new speed event (value in m/s)
    pub fn new(value: f64, name: impl Into<String>, direction: EventDirection) -> Self {
        let value_fn = |_t: Epoch, state: &DVector<f64>, _params: Option<&DVector<f64>>| {
            let v = Vector3::new(state[3], state[4], state[5]);
            v.norm()
        };

        Self {
            inner: DValueEvent::new(name, value_fn, value, direction),
        }
    }

    /// Set instance number for display name
    pub fn with_instance(mut self, instance: usize) -> Self {
        self.inner = self.inner.with_instance(instance);
        self
    }

    /// Set custom tolerances for event detection
    pub fn with_tolerances(mut self, time_tol: f64, value_tol: f64) -> Self {
        self.inner = self.inner.with_tolerances(time_tol, value_tol);
        self
    }

    /// Set custom step reduction factor for bisection search
    pub fn with_step_reduction_factor(mut self, factor: f64) -> Self {
        self.inner = self.inner.with_step_reduction_factor(factor);
        self
    }

    /// Set event callback
    pub fn with_callback(mut self, callback: DEventCallback) -> Self {
        self.inner = self.inner.with_callback(callback);
        self
    }

    /// Mark as terminal event
    pub fn set_terminal(mut self) -> Self {
        self.inner = self.inner.set_terminal();
        self
    }
}

impl_devent_detector_delegate!(DSpeedEvent, inner);

/// Longitude event detector (static-sized)
///
/// Detects when geodetic longitude crosses a value value.
/// Requires EOP initialization for ECI->ECEF transformation.
///
/// # Examples
/// ```
/// use brahe::events::{SLongitudeEvent, EventDirection};
/// use brahe::AngleFormat;
///
/// // Detect when crossing prime meridian (0 degrees longitude)
/// let event = SLongitudeEvent::<6, 0>::new(
///     0.0,
///     "Prime Meridian",
///     EventDirection::Any,
///     AngleFormat::Degrees
/// );
/// ```
pub struct SLongitudeEvent<const S: usize, const P: usize> {
    inner: SValueEvent<S, P>,
}

impl<const S: usize, const P: usize> SLongitudeEvent<S, P> {
    /// Create new longitude event
    ///
    /// # Arguments
    /// * `value` - Longitude value
    /// * `name` - Event name
    /// * `direction` - Detection direction
    /// * `angle_format` - Whether value is in degrees or radians
    pub fn new(
        value: f64,
        name: impl Into<String>,
        direction: EventDirection,
        angle_format: AngleFormat,
    ) -> Self {
        let value_rad = match angle_format {
            AngleFormat::Degrees => value.to_radians(),
            AngleFormat::Radians => value,
        };

        let value_fn = |t: Epoch, state: &SVector<f64, S>, _params: Option<&SVector<f64, P>>| {
            let r_eci = state.fixed_rows::<3>(0).into_owned();
            let r_ecef = position_eci_to_ecef(t, r_eci);
            let geodetic = position_ecef_to_geodetic(r_ecef, AngleFormat::Radians);
            geodetic[0] // Longitude
        };

        Self {
            inner: SValueEvent::new(name, value_fn, value_rad, direction),
        }
    }

    /// Set instance number for display name
    pub fn with_instance(mut self, instance: usize) -> Self {
        self.inner = self.inner.with_instance(instance);
        self
    }

    /// Set custom tolerances for event detection
    pub fn with_tolerances(mut self, time_tol: f64, value_tol: f64) -> Self {
        self.inner = self.inner.with_tolerances(time_tol, value_tol);
        self
    }

    /// Set custom step reduction factor for bisection search
    pub fn with_step_reduction_factor(mut self, factor: f64) -> Self {
        self.inner = self.inner.with_step_reduction_factor(factor);
        self
    }

    /// Set event callback
    pub fn with_callback(mut self, callback: SEventCallback<S, P>) -> Self {
        self.inner = self.inner.with_callback(callback);
        self
    }

    /// Mark as terminal event (stops propagation)
    pub fn set_terminal(mut self) -> Self {
        self.inner = self.inner.set_terminal();
        self
    }
}

impl_sevent_detector_delegate!(SLongitudeEvent<S, P>, inner, S, P);

/// Dynamic-sized longitude event detector
pub struct DLongitudeEvent {
    inner: DValueEvent,
}

impl DLongitudeEvent {
    /// Create new longitude event
    ///
    /// # Arguments
    /// * `value` - Longitude value
    /// * `name` - Event name
    /// * `direction` - Detection direction
    /// * `angle_format` - Whether value is in degrees or radians
    pub fn new(
        value: f64,
        name: impl Into<String>,
        direction: EventDirection,
        angle_format: AngleFormat,
    ) -> Self {
        let value_rad = match angle_format {
            AngleFormat::Degrees => value.to_radians(),
            AngleFormat::Radians => value,
        };

        let value_fn = |t: Epoch, state: &DVector<f64>, _params: Option<&DVector<f64>>| {
            let r_eci = Vector3::new(state[0], state[1], state[2]);
            let r_ecef = position_eci_to_ecef(t, r_eci);
            let geodetic = position_ecef_to_geodetic(r_ecef, AngleFormat::Radians);
            geodetic[0]
        };

        Self {
            inner: DValueEvent::new(name, value_fn, value_rad, direction),
        }
    }

    /// Set instance number for display name
    pub fn with_instance(mut self, instance: usize) -> Self {
        self.inner = self.inner.with_instance(instance);
        self
    }

    /// Set custom tolerances for event detection
    pub fn with_tolerances(mut self, time_tol: f64, value_tol: f64) -> Self {
        self.inner = self.inner.with_tolerances(time_tol, value_tol);
        self
    }

    /// Set custom step reduction factor for bisection search
    pub fn with_step_reduction_factor(mut self, factor: f64) -> Self {
        self.inner = self.inner.with_step_reduction_factor(factor);
        self
    }

    /// Set event callback
    pub fn with_callback(mut self, callback: DEventCallback) -> Self {
        self.inner = self.inner.with_callback(callback);
        self
    }

    /// Mark as terminal event
    pub fn set_terminal(mut self) -> Self {
        self.inner = self.inner.set_terminal();
        self
    }
}

impl_devent_detector_delegate!(DLongitudeEvent, inner);

/// Latitude event detector (static-sized)
///
/// Detects when geodetic latitude crosses a value value.
/// Requires EOP initialization for ECI->ECEF transformation.
///
/// # Examples
/// ```
/// use brahe::events::{SLatitudeEvent, EventDirection};
/// use brahe::AngleFormat;
///
/// // Detect when crossing equator (0 degrees latitude)
/// let event = SLatitudeEvent::<6, 0>::new(
///     0.0,
///     "Equator Crossing",
///     EventDirection::Any,
///     AngleFormat::Degrees
/// );
/// ```
pub struct SLatitudeEvent<const S: usize, const P: usize> {
    inner: SValueEvent<S, P>,
}

impl<const S: usize, const P: usize> SLatitudeEvent<S, P> {
    /// Create new latitude event
    ///
    /// # Arguments
    /// * `value` - Latitude value
    /// * `name` - Event name
    /// * `direction` - Detection direction
    /// * `angle_format` - Whether value is in degrees or radians
    pub fn new(
        value: f64,
        name: impl Into<String>,
        direction: EventDirection,
        angle_format: AngleFormat,
    ) -> Self {
        let value_rad = match angle_format {
            AngleFormat::Degrees => value.to_radians(),
            AngleFormat::Radians => value,
        };

        let value_fn = |t: Epoch, state: &SVector<f64, S>, _params: Option<&SVector<f64, P>>| {
            let r_eci = state.fixed_rows::<3>(0).into_owned();
            let r_ecef = position_eci_to_ecef(t, r_eci);
            let geodetic = position_ecef_to_geodetic(r_ecef, AngleFormat::Radians);
            geodetic[1] // Latitude
        };

        Self {
            inner: SValueEvent::new(name, value_fn, value_rad, direction),
        }
    }

    /// Set instance number for display name
    pub fn with_instance(mut self, instance: usize) -> Self {
        self.inner = self.inner.with_instance(instance);
        self
    }

    /// Set custom tolerances for event detection
    pub fn with_tolerances(mut self, time_tol: f64, value_tol: f64) -> Self {
        self.inner = self.inner.with_tolerances(time_tol, value_tol);
        self
    }

    /// Set custom step reduction factor for bisection search
    pub fn with_step_reduction_factor(mut self, factor: f64) -> Self {
        self.inner = self.inner.with_step_reduction_factor(factor);
        self
    }

    /// Set event callback
    pub fn with_callback(mut self, callback: SEventCallback<S, P>) -> Self {
        self.inner = self.inner.with_callback(callback);
        self
    }

    /// Mark as terminal event (stops propagation)
    pub fn set_terminal(mut self) -> Self {
        self.inner = self.inner.set_terminal();
        self
    }
}

impl_sevent_detector_delegate!(SLatitudeEvent<S, P>, inner, S, P);

/// Dynamic-sized latitude event detector
pub struct DLatitudeEvent {
    inner: DValueEvent,
}

impl DLatitudeEvent {
    /// Create new latitude event
    ///
    /// # Arguments
    /// * `value` - Latitude value
    /// * `name` - Event name
    /// * `direction` - Detection direction
    /// * `angle_format` - Whether value is in degrees or radians
    pub fn new(
        value: f64,
        name: impl Into<String>,
        direction: EventDirection,
        angle_format: AngleFormat,
    ) -> Self {
        let value_rad = match angle_format {
            AngleFormat::Degrees => value.to_radians(),
            AngleFormat::Radians => value,
        };

        let value_fn = |t: Epoch, state: &DVector<f64>, _params: Option<&DVector<f64>>| {
            let r_eci = Vector3::new(state[0], state[1], state[2]);
            let r_ecef = position_eci_to_ecef(t, r_eci);
            let geodetic = position_ecef_to_geodetic(r_ecef, AngleFormat::Radians);
            geodetic[1]
        };

        Self {
            inner: DValueEvent::new(name, value_fn, value_rad, direction),
        }
    }

    /// Set instance number for display name
    pub fn with_instance(mut self, instance: usize) -> Self {
        self.inner = self.inner.with_instance(instance);
        self
    }

    /// Set custom tolerances for event detection
    pub fn with_tolerances(mut self, time_tol: f64, value_tol: f64) -> Self {
        self.inner = self.inner.with_tolerances(time_tol, value_tol);
        self
    }

    /// Set custom step reduction factor for bisection search
    pub fn with_step_reduction_factor(mut self, factor: f64) -> Self {
        self.inner = self.inner.with_step_reduction_factor(factor);
        self
    }

    /// Set event callback
    pub fn with_callback(mut self, callback: DEventCallback) -> Self {
        self.inner = self.inner.with_callback(callback);
        self
    }

    /// Mark as terminal event
    pub fn set_terminal(mut self) -> Self {
        self.inner = self.inner.set_terminal();
        self
    }
}

impl_devent_detector_delegate!(DLatitudeEvent, inner);

// =============================================================================
// Eclipse/Shadow Events
// =============================================================================

/// Helper function to get sun position based on ephemeris source
fn get_sun_position(t: Epoch, source: Option<EphemerisSource>) -> Vector3<f64> {
    match source {
        None | Some(EphemerisSource::LowPrecision) => sun_position(t),
        Some(source) => {
            let kernel =
                SPKKernel::try_from(source).expect("DE ephemeris source should map to a DE kernel");
            sun_position_de(t, kernel).unwrap_or_else(|_| sun_position(t))
        }
    }
}

/// Umbra event detector (static-sized)
///
/// Detects when spacecraft enters/exits Earth's umbra (full shadow).
/// Uses the conical shadow model.
///
/// # Examples
/// ```
/// use brahe::events::{SUmbraEvent, EdgeType};
/// use brahe::propagators::EphemerisSource;
///
/// // Detect umbra entry with low-precision ephemeris (default)
/// let event = SUmbraEvent::<6, 0>::new("Umbra Entry", EdgeType::RisingEdge, None);
///
/// // Detect umbra exit with high-precision DE440s ephemeris
/// let event = SUmbraEvent::<6, 0>::new("Umbra Exit", EdgeType::FallingEdge, Some(EphemerisSource::DE440s));
/// ```
pub struct SUmbraEvent<const S: usize, const P: usize> {
    inner: SBinaryEvent<S, P>,
}

impl<const S: usize, const P: usize> SUmbraEvent<S, P> {
    /// Create new umbra event
    ///
    /// # Arguments
    /// * `name` - Event name
    /// * `edge` - Edge type (RisingEdge = entering umbra, FallingEdge = exiting umbra)
    /// * `ephemeris_source` - Source for sun position (None for low-precision analytical)
    pub fn new(
        name: impl Into<String>,
        edge: EdgeType,
        ephemeris_source: Option<EphemerisSource>,
    ) -> Self {
        let condition_fn =
            move |t: Epoch, state: &SVector<f64, S>, _params: Option<&SVector<f64, P>>| {
                let r_eci = Vector3::new(state[0], state[1], state[2]);
                let r_sun = get_sun_position(t, ephemeris_source);
                let illumination = eclipse_conical(r_eci, r_sun);
                // In umbra when illumination is exactly 0.0
                illumination == 0.0
            };

        Self {
            inner: SBinaryEvent::new(name, condition_fn, edge),
        }
    }

    /// Set instance number for display name
    pub fn with_instance(mut self, instance: usize) -> Self {
        self.inner = self.inner.with_instance(instance);
        self
    }

    /// Set custom tolerances for event detection
    pub fn with_tolerances(mut self, time_tol: f64, value_tol: f64) -> Self {
        self.inner = self.inner.with_tolerances(time_tol, value_tol);
        self
    }

    /// Set custom step reduction factor for bisection search
    pub fn with_step_reduction_factor(mut self, factor: f64) -> Self {
        self.inner = self.inner.with_step_reduction_factor(factor);
        self
    }

    /// Set event callback
    pub fn with_callback(mut self, callback: SEventCallback<S, P>) -> Self {
        self.inner = self.inner.with_callback(callback);
        self
    }

    /// Mark as terminal event (stops propagation)
    pub fn set_terminal(mut self) -> Self {
        self.inner = self.inner.set_terminal();
        self
    }
}

impl_sevent_detector_delegate!(SUmbraEvent<S, P>, inner, S, P);

/// Dynamic-sized umbra event detector
pub struct DUmbraEvent {
    inner: DBinaryEvent,
}

impl DUmbraEvent {
    /// Create new umbra event
    ///
    /// # Arguments
    /// * `name` - Event name
    /// * `edge` - Edge type (RisingEdge = entering umbra, FallingEdge = exiting umbra)
    /// * `ephemeris_source` - Source for sun position (None for low-precision analytical)
    pub fn new(
        name: impl Into<String>,
        edge: EdgeType,
        ephemeris_source: Option<EphemerisSource>,
    ) -> Self {
        let condition_fn = move |t: Epoch, state: &DVector<f64>, _params: Option<&DVector<f64>>| {
            let r_eci = Vector3::new(state[0], state[1], state[2]);
            let r_sun = get_sun_position(t, ephemeris_source);
            let illumination = eclipse_conical(r_eci, r_sun);
            illumination == 0.0
        };

        Self {
            inner: DBinaryEvent::new(name, condition_fn, edge),
        }
    }

    /// Set instance number for display name
    pub fn with_instance(mut self, instance: usize) -> Self {
        self.inner = self.inner.with_instance(instance);
        self
    }

    /// Set custom tolerances for event detection
    pub fn with_tolerances(mut self, time_tol: f64, value_tol: f64) -> Self {
        self.inner = self.inner.with_tolerances(time_tol, value_tol);
        self
    }

    /// Set custom step reduction factor for bisection search
    pub fn with_step_reduction_factor(mut self, factor: f64) -> Self {
        self.inner = self.inner.with_step_reduction_factor(factor);
        self
    }

    /// Set event callback
    pub fn with_callback(mut self, callback: DEventCallback) -> Self {
        self.inner = self.inner.with_callback(callback);
        self
    }

    /// Mark as terminal event
    pub fn set_terminal(mut self) -> Self {
        self.inner = self.inner.set_terminal();
        self
    }
}

impl_devent_detector_delegate!(DUmbraEvent, inner);

/// Penumbra event detector (static-sized)
///
/// Detects when spacecraft enters/exits Earth's penumbra (partial shadow).
/// Uses the conical shadow model.
///
/// # Examples
/// ```
/// use brahe::events::{SPenumbraEvent, EdgeType};
/// use brahe::propagators::EphemerisSource;
///
/// // Detect penumbra entry with low-precision ephemeris
/// let event = SPenumbraEvent::<6, 0>::new("Penumbra Entry", EdgeType::RisingEdge, None);
///
/// // Detect penumbra with high-precision DE440s ephemeris
/// let event = SPenumbraEvent::<6, 0>::new("Penumbra Entry", EdgeType::RisingEdge, Some(EphemerisSource::DE440s));
/// ```
pub struct SPenumbraEvent<const S: usize, const P: usize> {
    inner: SBinaryEvent<S, P>,
}

impl<const S: usize, const P: usize> SPenumbraEvent<S, P> {
    /// Create new penumbra event
    ///
    /// # Arguments
    /// * `name` - Event name
    /// * `edge` - Edge type (RisingEdge = entering penumbra, FallingEdge = exiting penumbra)
    /// * `ephemeris_source` - Source for sun position (None for low-precision analytical)
    pub fn new(
        name: impl Into<String>,
        edge: EdgeType,
        ephemeris_source: Option<EphemerisSource>,
    ) -> Self {
        let condition_fn =
            move |t: Epoch, state: &SVector<f64, S>, _params: Option<&SVector<f64, P>>| {
                let r_eci = Vector3::new(state[0], state[1], state[2]);
                let r_sun = get_sun_position(t, ephemeris_source);
                let illumination = eclipse_conical(r_eci, r_sun);
                // In penumbra when 0 < illumination < 1
                illumination > 0.0 && illumination < 1.0
            };

        Self {
            inner: SBinaryEvent::new(name, condition_fn, edge),
        }
    }

    /// Set instance number for display name
    pub fn with_instance(mut self, instance: usize) -> Self {
        self.inner = self.inner.with_instance(instance);
        self
    }

    /// Set custom tolerances for event detection
    pub fn with_tolerances(mut self, time_tol: f64, value_tol: f64) -> Self {
        self.inner = self.inner.with_tolerances(time_tol, value_tol);
        self
    }

    /// Set custom step reduction factor for bisection search
    pub fn with_step_reduction_factor(mut self, factor: f64) -> Self {
        self.inner = self.inner.with_step_reduction_factor(factor);
        self
    }

    /// Set event callback
    pub fn with_callback(mut self, callback: SEventCallback<S, P>) -> Self {
        self.inner = self.inner.with_callback(callback);
        self
    }

    /// Mark as terminal event (stops propagation)
    pub fn set_terminal(mut self) -> Self {
        self.inner = self.inner.set_terminal();
        self
    }
}

impl_sevent_detector_delegate!(SPenumbraEvent<S, P>, inner, S, P);

/// Dynamic-sized penumbra event detector
pub struct DPenumbraEvent {
    inner: DBinaryEvent,
}

impl DPenumbraEvent {
    /// Create new penumbra event
    ///
    /// # Arguments
    /// * `name` - Event name
    /// * `edge` - Edge type (RisingEdge = entering penumbra, FallingEdge = exiting penumbra)
    /// * `ephemeris_source` - Source for sun position (None for low-precision analytical)
    pub fn new(
        name: impl Into<String>,
        edge: EdgeType,
        ephemeris_source: Option<EphemerisSource>,
    ) -> Self {
        let condition_fn = move |t: Epoch, state: &DVector<f64>, _params: Option<&DVector<f64>>| {
            let r_eci = Vector3::new(state[0], state[1], state[2]);
            let r_sun = get_sun_position(t, ephemeris_source);
            let illumination = eclipse_conical(r_eci, r_sun);
            illumination > 0.0 && illumination < 1.0
        };

        Self {
            inner: DBinaryEvent::new(name, condition_fn, edge),
        }
    }

    /// Set instance number for display name
    pub fn with_instance(mut self, instance: usize) -> Self {
        self.inner = self.inner.with_instance(instance);
        self
    }

    /// Set custom tolerances for event detection
    pub fn with_tolerances(mut self, time_tol: f64, value_tol: f64) -> Self {
        self.inner = self.inner.with_tolerances(time_tol, value_tol);
        self
    }

    /// Set custom step reduction factor for bisection search
    pub fn with_step_reduction_factor(mut self, factor: f64) -> Self {
        self.inner = self.inner.with_step_reduction_factor(factor);
        self
    }

    /// Set event callback
    pub fn with_callback(mut self, callback: DEventCallback) -> Self {
        self.inner = self.inner.with_callback(callback);
        self
    }

    /// Mark as terminal event
    pub fn set_terminal(mut self) -> Self {
        self.inner = self.inner.set_terminal();
        self
    }
}

impl_devent_detector_delegate!(DPenumbraEvent, inner);

/// Eclipse event detector (static-sized)
///
/// Detects when spacecraft enters/exits eclipse (either umbra or penumbra).
/// Uses the conical shadow model.
///
/// # Examples
/// ```
/// use brahe::events::{SEclipseEvent, EdgeType};
/// use brahe::propagators::EphemerisSource;
///
/// // Detect any eclipse entry with low-precision ephemeris
/// let event = SEclipseEvent::<6, 0>::new("Eclipse Entry", EdgeType::RisingEdge, None);
///
/// // Detect eclipse exit with high-precision DE440s ephemeris
/// let event = SEclipseEvent::<6, 0>::new("Eclipse Exit", EdgeType::FallingEdge, Some(EphemerisSource::DE440s));
/// ```
pub struct SEclipseEvent<const S: usize, const P: usize> {
    inner: SBinaryEvent<S, P>,
}

impl<const S: usize, const P: usize> SEclipseEvent<S, P> {
    /// Create new eclipse event
    ///
    /// # Arguments
    /// * `name` - Event name
    /// * `edge` - Edge type (RisingEdge = entering eclipse, FallingEdge = exiting eclipse)
    /// * `ephemeris_source` - Source for sun position (None for low-precision analytical)
    pub fn new(
        name: impl Into<String>,
        edge: EdgeType,
        ephemeris_source: Option<EphemerisSource>,
    ) -> Self {
        let condition_fn =
            move |t: Epoch, state: &SVector<f64, S>, _params: Option<&SVector<f64, P>>| {
                let r_eci = Vector3::new(state[0], state[1], state[2]);
                let r_sun = get_sun_position(t, ephemeris_source);
                let illumination = eclipse_conical(r_eci, r_sun);
                // In eclipse when illumination < 1.0 (either umbra or penumbra)
                illumination < 1.0
            };

        Self {
            inner: SBinaryEvent::new(name, condition_fn, edge),
        }
    }

    /// Set instance number for display name
    pub fn with_instance(mut self, instance: usize) -> Self {
        self.inner = self.inner.with_instance(instance);
        self
    }

    /// Set custom tolerances for event detection
    pub fn with_tolerances(mut self, time_tol: f64, value_tol: f64) -> Self {
        self.inner = self.inner.with_tolerances(time_tol, value_tol);
        self
    }

    /// Set custom step reduction factor for bisection search
    pub fn with_step_reduction_factor(mut self, factor: f64) -> Self {
        self.inner = self.inner.with_step_reduction_factor(factor);
        self
    }

    /// Set event callback
    pub fn with_callback(mut self, callback: SEventCallback<S, P>) -> Self {
        self.inner = self.inner.with_callback(callback);
        self
    }

    /// Mark as terminal event (stops propagation)
    pub fn set_terminal(mut self) -> Self {
        self.inner = self.inner.set_terminal();
        self
    }
}

impl_sevent_detector_delegate!(SEclipseEvent<S, P>, inner, S, P);

/// Dynamic-sized eclipse event detector
pub struct DEclipseEvent {
    inner: DBinaryEvent,
}

impl DEclipseEvent {
    /// Create new eclipse event
    ///
    /// # Arguments
    /// * `name` - Event name
    /// * `edge` - Edge type (RisingEdge = entering eclipse, FallingEdge = exiting eclipse)
    /// * `ephemeris_source` - Source for sun position (None for low-precision analytical)
    pub fn new(
        name: impl Into<String>,
        edge: EdgeType,
        ephemeris_source: Option<EphemerisSource>,
    ) -> Self {
        let condition_fn = move |t: Epoch, state: &DVector<f64>, _params: Option<&DVector<f64>>| {
            let r_eci = Vector3::new(state[0], state[1], state[2]);
            let r_sun = get_sun_position(t, ephemeris_source);
            let illumination = eclipse_conical(r_eci, r_sun);
            illumination < 1.0
        };

        Self {
            inner: DBinaryEvent::new(name, condition_fn, edge),
        }
    }

    /// Set instance number for display name
    pub fn with_instance(mut self, instance: usize) -> Self {
        self.inner = self.inner.with_instance(instance);
        self
    }

    /// Set custom tolerances for event detection
    pub fn with_tolerances(mut self, time_tol: f64, value_tol: f64) -> Self {
        self.inner = self.inner.with_tolerances(time_tol, value_tol);
        self
    }

    /// Set custom step reduction factor for bisection search
    pub fn with_step_reduction_factor(mut self, factor: f64) -> Self {
        self.inner = self.inner.with_step_reduction_factor(factor);
        self
    }

    /// Set event callback
    pub fn with_callback(mut self, callback: DEventCallback) -> Self {
        self.inner = self.inner.with_callback(callback);
        self
    }

    /// Mark as terminal event
    pub fn set_terminal(mut self) -> Self {
        self.inner = self.inner.set_terminal();
        self
    }
}

impl_devent_detector_delegate!(DEclipseEvent, inner);

/// Sunlit event detector (static-sized)
///
/// Detects when spacecraft enters/exits sunlight (fully illuminated).
/// Uses the conical shadow model.
///
/// # Examples
/// ```
/// use brahe::events::{SSunlitEvent, EdgeType};
/// use brahe::propagators::EphemerisSource;
///
/// // Detect entering sunlight with low-precision ephemeris
/// let event = SSunlitEvent::<6, 0>::new("Enter Sunlight", EdgeType::RisingEdge, None);
///
/// // Detect leaving sunlight with high-precision DE440s ephemeris
/// let event = SSunlitEvent::<6, 0>::new("Leave Sunlight", EdgeType::FallingEdge, Some(EphemerisSource::DE440s));
/// ```
pub struct SSunlitEvent<const S: usize, const P: usize> {
    inner: SBinaryEvent<S, P>,
}

impl<const S: usize, const P: usize> SSunlitEvent<S, P> {
    /// Create new sunlit event
    ///
    /// # Arguments
    /// * `name` - Event name
    /// * `edge` - Edge type (RisingEdge = entering sunlight, FallingEdge = leaving sunlight)
    /// * `ephemeris_source` - Source for sun position (None for low-precision analytical)
    pub fn new(
        name: impl Into<String>,
        edge: EdgeType,
        ephemeris_source: Option<EphemerisSource>,
    ) -> Self {
        let condition_fn =
            move |t: Epoch, state: &SVector<f64, S>, _params: Option<&SVector<f64, P>>| {
                let r_eci = Vector3::new(state[0], state[1], state[2]);
                let r_sun = get_sun_position(t, ephemeris_source);
                let illumination = eclipse_conical(r_eci, r_sun);
                // Sunlit when illumination == 1.0
                illumination == 1.0
            };

        Self {
            inner: SBinaryEvent::new(name, condition_fn, edge),
        }
    }

    /// Set instance number for display name
    pub fn with_instance(mut self, instance: usize) -> Self {
        self.inner = self.inner.with_instance(instance);
        self
    }

    /// Set custom tolerances for event detection
    pub fn with_tolerances(mut self, time_tol: f64, value_tol: f64) -> Self {
        self.inner = self.inner.with_tolerances(time_tol, value_tol);
        self
    }

    /// Set custom step reduction factor for bisection search
    pub fn with_step_reduction_factor(mut self, factor: f64) -> Self {
        self.inner = self.inner.with_step_reduction_factor(factor);
        self
    }

    /// Set event callback
    pub fn with_callback(mut self, callback: SEventCallback<S, P>) -> Self {
        self.inner = self.inner.with_callback(callback);
        self
    }

    /// Mark as terminal event (stops propagation)
    pub fn set_terminal(mut self) -> Self {
        self.inner = self.inner.set_terminal();
        self
    }
}

impl_sevent_detector_delegate!(SSunlitEvent<S, P>, inner, S, P);

/// Dynamic-sized sunlit event detector
pub struct DSunlitEvent {
    inner: DBinaryEvent,
}

impl DSunlitEvent {
    /// Create new sunlit event
    ///
    /// # Arguments
    /// * `name` - Event name
    /// * `edge` - Edge type (RisingEdge = entering sunlight, FallingEdge = leaving sunlight)
    /// * `ephemeris_source` - Source for sun position (None for low-precision analytical)
    pub fn new(
        name: impl Into<String>,
        edge: EdgeType,
        ephemeris_source: Option<EphemerisSource>,
    ) -> Self {
        let condition_fn = move |t: Epoch, state: &DVector<f64>, _params: Option<&DVector<f64>>| {
            let r_eci = Vector3::new(state[0], state[1], state[2]);
            let r_sun = get_sun_position(t, ephemeris_source);
            let illumination = eclipse_conical(r_eci, r_sun);
            illumination == 1.0
        };

        Self {
            inner: DBinaryEvent::new(name, condition_fn, edge),
        }
    }

    /// Set instance number for display name
    pub fn with_instance(mut self, instance: usize) -> Self {
        self.inner = self.inner.with_instance(instance);
        self
    }

    /// Set custom tolerances for event detection
    pub fn with_tolerances(mut self, time_tol: f64, value_tol: f64) -> Self {
        self.inner = self.inner.with_tolerances(time_tol, value_tol);
        self
    }

    /// Set custom step reduction factor for bisection search
    pub fn with_step_reduction_factor(mut self, factor: f64) -> Self {
        self.inner = self.inner.with_step_reduction_factor(factor);
        self
    }

    /// Set event callback
    pub fn with_callback(mut self, callback: DEventCallback) -> Self {
        self.inner = self.inner.with_callback(callback);
        self
    }

    /// Mark as terminal event
    pub fn set_terminal(mut self) -> Self {
        self.inner = self.inner.set_terminal();
        self
    }
}

impl_devent_detector_delegate!(DSunlitEvent, inner);

// =============================================================================
// AOI (Area of Interest) Event Detectors
// =============================================================================

/// Prepare AOI vertices for point-in-polygon testing.
///
/// Converts polygon vertices to (lon, lat) pairs in radians and pre-processes
/// for efficient point-in-polygon testing. This function is called once during
/// event detector construction.
fn prepare_aoi_vertices(
    vertices: &[(f64, f64)],
    angle_format: AngleFormat,
) -> Arc<Vec<(f64, f64)>> {
    let radians_vertices: Vec<(f64, f64)> = match angle_format {
        AngleFormat::Degrees => vertices
            .iter()
            .map(|(lon, lat)| (lon.to_radians(), lat.to_radians()))
            .collect(),
        AngleFormat::Radians => vertices.to_vec(),
    };

    Arc::new(radians_vertices)
}

/// Extract polygon vertices from PolygonLocation as (lon, lat) tuples in radians.
fn polygon_location_to_vertices(polygon: &PolygonLocation) -> Vec<(f64, f64)> {
    polygon
        .vertices()
        .iter()
        .map(|v| {
            // PolygonLocation stores vertices as [lon, lat, alt] in degrees
            (v.x.to_radians(), v.y.to_radians())
        })
        .collect()
}

/// AOI Entry event detector (static-sized).
///
/// Detects when a satellite's sub-satellite point enters a polygonal Area of Interest.
/// The sub-satellite point is computed by transforming the ECI position to geodetic
/// coordinates.
///
/// # Coordinate Handling
/// - Input coordinates can be in degrees or radians (specified via `AngleFormat`)
/// - Polygons crossing the anti-meridian (±180° longitude) are automatically handled
/// - Polygon should be closed (first vertex == last vertex), but auto-closes if not
///
/// # Event Triggering
/// - Triggers on `RisingEdge`: when point transitions from outside to inside the AOI
/// - Returns `true` when satellite is inside the AOI
///
/// # Examples
/// ```
/// use brahe::events::SAOIEntryEvent;
/// use brahe::constants::AngleFormat;
///
/// // Define AOI as [(lon, lat)] pairs in degrees
/// let aoi_vertices = vec![
///     (10.0, 50.0),
///     (20.0, 50.0),
///     (20.0, 55.0),
///     (10.0, 55.0),
///     (10.0, 50.0),  // Closed polygon
/// ];
///
/// let event = SAOIEntryEvent::<6, 0>::from_coordinates(
///     &aoi_vertices,
///     "Europe AOI Entry",
///     AngleFormat::Degrees,
/// );
/// ```
pub struct SAOIEntryEvent<const S: usize, const P: usize> {
    inner: SBinaryEvent<S, P>,
}

impl<const S: usize, const P: usize> SAOIEntryEvent<S, P> {
    /// Create an AOI entry event from a PolygonLocation.
    ///
    /// # Arguments
    /// * `polygon` - PolygonLocation defining the Area of Interest
    /// * `name` - Event name for identification
    ///
    /// # Example
    /// ```
    /// use brahe::events::SAOIEntryEvent;
    /// use brahe::access::PolygonLocation;
    /// use nalgebra::Vector3;
    ///
    /// let vertices = vec![
    ///     Vector3::new(10.0, 50.0, 0.0),
    ///     Vector3::new(20.0, 50.0, 0.0),
    ///     Vector3::new(20.0, 55.0, 0.0),
    ///     Vector3::new(10.0, 55.0, 0.0),
    ///     Vector3::new(10.0, 50.0, 0.0),
    /// ];
    /// let polygon = PolygonLocation::new(vertices).unwrap();
    ///
    /// let event = SAOIEntryEvent::<6, 0>::from_polygon(&polygon, "AOI Entry");
    /// ```
    pub fn from_polygon(polygon: &PolygonLocation, name: impl Into<String>) -> Self {
        let vertices = polygon_location_to_vertices(polygon);
        let vertices_arc = Arc::new(vertices);

        let condition_fn =
            move |t: Epoch, state: &SVector<f64, S>, _params: Option<&SVector<f64, P>>| {
                // Extract ECI position from state
                let r_eci = Vector3::new(state[0], state[1], state[2]);

                // Transform ECI -> ECEF -> Geodetic
                let r_ecef = position_eci_to_ecef(t, r_eci);
                let geodetic = position_ecef_to_geodetic(r_ecef, AngleFormat::Radians);

                // geodetic is [lon, lat, alt] in radians
                let lon = geodetic[0];
                let lat = geodetic[1];

                // Test if sub-satellite point is inside AOI
                point_in_polygon(lon, lat, &vertices_arc)
            };

        Self {
            inner: SBinaryEvent::new(name, condition_fn, EdgeType::RisingEdge),
        }
    }

    /// Create an AOI entry event from raw coordinate pairs.
    ///
    /// # Arguments
    /// * `vertices` - AOI vertices as (longitude, latitude) pairs
    /// * `name` - Event name for identification
    /// * `angle_format` - Whether coordinates are in degrees or radians
    ///
    /// # Example
    /// ```
    /// use brahe::events::SAOIEntryEvent;
    /// use brahe::constants::AngleFormat;
    ///
    /// let vertices = vec![
    ///     (10.0, 50.0),
    ///     (20.0, 50.0),
    ///     (20.0, 55.0),
    ///     (10.0, 55.0),
    ///     (10.0, 50.0),
    /// ];
    ///
    /// let event = SAOIEntryEvent::<6, 0>::from_coordinates(
    ///     &vertices,
    ///     "AOI Entry",
    ///     AngleFormat::Degrees,
    /// );
    /// ```
    pub fn from_coordinates(
        vertices: &[(f64, f64)],
        name: impl Into<String>,
        angle_format: AngleFormat,
    ) -> Self {
        let vertices_arc = prepare_aoi_vertices(vertices, angle_format);

        let condition_fn =
            move |t: Epoch, state: &SVector<f64, S>, _params: Option<&SVector<f64, P>>| {
                let r_eci = Vector3::new(state[0], state[1], state[2]);
                let r_ecef = position_eci_to_ecef(t, r_eci);
                let geodetic = position_ecef_to_geodetic(r_ecef, AngleFormat::Radians);

                let lon = geodetic[0];
                let lat = geodetic[1];

                point_in_polygon(lon, lat, &vertices_arc)
            };

        Self {
            inner: SBinaryEvent::new(name, condition_fn, EdgeType::RisingEdge),
        }
    }

    /// Set instance number for display name.
    pub fn with_instance(mut self, instance: usize) -> Self {
        self.inner = self.inner.with_instance(instance);
        self
    }

    /// Set custom tolerances for event detection.
    pub fn with_tolerances(mut self, time_tol: f64, value_tol: f64) -> Self {
        self.inner = self.inner.with_tolerances(time_tol, value_tol);
        self
    }

    /// Set custom step reduction factor for bisection search.
    pub fn with_step_reduction_factor(mut self, factor: f64) -> Self {
        self.inner = self.inner.with_step_reduction_factor(factor);
        self
    }

    /// Set event callback.
    pub fn with_callback(mut self, callback: SEventCallback<S, P>) -> Self {
        self.inner = self.inner.with_callback(callback);
        self
    }

    /// Mark as terminal event (stops propagation).
    pub fn set_terminal(mut self) -> Self {
        self.inner = self.inner.set_terminal();
        self
    }
}

impl_sevent_detector_delegate!(SAOIEntryEvent<S, P>, inner, S, P);

/// Dynamic-sized AOI Entry event detector.
///
/// See [`SAOIEntryEvent`] for detailed documentation.
pub struct DAOIEntryEvent {
    inner: DBinaryEvent,
}

impl DAOIEntryEvent {
    /// Create an AOI entry event from a PolygonLocation.
    pub fn from_polygon(polygon: &PolygonLocation, name: impl Into<String>) -> Self {
        let vertices = polygon_location_to_vertices(polygon);
        let vertices_arc = Arc::new(vertices);

        let condition_fn = move |t: Epoch, state: &DVector<f64>, _params: Option<&DVector<f64>>| {
            let r_eci = Vector3::new(state[0], state[1], state[2]);
            let r_ecef = position_eci_to_ecef(t, r_eci);
            let geodetic = position_ecef_to_geodetic(r_ecef, AngleFormat::Radians);

            let lon = geodetic[0];
            let lat = geodetic[1];

            point_in_polygon(lon, lat, &vertices_arc)
        };

        Self {
            inner: DBinaryEvent::new(name, condition_fn, EdgeType::RisingEdge),
        }
    }

    /// Create an AOI entry event from raw coordinate pairs.
    pub fn from_coordinates(
        vertices: &[(f64, f64)],
        name: impl Into<String>,
        angle_format: AngleFormat,
    ) -> Self {
        let vertices_arc = prepare_aoi_vertices(vertices, angle_format);

        let condition_fn = move |t: Epoch, state: &DVector<f64>, _params: Option<&DVector<f64>>| {
            let r_eci = Vector3::new(state[0], state[1], state[2]);
            let r_ecef = position_eci_to_ecef(t, r_eci);
            let geodetic = position_ecef_to_geodetic(r_ecef, AngleFormat::Radians);

            let lon = geodetic[0];
            let lat = geodetic[1];

            point_in_polygon(lon, lat, &vertices_arc)
        };

        Self {
            inner: DBinaryEvent::new(name, condition_fn, EdgeType::RisingEdge),
        }
    }

    /// Set instance number for display name.
    pub fn with_instance(mut self, instance: usize) -> Self {
        self.inner = self.inner.with_instance(instance);
        self
    }

    /// Set custom tolerances for event detection.
    pub fn with_tolerances(mut self, time_tol: f64, value_tol: f64) -> Self {
        self.inner = self.inner.with_tolerances(time_tol, value_tol);
        self
    }

    /// Set custom step reduction factor for bisection search.
    pub fn with_step_reduction_factor(mut self, factor: f64) -> Self {
        self.inner = self.inner.with_step_reduction_factor(factor);
        self
    }

    /// Set event callback.
    pub fn with_callback(mut self, callback: DEventCallback) -> Self {
        self.inner = self.inner.with_callback(callback);
        self
    }

    /// Mark as terminal event.
    pub fn set_terminal(mut self) -> Self {
        self.inner = self.inner.set_terminal();
        self
    }
}

impl_devent_detector_delegate!(DAOIEntryEvent, inner);

/// AOI Exit event detector (static-sized).
///
/// Detects when a satellite's sub-satellite point exits a polygonal Area of Interest.
/// The sub-satellite point is computed by transforming the ECI position to geodetic
/// coordinates.
///
/// # Coordinate Handling
/// - Input coordinates can be in degrees or radians (specified via `AngleFormat`)
/// - Polygons crossing the anti-meridian (±180° longitude) are automatically handled
/// - Polygon should be closed (first vertex == last vertex), but auto-closes if not
///
/// # Event Triggering
/// - Triggers on `FallingEdge`: when point transitions from inside to outside the AOI
/// - Returns `true` when satellite is inside the AOI
///
/// # Examples
/// ```
/// use brahe::events::SAOIExitEvent;
/// use brahe::constants::AngleFormat;
///
/// let aoi_vertices = vec![
///     (10.0, 50.0),
///     (20.0, 50.0),
///     (20.0, 55.0),
///     (10.0, 55.0),
///     (10.0, 50.0),
/// ];
///
/// let event = SAOIExitEvent::<6, 0>::from_coordinates(
///     &aoi_vertices,
///     "Europe AOI Exit",
///     AngleFormat::Degrees,
/// );
/// ```
pub struct SAOIExitEvent<const S: usize, const P: usize> {
    inner: SBinaryEvent<S, P>,
}

impl<const S: usize, const P: usize> SAOIExitEvent<S, P> {
    /// Create an AOI exit event from a PolygonLocation.
    ///
    /// # Arguments
    /// * `polygon` - PolygonLocation defining the Area of Interest
    /// * `name` - Event name for identification
    pub fn from_polygon(polygon: &PolygonLocation, name: impl Into<String>) -> Self {
        let vertices = polygon_location_to_vertices(polygon);
        let vertices_arc = Arc::new(vertices);

        let condition_fn =
            move |t: Epoch, state: &SVector<f64, S>, _params: Option<&SVector<f64, P>>| {
                let r_eci = Vector3::new(state[0], state[1], state[2]);
                let r_ecef = position_eci_to_ecef(t, r_eci);
                let geodetic = position_ecef_to_geodetic(r_ecef, AngleFormat::Radians);

                let lon = geodetic[0];
                let lat = geodetic[1];

                point_in_polygon(lon, lat, &vertices_arc)
            };

        Self {
            inner: SBinaryEvent::new(name, condition_fn, EdgeType::FallingEdge),
        }
    }

    /// Create an AOI exit event from raw coordinate pairs.
    ///
    /// # Arguments
    /// * `vertices` - AOI vertices as (longitude, latitude) pairs
    /// * `name` - Event name for identification
    /// * `angle_format` - Whether coordinates are in degrees or radians
    pub fn from_coordinates(
        vertices: &[(f64, f64)],
        name: impl Into<String>,
        angle_format: AngleFormat,
    ) -> Self {
        let vertices_arc = prepare_aoi_vertices(vertices, angle_format);

        let condition_fn =
            move |t: Epoch, state: &SVector<f64, S>, _params: Option<&SVector<f64, P>>| {
                let r_eci = Vector3::new(state[0], state[1], state[2]);
                let r_ecef = position_eci_to_ecef(t, r_eci);
                let geodetic = position_ecef_to_geodetic(r_ecef, AngleFormat::Radians);

                let lon = geodetic[0];
                let lat = geodetic[1];

                point_in_polygon(lon, lat, &vertices_arc)
            };

        Self {
            inner: SBinaryEvent::new(name, condition_fn, EdgeType::FallingEdge),
        }
    }

    /// Set instance number for display name.
    pub fn with_instance(mut self, instance: usize) -> Self {
        self.inner = self.inner.with_instance(instance);
        self
    }

    /// Set custom tolerances for event detection.
    pub fn with_tolerances(mut self, time_tol: f64, value_tol: f64) -> Self {
        self.inner = self.inner.with_tolerances(time_tol, value_tol);
        self
    }

    /// Set custom step reduction factor for bisection search.
    pub fn with_step_reduction_factor(mut self, factor: f64) -> Self {
        self.inner = self.inner.with_step_reduction_factor(factor);
        self
    }

    /// Set event callback.
    pub fn with_callback(mut self, callback: SEventCallback<S, P>) -> Self {
        self.inner = self.inner.with_callback(callback);
        self
    }

    /// Mark as terminal event (stops propagation).
    pub fn set_terminal(mut self) -> Self {
        self.inner = self.inner.set_terminal();
        self
    }
}

impl_sevent_detector_delegate!(SAOIExitEvent<S, P>, inner, S, P);

/// Dynamic-sized AOI Exit event detector.
///
/// See [`SAOIExitEvent`] for detailed documentation.
pub struct DAOIExitEvent {
    inner: DBinaryEvent,
}

impl DAOIExitEvent {
    /// Create an AOI exit event from a PolygonLocation.
    pub fn from_polygon(polygon: &PolygonLocation, name: impl Into<String>) -> Self {
        let vertices = polygon_location_to_vertices(polygon);
        let vertices_arc = Arc::new(vertices);

        let condition_fn = move |t: Epoch, state: &DVector<f64>, _params: Option<&DVector<f64>>| {
            let r_eci = Vector3::new(state[0], state[1], state[2]);
            let r_ecef = position_eci_to_ecef(t, r_eci);
            let geodetic = position_ecef_to_geodetic(r_ecef, AngleFormat::Radians);

            let lon = geodetic[0];
            let lat = geodetic[1];

            point_in_polygon(lon, lat, &vertices_arc)
        };

        Self {
            inner: DBinaryEvent::new(name, condition_fn, EdgeType::FallingEdge),
        }
    }

    /// Create an AOI exit event from raw coordinate pairs.
    pub fn from_coordinates(
        vertices: &[(f64, f64)],
        name: impl Into<String>,
        angle_format: AngleFormat,
    ) -> Self {
        let vertices_arc = prepare_aoi_vertices(vertices, angle_format);

        let condition_fn = move |t: Epoch, state: &DVector<f64>, _params: Option<&DVector<f64>>| {
            let r_eci = Vector3::new(state[0], state[1], state[2]);
            let r_ecef = position_eci_to_ecef(t, r_eci);
            let geodetic = position_ecef_to_geodetic(r_ecef, AngleFormat::Radians);

            let lon = geodetic[0];
            let lat = geodetic[1];

            point_in_polygon(lon, lat, &vertices_arc)
        };

        Self {
            inner: DBinaryEvent::new(name, condition_fn, EdgeType::FallingEdge),
        }
    }

    /// Set instance number for display name.
    pub fn with_instance(mut self, instance: usize) -> Self {
        self.inner = self.inner.with_instance(instance);
        self
    }

    /// Set custom tolerances for event detection.
    pub fn with_tolerances(mut self, time_tol: f64, value_tol: f64) -> Self {
        self.inner = self.inner.with_tolerances(time_tol, value_tol);
        self
    }

    /// Set custom step reduction factor for bisection search.
    pub fn with_step_reduction_factor(mut self, factor: f64) -> Self {
        self.inner = self.inner.with_step_reduction_factor(factor);
        self
    }

    /// Set event callback.
    pub fn with_callback(mut self, callback: DEventCallback) -> Self {
        self.inner = self.inner.with_callback(callback);
        self
    }

    /// Mark as terminal event.
    pub fn set_terminal(mut self) -> Self {
        self.inner = self.inner.set_terminal();
        self
    }
}

impl_devent_detector_delegate!(DAOIExitEvent, inner);

#[cfg(test)]
#[allow(non_snake_case)]
#[cfg_attr(coverage_nightly, coverage(off))]
mod tests {
    use super::*;
    use crate::constants::R_EARTH;
    use crate::time::TimeSystem;
    use crate::utils::testing::setup_global_test_eop;
    use nalgebra::{DVector, Vector3, Vector6};
    use serial_test::serial;
    use std::sync::Arc;
    use std::sync::atomic::{AtomicBool, Ordering};

    // =========================================================================
    // SAltitudeEvent Tests
    // =========================================================================

    #[test]
    fn test_SAltitudeEvent_new() {
        setup_global_test_eop();

        let event = SAltitudeEvent::<6, 0>::new(500e3, "Low Altitude", EventDirection::Decreasing);

        assert_eq!(event.name(), "Low Altitude");
        assert_eq!(event.target_value(), 500e3);
        assert_eq!(event.direction(), EventDirection::Decreasing);
        assert_eq!(event.action(), EventAction::Continue);
    }

    #[test]
    #[serial]
    fn test_SAltitudeEvent_evaluate() {
        setup_global_test_eop();

        let value = 500e3;
        let event = SAltitudeEvent::<6, 0>::new(value, "Altitude Test", EventDirection::Any);

        let epoch = Epoch::from_jd(2451545.0, TimeSystem::UTC);

        // Low altitude state (below value) - ECI position along x-axis
        let state_low = Vector6::new(6000e3, 0.0, 0.0, 0.0, 7.5e3, 0.0);
        let val_low = event.evaluate(epoch, &state_low, None);
        // Should be negative (altitude < value)
        assert!(val_low < 0.0);

        // High altitude state (above value) - ECI position along x-axis
        let state_high = Vector6::new(R_EARTH + 1000e3, 0.0, 0.0, 0.0, 7.0e3, 0.0);
        let val_high = event.evaluate(epoch, &state_high, None);
        // Should be positive (altitude > value)
        assert!(val_high > 0.0);
    }

    #[test]
    fn test_SAltitudeEvent_target_value() {
        setup_global_test_eop();

        let event = SAltitudeEvent::<6, 0>::new(350e3, "Reentry", EventDirection::Decreasing);
        assert_eq!(event.target_value(), 350e3);
    }

    #[test]
    fn test_SAltitudeEvent_name() {
        setup_global_test_eop();

        let event = SAltitudeEvent::<6, 0>::new(500e3, "Test Name", EventDirection::Any);
        assert_eq!(event.name(), "Test Name");
    }

    #[test]
    fn test_SAltitudeEvent_with_instance() {
        setup_global_test_eop();

        let event = SAltitudeEvent::<6, 0>::new(500e3, "Altitude Check", EventDirection::Any)
            .with_instance(3);

        assert_eq!(event.name(), "Altitude Check 3");
    }

    #[test]
    fn test_SAltitudeEvent_with_tolerances() {
        setup_global_test_eop();

        // Default tolerances from SValueEvent
        let event = SAltitudeEvent::<6, 0>::new(500e3, "Test", EventDirection::Any);
        assert_eq!(event.time_tolerance(), 1e-6);
        assert_eq!(event.value_tolerance(), 1e-9);

        // Custom tolerances
        let event = SAltitudeEvent::<6, 0>::new(500e3, "Test", EventDirection::Any)
            .with_tolerances(1e-4, 1e-6);

        assert_eq!(event.time_tolerance(), 1e-4);
        assert_eq!(event.value_tolerance(), 1e-6);
    }

    #[test]
    fn test_SAltitudeEvent_time_tolerance() {
        setup_global_test_eop();

        let event = SAltitudeEvent::<6, 0>::new(500e3, "Test", EventDirection::Any)
            .with_tolerances(5e-5, 1e-9);

        assert_eq!(event.time_tolerance(), 5e-5);
    }

    #[test]
    fn test_SAltitudeEvent_value_tolerance() {
        setup_global_test_eop();

        let event = SAltitudeEvent::<6, 0>::new(500e3, "Test", EventDirection::Any)
            .with_tolerances(1e-6, 5e-8);

        assert_eq!(event.value_tolerance(), 5e-8);
    }

    #[test]
    fn test_SAltitudeEvent_with_step_reduction_factor() {
        setup_global_test_eop();

        // Default
        let event = SAltitudeEvent::<6, 0>::new(500e3, "Test", EventDirection::Any);
        assert_eq!(event.step_reduction_factor(), 0.2);

        // Custom
        let event = SAltitudeEvent::<6, 0>::new(500e3, "Test", EventDirection::Any)
            .with_step_reduction_factor(0.15);
        assert_eq!(event.step_reduction_factor(), 0.15);
    }

    #[test]
    fn test_SAltitudeEvent_step_reduction_factor() {
        setup_global_test_eop();

        let event = SAltitudeEvent::<6, 0>::new(500e3, "Test", EventDirection::Any)
            .with_step_reduction_factor(0.1);

        assert_eq!(event.step_reduction_factor(), 0.1);
    }

    #[test]
    fn test_SAltitudeEvent_with_callback() {
        setup_global_test_eop();

        let called = Arc::new(AtomicBool::new(false));
        let called_clone = called.clone();

        let callback: SEventCallback<6, 0> = Box::new(move |_t, _state, _params| {
            called_clone.store(true, Ordering::SeqCst);
            (None, None, EventAction::Continue)
        });

        let event =
            SAltitudeEvent::<6, 0>::new(500e3, "Test", EventDirection::Any).with_callback(callback);

        // Callback should exist
        assert!(event.callback().is_some());

        // Execute callback
        let epoch = Epoch::from_jd(2451545.0, TimeSystem::UTC);
        let state = Vector6::zeros();
        if let Some(cb) = event.callback() {
            cb(epoch, &state, None);
        }
        assert!(called.load(Ordering::SeqCst));
    }

    #[test]
    fn test_SAltitudeEvent_callback_none() {
        setup_global_test_eop();

        let event = SAltitudeEvent::<6, 0>::new(500e3, "Test", EventDirection::Any);
        assert!(event.callback().is_none());
    }

    #[test]
    fn test_SAltitudeEvent_set_terminal() {
        setup_global_test_eop();

        let event = SAltitudeEvent::<6, 0>::new(500e3, "Test", EventDirection::Any);
        assert_eq!(event.action(), EventAction::Continue);

        let event = SAltitudeEvent::<6, 0>::new(500e3, "Test", EventDirection::Any).set_terminal();
        assert_eq!(event.action(), EventAction::Stop);
    }

    #[test]
    fn test_SAltitudeEvent_action_continue() {
        setup_global_test_eop();

        let event = SAltitudeEvent::<6, 0>::new(500e3, "Test", EventDirection::Any);
        assert_eq!(event.action(), EventAction::Continue);
    }

    #[test]
    fn test_SAltitudeEvent_action_stop() {
        setup_global_test_eop();

        let event = SAltitudeEvent::<6, 0>::new(500e3, "Test", EventDirection::Any).set_terminal();
        assert_eq!(event.action(), EventAction::Stop);
    }

    #[test]
    fn test_SAltitudeEvent_direction_increasing() {
        setup_global_test_eop();

        let event = SAltitudeEvent::<6, 0>::new(500e3, "Ascending", EventDirection::Increasing);
        assert_eq!(event.direction(), EventDirection::Increasing);
    }

    #[test]
    fn test_SAltitudeEvent_direction_decreasing() {
        setup_global_test_eop();

        let event = SAltitudeEvent::<6, 0>::new(500e3, "Descending", EventDirection::Decreasing);
        assert_eq!(event.direction(), EventDirection::Decreasing);
    }

    #[test]
    fn test_SAltitudeEvent_direction_any() {
        setup_global_test_eop();

        let event = SAltitudeEvent::<6, 0>::new(500e3, "Any", EventDirection::Any);
        assert_eq!(event.direction(), EventDirection::Any);
    }

    #[test]
    #[serial]
    fn test_SAltitudeEvent_7d_state() {
        setup_global_test_eop();

        // Test with 7D state (e.g., with mass)
        let event = SAltitudeEvent::<7, 4>::new(500e3, "7D State", EventDirection::Decreasing);

        let epoch = Epoch::from_jd(2451545.0, TimeSystem::UTC);
        let state: SVector<f64, 7> =
            SVector::from([R_EARTH + 600e3, 0.0, 0.0, 0.0, 7.5e3, 0.0, 1000.0]);

        let altitude = event.evaluate(epoch, &state, None);
        // Altitude should be approximately 600 km
        assert!(altitude > 500e3);
    }

    #[test]
    fn test_SAltitudeEvent_builder_chaining() {
        setup_global_test_eop();

        let callback: SEventCallback<6, 0> =
            Box::new(|_t, _state, _params| (None, None, EventAction::Stop));

        let event = SAltitudeEvent::<6, 0>::new(200e3, "Reentry", EventDirection::Decreasing)
            .with_instance(1)
            .with_tolerances(1e-5, 1e-8)
            .with_step_reduction_factor(0.1)
            .with_callback(callback)
            .set_terminal();

        assert_eq!(event.name(), "Reentry 1");
        assert_eq!(event.time_tolerance(), 1e-5);
        assert_eq!(event.value_tolerance(), 1e-8);
        assert_eq!(event.step_reduction_factor(), 0.1);
        assert!(event.callback().is_some());
        assert_eq!(event.action(), EventAction::Stop);
    }

    // =========================================================================
    // DAltitudeEvent Tests
    // =========================================================================

    #[test]
    fn test_DAltitudeEvent_new() {
        setup_global_test_eop();

        let event = DAltitudeEvent::new(500e3, "Low Altitude", EventDirection::Decreasing);

        assert_eq!(event.name(), "Low Altitude");
        assert_eq!(event.target_value(), 500e3);
        assert_eq!(event.direction(), EventDirection::Decreasing);
        assert_eq!(event.action(), EventAction::Continue);
    }

    #[test]
    #[serial]
    fn test_DAltitudeEvent_evaluate() {
        setup_global_test_eop();

        let value = 500e3;
        let event = DAltitudeEvent::new(value, "Altitude Test", EventDirection::Any);

        let epoch = Epoch::from_jd(2451545.0, TimeSystem::UTC);

        // Low altitude state (below value)
        let state_low = DVector::from_vec(vec![6000e3, 0.0, 0.0, 0.0, 7.5e3, 0.0]);
        let val_low = event.evaluate(epoch, &state_low, None);
        assert!(val_low < 0.0);

        // High altitude state (above value)
        let state_high = DVector::from_vec(vec![R_EARTH + 1000e3, 0.0, 0.0, 0.0, 7.0e3, 0.0]);
        let val_high = event.evaluate(epoch, &state_high, None);
        assert!(val_high > 0.0);
    }

    #[test]
    fn test_DAltitudeEvent_target_value() {
        setup_global_test_eop();

        let event = DAltitudeEvent::new(350e3, "Reentry", EventDirection::Decreasing);
        assert_eq!(event.target_value(), 350e3);
    }

    #[test]
    fn test_DAltitudeEvent_name() {
        setup_global_test_eop();

        let event = DAltitudeEvent::new(500e3, "Test Name", EventDirection::Any);
        assert_eq!(event.name(), "Test Name");
    }

    #[test]
    fn test_DAltitudeEvent_with_instance() {
        setup_global_test_eop();

        let event =
            DAltitudeEvent::new(500e3, "Altitude Check", EventDirection::Any).with_instance(2);

        assert_eq!(event.name(), "Altitude Check 2");
    }

    #[test]
    fn test_DAltitudeEvent_with_tolerances() {
        setup_global_test_eop();

        let event =
            DAltitudeEvent::new(500e3, "Test", EventDirection::Any).with_tolerances(1e-4, 1e-7);

        assert_eq!(event.time_tolerance(), 1e-4);
        assert_eq!(event.value_tolerance(), 1e-7);
    }

    #[test]
    fn test_DAltitudeEvent_time_tolerance() {
        setup_global_test_eop();

        let event =
            DAltitudeEvent::new(500e3, "Test", EventDirection::Any).with_tolerances(2e-5, 1e-9);

        assert_eq!(event.time_tolerance(), 2e-5);
    }

    #[test]
    fn test_DAltitudeEvent_value_tolerance() {
        setup_global_test_eop();

        let event =
            DAltitudeEvent::new(500e3, "Test", EventDirection::Any).with_tolerances(1e-6, 3e-8);

        assert_eq!(event.value_tolerance(), 3e-8);
    }

    #[test]
    fn test_DAltitudeEvent_with_step_reduction_factor() {
        setup_global_test_eop();

        // Default
        let event = DAltitudeEvent::new(500e3, "Test", EventDirection::Any);
        assert_eq!(event.step_reduction_factor(), 0.2);

        // Custom
        let event = DAltitudeEvent::new(500e3, "Test", EventDirection::Any)
            .with_step_reduction_factor(0.25);
        assert_eq!(event.step_reduction_factor(), 0.25);
    }

    #[test]
    fn test_DAltitudeEvent_step_reduction_factor() {
        setup_global_test_eop();

        let event = DAltitudeEvent::new(500e3, "Test", EventDirection::Any)
            .with_step_reduction_factor(0.12);

        assert_eq!(event.step_reduction_factor(), 0.12);
    }

    #[test]
    fn test_DAltitudeEvent_with_callback() {
        setup_global_test_eop();

        let called = Arc::new(AtomicBool::new(false));
        let called_clone = called.clone();

        let callback: DEventCallback = Box::new(move |_t, _state, _params| {
            called_clone.store(true, Ordering::SeqCst);
            (None, None, EventAction::Continue)
        });

        let event = DAltitudeEvent::new(500e3, "Test", EventDirection::Any).with_callback(callback);

        assert!(event.callback().is_some());
    }

    #[test]
    fn test_DAltitudeEvent_callback_none() {
        setup_global_test_eop();

        let event = DAltitudeEvent::new(500e3, "Test", EventDirection::Any);
        assert!(event.callback().is_none());
    }

    #[test]
    fn test_DAltitudeEvent_set_terminal() {
        setup_global_test_eop();

        let event = DAltitudeEvent::new(500e3, "Test", EventDirection::Any);
        assert_eq!(event.action(), EventAction::Continue);

        let event = DAltitudeEvent::new(500e3, "Test", EventDirection::Any).set_terminal();
        assert_eq!(event.action(), EventAction::Stop);
    }

    #[test]
    fn test_DAltitudeEvent_action_continue() {
        setup_global_test_eop();

        let event = DAltitudeEvent::new(500e3, "Test", EventDirection::Any);
        assert_eq!(event.action(), EventAction::Continue);
    }

    #[test]
    fn test_DAltitudeEvent_action_stop() {
        setup_global_test_eop();

        let event = DAltitudeEvent::new(500e3, "Test", EventDirection::Any).set_terminal();
        assert_eq!(event.action(), EventAction::Stop);
    }

    #[test]
    fn test_DAltitudeEvent_direction_increasing() {
        setup_global_test_eop();

        let event = DAltitudeEvent::new(500e3, "Ascending", EventDirection::Increasing);
        assert_eq!(event.direction(), EventDirection::Increasing);
    }

    #[test]
    fn test_DAltitudeEvent_direction_decreasing() {
        setup_global_test_eop();

        let event = DAltitudeEvent::new(500e3, "Descending", EventDirection::Decreasing);
        assert_eq!(event.direction(), EventDirection::Decreasing);
    }

    #[test]
    fn test_DAltitudeEvent_direction_any() {
        setup_global_test_eop();

        let event = DAltitudeEvent::new(500e3, "Any", EventDirection::Any);
        assert_eq!(event.direction(), EventDirection::Any);
    }

    #[test]
    fn test_DAltitudeEvent_builder_chaining() {
        setup_global_test_eop();

        let callback: DEventCallback =
            Box::new(|_t, _state, _params| (None, None, EventAction::Stop));

        let event = DAltitudeEvent::new(200e3, "Reentry", EventDirection::Decreasing)
            .with_instance(1)
            .with_tolerances(1e-5, 1e-8)
            .with_step_reduction_factor(0.1)
            .with_callback(callback)
            .set_terminal();

        assert_eq!(event.name(), "Reentry 1");
        assert_eq!(event.time_tolerance(), 1e-5);
        assert_eq!(event.value_tolerance(), 1e-8);
        assert_eq!(event.step_reduction_factor(), 0.1);
        assert!(event.callback().is_some());
        assert_eq!(event.action(), EventAction::Stop);
    }

    // =========================================================================
    // Orbital Element Event Tests
    // =========================================================================

    #[test]
    fn test_SSemiMajorAxisEvent_new() {
        let event =
            SSemiMajorAxisEvent::<6, 0>::new(7000e3, "SMA Check", EventDirection::Increasing);
        assert_eq!(event.name(), "SMA Check");
        assert_eq!(event.target_value(), 7000e3);
        assert_eq!(event.direction(), EventDirection::Increasing);
    }

    #[test]
    fn test_SSemiMajorAxisEvent_evaluate() {
        let event = SSemiMajorAxisEvent::<6, 0>::new(7000e3, "SMA Check", EventDirection::Any);
        let epoch = Epoch::from_jd(2451545.0, TimeSystem::UTC);

        // Circular orbit: v = sqrt(GM/r) for r = 7000 km
        // GM_EARTH = 3.986004418e14 m^3/s^2
        // v = sqrt(3.986e14 / 7e6) = 7546.05 m/s
        let state = Vector6::new(7000e3, 0.0, 0.0, 0.0, 7546.05, 0.0);
        let val = event.evaluate(epoch, &state, None);
        // evaluate() returns raw SMA value
        // For a circular orbit at r=7000km, SMA should be very close to 7000 km
        assert!((val - 7000e3).abs() < 1000.0); // Within 1 km of expected SMA
    }

    #[test]
    fn test_DSemiMajorAxisEvent_new() {
        let event = DSemiMajorAxisEvent::new(7000e3, "SMA Check", EventDirection::Decreasing);
        assert_eq!(event.name(), "SMA Check");
        assert_eq!(event.target_value(), 7000e3);
        assert_eq!(event.direction(), EventDirection::Decreasing);
    }

    #[test]
    fn test_SSemiMajorAxisEvent_target_value() {
        let event = SSemiMajorAxisEvent::<6, 0>::new(8000e3, "SMA", EventDirection::Any);
        assert_eq!(event.target_value(), 8000e3);
    }

    #[test]
    fn test_SSemiMajorAxisEvent_name() {
        let event = SSemiMajorAxisEvent::<6, 0>::new(7000e3, "Test SMA Name", EventDirection::Any);
        assert_eq!(event.name(), "Test SMA Name");
    }

    #[test]
    fn test_SSemiMajorAxisEvent_with_instance() {
        let event = SSemiMajorAxisEvent::<6, 0>::new(7000e3, "SMA Check", EventDirection::Any)
            .with_instance(3);
        assert_eq!(event.name(), "SMA Check 3");
    }

    #[test]
    fn test_SSemiMajorAxisEvent_with_tolerances() {
        let event = SSemiMajorAxisEvent::<6, 0>::new(7000e3, "Test", EventDirection::Any);
        assert_eq!(event.time_tolerance(), 1e-6);
        assert_eq!(event.value_tolerance(), 1e-9);

        let event = SSemiMajorAxisEvent::<6, 0>::new(7000e3, "Test", EventDirection::Any)
            .with_tolerances(1e-4, 1e-6);
        assert_eq!(event.time_tolerance(), 1e-4);
        assert_eq!(event.value_tolerance(), 1e-6);
    }

    #[test]
    fn test_SSemiMajorAxisEvent_time_tolerance() {
        let event = SSemiMajorAxisEvent::<6, 0>::new(7000e3, "Test", EventDirection::Any)
            .with_tolerances(5e-5, 1e-9);
        assert_eq!(event.time_tolerance(), 5e-5);
    }

    #[test]
    fn test_SSemiMajorAxisEvent_value_tolerance() {
        let event = SSemiMajorAxisEvent::<6, 0>::new(7000e3, "Test", EventDirection::Any)
            .with_tolerances(1e-6, 5e-8);
        assert_eq!(event.value_tolerance(), 5e-8);
    }

    #[test]
    fn test_SSemiMajorAxisEvent_with_step_reduction_factor() {
        let event = SSemiMajorAxisEvent::<6, 0>::new(7000e3, "Test", EventDirection::Any);
        assert_eq!(event.step_reduction_factor(), 0.2);

        let event = SSemiMajorAxisEvent::<6, 0>::new(7000e3, "Test", EventDirection::Any)
            .with_step_reduction_factor(0.15);
        assert_eq!(event.step_reduction_factor(), 0.15);
    }

    #[test]
    fn test_SSemiMajorAxisEvent_step_reduction_factor() {
        let event = SSemiMajorAxisEvent::<6, 0>::new(7000e3, "Test", EventDirection::Any)
            .with_step_reduction_factor(0.1);
        assert_eq!(event.step_reduction_factor(), 0.1);
    }

    #[test]
    fn test_SSemiMajorAxisEvent_with_callback() {
        let called = Arc::new(AtomicBool::new(false));
        let called_clone = called.clone();

        let callback: SEventCallback<6, 0> = Box::new(move |_t, _state, _params| {
            called_clone.store(true, Ordering::SeqCst);
            (None, None, EventAction::Continue)
        });

        let event = SSemiMajorAxisEvent::<6, 0>::new(7000e3, "Test", EventDirection::Any)
            .with_callback(callback);
        assert!(event.callback().is_some());

        let epoch = Epoch::from_jd(2451545.0, TimeSystem::UTC);
        let state = Vector6::zeros();
        if let Some(cb) = event.callback() {
            cb(epoch, &state, None);
        }
        assert!(called.load(Ordering::SeqCst));
    }

    #[test]
    fn test_SSemiMajorAxisEvent_callback_none() {
        let event = SSemiMajorAxisEvent::<6, 0>::new(7000e3, "Test", EventDirection::Any);
        assert!(event.callback().is_none());
    }

    #[test]
    fn test_SSemiMajorAxisEvent_set_terminal() {
        let event = SSemiMajorAxisEvent::<6, 0>::new(7000e3, "Test", EventDirection::Any);
        assert_eq!(event.action(), EventAction::Continue);

        let event =
            SSemiMajorAxisEvent::<6, 0>::new(7000e3, "Test", EventDirection::Any).set_terminal();
        assert_eq!(event.action(), EventAction::Stop);
    }

    #[test]
    fn test_SSemiMajorAxisEvent_action_continue() {
        let event = SSemiMajorAxisEvent::<6, 0>::new(7000e3, "Test", EventDirection::Any);
        assert_eq!(event.action(), EventAction::Continue);
    }

    #[test]
    fn test_SSemiMajorAxisEvent_action_stop() {
        let event =
            SSemiMajorAxisEvent::<6, 0>::new(7000e3, "Test", EventDirection::Any).set_terminal();
        assert_eq!(event.action(), EventAction::Stop);
    }

    #[test]
    fn test_SSemiMajorAxisEvent_direction_increasing() {
        let event =
            SSemiMajorAxisEvent::<6, 0>::new(7000e3, "Ascending", EventDirection::Increasing);
        assert_eq!(event.direction(), EventDirection::Increasing);
    }

    #[test]
    fn test_SSemiMajorAxisEvent_direction_decreasing() {
        let event =
            SSemiMajorAxisEvent::<6, 0>::new(7000e3, "Descending", EventDirection::Decreasing);
        assert_eq!(event.direction(), EventDirection::Decreasing);
    }

    #[test]
    fn test_SSemiMajorAxisEvent_direction_any() {
        let event = SSemiMajorAxisEvent::<6, 0>::new(7000e3, "Any", EventDirection::Any);
        assert_eq!(event.direction(), EventDirection::Any);
    }

    #[test]
    fn test_SSemiMajorAxisEvent_builder_chaining() {
        let callback: SEventCallback<6, 0> =
            Box::new(|_t, _state, _params| (None, None, EventAction::Stop));

        let event = SSemiMajorAxisEvent::<6, 0>::new(7000e3, "SMA", EventDirection::Increasing)
            .with_instance(1)
            .with_tolerances(1e-5, 1e-8)
            .with_step_reduction_factor(0.1)
            .with_callback(callback)
            .set_terminal();

        assert_eq!(event.name(), "SMA 1");
        assert_eq!(event.time_tolerance(), 1e-5);
        assert_eq!(event.value_tolerance(), 1e-8);
        assert_eq!(event.step_reduction_factor(), 0.1);
        assert!(event.callback().is_some());
        assert_eq!(event.action(), EventAction::Stop);
    }

    // DSemiMajorAxisEvent additional tests

    #[test]
    fn test_DSemiMajorAxisEvent_evaluate() {
        let event = DSemiMajorAxisEvent::new(7000e3, "SMA Test", EventDirection::Any);
        let epoch = Epoch::from_jd(2451545.0, TimeSystem::UTC);

        let state = DVector::from_vec(vec![7000e3, 0.0, 0.0, 0.0, 7546.05, 0.0]);
        let val = event.evaluate(epoch, &state, None);
        assert!((val - 7000e3).abs() < 1000.0);
    }

    #[test]
    fn test_DSemiMajorAxisEvent_target_value() {
        let event = DSemiMajorAxisEvent::new(8000e3, "SMA", EventDirection::Any);
        assert_eq!(event.target_value(), 8000e3);
    }

    #[test]
    fn test_DSemiMajorAxisEvent_name() {
        let event = DSemiMajorAxisEvent::new(7000e3, "Test SMA Name", EventDirection::Any);
        assert_eq!(event.name(), "Test SMA Name");
    }

    #[test]
    fn test_DSemiMajorAxisEvent_with_instance() {
        let event =
            DSemiMajorAxisEvent::new(7000e3, "SMA Check", EventDirection::Any).with_instance(2);
        assert_eq!(event.name(), "SMA Check 2");
    }

    #[test]
    fn test_DSemiMajorAxisEvent_with_tolerances() {
        let event = DSemiMajorAxisEvent::new(7000e3, "Test", EventDirection::Any)
            .with_tolerances(1e-4, 1e-7);
        assert_eq!(event.time_tolerance(), 1e-4);
        assert_eq!(event.value_tolerance(), 1e-7);
    }

    #[test]
    fn test_DSemiMajorAxisEvent_time_tolerance() {
        let event = DSemiMajorAxisEvent::new(7000e3, "Test", EventDirection::Any)
            .with_tolerances(2e-5, 1e-9);
        assert_eq!(event.time_tolerance(), 2e-5);
    }

    #[test]
    fn test_DSemiMajorAxisEvent_value_tolerance() {
        let event = DSemiMajorAxisEvent::new(7000e3, "Test", EventDirection::Any)
            .with_tolerances(1e-6, 3e-8);
        assert_eq!(event.value_tolerance(), 3e-8);
    }

    #[test]
    fn test_DSemiMajorAxisEvent_with_step_reduction_factor() {
        let event = DSemiMajorAxisEvent::new(7000e3, "Test", EventDirection::Any);
        assert_eq!(event.step_reduction_factor(), 0.2);

        let event = DSemiMajorAxisEvent::new(7000e3, "Test", EventDirection::Any)
            .with_step_reduction_factor(0.25);
        assert_eq!(event.step_reduction_factor(), 0.25);
    }

    #[test]
    fn test_DSemiMajorAxisEvent_step_reduction_factor() {
        let event = DSemiMajorAxisEvent::new(7000e3, "Test", EventDirection::Any)
            .with_step_reduction_factor(0.12);
        assert_eq!(event.step_reduction_factor(), 0.12);
    }

    #[test]
    fn test_DSemiMajorAxisEvent_with_callback() {
        let called = Arc::new(AtomicBool::new(false));
        let called_clone = called.clone();

        let callback: DEventCallback = Box::new(move |_t, _state, _params| {
            called_clone.store(true, Ordering::SeqCst);
            (None, None, EventAction::Continue)
        });

        let event =
            DSemiMajorAxisEvent::new(7000e3, "Test", EventDirection::Any).with_callback(callback);
        assert!(event.callback().is_some());
    }

    #[test]
    fn test_DSemiMajorAxisEvent_callback_none() {
        let event = DSemiMajorAxisEvent::new(7000e3, "Test", EventDirection::Any);
        assert!(event.callback().is_none());
    }

    #[test]
    fn test_DSemiMajorAxisEvent_set_terminal() {
        let event = DSemiMajorAxisEvent::new(7000e3, "Test", EventDirection::Any);
        assert_eq!(event.action(), EventAction::Continue);

        let event = DSemiMajorAxisEvent::new(7000e3, "Test", EventDirection::Any).set_terminal();
        assert_eq!(event.action(), EventAction::Stop);
    }

    #[test]
    fn test_DSemiMajorAxisEvent_action_continue() {
        let event = DSemiMajorAxisEvent::new(7000e3, "Test", EventDirection::Any);
        assert_eq!(event.action(), EventAction::Continue);
    }

    #[test]
    fn test_DSemiMajorAxisEvent_action_stop() {
        let event = DSemiMajorAxisEvent::new(7000e3, "Test", EventDirection::Any).set_terminal();
        assert_eq!(event.action(), EventAction::Stop);
    }

    #[test]
    fn test_DSemiMajorAxisEvent_direction_increasing() {
        let event = DSemiMajorAxisEvent::new(7000e3, "Ascending", EventDirection::Increasing);
        assert_eq!(event.direction(), EventDirection::Increasing);
    }

    #[test]
    fn test_DSemiMajorAxisEvent_direction_decreasing() {
        let event = DSemiMajorAxisEvent::new(7000e3, "Descending", EventDirection::Decreasing);
        assert_eq!(event.direction(), EventDirection::Decreasing);
    }

    #[test]
    fn test_DSemiMajorAxisEvent_direction_any() {
        let event = DSemiMajorAxisEvent::new(7000e3, "Any", EventDirection::Any);
        assert_eq!(event.direction(), EventDirection::Any);
    }

    #[test]
    fn test_DSemiMajorAxisEvent_builder_chaining() {
        let callback: DEventCallback =
            Box::new(|_t, _state, _params| (None, None, EventAction::Stop));

        let event = DSemiMajorAxisEvent::new(7000e3, "SMA", EventDirection::Decreasing)
            .with_instance(1)
            .with_tolerances(1e-5, 1e-8)
            .with_step_reduction_factor(0.1)
            .with_callback(callback)
            .set_terminal();

        assert_eq!(event.name(), "SMA 1");
        assert_eq!(event.time_tolerance(), 1e-5);
        assert_eq!(event.value_tolerance(), 1e-8);
        assert_eq!(event.step_reduction_factor(), 0.1);
        assert!(event.callback().is_some());
        assert_eq!(event.action(), EventAction::Stop);
    }

    #[test]
    fn test_SEccentricityEvent_new() {
        let event = SEccentricityEvent::<6, 0>::new(0.1, "Ecc value", EventDirection::Increasing);
        assert_eq!(event.name(), "Ecc value");
        assert_eq!(event.target_value(), 0.1);
        assert_eq!(event.direction(), EventDirection::Increasing);
    }

    #[test]
    fn test_SEccentricityEvent_evaluate() {
        let event = SEccentricityEvent::<6, 0>::new(0.1, "Ecc Test", EventDirection::Any);
        let epoch = Epoch::from_jd(2451545.0, TimeSystem::UTC);

        // Circular orbit: v = sqrt(GM/r) at r = 7000 km gives e ≈ 0
        let state = Vector6::new(7000e3, 0.0, 0.0, 0.0, 7546.05, 0.0);
        let val = event.evaluate(epoch, &state, None);
        // evaluate() returns raw eccentricity value
        // For a circular orbit, eccentricity should be very small (close to 0)
        assert!(val >= 0.0); // Eccentricity is always non-negative
        assert!(val < 0.01); // Should be nearly circular
    }

    #[test]
    fn test_DEccentricityEvent_new() {
        let event = DEccentricityEvent::new(0.5, "Ecc Check", EventDirection::Any);
        assert_eq!(event.name(), "Ecc Check");
        assert_eq!(event.target_value(), 0.5);
    }

    #[test]
    fn test_SEccentricityEvent_target_value() {
        let event = SEccentricityEvent::<6, 0>::new(0.25, "Ecc", EventDirection::Any);
        assert_eq!(event.target_value(), 0.25);
    }

    #[test]
    fn test_SEccentricityEvent_name() {
        let event = SEccentricityEvent::<6, 0>::new(0.1, "Test Ecc Name", EventDirection::Any);
        assert_eq!(event.name(), "Test Ecc Name");
    }

    #[test]
    fn test_SEccentricityEvent_with_instance() {
        let event =
            SEccentricityEvent::<6, 0>::new(0.1, "Ecc Check", EventDirection::Any).with_instance(3);
        assert_eq!(event.name(), "Ecc Check 3");
    }

    #[test]
    fn test_SEccentricityEvent_with_tolerances() {
        let event = SEccentricityEvent::<6, 0>::new(0.1, "Test", EventDirection::Any);
        assert_eq!(event.time_tolerance(), 1e-6);
        assert_eq!(event.value_tolerance(), 1e-9);

        let event = SEccentricityEvent::<6, 0>::new(0.1, "Test", EventDirection::Any)
            .with_tolerances(1e-4, 1e-6);
        assert_eq!(event.time_tolerance(), 1e-4);
        assert_eq!(event.value_tolerance(), 1e-6);
    }

    #[test]
    fn test_SEccentricityEvent_time_tolerance() {
        let event = SEccentricityEvent::<6, 0>::new(0.1, "Test", EventDirection::Any)
            .with_tolerances(5e-5, 1e-9);
        assert_eq!(event.time_tolerance(), 5e-5);
    }

    #[test]
    fn test_SEccentricityEvent_value_tolerance() {
        let event = SEccentricityEvent::<6, 0>::new(0.1, "Test", EventDirection::Any)
            .with_tolerances(1e-6, 5e-8);
        assert_eq!(event.value_tolerance(), 5e-8);
    }

    #[test]
    fn test_SEccentricityEvent_with_step_reduction_factor() {
        let event = SEccentricityEvent::<6, 0>::new(0.1, "Test", EventDirection::Any);
        assert_eq!(event.step_reduction_factor(), 0.2);

        let event = SEccentricityEvent::<6, 0>::new(0.1, "Test", EventDirection::Any)
            .with_step_reduction_factor(0.15);
        assert_eq!(event.step_reduction_factor(), 0.15);
    }

    #[test]
    fn test_SEccentricityEvent_step_reduction_factor() {
        let event = SEccentricityEvent::<6, 0>::new(0.1, "Test", EventDirection::Any)
            .with_step_reduction_factor(0.1);
        assert_eq!(event.step_reduction_factor(), 0.1);
    }

    #[test]
    fn test_SEccentricityEvent_with_callback() {
        let called = Arc::new(AtomicBool::new(false));
        let called_clone = called.clone();

        let callback: SEventCallback<6, 0> = Box::new(move |_t, _state, _params| {
            called_clone.store(true, Ordering::SeqCst);
            (None, None, EventAction::Continue)
        });

        let event = SEccentricityEvent::<6, 0>::new(0.1, "Test", EventDirection::Any)
            .with_callback(callback);
        assert!(event.callback().is_some());

        let epoch = Epoch::from_jd(2451545.0, TimeSystem::UTC);
        let state = Vector6::zeros();
        if let Some(cb) = event.callback() {
            cb(epoch, &state, None);
        }
        assert!(called.load(Ordering::SeqCst));
    }

    #[test]
    fn test_SEccentricityEvent_callback_none() {
        let event = SEccentricityEvent::<6, 0>::new(0.1, "Test", EventDirection::Any);
        assert!(event.callback().is_none());
    }

    #[test]
    fn test_SEccentricityEvent_set_terminal() {
        let event = SEccentricityEvent::<6, 0>::new(0.1, "Test", EventDirection::Any);
        assert_eq!(event.action(), EventAction::Continue);

        let event =
            SEccentricityEvent::<6, 0>::new(0.1, "Test", EventDirection::Any).set_terminal();
        assert_eq!(event.action(), EventAction::Stop);
    }

    #[test]
    fn test_SEccentricityEvent_action_continue() {
        let event = SEccentricityEvent::<6, 0>::new(0.1, "Test", EventDirection::Any);
        assert_eq!(event.action(), EventAction::Continue);
    }

    #[test]
    fn test_SEccentricityEvent_action_stop() {
        let event =
            SEccentricityEvent::<6, 0>::new(0.1, "Test", EventDirection::Any).set_terminal();
        assert_eq!(event.action(), EventAction::Stop);
    }

    #[test]
    fn test_SEccentricityEvent_direction_increasing() {
        let event = SEccentricityEvent::<6, 0>::new(0.1, "Ascending", EventDirection::Increasing);
        assert_eq!(event.direction(), EventDirection::Increasing);
    }

    #[test]
    fn test_SEccentricityEvent_direction_decreasing() {
        let event = SEccentricityEvent::<6, 0>::new(0.1, "Descending", EventDirection::Decreasing);
        assert_eq!(event.direction(), EventDirection::Decreasing);
    }

    #[test]
    fn test_SEccentricityEvent_direction_any() {
        let event = SEccentricityEvent::<6, 0>::new(0.1, "Any", EventDirection::Any);
        assert_eq!(event.direction(), EventDirection::Any);
    }

    #[test]
    fn test_SEccentricityEvent_builder_chaining() {
        let callback: SEventCallback<6, 0> =
            Box::new(|_t, _state, _params| (None, None, EventAction::Stop));

        let event = SEccentricityEvent::<6, 0>::new(0.1, "Ecc", EventDirection::Increasing)
            .with_instance(1)
            .with_tolerances(1e-5, 1e-8)
            .with_step_reduction_factor(0.1)
            .with_callback(callback)
            .set_terminal();

        assert_eq!(event.name(), "Ecc 1");
        assert_eq!(event.time_tolerance(), 1e-5);
        assert_eq!(event.value_tolerance(), 1e-8);
        assert_eq!(event.step_reduction_factor(), 0.1);
        assert!(event.callback().is_some());
        assert_eq!(event.action(), EventAction::Stop);
    }

    // DEccentricityEvent additional tests

    #[test]
    fn test_DEccentricityEvent_evaluate() {
        let event = DEccentricityEvent::new(0.1, "Ecc Test", EventDirection::Any);
        let epoch = Epoch::from_jd(2451545.0, TimeSystem::UTC);

        let state = DVector::from_vec(vec![7000e3, 0.0, 0.0, 0.0, 7546.05, 0.0]);
        let val = event.evaluate(epoch, &state, None);
        assert!(val >= 0.0);
        assert!(val < 0.01);
    }

    #[test]
    fn test_DEccentricityEvent_target_value() {
        let event = DEccentricityEvent::new(0.25, "Ecc", EventDirection::Any);
        assert_eq!(event.target_value(), 0.25);
    }

    #[test]
    fn test_DEccentricityEvent_name() {
        let event = DEccentricityEvent::new(0.1, "Test Ecc Name", EventDirection::Any);
        assert_eq!(event.name(), "Test Ecc Name");
    }

    #[test]
    fn test_DEccentricityEvent_with_instance() {
        let event = DEccentricityEvent::new(0.1, "Ecc Check", EventDirection::Any).with_instance(2);
        assert_eq!(event.name(), "Ecc Check 2");
    }

    #[test]
    fn test_DEccentricityEvent_with_tolerances() {
        let event =
            DEccentricityEvent::new(0.1, "Test", EventDirection::Any).with_tolerances(1e-4, 1e-7);
        assert_eq!(event.time_tolerance(), 1e-4);
        assert_eq!(event.value_tolerance(), 1e-7);
    }

    #[test]
    fn test_DEccentricityEvent_time_tolerance() {
        let event =
            DEccentricityEvent::new(0.1, "Test", EventDirection::Any).with_tolerances(2e-5, 1e-9);
        assert_eq!(event.time_tolerance(), 2e-5);
    }

    #[test]
    fn test_DEccentricityEvent_value_tolerance() {
        let event =
            DEccentricityEvent::new(0.1, "Test", EventDirection::Any).with_tolerances(1e-6, 3e-8);
        assert_eq!(event.value_tolerance(), 3e-8);
    }

    #[test]
    fn test_DEccentricityEvent_with_step_reduction_factor() {
        let event = DEccentricityEvent::new(0.1, "Test", EventDirection::Any);
        assert_eq!(event.step_reduction_factor(), 0.2);

        let event = DEccentricityEvent::new(0.1, "Test", EventDirection::Any)
            .with_step_reduction_factor(0.25);
        assert_eq!(event.step_reduction_factor(), 0.25);
    }

    #[test]
    fn test_DEccentricityEvent_step_reduction_factor() {
        let event = DEccentricityEvent::new(0.1, "Test", EventDirection::Any)
            .with_step_reduction_factor(0.12);
        assert_eq!(event.step_reduction_factor(), 0.12);
    }

    #[test]
    fn test_DEccentricityEvent_with_callback() {
        let called = Arc::new(AtomicBool::new(false));
        let called_clone = called.clone();

        let callback: DEventCallback = Box::new(move |_t, _state, _params| {
            called_clone.store(true, Ordering::SeqCst);
            (None, None, EventAction::Continue)
        });

        let event =
            DEccentricityEvent::new(0.1, "Test", EventDirection::Any).with_callback(callback);
        assert!(event.callback().is_some());
    }

    #[test]
    fn test_DEccentricityEvent_callback_none() {
        let event = DEccentricityEvent::new(0.1, "Test", EventDirection::Any);
        assert!(event.callback().is_none());
    }

    #[test]
    fn test_DEccentricityEvent_set_terminal() {
        let event = DEccentricityEvent::new(0.1, "Test", EventDirection::Any);
        assert_eq!(event.action(), EventAction::Continue);

        let event = DEccentricityEvent::new(0.1, "Test", EventDirection::Any).set_terminal();
        assert_eq!(event.action(), EventAction::Stop);
    }

    #[test]
    fn test_DEccentricityEvent_action_continue() {
        let event = DEccentricityEvent::new(0.1, "Test", EventDirection::Any);
        assert_eq!(event.action(), EventAction::Continue);
    }

    #[test]
    fn test_DEccentricityEvent_action_stop() {
        let event = DEccentricityEvent::new(0.1, "Test", EventDirection::Any).set_terminal();
        assert_eq!(event.action(), EventAction::Stop);
    }

    #[test]
    fn test_DEccentricityEvent_direction_increasing() {
        let event = DEccentricityEvent::new(0.1, "Ascending", EventDirection::Increasing);
        assert_eq!(event.direction(), EventDirection::Increasing);
    }

    #[test]
    fn test_DEccentricityEvent_direction_decreasing() {
        let event = DEccentricityEvent::new(0.1, "Descending", EventDirection::Decreasing);
        assert_eq!(event.direction(), EventDirection::Decreasing);
    }

    #[test]
    fn test_DEccentricityEvent_direction_any() {
        let event = DEccentricityEvent::new(0.1, "Any", EventDirection::Any);
        assert_eq!(event.direction(), EventDirection::Any);
    }

    #[test]
    fn test_DEccentricityEvent_builder_chaining() {
        let callback: DEventCallback =
            Box::new(|_t, _state, _params| (None, None, EventAction::Stop));

        let event = DEccentricityEvent::new(0.1, "Ecc", EventDirection::Decreasing)
            .with_instance(1)
            .with_tolerances(1e-5, 1e-8)
            .with_step_reduction_factor(0.1)
            .with_callback(callback)
            .set_terminal();

        assert_eq!(event.name(), "Ecc 1");
        assert_eq!(event.time_tolerance(), 1e-5);
        assert_eq!(event.value_tolerance(), 1e-8);
        assert_eq!(event.step_reduction_factor(), 0.1);
        assert!(event.callback().is_some());
        assert_eq!(event.action(), EventAction::Stop);
    }

    #[test]
    fn test_SInclinationEvent_new() {
        let inc_rad = std::f64::consts::PI / 4.0; // 45 degrees
        let event = SInclinationEvent::<6, 0>::new(
            inc_rad,
            "Inc value",
            EventDirection::Increasing,
            AngleFormat::Radians,
        );
        assert_eq!(event.name(), "Inc value");
        assert_eq!(event.target_value(), inc_rad);
    }

    #[test]
    fn test_DInclinationEvent_new() {
        let event =
            DInclinationEvent::new(45.0, "Inc Check", EventDirection::Any, AngleFormat::Degrees);
        assert_eq!(event.name(), "Inc Check");
    }

    #[test]
    fn test_SInclinationEvent_new_degrees() {
        let event = SInclinationEvent::<6, 0>::new(
            45.0,
            "Inc Deg",
            EventDirection::Increasing,
            AngleFormat::Degrees,
        );
        assert_eq!(event.name(), "Inc Deg");
        assert!((event.target_value() - std::f64::consts::PI / 4.0).abs() < 1e-10);
    }

    #[test]
    fn test_SInclinationEvent_target_value() {
        let event =
            SInclinationEvent::<6, 0>::new(0.5, "Inc", EventDirection::Any, AngleFormat::Radians);
        assert_eq!(event.target_value(), 0.5);
    }

    #[test]
    fn test_SInclinationEvent_name() {
        let event = SInclinationEvent::<6, 0>::new(
            45.0,
            "Test Inc Name",
            EventDirection::Any,
            AngleFormat::Degrees,
        );
        assert_eq!(event.name(), "Test Inc Name");
    }

    #[test]
    fn test_SInclinationEvent_with_instance() {
        let event = SInclinationEvent::<6, 0>::new(
            45.0,
            "Inc Check",
            EventDirection::Any,
            AngleFormat::Degrees,
        )
        .with_instance(3);
        assert_eq!(event.name(), "Inc Check 3");
    }

    #[test]
    fn test_SInclinationEvent_with_tolerances() {
        let event =
            SInclinationEvent::<6, 0>::new(45.0, "Test", EventDirection::Any, AngleFormat::Degrees);
        assert_eq!(event.time_tolerance(), 1e-6);
        assert_eq!(event.value_tolerance(), 1e-9);

        let event =
            SInclinationEvent::<6, 0>::new(45.0, "Test", EventDirection::Any, AngleFormat::Degrees)
                .with_tolerances(1e-4, 1e-6);
        assert_eq!(event.time_tolerance(), 1e-4);
        assert_eq!(event.value_tolerance(), 1e-6);
    }

    #[test]
    fn test_SInclinationEvent_time_tolerance() {
        let event =
            SInclinationEvent::<6, 0>::new(45.0, "Test", EventDirection::Any, AngleFormat::Degrees)
                .with_tolerances(5e-5, 1e-9);
        assert_eq!(event.time_tolerance(), 5e-5);
    }

    #[test]
    fn test_SInclinationEvent_value_tolerance() {
        let event =
            SInclinationEvent::<6, 0>::new(45.0, "Test", EventDirection::Any, AngleFormat::Degrees)
                .with_tolerances(1e-6, 5e-8);
        assert_eq!(event.value_tolerance(), 5e-8);
    }

    #[test]
    fn test_SInclinationEvent_with_step_reduction_factor() {
        let event =
            SInclinationEvent::<6, 0>::new(45.0, "Test", EventDirection::Any, AngleFormat::Degrees);
        assert_eq!(event.step_reduction_factor(), 0.2);

        let event =
            SInclinationEvent::<6, 0>::new(45.0, "Test", EventDirection::Any, AngleFormat::Degrees)
                .with_step_reduction_factor(0.15);
        assert_eq!(event.step_reduction_factor(), 0.15);
    }

    #[test]
    fn test_SInclinationEvent_step_reduction_factor() {
        let event =
            SInclinationEvent::<6, 0>::new(45.0, "Test", EventDirection::Any, AngleFormat::Degrees)
                .with_step_reduction_factor(0.1);
        assert_eq!(event.step_reduction_factor(), 0.1);
    }

    #[test]
    fn test_SInclinationEvent_with_callback() {
        let called = Arc::new(AtomicBool::new(false));
        let called_clone = called.clone();

        let callback: SEventCallback<6, 0> = Box::new(move |_t, _state, _params| {
            called_clone.store(true, Ordering::SeqCst);
            (None, None, EventAction::Continue)
        });

        let event =
            SInclinationEvent::<6, 0>::new(45.0, "Test", EventDirection::Any, AngleFormat::Degrees)
                .with_callback(callback);
        assert!(event.callback().is_some());

        let epoch = Epoch::from_jd(2451545.0, TimeSystem::UTC);
        let state = Vector6::zeros();
        if let Some(cb) = event.callback() {
            cb(epoch, &state, None);
        }
        assert!(called.load(Ordering::SeqCst));
    }

    #[test]
    fn test_SInclinationEvent_callback_none() {
        let event =
            SInclinationEvent::<6, 0>::new(45.0, "Test", EventDirection::Any, AngleFormat::Degrees);
        assert!(event.callback().is_none());
    }

    #[test]
    fn test_SInclinationEvent_set_terminal() {
        let event =
            SInclinationEvent::<6, 0>::new(45.0, "Test", EventDirection::Any, AngleFormat::Degrees);
        assert_eq!(event.action(), EventAction::Continue);

        let event =
            SInclinationEvent::<6, 0>::new(45.0, "Test", EventDirection::Any, AngleFormat::Degrees)
                .set_terminal();
        assert_eq!(event.action(), EventAction::Stop);
    }

    #[test]
    fn test_SInclinationEvent_action_continue() {
        let event =
            SInclinationEvent::<6, 0>::new(45.0, "Test", EventDirection::Any, AngleFormat::Degrees);
        assert_eq!(event.action(), EventAction::Continue);
    }

    #[test]
    fn test_SInclinationEvent_action_stop() {
        let event =
            SInclinationEvent::<6, 0>::new(45.0, "Test", EventDirection::Any, AngleFormat::Degrees)
                .set_terminal();
        assert_eq!(event.action(), EventAction::Stop);
    }

    #[test]
    fn test_SInclinationEvent_direction_increasing() {
        let event = SInclinationEvent::<6, 0>::new(
            45.0,
            "Ascending",
            EventDirection::Increasing,
            AngleFormat::Degrees,
        );
        assert_eq!(event.direction(), EventDirection::Increasing);
    }

    #[test]
    fn test_SInclinationEvent_direction_decreasing() {
        let event = SInclinationEvent::<6, 0>::new(
            45.0,
            "Descending",
            EventDirection::Decreasing,
            AngleFormat::Degrees,
        );
        assert_eq!(event.direction(), EventDirection::Decreasing);
    }

    #[test]
    fn test_SInclinationEvent_direction_any() {
        let event =
            SInclinationEvent::<6, 0>::new(45.0, "Any", EventDirection::Any, AngleFormat::Degrees);
        assert_eq!(event.direction(), EventDirection::Any);
    }

    #[test]
    fn test_SInclinationEvent_builder_chaining() {
        let callback: SEventCallback<6, 0> =
            Box::new(|_t, _state, _params| (None, None, EventAction::Stop));

        let event = SInclinationEvent::<6, 0>::new(
            45.0,
            "Inc",
            EventDirection::Increasing,
            AngleFormat::Degrees,
        )
        .with_instance(1)
        .with_tolerances(1e-5, 1e-8)
        .with_step_reduction_factor(0.1)
        .with_callback(callback)
        .set_terminal();

        assert_eq!(event.name(), "Inc 1");
        assert_eq!(event.time_tolerance(), 1e-5);
        assert_eq!(event.value_tolerance(), 1e-8);
        assert_eq!(event.step_reduction_factor(), 0.1);
        assert!(event.callback().is_some());
        assert_eq!(event.action(), EventAction::Stop);
    }

    // DInclinationEvent additional tests

    #[test]
    fn test_DInclinationEvent_new_radians() {
        let event = DInclinationEvent::new(
            std::f64::consts::PI / 4.0,
            "Inc Rad",
            EventDirection::Any,
            AngleFormat::Radians,
        );
        assert_eq!(event.name(), "Inc Rad");
    }

    #[test]
    fn test_DInclinationEvent_target_value() {
        let event = DInclinationEvent::new(0.5, "Inc", EventDirection::Any, AngleFormat::Radians);
        assert_eq!(event.target_value(), 0.5);
    }

    #[test]
    fn test_DInclinationEvent_name() {
        let event = DInclinationEvent::new(
            45.0,
            "Test Inc Name",
            EventDirection::Any,
            AngleFormat::Degrees,
        );
        assert_eq!(event.name(), "Test Inc Name");
    }

    #[test]
    fn test_DInclinationEvent_with_instance() {
        let event =
            DInclinationEvent::new(45.0, "Inc Check", EventDirection::Any, AngleFormat::Degrees)
                .with_instance(2);
        assert_eq!(event.name(), "Inc Check 2");
    }

    #[test]
    fn test_DInclinationEvent_with_tolerances() {
        let event = DInclinationEvent::new(45.0, "Test", EventDirection::Any, AngleFormat::Degrees)
            .with_tolerances(1e-4, 1e-7);
        assert_eq!(event.time_tolerance(), 1e-4);
        assert_eq!(event.value_tolerance(), 1e-7);
    }

    #[test]
    fn test_DInclinationEvent_time_tolerance() {
        let event = DInclinationEvent::new(45.0, "Test", EventDirection::Any, AngleFormat::Degrees)
            .with_tolerances(2e-5, 1e-9);
        assert_eq!(event.time_tolerance(), 2e-5);
    }

    #[test]
    fn test_DInclinationEvent_value_tolerance() {
        let event = DInclinationEvent::new(45.0, "Test", EventDirection::Any, AngleFormat::Degrees)
            .with_tolerances(1e-6, 3e-8);
        assert_eq!(event.value_tolerance(), 3e-8);
    }

    #[test]
    fn test_DInclinationEvent_with_step_reduction_factor() {
        let event = DInclinationEvent::new(45.0, "Test", EventDirection::Any, AngleFormat::Degrees);
        assert_eq!(event.step_reduction_factor(), 0.2);

        let event = DInclinationEvent::new(45.0, "Test", EventDirection::Any, AngleFormat::Degrees)
            .with_step_reduction_factor(0.25);
        assert_eq!(event.step_reduction_factor(), 0.25);
    }

    #[test]
    fn test_DInclinationEvent_step_reduction_factor() {
        let event = DInclinationEvent::new(45.0, "Test", EventDirection::Any, AngleFormat::Degrees)
            .with_step_reduction_factor(0.12);
        assert_eq!(event.step_reduction_factor(), 0.12);
    }

    #[test]
    fn test_DInclinationEvent_with_callback() {
        let called = Arc::new(AtomicBool::new(false));
        let called_clone = called.clone();

        let callback: DEventCallback = Box::new(move |_t, _state, _params| {
            called_clone.store(true, Ordering::SeqCst);
            (None, None, EventAction::Continue)
        });

        let event = DInclinationEvent::new(45.0, "Test", EventDirection::Any, AngleFormat::Degrees)
            .with_callback(callback);
        assert!(event.callback().is_some());
    }

    #[test]
    fn test_DInclinationEvent_callback_none() {
        let event = DInclinationEvent::new(45.0, "Test", EventDirection::Any, AngleFormat::Degrees);
        assert!(event.callback().is_none());
    }

    #[test]
    fn test_DInclinationEvent_set_terminal() {
        let event = DInclinationEvent::new(45.0, "Test", EventDirection::Any, AngleFormat::Degrees);
        assert_eq!(event.action(), EventAction::Continue);

        let event = DInclinationEvent::new(45.0, "Test", EventDirection::Any, AngleFormat::Degrees)
            .set_terminal();
        assert_eq!(event.action(), EventAction::Stop);
    }

    #[test]
    fn test_DInclinationEvent_action_continue() {
        let event = DInclinationEvent::new(45.0, "Test", EventDirection::Any, AngleFormat::Degrees);
        assert_eq!(event.action(), EventAction::Continue);
    }

    #[test]
    fn test_DInclinationEvent_action_stop() {
        let event = DInclinationEvent::new(45.0, "Test", EventDirection::Any, AngleFormat::Degrees)
            .set_terminal();
        assert_eq!(event.action(), EventAction::Stop);
    }

    #[test]
    fn test_DInclinationEvent_direction_increasing() {
        let event = DInclinationEvent::new(
            45.0,
            "Ascending",
            EventDirection::Increasing,
            AngleFormat::Degrees,
        );
        assert_eq!(event.direction(), EventDirection::Increasing);
    }

    #[test]
    fn test_DInclinationEvent_direction_decreasing() {
        let event = DInclinationEvent::new(
            45.0,
            "Descending",
            EventDirection::Decreasing,
            AngleFormat::Degrees,
        );
        assert_eq!(event.direction(), EventDirection::Decreasing);
    }

    #[test]
    fn test_DInclinationEvent_direction_any() {
        let event = DInclinationEvent::new(45.0, "Any", EventDirection::Any, AngleFormat::Degrees);
        assert_eq!(event.direction(), EventDirection::Any);
    }

    #[test]
    fn test_DInclinationEvent_builder_chaining() {
        let callback: DEventCallback =
            Box::new(|_t, _state, _params| (None, None, EventAction::Stop));

        let event = DInclinationEvent::new(
            45.0,
            "Inc",
            EventDirection::Decreasing,
            AngleFormat::Degrees,
        )
        .with_instance(1)
        .with_tolerances(1e-5, 1e-8)
        .with_step_reduction_factor(0.1)
        .with_callback(callback)
        .set_terminal();

        assert_eq!(event.name(), "Inc 1");
        assert_eq!(event.time_tolerance(), 1e-5);
        assert_eq!(event.value_tolerance(), 1e-8);
        assert_eq!(event.step_reduction_factor(), 0.1);
        assert!(event.callback().is_some());
        assert_eq!(event.action(), EventAction::Stop);
    }

    #[test]
    fn test_SArgumentOfPerigeeEvent_new() {
        let event = SArgumentOfPerigeeEvent::<6, 0>::new(
            std::f64::consts::PI / 2.0,
            "AoP",
            EventDirection::Increasing,
            AngleFormat::Radians,
        );
        assert_eq!(event.name(), "AoP");
    }

    #[test]
    fn test_DArgumentOfPerigeeEvent_new() {
        let event =
            DArgumentOfPerigeeEvent::new(90.0, "AoP", EventDirection::Any, AngleFormat::Degrees);
        assert_eq!(event.name(), "AoP");
    }

    #[test]
    fn test_SArgumentOfPerigeeEvent_target_value() {
        let event = SArgumentOfPerigeeEvent::<6, 0>::new(
            0.5,
            "AoP",
            EventDirection::Any,
            AngleFormat::Radians,
        );
        assert_eq!(event.target_value(), 0.5);
    }

    #[test]
    fn test_SArgumentOfPerigeeEvent_with_instance() {
        let event = SArgumentOfPerigeeEvent::<6, 0>::new(
            90.0,
            "AoP Check",
            EventDirection::Any,
            AngleFormat::Degrees,
        )
        .with_instance(3);
        assert_eq!(event.name(), "AoP Check 3");
    }

    #[test]
    fn test_SArgumentOfPerigeeEvent_with_tolerances() {
        let event = SArgumentOfPerigeeEvent::<6, 0>::new(
            90.0,
            "Test",
            EventDirection::Any,
            AngleFormat::Degrees,
        )
        .with_tolerances(1e-4, 1e-6);
        assert_eq!(event.time_tolerance(), 1e-4);
        assert_eq!(event.value_tolerance(), 1e-6);
    }

    #[test]
    fn test_SArgumentOfPerigeeEvent_with_step_reduction_factor() {
        let event = SArgumentOfPerigeeEvent::<6, 0>::new(
            90.0,
            "Test",
            EventDirection::Any,
            AngleFormat::Degrees,
        )
        .with_step_reduction_factor(0.15);
        assert_eq!(event.step_reduction_factor(), 0.15);
    }

    #[test]
    fn test_SArgumentOfPerigeeEvent_with_callback() {
        let callback: SEventCallback<6, 0> =
            Box::new(move |_t, _state, _params| (None, None, EventAction::Continue));
        let event = SArgumentOfPerigeeEvent::<6, 0>::new(
            90.0,
            "Test",
            EventDirection::Any,
            AngleFormat::Degrees,
        )
        .with_callback(callback);
        assert!(event.callback().is_some());
    }

    #[test]
    fn test_SArgumentOfPerigeeEvent_callback_none() {
        let event = SArgumentOfPerigeeEvent::<6, 0>::new(
            90.0,
            "Test",
            EventDirection::Any,
            AngleFormat::Degrees,
        );
        assert!(event.callback().is_none());
    }

    #[test]
    fn test_SArgumentOfPerigeeEvent_set_terminal() {
        let event = SArgumentOfPerigeeEvent::<6, 0>::new(
            90.0,
            "Test",
            EventDirection::Any,
            AngleFormat::Degrees,
        )
        .set_terminal();
        assert_eq!(event.action(), EventAction::Stop);
    }

    #[test]
    fn test_SArgumentOfPerigeeEvent_direction_increasing() {
        let event = SArgumentOfPerigeeEvent::<6, 0>::new(
            90.0,
            "AoP",
            EventDirection::Increasing,
            AngleFormat::Degrees,
        );
        assert_eq!(event.direction(), EventDirection::Increasing);
    }

    #[test]
    fn test_SArgumentOfPerigeeEvent_direction_decreasing() {
        let event = SArgumentOfPerigeeEvent::<6, 0>::new(
            90.0,
            "AoP",
            EventDirection::Decreasing,
            AngleFormat::Degrees,
        );
        assert_eq!(event.direction(), EventDirection::Decreasing);
    }

    #[test]
    fn test_SArgumentOfPerigeeEvent_direction_any() {
        let event = SArgumentOfPerigeeEvent::<6, 0>::new(
            90.0,
            "AoP",
            EventDirection::Any,
            AngleFormat::Degrees,
        );
        assert_eq!(event.direction(), EventDirection::Any);
    }

    #[test]
    fn test_SArgumentOfPerigeeEvent_builder_chaining() {
        let callback: SEventCallback<6, 0> =
            Box::new(|_t, _state, _params| (None, None, EventAction::Stop));
        let event = SArgumentOfPerigeeEvent::<6, 0>::new(
            90.0,
            "AoP",
            EventDirection::Increasing,
            AngleFormat::Degrees,
        )
        .with_instance(1)
        .with_tolerances(1e-5, 1e-8)
        .with_step_reduction_factor(0.1)
        .with_callback(callback)
        .set_terminal();

        assert_eq!(event.name(), "AoP 1");
        assert_eq!(event.time_tolerance(), 1e-5);
        assert_eq!(event.step_reduction_factor(), 0.1);
        assert!(event.callback().is_some());
        assert_eq!(event.action(), EventAction::Stop);
    }

    // DArgumentOfPerigeeEvent additional tests

    #[test]
    fn test_DArgumentOfPerigeeEvent_target_value() {
        let event =
            DArgumentOfPerigeeEvent::new(0.5, "AoP", EventDirection::Any, AngleFormat::Radians);
        assert_eq!(event.target_value(), 0.5);
    }

    #[test]
    fn test_DArgumentOfPerigeeEvent_with_instance() {
        let event = DArgumentOfPerigeeEvent::new(
            90.0,
            "AoP Check",
            EventDirection::Any,
            AngleFormat::Degrees,
        )
        .with_instance(2);
        assert_eq!(event.name(), "AoP Check 2");
    }

    #[test]
    fn test_DArgumentOfPerigeeEvent_with_tolerances() {
        let event =
            DArgumentOfPerigeeEvent::new(90.0, "Test", EventDirection::Any, AngleFormat::Degrees)
                .with_tolerances(1e-4, 1e-7);
        assert_eq!(event.time_tolerance(), 1e-4);
        assert_eq!(event.value_tolerance(), 1e-7);
    }

    #[test]
    fn test_DArgumentOfPerigeeEvent_with_step_reduction_factor() {
        let event =
            DArgumentOfPerigeeEvent::new(90.0, "Test", EventDirection::Any, AngleFormat::Degrees)
                .with_step_reduction_factor(0.25);
        assert_eq!(event.step_reduction_factor(), 0.25);
    }

    #[test]
    fn test_DArgumentOfPerigeeEvent_with_callback() {
        let callback: DEventCallback =
            Box::new(move |_t, _state, _params| (None, None, EventAction::Continue));
        let event =
            DArgumentOfPerigeeEvent::new(90.0, "Test", EventDirection::Any, AngleFormat::Degrees)
                .with_callback(callback);
        assert!(event.callback().is_some());
    }

    #[test]
    fn test_DArgumentOfPerigeeEvent_callback_none() {
        let event =
            DArgumentOfPerigeeEvent::new(90.0, "Test", EventDirection::Any, AngleFormat::Degrees);
        assert!(event.callback().is_none());
    }

    #[test]
    fn test_DArgumentOfPerigeeEvent_set_terminal() {
        let event =
            DArgumentOfPerigeeEvent::new(90.0, "Test", EventDirection::Any, AngleFormat::Degrees)
                .set_terminal();
        assert_eq!(event.action(), EventAction::Stop);
    }

    #[test]
    fn test_DArgumentOfPerigeeEvent_direction_increasing() {
        let event = DArgumentOfPerigeeEvent::new(
            90.0,
            "AoP",
            EventDirection::Increasing,
            AngleFormat::Degrees,
        );
        assert_eq!(event.direction(), EventDirection::Increasing);
    }

    #[test]
    fn test_DArgumentOfPerigeeEvent_direction_decreasing() {
        let event = DArgumentOfPerigeeEvent::new(
            90.0,
            "AoP",
            EventDirection::Decreasing,
            AngleFormat::Degrees,
        );
        assert_eq!(event.direction(), EventDirection::Decreasing);
    }

    #[test]
    fn test_DArgumentOfPerigeeEvent_direction_any() {
        let event =
            DArgumentOfPerigeeEvent::new(90.0, "AoP", EventDirection::Any, AngleFormat::Degrees);
        assert_eq!(event.direction(), EventDirection::Any);
    }

    #[test]
    fn test_DArgumentOfPerigeeEvent_builder_chaining() {
        let callback: DEventCallback =
            Box::new(|_t, _state, _params| (None, None, EventAction::Stop));
        let event = DArgumentOfPerigeeEvent::new(
            90.0,
            "AoP",
            EventDirection::Decreasing,
            AngleFormat::Degrees,
        )
        .with_instance(1)
        .with_tolerances(1e-5, 1e-8)
        .with_step_reduction_factor(0.1)
        .with_callback(callback)
        .set_terminal();

        assert_eq!(event.name(), "AoP 1");
        assert_eq!(event.time_tolerance(), 1e-5);
        assert_eq!(event.step_reduction_factor(), 0.1);
        assert!(event.callback().is_some());
        assert_eq!(event.action(), EventAction::Stop);
    }

    #[test]
    fn test_SMeanAnomalyEvent_new() {
        let event = SMeanAnomalyEvent::<6, 0>::new(
            0.0,
            "Periapsis",
            EventDirection::Increasing,
            AngleFormat::Radians,
        );
        assert_eq!(event.name(), "Periapsis");
        assert_eq!(event.target_value(), 0.0);
    }

    #[test]
    fn test_DMeanAnomalyEvent_new() {
        let event =
            DMeanAnomalyEvent::new(180.0, "Apoapsis", EventDirection::Any, AngleFormat::Degrees);
        assert_eq!(event.name(), "Apoapsis");
    }

    #[test]
    fn test_SMeanAnomalyEvent_with_instance() {
        let event =
            SMeanAnomalyEvent::<6, 0>::new(0.0, "MA", EventDirection::Any, AngleFormat::Radians)
                .with_instance(3);
        assert_eq!(event.name(), "MA 3");
    }

    #[test]
    fn test_SMeanAnomalyEvent_with_tolerances() {
        let event =
            SMeanAnomalyEvent::<6, 0>::new(0.0, "Test", EventDirection::Any, AngleFormat::Radians)
                .with_tolerances(1e-4, 1e-6);
        assert_eq!(event.time_tolerance(), 1e-4);
        assert_eq!(event.value_tolerance(), 1e-6);
    }

    #[test]
    fn test_SMeanAnomalyEvent_with_step_reduction_factor() {
        let event =
            SMeanAnomalyEvent::<6, 0>::new(0.0, "Test", EventDirection::Any, AngleFormat::Radians)
                .with_step_reduction_factor(0.15);
        assert_eq!(event.step_reduction_factor(), 0.15);
    }

    #[test]
    fn test_SMeanAnomalyEvent_with_callback() {
        let callback: SEventCallback<6, 0> =
            Box::new(move |_t, _state, _params| (None, None, EventAction::Continue));
        let event =
            SMeanAnomalyEvent::<6, 0>::new(0.0, "Test", EventDirection::Any, AngleFormat::Radians)
                .with_callback(callback);
        assert!(event.callback().is_some());
    }

    #[test]
    fn test_SMeanAnomalyEvent_callback_none() {
        let event =
            SMeanAnomalyEvent::<6, 0>::new(0.0, "Test", EventDirection::Any, AngleFormat::Radians);
        assert!(event.callback().is_none());
    }

    #[test]
    fn test_SMeanAnomalyEvent_set_terminal() {
        let event =
            SMeanAnomalyEvent::<6, 0>::new(0.0, "Test", EventDirection::Any, AngleFormat::Radians)
                .set_terminal();
        assert_eq!(event.action(), EventAction::Stop);
    }

    #[test]
    fn test_SMeanAnomalyEvent_direction_increasing() {
        let event = SMeanAnomalyEvent::<6, 0>::new(
            0.0,
            "MA",
            EventDirection::Increasing,
            AngleFormat::Radians,
        );
        assert_eq!(event.direction(), EventDirection::Increasing);
    }

    #[test]
    fn test_SMeanAnomalyEvent_direction_decreasing() {
        let event = SMeanAnomalyEvent::<6, 0>::new(
            0.0,
            "MA",
            EventDirection::Decreasing,
            AngleFormat::Radians,
        );
        assert_eq!(event.direction(), EventDirection::Decreasing);
    }

    #[test]
    fn test_SMeanAnomalyEvent_direction_any() {
        let event =
            SMeanAnomalyEvent::<6, 0>::new(0.0, "MA", EventDirection::Any, AngleFormat::Radians);
        assert_eq!(event.direction(), EventDirection::Any);
    }

    #[test]
    fn test_SMeanAnomalyEvent_builder_chaining() {
        let callback: SEventCallback<6, 0> =
            Box::new(|_t, _state, _params| (None, None, EventAction::Stop));
        let event = SMeanAnomalyEvent::<6, 0>::new(
            0.0,
            "MA",
            EventDirection::Increasing,
            AngleFormat::Radians,
        )
        .with_instance(1)
        .with_tolerances(1e-5, 1e-8)
        .with_step_reduction_factor(0.1)
        .with_callback(callback)
        .set_terminal();

        assert_eq!(event.name(), "MA 1");
        assert_eq!(event.time_tolerance(), 1e-5);
        assert_eq!(event.step_reduction_factor(), 0.1);
        assert!(event.callback().is_some());
        assert_eq!(event.action(), EventAction::Stop);
    }

    // DMeanAnomalyEvent additional tests

    #[test]
    fn test_DMeanAnomalyEvent_with_instance() {
        let event = DMeanAnomalyEvent::new(180.0, "MA", EventDirection::Any, AngleFormat::Degrees)
            .with_instance(2);
        assert_eq!(event.name(), "MA 2");
    }

    #[test]
    fn test_DMeanAnomalyEvent_with_tolerances() {
        let event =
            DMeanAnomalyEvent::new(180.0, "Test", EventDirection::Any, AngleFormat::Degrees)
                .with_tolerances(1e-4, 1e-7);
        assert_eq!(event.time_tolerance(), 1e-4);
        assert_eq!(event.value_tolerance(), 1e-7);
    }

    #[test]
    fn test_DMeanAnomalyEvent_with_step_reduction_factor() {
        let event =
            DMeanAnomalyEvent::new(180.0, "Test", EventDirection::Any, AngleFormat::Degrees)
                .with_step_reduction_factor(0.25);
        assert_eq!(event.step_reduction_factor(), 0.25);
    }

    #[test]
    fn test_DMeanAnomalyEvent_with_callback() {
        let callback: DEventCallback =
            Box::new(move |_t, _state, _params| (None, None, EventAction::Continue));
        let event =
            DMeanAnomalyEvent::new(180.0, "Test", EventDirection::Any, AngleFormat::Degrees)
                .with_callback(callback);
        assert!(event.callback().is_some());
    }

    #[test]
    fn test_DMeanAnomalyEvent_callback_none() {
        let event =
            DMeanAnomalyEvent::new(180.0, "Test", EventDirection::Any, AngleFormat::Degrees);
        assert!(event.callback().is_none());
    }

    #[test]
    fn test_DMeanAnomalyEvent_set_terminal() {
        let event =
            DMeanAnomalyEvent::new(180.0, "Test", EventDirection::Any, AngleFormat::Degrees)
                .set_terminal();
        assert_eq!(event.action(), EventAction::Stop);
    }

    #[test]
    fn test_DMeanAnomalyEvent_direction_increasing() {
        let event = DMeanAnomalyEvent::new(
            180.0,
            "MA",
            EventDirection::Increasing,
            AngleFormat::Degrees,
        );
        assert_eq!(event.direction(), EventDirection::Increasing);
    }

    #[test]
    fn test_DMeanAnomalyEvent_direction_decreasing() {
        let event = DMeanAnomalyEvent::new(
            180.0,
            "MA",
            EventDirection::Decreasing,
            AngleFormat::Degrees,
        );
        assert_eq!(event.direction(), EventDirection::Decreasing);
    }

    #[test]
    fn test_DMeanAnomalyEvent_direction_any() {
        let event = DMeanAnomalyEvent::new(180.0, "MA", EventDirection::Any, AngleFormat::Degrees);
        assert_eq!(event.direction(), EventDirection::Any);
    }

    #[test]
    fn test_DMeanAnomalyEvent_builder_chaining() {
        let callback: DEventCallback =
            Box::new(|_t, _state, _params| (None, None, EventAction::Stop));
        let event = DMeanAnomalyEvent::new(
            180.0,
            "MA",
            EventDirection::Decreasing,
            AngleFormat::Degrees,
        )
        .with_instance(1)
        .with_tolerances(1e-5, 1e-8)
        .with_step_reduction_factor(0.1)
        .with_callback(callback)
        .set_terminal();

        assert_eq!(event.name(), "MA 1");
        assert_eq!(event.time_tolerance(), 1e-5);
        assert_eq!(event.step_reduction_factor(), 0.1);
        assert!(event.callback().is_some());
        assert_eq!(event.action(), EventAction::Stop);
    }

    #[test]
    fn test_SEccentricAnomalyEvent_new() {
        let event = SEccentricAnomalyEvent::<6, 0>::new(
            0.0,
            "Periapsis EA",
            EventDirection::Increasing,
            AngleFormat::Radians,
        );
        assert_eq!(event.name(), "Periapsis EA");
    }

    #[test]
    fn test_DEccentricAnomalyEvent_new() {
        let event =
            DEccentricAnomalyEvent::new(0.0, "EA Check", EventDirection::Any, AngleFormat::Degrees);
        assert_eq!(event.name(), "EA Check");
    }

    #[test]
    fn test_SEccentricAnomalyEvent_with_instance() {
        let event = SEccentricAnomalyEvent::<6, 0>::new(
            0.0,
            "EA",
            EventDirection::Any,
            AngleFormat::Radians,
        )
        .with_instance(3);
        assert_eq!(event.name(), "EA 3");
    }

    #[test]
    fn test_SEccentricAnomalyEvent_with_tolerances() {
        let event = SEccentricAnomalyEvent::<6, 0>::new(
            0.0,
            "Test",
            EventDirection::Any,
            AngleFormat::Radians,
        )
        .with_tolerances(1e-4, 1e-6);
        assert_eq!(event.time_tolerance(), 1e-4);
        assert_eq!(event.value_tolerance(), 1e-6);
    }

    #[test]
    fn test_SEccentricAnomalyEvent_with_step_reduction_factor() {
        let event = SEccentricAnomalyEvent::<6, 0>::new(
            0.0,
            "Test",
            EventDirection::Any,
            AngleFormat::Radians,
        )
        .with_step_reduction_factor(0.15);
        assert_eq!(event.step_reduction_factor(), 0.15);
    }

    #[test]
    fn test_SEccentricAnomalyEvent_with_callback() {
        let callback: SEventCallback<6, 0> =
            Box::new(move |_t, _state, _params| (None, None, EventAction::Continue));
        let event = SEccentricAnomalyEvent::<6, 0>::new(
            0.0,
            "Test",
            EventDirection::Any,
            AngleFormat::Radians,
        )
        .with_callback(callback);
        assert!(event.callback().is_some());
    }

    #[test]
    fn test_SEccentricAnomalyEvent_callback_none() {
        let event = SEccentricAnomalyEvent::<6, 0>::new(
            0.0,
            "Test",
            EventDirection::Any,
            AngleFormat::Radians,
        );
        assert!(event.callback().is_none());
    }

    #[test]
    fn test_SEccentricAnomalyEvent_set_terminal() {
        let event = SEccentricAnomalyEvent::<6, 0>::new(
            0.0,
            "Test",
            EventDirection::Any,
            AngleFormat::Radians,
        )
        .set_terminal();
        assert_eq!(event.action(), EventAction::Stop);
    }

    #[test]
    fn test_SEccentricAnomalyEvent_direction_increasing() {
        let event = SEccentricAnomalyEvent::<6, 0>::new(
            0.0,
            "EA",
            EventDirection::Increasing,
            AngleFormat::Radians,
        );
        assert_eq!(event.direction(), EventDirection::Increasing);
    }

    #[test]
    fn test_SEccentricAnomalyEvent_direction_decreasing() {
        let event = SEccentricAnomalyEvent::<6, 0>::new(
            0.0,
            "EA",
            EventDirection::Decreasing,
            AngleFormat::Radians,
        );
        assert_eq!(event.direction(), EventDirection::Decreasing);
    }

    #[test]
    fn test_SEccentricAnomalyEvent_direction_any() {
        let event = SEccentricAnomalyEvent::<6, 0>::new(
            0.0,
            "EA",
            EventDirection::Any,
            AngleFormat::Radians,
        );
        assert_eq!(event.direction(), EventDirection::Any);
    }

    #[test]
    fn test_SEccentricAnomalyEvent_builder_chaining() {
        let callback: SEventCallback<6, 0> =
            Box::new(|_t, _state, _params| (None, None, EventAction::Stop));
        let event = SEccentricAnomalyEvent::<6, 0>::new(
            0.0,
            "EA",
            EventDirection::Increasing,
            AngleFormat::Radians,
        )
        .with_instance(1)
        .with_tolerances(1e-5, 1e-8)
        .with_step_reduction_factor(0.1)
        .with_callback(callback)
        .set_terminal();

        assert_eq!(event.name(), "EA 1");
        assert_eq!(event.time_tolerance(), 1e-5);
        assert_eq!(event.step_reduction_factor(), 0.1);
        assert!(event.callback().is_some());
        assert_eq!(event.action(), EventAction::Stop);
    }

    // DEccentricAnomalyEvent additional tests

    #[test]
    fn test_DEccentricAnomalyEvent_with_instance() {
        let event =
            DEccentricAnomalyEvent::new(0.0, "EA", EventDirection::Any, AngleFormat::Degrees)
                .with_instance(2);
        assert_eq!(event.name(), "EA 2");
    }

    #[test]
    fn test_DEccentricAnomalyEvent_with_tolerances() {
        let event =
            DEccentricAnomalyEvent::new(0.0, "Test", EventDirection::Any, AngleFormat::Degrees)
                .with_tolerances(1e-4, 1e-7);
        assert_eq!(event.time_tolerance(), 1e-4);
        assert_eq!(event.value_tolerance(), 1e-7);
    }

    #[test]
    fn test_DEccentricAnomalyEvent_with_step_reduction_factor() {
        let event =
            DEccentricAnomalyEvent::new(0.0, "Test", EventDirection::Any, AngleFormat::Degrees)
                .with_step_reduction_factor(0.25);
        assert_eq!(event.step_reduction_factor(), 0.25);
    }

    #[test]
    fn test_DEccentricAnomalyEvent_with_callback() {
        let callback: DEventCallback =
            Box::new(move |_t, _state, _params| (None, None, EventAction::Continue));
        let event =
            DEccentricAnomalyEvent::new(0.0, "Test", EventDirection::Any, AngleFormat::Degrees)
                .with_callback(callback);
        assert!(event.callback().is_some());
    }

    #[test]
    fn test_DEccentricAnomalyEvent_callback_none() {
        let event =
            DEccentricAnomalyEvent::new(0.0, "Test", EventDirection::Any, AngleFormat::Degrees);
        assert!(event.callback().is_none());
    }

    #[test]
    fn test_DEccentricAnomalyEvent_set_terminal() {
        let event =
            DEccentricAnomalyEvent::new(0.0, "Test", EventDirection::Any, AngleFormat::Degrees)
                .set_terminal();
        assert_eq!(event.action(), EventAction::Stop);
    }

    #[test]
    fn test_DEccentricAnomalyEvent_direction_increasing() {
        let event = DEccentricAnomalyEvent::new(
            0.0,
            "EA",
            EventDirection::Increasing,
            AngleFormat::Degrees,
        );
        assert_eq!(event.direction(), EventDirection::Increasing);
    }

    #[test]
    fn test_DEccentricAnomalyEvent_direction_decreasing() {
        let event = DEccentricAnomalyEvent::new(
            0.0,
            "EA",
            EventDirection::Decreasing,
            AngleFormat::Degrees,
        );
        assert_eq!(event.direction(), EventDirection::Decreasing);
    }

    #[test]
    fn test_DEccentricAnomalyEvent_direction_any() {
        let event =
            DEccentricAnomalyEvent::new(0.0, "EA", EventDirection::Any, AngleFormat::Degrees);
        assert_eq!(event.direction(), EventDirection::Any);
    }

    #[test]
    fn test_DEccentricAnomalyEvent_builder_chaining() {
        let callback: DEventCallback =
            Box::new(|_t, _state, _params| (None, None, EventAction::Stop));
        let event = DEccentricAnomalyEvent::new(
            0.0,
            "EA",
            EventDirection::Decreasing,
            AngleFormat::Degrees,
        )
        .with_instance(1)
        .with_tolerances(1e-5, 1e-8)
        .with_step_reduction_factor(0.1)
        .with_callback(callback)
        .set_terminal();

        assert_eq!(event.name(), "EA 1");
        assert_eq!(event.time_tolerance(), 1e-5);
        assert_eq!(event.step_reduction_factor(), 0.1);
        assert!(event.callback().is_some());
        assert_eq!(event.action(), EventAction::Stop);
    }

    #[test]
    fn test_STrueAnomalyEvent_new() {
        let event = STrueAnomalyEvent::<6, 0>::new(
            0.0,
            "Periapsis TA",
            EventDirection::Increasing,
            AngleFormat::Radians,
        );
        assert_eq!(event.name(), "Periapsis TA");
    }

    #[test]
    fn test_DTrueAnomalyEvent_new() {
        let event = DTrueAnomalyEvent::new(
            180.0,
            "Apoapsis TA",
            EventDirection::Any,
            AngleFormat::Degrees,
        );
        assert_eq!(event.name(), "Apoapsis TA");
    }

    #[test]
    fn test_STrueAnomalyEvent_with_instance() {
        let event =
            STrueAnomalyEvent::<6, 0>::new(0.0, "TA", EventDirection::Any, AngleFormat::Radians)
                .with_instance(3);
        assert_eq!(event.name(), "TA 3");
    }

    #[test]
    fn test_STrueAnomalyEvent_with_tolerances() {
        let event =
            STrueAnomalyEvent::<6, 0>::new(0.0, "Test", EventDirection::Any, AngleFormat::Radians)
                .with_tolerances(1e-4, 1e-6);
        assert_eq!(event.time_tolerance(), 1e-4);
        assert_eq!(event.value_tolerance(), 1e-6);
    }

    #[test]
    fn test_STrueAnomalyEvent_with_step_reduction_factor() {
        let event =
            STrueAnomalyEvent::<6, 0>::new(0.0, "Test", EventDirection::Any, AngleFormat::Radians)
                .with_step_reduction_factor(0.15);
        assert_eq!(event.step_reduction_factor(), 0.15);
    }

    #[test]
    fn test_STrueAnomalyEvent_with_callback() {
        let callback: SEventCallback<6, 0> =
            Box::new(move |_t, _state, _params| (None, None, EventAction::Continue));
        let event =
            STrueAnomalyEvent::<6, 0>::new(0.0, "Test", EventDirection::Any, AngleFormat::Radians)
                .with_callback(callback);
        assert!(event.callback().is_some());
    }

    #[test]
    fn test_STrueAnomalyEvent_callback_none() {
        let event =
            STrueAnomalyEvent::<6, 0>::new(0.0, "Test", EventDirection::Any, AngleFormat::Radians);
        assert!(event.callback().is_none());
    }

    #[test]
    fn test_STrueAnomalyEvent_set_terminal() {
        let event =
            STrueAnomalyEvent::<6, 0>::new(0.0, "Test", EventDirection::Any, AngleFormat::Radians)
                .set_terminal();
        assert_eq!(event.action(), EventAction::Stop);
    }

    #[test]
    fn test_STrueAnomalyEvent_direction_increasing() {
        let event = STrueAnomalyEvent::<6, 0>::new(
            0.0,
            "TA",
            EventDirection::Increasing,
            AngleFormat::Radians,
        );
        assert_eq!(event.direction(), EventDirection::Increasing);
    }

    #[test]
    fn test_STrueAnomalyEvent_direction_decreasing() {
        let event = STrueAnomalyEvent::<6, 0>::new(
            0.0,
            "TA",
            EventDirection::Decreasing,
            AngleFormat::Radians,
        );
        assert_eq!(event.direction(), EventDirection::Decreasing);
    }

    #[test]
    fn test_STrueAnomalyEvent_direction_any() {
        let event =
            STrueAnomalyEvent::<6, 0>::new(0.0, "TA", EventDirection::Any, AngleFormat::Radians);
        assert_eq!(event.direction(), EventDirection::Any);
    }

    #[test]
    fn test_STrueAnomalyEvent_builder_chaining() {
        let callback: SEventCallback<6, 0> =
            Box::new(|_t, _state, _params| (None, None, EventAction::Stop));
        let event = STrueAnomalyEvent::<6, 0>::new(
            0.0,
            "TA",
            EventDirection::Increasing,
            AngleFormat::Radians,
        )
        .with_instance(1)
        .with_tolerances(1e-5, 1e-8)
        .with_step_reduction_factor(0.1)
        .with_callback(callback)
        .set_terminal();

        assert_eq!(event.name(), "TA 1");
        assert_eq!(event.time_tolerance(), 1e-5);
        assert_eq!(event.step_reduction_factor(), 0.1);
        assert!(event.callback().is_some());
        assert_eq!(event.action(), EventAction::Stop);
    }

    // DTrueAnomalyEvent additional tests

    #[test]
    fn test_DTrueAnomalyEvent_with_instance() {
        let event = DTrueAnomalyEvent::new(180.0, "TA", EventDirection::Any, AngleFormat::Degrees)
            .with_instance(2);
        assert_eq!(event.name(), "TA 2");
    }

    #[test]
    fn test_DTrueAnomalyEvent_with_tolerances() {
        let event =
            DTrueAnomalyEvent::new(180.0, "Test", EventDirection::Any, AngleFormat::Degrees)
                .with_tolerances(1e-4, 1e-7);
        assert_eq!(event.time_tolerance(), 1e-4);
        assert_eq!(event.value_tolerance(), 1e-7);
    }

    #[test]
    fn test_DTrueAnomalyEvent_with_step_reduction_factor() {
        let event =
            DTrueAnomalyEvent::new(180.0, "Test", EventDirection::Any, AngleFormat::Degrees)
                .with_step_reduction_factor(0.25);
        assert_eq!(event.step_reduction_factor(), 0.25);
    }

    #[test]
    fn test_DTrueAnomalyEvent_with_callback() {
        let callback: DEventCallback =
            Box::new(move |_t, _state, _params| (None, None, EventAction::Continue));
        let event =
            DTrueAnomalyEvent::new(180.0, "Test", EventDirection::Any, AngleFormat::Degrees)
                .with_callback(callback);
        assert!(event.callback().is_some());
    }

    #[test]
    fn test_DTrueAnomalyEvent_callback_none() {
        let event =
            DTrueAnomalyEvent::new(180.0, "Test", EventDirection::Any, AngleFormat::Degrees);
        assert!(event.callback().is_none());
    }

    #[test]
    fn test_DTrueAnomalyEvent_set_terminal() {
        let event =
            DTrueAnomalyEvent::new(180.0, "Test", EventDirection::Any, AngleFormat::Degrees)
                .set_terminal();
        assert_eq!(event.action(), EventAction::Stop);
    }

    #[test]
    fn test_DTrueAnomalyEvent_direction_increasing() {
        let event = DTrueAnomalyEvent::new(
            180.0,
            "TA",
            EventDirection::Increasing,
            AngleFormat::Degrees,
        );
        assert_eq!(event.direction(), EventDirection::Increasing);
    }

    #[test]
    fn test_DTrueAnomalyEvent_direction_decreasing() {
        let event = DTrueAnomalyEvent::new(
            180.0,
            "TA",
            EventDirection::Decreasing,
            AngleFormat::Degrees,
        );
        assert_eq!(event.direction(), EventDirection::Decreasing);
    }

    #[test]
    fn test_DTrueAnomalyEvent_direction_any() {
        let event = DTrueAnomalyEvent::new(180.0, "TA", EventDirection::Any, AngleFormat::Degrees);
        assert_eq!(event.direction(), EventDirection::Any);
    }

    #[test]
    fn test_DTrueAnomalyEvent_builder_chaining() {
        let callback: DEventCallback =
            Box::new(|_t, _state, _params| (None, None, EventAction::Stop));
        let event = DTrueAnomalyEvent::new(
            180.0,
            "TA",
            EventDirection::Decreasing,
            AngleFormat::Degrees,
        )
        .with_instance(1)
        .with_tolerances(1e-5, 1e-8)
        .with_step_reduction_factor(0.1)
        .with_callback(callback)
        .set_terminal();

        assert_eq!(event.name(), "TA 1");
        assert_eq!(event.time_tolerance(), 1e-5);
        assert_eq!(event.step_reduction_factor(), 0.1);
        assert!(event.callback().is_some());
        assert_eq!(event.action(), EventAction::Stop);
    }

    #[test]
    fn test_SArgumentOfLatitudeEvent_new() {
        let event = SArgumentOfLatitudeEvent::<6, 0>::new(
            0.5,
            "AoL Check",
            EventDirection::Increasing,
            AngleFormat::Radians,
        );
        assert_eq!(event.name(), "AoL Check");
        assert_eq!(event.target_value(), 0.5);
    }

    #[test]
    fn test_DArgumentOfLatitudeEvent_new() {
        let event =
            DArgumentOfLatitudeEvent::new(30.0, "AoL", EventDirection::Any, AngleFormat::Degrees);
        assert_eq!(event.name(), "AoL");
    }

    #[test]
    fn test_SArgumentOfLatitudeEvent_with_instance() {
        let event = SArgumentOfLatitudeEvent::<6, 0>::new(
            0.5,
            "AoL",
            EventDirection::Any,
            AngleFormat::Radians,
        )
        .with_instance(3);
        assert_eq!(event.name(), "AoL 3");
    }

    #[test]
    fn test_SArgumentOfLatitudeEvent_with_tolerances() {
        let event = SArgumentOfLatitudeEvent::<6, 0>::new(
            0.5,
            "Test",
            EventDirection::Any,
            AngleFormat::Radians,
        )
        .with_tolerances(1e-4, 1e-6);
        assert_eq!(event.time_tolerance(), 1e-4);
        assert_eq!(event.value_tolerance(), 1e-6);
    }

    #[test]
    fn test_SArgumentOfLatitudeEvent_with_step_reduction_factor() {
        let event = SArgumentOfLatitudeEvent::<6, 0>::new(
            0.5,
            "Test",
            EventDirection::Any,
            AngleFormat::Radians,
        )
        .with_step_reduction_factor(0.15);
        assert_eq!(event.step_reduction_factor(), 0.15);
    }

    #[test]
    fn test_SArgumentOfLatitudeEvent_with_callback() {
        let callback: SEventCallback<6, 0> =
            Box::new(move |_t, _state, _params| (None, None, EventAction::Continue));
        let event = SArgumentOfLatitudeEvent::<6, 0>::new(
            0.5,
            "Test",
            EventDirection::Any,
            AngleFormat::Radians,
        )
        .with_callback(callback);
        assert!(event.callback().is_some());
    }

    #[test]
    fn test_SArgumentOfLatitudeEvent_callback_none() {
        let event = SArgumentOfLatitudeEvent::<6, 0>::new(
            0.5,
            "Test",
            EventDirection::Any,
            AngleFormat::Radians,
        );
        assert!(event.callback().is_none());
    }

    #[test]
    fn test_SArgumentOfLatitudeEvent_set_terminal() {
        let event = SArgumentOfLatitudeEvent::<6, 0>::new(
            0.5,
            "Test",
            EventDirection::Any,
            AngleFormat::Radians,
        )
        .set_terminal();
        assert_eq!(event.action(), EventAction::Stop);
    }

    #[test]
    fn test_SArgumentOfLatitudeEvent_direction_increasing() {
        let event = SArgumentOfLatitudeEvent::<6, 0>::new(
            0.5,
            "AoL",
            EventDirection::Increasing,
            AngleFormat::Radians,
        );
        assert_eq!(event.direction(), EventDirection::Increasing);
    }

    #[test]
    fn test_SArgumentOfLatitudeEvent_direction_decreasing() {
        let event = SArgumentOfLatitudeEvent::<6, 0>::new(
            0.5,
            "AoL",
            EventDirection::Decreasing,
            AngleFormat::Radians,
        );
        assert_eq!(event.direction(), EventDirection::Decreasing);
    }

    #[test]
    fn test_SArgumentOfLatitudeEvent_direction_any() {
        let event = SArgumentOfLatitudeEvent::<6, 0>::new(
            0.5,
            "AoL",
            EventDirection::Any,
            AngleFormat::Radians,
        );
        assert_eq!(event.direction(), EventDirection::Any);
    }

    #[test]
    fn test_SArgumentOfLatitudeEvent_builder_chaining() {
        let callback: SEventCallback<6, 0> =
            Box::new(|_t, _state, _params| (None, None, EventAction::Stop));
        let event = SArgumentOfLatitudeEvent::<6, 0>::new(
            0.5,
            "AoL",
            EventDirection::Increasing,
            AngleFormat::Radians,
        )
        .with_instance(1)
        .with_tolerances(1e-5, 1e-8)
        .with_step_reduction_factor(0.1)
        .with_callback(callback)
        .set_terminal();

        assert_eq!(event.name(), "AoL 1");
        assert_eq!(event.time_tolerance(), 1e-5);
        assert_eq!(event.step_reduction_factor(), 0.1);
        assert!(event.callback().is_some());
        assert_eq!(event.action(), EventAction::Stop);
    }

    // DArgumentOfLatitudeEvent additional tests

    #[test]
    fn test_DArgumentOfLatitudeEvent_with_instance() {
        let event =
            DArgumentOfLatitudeEvent::new(30.0, "AoL", EventDirection::Any, AngleFormat::Degrees)
                .with_instance(2);
        assert_eq!(event.name(), "AoL 2");
    }

    #[test]
    fn test_DArgumentOfLatitudeEvent_with_tolerances() {
        let event =
            DArgumentOfLatitudeEvent::new(30.0, "Test", EventDirection::Any, AngleFormat::Degrees)
                .with_tolerances(1e-4, 1e-7);
        assert_eq!(event.time_tolerance(), 1e-4);
        assert_eq!(event.value_tolerance(), 1e-7);
    }

    #[test]
    fn test_DArgumentOfLatitudeEvent_with_step_reduction_factor() {
        let event =
            DArgumentOfLatitudeEvent::new(30.0, "Test", EventDirection::Any, AngleFormat::Degrees)
                .with_step_reduction_factor(0.25);
        assert_eq!(event.step_reduction_factor(), 0.25);
    }

    #[test]
    fn test_DArgumentOfLatitudeEvent_with_callback() {
        let callback: DEventCallback =
            Box::new(move |_t, _state, _params| (None, None, EventAction::Continue));
        let event =
            DArgumentOfLatitudeEvent::new(30.0, "Test", EventDirection::Any, AngleFormat::Degrees)
                .with_callback(callback);
        assert!(event.callback().is_some());
    }

    #[test]
    fn test_DArgumentOfLatitudeEvent_callback_none() {
        let event =
            DArgumentOfLatitudeEvent::new(30.0, "Test", EventDirection::Any, AngleFormat::Degrees);
        assert!(event.callback().is_none());
    }

    #[test]
    fn test_DArgumentOfLatitudeEvent_set_terminal() {
        let event =
            DArgumentOfLatitudeEvent::new(30.0, "Test", EventDirection::Any, AngleFormat::Degrees)
                .set_terminal();
        assert_eq!(event.action(), EventAction::Stop);
    }

    #[test]
    fn test_DArgumentOfLatitudeEvent_direction_increasing() {
        let event = DArgumentOfLatitudeEvent::new(
            30.0,
            "AoL",
            EventDirection::Increasing,
            AngleFormat::Degrees,
        );
        assert_eq!(event.direction(), EventDirection::Increasing);
    }

    #[test]
    fn test_DArgumentOfLatitudeEvent_direction_decreasing() {
        let event = DArgumentOfLatitudeEvent::new(
            30.0,
            "AoL",
            EventDirection::Decreasing,
            AngleFormat::Degrees,
        );
        assert_eq!(event.direction(), EventDirection::Decreasing);
    }

    #[test]
    fn test_DArgumentOfLatitudeEvent_direction_any() {
        let event =
            DArgumentOfLatitudeEvent::new(30.0, "AoL", EventDirection::Any, AngleFormat::Degrees);
        assert_eq!(event.direction(), EventDirection::Any);
    }

    #[test]
    fn test_DArgumentOfLatitudeEvent_builder_chaining() {
        let callback: DEventCallback =
            Box::new(|_t, _state, _params| (None, None, EventAction::Stop));
        let event = DArgumentOfLatitudeEvent::new(
            30.0,
            "AoL",
            EventDirection::Decreasing,
            AngleFormat::Degrees,
        )
        .with_instance(1)
        .with_tolerances(1e-5, 1e-8)
        .with_step_reduction_factor(0.1)
        .with_callback(callback)
        .set_terminal();

        assert_eq!(event.name(), "AoL 1");
        assert_eq!(event.time_tolerance(), 1e-5);
        assert_eq!(event.step_reduction_factor(), 0.1);
        assert!(event.callback().is_some());
        assert_eq!(event.action(), EventAction::Stop);
    }

    // =========================================================================
    // Node Crossing Event Tests
    // =========================================================================

    #[test]
    fn test_SAscendingNodeEvent_new() {
        let event = SAscendingNodeEvent::<6, 0>::new("Ascending Node");
        assert_eq!(event.name(), "Ascending Node");
        assert_eq!(event.target_value(), 0.0); // Zero crossing
        assert_eq!(event.direction(), EventDirection::Increasing);
    }

    #[test]
    fn test_SAscendingNodeEvent_builder_chaining() {
        let event = SAscendingNodeEvent::<6, 0>::new("Asc Node")
            .with_instance(2)
            .with_tolerances(1e-5, 1e-8)
            .with_step_reduction_factor(0.15);
        assert_eq!(event.name(), "Asc Node 2");
        assert_eq!(event.time_tolerance(), 1e-5);
        assert_eq!(event.step_reduction_factor(), 0.15);
    }

    #[test]
    fn test_DAscendingNodeEvent_new() {
        let event = DAscendingNodeEvent::new("Ascending");
        assert_eq!(event.name(), "Ascending");
        assert_eq!(event.target_value(), 0.0);
    }

    #[test]
    fn test_SDescendingNodeEvent_new() {
        let event = SDescendingNodeEvent::<6, 0>::new("Descending Node");
        assert_eq!(event.name(), "Descending Node");
        assert_eq!(event.target_value(), 0.0); // Zero crossing (shifted)
        assert_eq!(event.direction(), EventDirection::Increasing);
    }

    #[test]
    fn test_DDescendingNodeEvent_new() {
        let event = DDescendingNodeEvent::new("Descending");
        assert_eq!(event.name(), "Descending");
    }

    // SAscendingNodeEvent additional tests

    #[test]
    fn test_SAscendingNodeEvent_with_instance() {
        let event = SAscendingNodeEvent::<6, 0>::new("Asc Node").with_instance(3);
        assert_eq!(event.name(), "Asc Node 3");
    }

    #[test]
    fn test_SAscendingNodeEvent_with_tolerances() {
        let event = SAscendingNodeEvent::<6, 0>::new("Test").with_tolerances(1e-4, 1e-6);
        assert_eq!(event.time_tolerance(), 1e-4);
        assert_eq!(event.value_tolerance(), 1e-6);
    }

    #[test]
    fn test_SAscendingNodeEvent_with_step_reduction_factor() {
        let event = SAscendingNodeEvent::<6, 0>::new("Test").with_step_reduction_factor(0.15);
        assert_eq!(event.step_reduction_factor(), 0.15);
    }

    #[test]
    fn test_SAscendingNodeEvent_with_callback() {
        let callback: SEventCallback<6, 0> =
            Box::new(move |_t, _state, _params| (None, None, EventAction::Continue));
        let event = SAscendingNodeEvent::<6, 0>::new("Test").with_callback(callback);
        assert!(event.callback().is_some());
    }

    #[test]
    fn test_SAscendingNodeEvent_callback_none() {
        let event = SAscendingNodeEvent::<6, 0>::new("Test");
        assert!(event.callback().is_none());
    }

    #[test]
    fn test_SAscendingNodeEvent_set_terminal() {
        let event = SAscendingNodeEvent::<6, 0>::new("Test").set_terminal();
        assert_eq!(event.action(), EventAction::Stop);
    }

    // DAscendingNodeEvent additional tests

    #[test]
    fn test_DAscendingNodeEvent_with_instance() {
        let event = DAscendingNodeEvent::new("Asc Node").with_instance(2);
        assert_eq!(event.name(), "Asc Node 2");
    }

    #[test]
    fn test_DAscendingNodeEvent_with_tolerances() {
        let event = DAscendingNodeEvent::new("Test").with_tolerances(1e-4, 1e-7);
        assert_eq!(event.time_tolerance(), 1e-4);
        assert_eq!(event.value_tolerance(), 1e-7);
    }

    #[test]
    fn test_DAscendingNodeEvent_with_step_reduction_factor() {
        let event = DAscendingNodeEvent::new("Test").with_step_reduction_factor(0.25);
        assert_eq!(event.step_reduction_factor(), 0.25);
    }

    #[test]
    fn test_DAscendingNodeEvent_with_callback() {
        let callback: DEventCallback =
            Box::new(move |_t, _state, _params| (None, None, EventAction::Continue));
        let event = DAscendingNodeEvent::new("Test").with_callback(callback);
        assert!(event.callback().is_some());
    }

    #[test]
    fn test_DAscendingNodeEvent_callback_none() {
        let event = DAscendingNodeEvent::new("Test");
        assert!(event.callback().is_none());
    }

    #[test]
    fn test_DAscendingNodeEvent_set_terminal() {
        let event = DAscendingNodeEvent::new("Test").set_terminal();
        assert_eq!(event.action(), EventAction::Stop);
    }

    #[test]
    fn test_DAscendingNodeEvent_builder_chaining() {
        let callback: DEventCallback =
            Box::new(|_t, _state, _params| (None, None, EventAction::Stop));
        let event = DAscendingNodeEvent::new("Asc Node")
            .with_instance(1)
            .with_tolerances(1e-5, 1e-8)
            .with_step_reduction_factor(0.1)
            .with_callback(callback)
            .set_terminal();

        assert_eq!(event.name(), "Asc Node 1");
        assert_eq!(event.time_tolerance(), 1e-5);
        assert_eq!(event.step_reduction_factor(), 0.1);
        assert!(event.callback().is_some());
        assert_eq!(event.action(), EventAction::Stop);
    }

    // SDescendingNodeEvent additional tests

    #[test]
    fn test_SDescendingNodeEvent_with_instance() {
        let event = SDescendingNodeEvent::<6, 0>::new("Desc Node").with_instance(3);
        assert_eq!(event.name(), "Desc Node 3");
    }

    #[test]
    fn test_SDescendingNodeEvent_with_tolerances() {
        let event = SDescendingNodeEvent::<6, 0>::new("Test").with_tolerances(1e-4, 1e-6);
        assert_eq!(event.time_tolerance(), 1e-4);
        assert_eq!(event.value_tolerance(), 1e-6);
    }

    #[test]
    fn test_SDescendingNodeEvent_with_step_reduction_factor() {
        let event = SDescendingNodeEvent::<6, 0>::new("Test").with_step_reduction_factor(0.15);
        assert_eq!(event.step_reduction_factor(), 0.15);
    }

    #[test]
    fn test_SDescendingNodeEvent_with_callback() {
        let callback: SEventCallback<6, 0> =
            Box::new(move |_t, _state, _params| (None, None, EventAction::Continue));
        let event = SDescendingNodeEvent::<6, 0>::new("Test").with_callback(callback);
        assert!(event.callback().is_some());
    }

    #[test]
    fn test_SDescendingNodeEvent_callback_none() {
        let event = SDescendingNodeEvent::<6, 0>::new("Test");
        assert!(event.callback().is_none());
    }

    #[test]
    fn test_SDescendingNodeEvent_set_terminal() {
        let event = SDescendingNodeEvent::<6, 0>::new("Test").set_terminal();
        assert_eq!(event.action(), EventAction::Stop);
    }

    #[test]
    fn test_SDescendingNodeEvent_builder_chaining() {
        let callback: SEventCallback<6, 0> =
            Box::new(|_t, _state, _params| (None, None, EventAction::Stop));
        let event = SDescendingNodeEvent::<6, 0>::new("Desc Node")
            .with_instance(1)
            .with_tolerances(1e-5, 1e-8)
            .with_step_reduction_factor(0.1)
            .with_callback(callback)
            .set_terminal();

        assert_eq!(event.name(), "Desc Node 1");
        assert_eq!(event.time_tolerance(), 1e-5);
        assert_eq!(event.step_reduction_factor(), 0.1);
        assert!(event.callback().is_some());
        assert_eq!(event.action(), EventAction::Stop);
    }

    // DDescendingNodeEvent additional tests

    #[test]
    fn test_DDescendingNodeEvent_with_instance() {
        let event = DDescendingNodeEvent::new("Desc Node").with_instance(2);
        assert_eq!(event.name(), "Desc Node 2");
    }

    #[test]
    fn test_DDescendingNodeEvent_with_tolerances() {
        let event = DDescendingNodeEvent::new("Test").with_tolerances(1e-4, 1e-7);
        assert_eq!(event.time_tolerance(), 1e-4);
        assert_eq!(event.value_tolerance(), 1e-7);
    }

    #[test]
    fn test_DDescendingNodeEvent_with_step_reduction_factor() {
        let event = DDescendingNodeEvent::new("Test").with_step_reduction_factor(0.25);
        assert_eq!(event.step_reduction_factor(), 0.25);
    }

    #[test]
    fn test_DDescendingNodeEvent_with_callback() {
        let callback: DEventCallback =
            Box::new(move |_t, _state, _params| (None, None, EventAction::Continue));
        let event = DDescendingNodeEvent::new("Test").with_callback(callback);
        assert!(event.callback().is_some());
    }

    #[test]
    fn test_DDescendingNodeEvent_callback_none() {
        let event = DDescendingNodeEvent::new("Test");
        assert!(event.callback().is_none());
    }

    #[test]
    fn test_DDescendingNodeEvent_set_terminal() {
        let event = DDescendingNodeEvent::new("Test").set_terminal();
        assert_eq!(event.action(), EventAction::Stop);
    }

    #[test]
    fn test_DDescendingNodeEvent_builder_chaining() {
        let callback: DEventCallback =
            Box::new(|_t, _state, _params| (None, None, EventAction::Stop));
        let event = DDescendingNodeEvent::new("Desc Node")
            .with_instance(1)
            .with_tolerances(1e-5, 1e-8)
            .with_step_reduction_factor(0.1)
            .with_callback(callback)
            .set_terminal();

        assert_eq!(event.name(), "Desc Node 1");
        assert_eq!(event.time_tolerance(), 1e-5);
        assert_eq!(event.step_reduction_factor(), 0.1);
        assert!(event.callback().is_some());
        assert_eq!(event.action(), EventAction::Stop);
    }

    // =========================================================================
    // State-Derived Event Tests
    // =========================================================================

    #[test]
    fn test_SSpeedEvent_new() {
        let event = SSpeedEvent::<6, 0>::new(7500.0, "Speed value", EventDirection::Increasing);
        assert_eq!(event.name(), "Speed value");
        assert_eq!(event.target_value(), 7500.0);
    }

    #[test]
    fn test_SSpeedEvent_evaluate() {
        let event = SSpeedEvent::<6, 0>::new(7000.0, "Speed", EventDirection::Any);
        let epoch = Epoch::from_jd(2451545.0, TimeSystem::UTC);

        // State with known velocity: [3000, 4000, 0] => |v| = 5000 m/s
        let state = Vector6::new(7000e3, 0.0, 0.0, 3000.0, 4000.0, 0.0);
        let val = event.evaluate(epoch, &state, None);
        // evaluate() returns raw speed value (5000), not relative to value
        assert!((val - 5000.0).abs() < 1e-6);
    }

    #[test]
    fn test_DSpeedEvent_new() {
        let event = DSpeedEvent::new(8000.0, "Speed", EventDirection::Decreasing);
        assert_eq!(event.name(), "Speed");
        assert_eq!(event.target_value(), 8000.0);
    }

    // SSpeedEvent additional tests

    #[test]
    fn test_SSpeedEvent_with_instance() {
        let event = SSpeedEvent::<6, 0>::new(7500.0, "Speed", EventDirection::Any).with_instance(3);
        assert_eq!(event.name(), "Speed 3");
    }

    #[test]
    fn test_SSpeedEvent_with_tolerances() {
        let event = SSpeedEvent::<6, 0>::new(7500.0, "Test", EventDirection::Any)
            .with_tolerances(1e-4, 1e-6);
        assert_eq!(event.time_tolerance(), 1e-4);
        assert_eq!(event.value_tolerance(), 1e-6);
    }

    #[test]
    fn test_SSpeedEvent_with_step_reduction_factor() {
        let event = SSpeedEvent::<6, 0>::new(7500.0, "Test", EventDirection::Any)
            .with_step_reduction_factor(0.15);
        assert_eq!(event.step_reduction_factor(), 0.15);
    }

    #[test]
    fn test_SSpeedEvent_with_callback() {
        let callback: SEventCallback<6, 0> =
            Box::new(move |_t, _state, _params| (None, None, EventAction::Continue));
        let event =
            SSpeedEvent::<6, 0>::new(7500.0, "Test", EventDirection::Any).with_callback(callback);
        assert!(event.callback().is_some());
    }

    #[test]
    fn test_SSpeedEvent_callback_none() {
        let event = SSpeedEvent::<6, 0>::new(7500.0, "Test", EventDirection::Any);
        assert!(event.callback().is_none());
    }

    #[test]
    fn test_SSpeedEvent_set_terminal() {
        let event = SSpeedEvent::<6, 0>::new(7500.0, "Test", EventDirection::Any).set_terminal();
        assert_eq!(event.action(), EventAction::Stop);
    }

    #[test]
    fn test_SSpeedEvent_direction_increasing() {
        let event = SSpeedEvent::<6, 0>::new(7500.0, "Speed", EventDirection::Increasing);
        assert_eq!(event.direction(), EventDirection::Increasing);
    }

    #[test]
    fn test_SSpeedEvent_direction_decreasing() {
        let event = SSpeedEvent::<6, 0>::new(7500.0, "Speed", EventDirection::Decreasing);
        assert_eq!(event.direction(), EventDirection::Decreasing);
    }

    #[test]
    fn test_SSpeedEvent_direction_any() {
        let event = SSpeedEvent::<6, 0>::new(7500.0, "Speed", EventDirection::Any);
        assert_eq!(event.direction(), EventDirection::Any);
    }

    #[test]
    fn test_SSpeedEvent_builder_chaining() {
        let callback: SEventCallback<6, 0> =
            Box::new(|_t, _state, _params| (None, None, EventAction::Stop));
        let event = SSpeedEvent::<6, 0>::new(7500.0, "Speed", EventDirection::Increasing)
            .with_instance(1)
            .with_tolerances(1e-5, 1e-8)
            .with_step_reduction_factor(0.1)
            .with_callback(callback)
            .set_terminal();

        assert_eq!(event.name(), "Speed 1");
        assert_eq!(event.time_tolerance(), 1e-5);
        assert_eq!(event.step_reduction_factor(), 0.1);
        assert!(event.callback().is_some());
        assert_eq!(event.action(), EventAction::Stop);
    }

    // DSpeedEvent additional tests

    #[test]
    fn test_DSpeedEvent_evaluate() {
        let event = DSpeedEvent::new(7000.0, "Speed", EventDirection::Any);
        let epoch = Epoch::from_jd(2451545.0, TimeSystem::UTC);

        let state = DVector::from_vec(vec![7000e3, 0.0, 0.0, 3000.0, 4000.0, 0.0]);
        let val = event.evaluate(epoch, &state, None);
        assert!((val - 5000.0).abs() < 1e-6);
    }

    #[test]
    fn test_DSpeedEvent_with_instance() {
        let event = DSpeedEvent::new(7500.0, "Speed", EventDirection::Any).with_instance(2);
        assert_eq!(event.name(), "Speed 2");
    }

    #[test]
    fn test_DSpeedEvent_with_tolerances() {
        let event =
            DSpeedEvent::new(7500.0, "Test", EventDirection::Any).with_tolerances(1e-4, 1e-7);
        assert_eq!(event.time_tolerance(), 1e-4);
        assert_eq!(event.value_tolerance(), 1e-7);
    }

    #[test]
    fn test_DSpeedEvent_with_step_reduction_factor() {
        let event =
            DSpeedEvent::new(7500.0, "Test", EventDirection::Any).with_step_reduction_factor(0.25);
        assert_eq!(event.step_reduction_factor(), 0.25);
    }

    #[test]
    fn test_DSpeedEvent_with_callback() {
        let callback: DEventCallback =
            Box::new(move |_t, _state, _params| (None, None, EventAction::Continue));
        let event = DSpeedEvent::new(7500.0, "Test", EventDirection::Any).with_callback(callback);
        assert!(event.callback().is_some());
    }

    #[test]
    fn test_DSpeedEvent_callback_none() {
        let event = DSpeedEvent::new(7500.0, "Test", EventDirection::Any);
        assert!(event.callback().is_none());
    }

    #[test]
    fn test_DSpeedEvent_set_terminal() {
        let event = DSpeedEvent::new(7500.0, "Test", EventDirection::Any).set_terminal();
        assert_eq!(event.action(), EventAction::Stop);
    }

    #[test]
    fn test_DSpeedEvent_direction_increasing() {
        let event = DSpeedEvent::new(7500.0, "Speed", EventDirection::Increasing);
        assert_eq!(event.direction(), EventDirection::Increasing);
    }

    #[test]
    fn test_DSpeedEvent_direction_decreasing() {
        let event = DSpeedEvent::new(7500.0, "Speed", EventDirection::Decreasing);
        assert_eq!(event.direction(), EventDirection::Decreasing);
    }

    #[test]
    fn test_DSpeedEvent_direction_any() {
        let event = DSpeedEvent::new(7500.0, "Speed", EventDirection::Any);
        assert_eq!(event.direction(), EventDirection::Any);
    }

    #[test]
    fn test_DSpeedEvent_builder_chaining() {
        let callback: DEventCallback =
            Box::new(|_t, _state, _params| (None, None, EventAction::Stop));
        let event = DSpeedEvent::new(7500.0, "Speed", EventDirection::Decreasing)
            .with_instance(1)
            .with_tolerances(1e-5, 1e-8)
            .with_step_reduction_factor(0.1)
            .with_callback(callback)
            .set_terminal();

        assert_eq!(event.name(), "Speed 1");
        assert_eq!(event.time_tolerance(), 1e-5);
        assert_eq!(event.step_reduction_factor(), 0.1);
        assert!(event.callback().is_some());
        assert_eq!(event.action(), EventAction::Stop);
    }

    #[test]
    fn test_SLongitudeEvent_new() {
        setup_global_test_eop();
        let event = SLongitudeEvent::<6, 0>::new(
            0.0,
            "Prime Meridian",
            EventDirection::Any,
            AngleFormat::Degrees,
        );
        assert_eq!(event.name(), "Prime Meridian");
        assert_eq!(event.target_value(), 0.0);
    }

    #[test]
    fn test_DLongitudeEvent_new() {
        setup_global_test_eop();
        let event = DLongitudeEvent::new(
            45.0,
            "Lon Check",
            EventDirection::Increasing,
            AngleFormat::Degrees,
        );
        assert_eq!(event.name(), "Lon Check");
    }

    #[test]
    #[serial]
    fn test_SLongitudeEvent_evaluate() {
        setup_global_test_eop();
        // State vector: ~500 km altitude orbit
        let state = Vector6::new(R_EARTH + 500e3, 0.0, 0.0, 0.0, 7612.0, 0.0);
        let event = SLongitudeEvent::<6, 0>::new(
            0.0,
            "Prime Meridian",
            EventDirection::Any,
            AngleFormat::Degrees,
        );
        // Evaluate at a specific time (longitude depends on satellite position and time)
        let epoch = Epoch::from_datetime(2024, 1, 1, 0, 0, 0.0, 0.0, TimeSystem::UTC);
        let result = event.evaluate(epoch, &state, None);
        // Result is longitude minus target (0.0 radians)
        assert!(result.is_finite());
    }

    #[test]
    fn test_SLongitudeEvent_target_value() {
        setup_global_test_eop();
        let event = SLongitudeEvent::<6, 0>::new(
            45.0,
            "Longitude",
            EventDirection::Any,
            AngleFormat::Degrees,
        );
        assert_eq!(event.target_value(), 45.0_f64.to_radians());
    }

    #[test]
    fn test_SLongitudeEvent_with_instance() {
        setup_global_test_eop();
        let event = SLongitudeEvent::<6, 0>::new(
            0.0,
            "Longitude",
            EventDirection::Any,
            AngleFormat::Degrees,
        )
        .with_instance(3);
        assert_eq!(event.name(), "Longitude 3");
    }

    #[test]
    fn test_SLongitudeEvent_with_tolerances() {
        setup_global_test_eop();
        let event = SLongitudeEvent::<6, 0>::new(
            0.0,
            "Longitude",
            EventDirection::Any,
            AngleFormat::Degrees,
        )
        .with_tolerances(1e-6, 1e-9);
        assert_eq!(event.time_tolerance(), 1e-6);
        assert_eq!(event.value_tolerance(), 1e-9);
    }

    #[test]
    fn test_SLongitudeEvent_time_tolerance() {
        setup_global_test_eop();
        let event = SLongitudeEvent::<6, 0>::new(
            0.0,
            "Longitude",
            EventDirection::Any,
            AngleFormat::Degrees,
        );
        // Default time tolerance
        assert!(event.time_tolerance() > 0.0);
    }

    #[test]
    fn test_SLongitudeEvent_value_tolerance() {
        setup_global_test_eop();
        let event = SLongitudeEvent::<6, 0>::new(
            0.0,
            "Longitude",
            EventDirection::Any,
            AngleFormat::Degrees,
        );
        // Default value tolerance
        assert!(event.value_tolerance() > 0.0);
    }

    #[test]
    fn test_SLongitudeEvent_with_step_reduction_factor() {
        setup_global_test_eop();
        let event = SLongitudeEvent::<6, 0>::new(
            0.0,
            "Longitude",
            EventDirection::Any,
            AngleFormat::Degrees,
        )
        .with_step_reduction_factor(0.25);
        assert_eq!(event.step_reduction_factor(), 0.25);
    }

    #[test]
    fn test_SLongitudeEvent_step_reduction_factor() {
        setup_global_test_eop();
        let event = SLongitudeEvent::<6, 0>::new(
            0.0,
            "Longitude",
            EventDirection::Any,
            AngleFormat::Degrees,
        );
        // Default step reduction factor
        assert!(event.step_reduction_factor() > 0.0);
        assert!(event.step_reduction_factor() < 1.0);
    }

    #[test]
    fn test_SLongitudeEvent_with_callback() {
        setup_global_test_eop();
        let callback: SEventCallback<6, 0> =
            Box::new(move |_t, _state, _params| (None, None, EventAction::Continue));
        let event = SLongitudeEvent::<6, 0>::new(
            0.0,
            "Longitude",
            EventDirection::Any,
            AngleFormat::Degrees,
        )
        .with_callback(callback);
        assert!(event.callback().is_some());
    }

    #[test]
    fn test_SLongitudeEvent_callback_none() {
        setup_global_test_eop();
        let event = SLongitudeEvent::<6, 0>::new(
            0.0,
            "Longitude",
            EventDirection::Any,
            AngleFormat::Degrees,
        );
        assert!(event.callback().is_none());
    }

    #[test]
    fn test_SLongitudeEvent_set_terminal() {
        setup_global_test_eop();
        let event = SLongitudeEvent::<6, 0>::new(
            0.0,
            "Longitude",
            EventDirection::Any,
            AngleFormat::Degrees,
        )
        .set_terminal();
        assert_eq!(event.action(), EventAction::Stop);
    }

    #[test]
    fn test_SLongitudeEvent_action_continue() {
        setup_global_test_eop();
        let event = SLongitudeEvent::<6, 0>::new(
            0.0,
            "Longitude",
            EventDirection::Any,
            AngleFormat::Degrees,
        );
        assert_eq!(event.action(), EventAction::Continue);
    }

    #[test]
    fn test_SLongitudeEvent_action_stop() {
        setup_global_test_eop();
        let event = SLongitudeEvent::<6, 0>::new(
            0.0,
            "Longitude",
            EventDirection::Any,
            AngleFormat::Degrees,
        )
        .set_terminal();
        assert_eq!(event.action(), EventAction::Stop);
    }

    #[test]
    fn test_SLongitudeEvent_direction_increasing() {
        setup_global_test_eop();
        let event = SLongitudeEvent::<6, 0>::new(
            0.0,
            "Longitude",
            EventDirection::Increasing,
            AngleFormat::Degrees,
        );
        assert_eq!(event.direction(), EventDirection::Increasing);
    }

    #[test]
    fn test_SLongitudeEvent_direction_decreasing() {
        setup_global_test_eop();
        let event = SLongitudeEvent::<6, 0>::new(
            0.0,
            "Longitude",
            EventDirection::Decreasing,
            AngleFormat::Degrees,
        );
        assert_eq!(event.direction(), EventDirection::Decreasing);
    }

    #[test]
    fn test_SLongitudeEvent_direction_any() {
        setup_global_test_eop();
        let event = SLongitudeEvent::<6, 0>::new(
            0.0,
            "Longitude",
            EventDirection::Any,
            AngleFormat::Degrees,
        );
        assert_eq!(event.direction(), EventDirection::Any);
    }

    #[test]
    fn test_SLongitudeEvent_angle_format_radians() {
        setup_global_test_eop();
        let event = SLongitudeEvent::<6, 0>::new(
            std::f64::consts::FRAC_PI_4,
            "Longitude",
            EventDirection::Any,
            AngleFormat::Radians,
        );
        assert_eq!(event.target_value(), std::f64::consts::FRAC_PI_4);
    }

    #[test]
    fn test_SLongitudeEvent_builder_chaining() {
        setup_global_test_eop();
        let callback: SEventCallback<6, 0> =
            Box::new(|_t, _state, _params| (None, None, EventAction::Stop));
        let event = SLongitudeEvent::<6, 0>::new(
            -74.0,
            "NYC Longitude",
            EventDirection::Increasing,
            AngleFormat::Degrees,
        )
        .with_instance(1)
        .with_tolerances(1e-5, 1e-8)
        .with_step_reduction_factor(0.1)
        .with_callback(callback)
        .set_terminal();

        assert_eq!(event.name(), "NYC Longitude 1");
        assert_eq!(event.time_tolerance(), 1e-5);
        assert_eq!(event.step_reduction_factor(), 0.1);
        assert!(event.callback().is_some());
        assert_eq!(event.action(), EventAction::Stop);
    }

    // DLongitudeEvent comprehensive tests

    #[test]
    #[serial]
    fn test_DLongitudeEvent_evaluate() {
        setup_global_test_eop();
        let state = DVector::from_vec(vec![R_EARTH + 500e3, 0.0, 0.0, 0.0, 7612.0, 0.0]);
        let event =
            DLongitudeEvent::new(0.0, "Longitude", EventDirection::Any, AngleFormat::Degrees);
        let epoch = Epoch::from_datetime(2024, 1, 1, 0, 0, 0.0, 0.0, TimeSystem::UTC);
        let result = event.evaluate(epoch, &state, None);
        assert!(result.is_finite());
    }

    #[test]
    fn test_DLongitudeEvent_target_value() {
        setup_global_test_eop();
        let event =
            DLongitudeEvent::new(90.0, "Longitude", EventDirection::Any, AngleFormat::Degrees);
        assert_eq!(event.target_value(), 90.0_f64.to_radians());
    }

    #[test]
    fn test_DLongitudeEvent_with_instance() {
        setup_global_test_eop();
        let event =
            DLongitudeEvent::new(0.0, "Longitude", EventDirection::Any, AngleFormat::Degrees)
                .with_instance(2);
        assert_eq!(event.name(), "Longitude 2");
    }

    #[test]
    fn test_DLongitudeEvent_with_tolerances() {
        setup_global_test_eop();
        let event =
            DLongitudeEvent::new(0.0, "Longitude", EventDirection::Any, AngleFormat::Degrees)
                .with_tolerances(1e-6, 1e-9);
        assert_eq!(event.time_tolerance(), 1e-6);
        assert_eq!(event.value_tolerance(), 1e-9);
    }

    #[test]
    fn test_DLongitudeEvent_time_tolerance() {
        setup_global_test_eop();
        let event =
            DLongitudeEvent::new(0.0, "Longitude", EventDirection::Any, AngleFormat::Degrees);
        assert!(event.time_tolerance() > 0.0);
    }

    #[test]
    fn test_DLongitudeEvent_value_tolerance() {
        setup_global_test_eop();
        let event =
            DLongitudeEvent::new(0.0, "Longitude", EventDirection::Any, AngleFormat::Degrees);
        assert!(event.value_tolerance() > 0.0);
    }

    #[test]
    fn test_DLongitudeEvent_with_step_reduction_factor() {
        setup_global_test_eop();
        let event =
            DLongitudeEvent::new(0.0, "Longitude", EventDirection::Any, AngleFormat::Degrees)
                .with_step_reduction_factor(0.15);
        assert_eq!(event.step_reduction_factor(), 0.15);
    }

    #[test]
    fn test_DLongitudeEvent_step_reduction_factor() {
        setup_global_test_eop();
        let event =
            DLongitudeEvent::new(0.0, "Longitude", EventDirection::Any, AngleFormat::Degrees);
        assert!(event.step_reduction_factor() > 0.0);
        assert!(event.step_reduction_factor() < 1.0);
    }

    #[test]
    fn test_DLongitudeEvent_with_callback() {
        setup_global_test_eop();
        let callback: DEventCallback =
            Box::new(move |_t, _state, _params| (None, None, EventAction::Continue));
        let event =
            DLongitudeEvent::new(0.0, "Longitude", EventDirection::Any, AngleFormat::Degrees)
                .with_callback(callback);
        assert!(event.callback().is_some());
    }

    #[test]
    fn test_DLongitudeEvent_callback_none() {
        setup_global_test_eop();
        let event =
            DLongitudeEvent::new(0.0, "Longitude", EventDirection::Any, AngleFormat::Degrees);
        assert!(event.callback().is_none());
    }

    #[test]
    fn test_DLongitudeEvent_set_terminal() {
        setup_global_test_eop();
        let event =
            DLongitudeEvent::new(0.0, "Longitude", EventDirection::Any, AngleFormat::Degrees)
                .set_terminal();
        assert_eq!(event.action(), EventAction::Stop);
    }

    #[test]
    fn test_DLongitudeEvent_action_continue() {
        setup_global_test_eop();
        let event =
            DLongitudeEvent::new(0.0, "Longitude", EventDirection::Any, AngleFormat::Degrees);
        assert_eq!(event.action(), EventAction::Continue);
    }

    #[test]
    fn test_DLongitudeEvent_action_stop() {
        setup_global_test_eop();
        let event =
            DLongitudeEvent::new(0.0, "Longitude", EventDirection::Any, AngleFormat::Degrees)
                .set_terminal();
        assert_eq!(event.action(), EventAction::Stop);
    }

    #[test]
    fn test_DLongitudeEvent_direction_increasing() {
        setup_global_test_eop();
        let event = DLongitudeEvent::new(
            0.0,
            "Longitude",
            EventDirection::Increasing,
            AngleFormat::Degrees,
        );
        assert_eq!(event.direction(), EventDirection::Increasing);
    }

    #[test]
    fn test_DLongitudeEvent_direction_decreasing() {
        setup_global_test_eop();
        let event = DLongitudeEvent::new(
            0.0,
            "Longitude",
            EventDirection::Decreasing,
            AngleFormat::Degrees,
        );
        assert_eq!(event.direction(), EventDirection::Decreasing);
    }

    #[test]
    fn test_DLongitudeEvent_direction_any() {
        setup_global_test_eop();
        let event =
            DLongitudeEvent::new(0.0, "Longitude", EventDirection::Any, AngleFormat::Degrees);
        assert_eq!(event.direction(), EventDirection::Any);
    }

    #[test]
    fn test_DLongitudeEvent_angle_format_radians() {
        setup_global_test_eop();
        let event = DLongitudeEvent::new(
            std::f64::consts::PI,
            "Longitude",
            EventDirection::Any,
            AngleFormat::Radians,
        );
        assert_eq!(event.target_value(), std::f64::consts::PI);
    }

    #[test]
    fn test_DLongitudeEvent_builder_chaining() {
        setup_global_test_eop();
        let callback: DEventCallback =
            Box::new(|_t, _state, _params| (None, None, EventAction::Stop));
        let event = DLongitudeEvent::new(
            139.0,
            "Tokyo Longitude",
            EventDirection::Decreasing,
            AngleFormat::Degrees,
        )
        .with_instance(5)
        .with_tolerances(1e-5, 1e-8)
        .with_step_reduction_factor(0.2)
        .with_callback(callback)
        .set_terminal();

        assert_eq!(event.name(), "Tokyo Longitude 5");
        assert_eq!(event.time_tolerance(), 1e-5);
        assert_eq!(event.step_reduction_factor(), 0.2);
        assert!(event.callback().is_some());
        assert_eq!(event.action(), EventAction::Stop);
    }

    #[test]
    fn test_SLatitudeEvent_new() {
        setup_global_test_eop();
        let event = SLatitudeEvent::<6, 0>::new(
            30.0,
            "Latitude Check",
            EventDirection::Any,
            AngleFormat::Degrees,
        );
        assert_eq!(event.name(), "Latitude Check");
        assert_eq!(event.target_value(), 30.0_f64.to_radians());
    }

    #[test]
    fn test_DLatitudeEvent_new() {
        setup_global_test_eop();
        let event = DLatitudeEvent::new(
            0.0,
            "Equator",
            EventDirection::Increasing,
            AngleFormat::Degrees,
        );
        assert_eq!(event.name(), "Equator");
    }

    #[test]
    #[serial]
    fn test_SLatitudeEvent_evaluate() {
        setup_global_test_eop();
        // Inclined orbit state vector
        let state = Vector6::new(6800e3, 0.0, 500e3, 0.0, 7000.0, 1000.0);
        let event = SLatitudeEvent::<6, 0>::new(
            0.0,
            "Equator Crossing",
            EventDirection::Any,
            AngleFormat::Degrees,
        );
        let epoch = Epoch::from_datetime(2024, 1, 1, 0, 0, 0.0, 0.0, TimeSystem::UTC);
        let result = event.evaluate(epoch, &state, None);
        assert!(result.is_finite());
    }

    #[test]
    fn test_SLatitudeEvent_target_value() {
        setup_global_test_eop();
        let event = SLatitudeEvent::<6, 0>::new(
            45.0,
            "Latitude",
            EventDirection::Any,
            AngleFormat::Degrees,
        );
        assert_eq!(event.target_value(), 45.0_f64.to_radians());
    }

    #[test]
    fn test_SLatitudeEvent_with_instance() {
        setup_global_test_eop();
        let event =
            SLatitudeEvent::<6, 0>::new(0.0, "Latitude", EventDirection::Any, AngleFormat::Degrees)
                .with_instance(4);
        assert_eq!(event.name(), "Latitude 4");
    }

    #[test]
    fn test_SLatitudeEvent_with_tolerances() {
        setup_global_test_eop();
        let event =
            SLatitudeEvent::<6, 0>::new(0.0, "Latitude", EventDirection::Any, AngleFormat::Degrees)
                .with_tolerances(1e-6, 1e-9);
        assert_eq!(event.time_tolerance(), 1e-6);
        assert_eq!(event.value_tolerance(), 1e-9);
    }

    #[test]
    fn test_SLatitudeEvent_time_tolerance() {
        setup_global_test_eop();
        let event =
            SLatitudeEvent::<6, 0>::new(0.0, "Latitude", EventDirection::Any, AngleFormat::Degrees);
        assert!(event.time_tolerance() > 0.0);
    }

    #[test]
    fn test_SLatitudeEvent_value_tolerance() {
        setup_global_test_eop();
        let event =
            SLatitudeEvent::<6, 0>::new(0.0, "Latitude", EventDirection::Any, AngleFormat::Degrees);
        assert!(event.value_tolerance() > 0.0);
    }

    #[test]
    fn test_SLatitudeEvent_with_step_reduction_factor() {
        setup_global_test_eop();
        let event =
            SLatitudeEvent::<6, 0>::new(0.0, "Latitude", EventDirection::Any, AngleFormat::Degrees)
                .with_step_reduction_factor(0.3);
        assert_eq!(event.step_reduction_factor(), 0.3);
    }

    #[test]
    fn test_SLatitudeEvent_step_reduction_factor() {
        setup_global_test_eop();
        let event =
            SLatitudeEvent::<6, 0>::new(0.0, "Latitude", EventDirection::Any, AngleFormat::Degrees);
        assert!(event.step_reduction_factor() > 0.0);
        assert!(event.step_reduction_factor() < 1.0);
    }

    #[test]
    fn test_SLatitudeEvent_with_callback() {
        setup_global_test_eop();
        let callback: SEventCallback<6, 0> =
            Box::new(move |_t, _state, _params| (None, None, EventAction::Continue));
        let event =
            SLatitudeEvent::<6, 0>::new(0.0, "Latitude", EventDirection::Any, AngleFormat::Degrees)
                .with_callback(callback);
        assert!(event.callback().is_some());
    }

    #[test]
    fn test_SLatitudeEvent_callback_none() {
        setup_global_test_eop();
        let event =
            SLatitudeEvent::<6, 0>::new(0.0, "Latitude", EventDirection::Any, AngleFormat::Degrees);
        assert!(event.callback().is_none());
    }

    #[test]
    fn test_SLatitudeEvent_set_terminal() {
        setup_global_test_eop();
        let event =
            SLatitudeEvent::<6, 0>::new(0.0, "Latitude", EventDirection::Any, AngleFormat::Degrees)
                .set_terminal();
        assert_eq!(event.action(), EventAction::Stop);
    }

    #[test]
    fn test_SLatitudeEvent_action_continue() {
        setup_global_test_eop();
        let event =
            SLatitudeEvent::<6, 0>::new(0.0, "Latitude", EventDirection::Any, AngleFormat::Degrees);
        assert_eq!(event.action(), EventAction::Continue);
    }

    #[test]
    fn test_SLatitudeEvent_action_stop() {
        setup_global_test_eop();
        let event =
            SLatitudeEvent::<6, 0>::new(0.0, "Latitude", EventDirection::Any, AngleFormat::Degrees)
                .set_terminal();
        assert_eq!(event.action(), EventAction::Stop);
    }

    #[test]
    fn test_SLatitudeEvent_direction_increasing() {
        setup_global_test_eop();
        let event = SLatitudeEvent::<6, 0>::new(
            0.0,
            "Latitude",
            EventDirection::Increasing,
            AngleFormat::Degrees,
        );
        assert_eq!(event.direction(), EventDirection::Increasing);
    }

    #[test]
    fn test_SLatitudeEvent_direction_decreasing() {
        setup_global_test_eop();
        let event = SLatitudeEvent::<6, 0>::new(
            0.0,
            "Latitude",
            EventDirection::Decreasing,
            AngleFormat::Degrees,
        );
        assert_eq!(event.direction(), EventDirection::Decreasing);
    }

    #[test]
    fn test_SLatitudeEvent_direction_any() {
        setup_global_test_eop();
        let event =
            SLatitudeEvent::<6, 0>::new(0.0, "Latitude", EventDirection::Any, AngleFormat::Degrees);
        assert_eq!(event.direction(), EventDirection::Any);
    }

    #[test]
    fn test_SLatitudeEvent_angle_format_radians() {
        setup_global_test_eop();
        let event = SLatitudeEvent::<6, 0>::new(
            std::f64::consts::FRAC_PI_6,
            "Latitude",
            EventDirection::Any,
            AngleFormat::Radians,
        );
        assert_eq!(event.target_value(), std::f64::consts::FRAC_PI_6);
    }

    #[test]
    fn test_SLatitudeEvent_builder_chaining() {
        setup_global_test_eop();
        let callback: SEventCallback<6, 0> =
            Box::new(|_t, _state, _params| (None, None, EventAction::Stop));
        let event = SLatitudeEvent::<6, 0>::new(
            40.7,
            "NYC Latitude",
            EventDirection::Increasing,
            AngleFormat::Degrees,
        )
        .with_instance(1)
        .with_tolerances(1e-5, 1e-8)
        .with_step_reduction_factor(0.1)
        .with_callback(callback)
        .set_terminal();

        assert_eq!(event.name(), "NYC Latitude 1");
        assert_eq!(event.time_tolerance(), 1e-5);
        assert_eq!(event.step_reduction_factor(), 0.1);
        assert!(event.callback().is_some());
        assert_eq!(event.action(), EventAction::Stop);
    }

    // DLatitudeEvent comprehensive tests

    #[test]
    #[serial]
    fn test_DLatitudeEvent_evaluate() {
        setup_global_test_eop();
        let state = DVector::from_vec(vec![6800e3, 0.0, 500e3, 0.0, 7000.0, 1000.0]);
        let event = DLatitudeEvent::new(0.0, "Latitude", EventDirection::Any, AngleFormat::Degrees);
        let epoch = Epoch::from_datetime(2024, 1, 1, 0, 0, 0.0, 0.0, TimeSystem::UTC);
        let result = event.evaluate(epoch, &state, None);
        assert!(result.is_finite());
    }

    #[test]
    fn test_DLatitudeEvent_target_value() {
        setup_global_test_eop();
        let event =
            DLatitudeEvent::new(-30.0, "Latitude", EventDirection::Any, AngleFormat::Degrees);
        assert_eq!(event.target_value(), (-30.0_f64).to_radians());
    }

    #[test]
    fn test_DLatitudeEvent_with_instance() {
        setup_global_test_eop();
        let event = DLatitudeEvent::new(0.0, "Latitude", EventDirection::Any, AngleFormat::Degrees)
            .with_instance(7);
        assert_eq!(event.name(), "Latitude 7");
    }

    #[test]
    fn test_DLatitudeEvent_with_tolerances() {
        setup_global_test_eop();
        let event = DLatitudeEvent::new(0.0, "Latitude", EventDirection::Any, AngleFormat::Degrees)
            .with_tolerances(1e-6, 1e-9);
        assert_eq!(event.time_tolerance(), 1e-6);
        assert_eq!(event.value_tolerance(), 1e-9);
    }

    #[test]
    fn test_DLatitudeEvent_time_tolerance() {
        setup_global_test_eop();
        let event = DLatitudeEvent::new(0.0, "Latitude", EventDirection::Any, AngleFormat::Degrees);
        assert!(event.time_tolerance() > 0.0);
    }

    #[test]
    fn test_DLatitudeEvent_value_tolerance() {
        setup_global_test_eop();
        let event = DLatitudeEvent::new(0.0, "Latitude", EventDirection::Any, AngleFormat::Degrees);
        assert!(event.value_tolerance() > 0.0);
    }

    #[test]
    fn test_DLatitudeEvent_with_step_reduction_factor() {
        setup_global_test_eop();
        let event = DLatitudeEvent::new(0.0, "Latitude", EventDirection::Any, AngleFormat::Degrees)
            .with_step_reduction_factor(0.35);
        assert_eq!(event.step_reduction_factor(), 0.35);
    }

    #[test]
    fn test_DLatitudeEvent_step_reduction_factor() {
        setup_global_test_eop();
        let event = DLatitudeEvent::new(0.0, "Latitude", EventDirection::Any, AngleFormat::Degrees);
        assert!(event.step_reduction_factor() > 0.0);
        assert!(event.step_reduction_factor() < 1.0);
    }

    #[test]
    fn test_DLatitudeEvent_with_callback() {
        setup_global_test_eop();
        let callback: DEventCallback =
            Box::new(move |_t, _state, _params| (None, None, EventAction::Continue));
        let event = DLatitudeEvent::new(0.0, "Latitude", EventDirection::Any, AngleFormat::Degrees)
            .with_callback(callback);
        assert!(event.callback().is_some());
    }

    #[test]
    fn test_DLatitudeEvent_callback_none() {
        setup_global_test_eop();
        let event = DLatitudeEvent::new(0.0, "Latitude", EventDirection::Any, AngleFormat::Degrees);
        assert!(event.callback().is_none());
    }

    #[test]
    fn test_DLatitudeEvent_set_terminal() {
        setup_global_test_eop();
        let event = DLatitudeEvent::new(0.0, "Latitude", EventDirection::Any, AngleFormat::Degrees)
            .set_terminal();
        assert_eq!(event.action(), EventAction::Stop);
    }

    #[test]
    fn test_DLatitudeEvent_action_continue() {
        setup_global_test_eop();
        let event = DLatitudeEvent::new(0.0, "Latitude", EventDirection::Any, AngleFormat::Degrees);
        assert_eq!(event.action(), EventAction::Continue);
    }

    #[test]
    fn test_DLatitudeEvent_action_stop() {
        setup_global_test_eop();
        let event = DLatitudeEvent::new(0.0, "Latitude", EventDirection::Any, AngleFormat::Degrees)
            .set_terminal();
        assert_eq!(event.action(), EventAction::Stop);
    }

    #[test]
    fn test_DLatitudeEvent_direction_increasing() {
        setup_global_test_eop();
        let event = DLatitudeEvent::new(
            0.0,
            "Latitude",
            EventDirection::Increasing,
            AngleFormat::Degrees,
        );
        assert_eq!(event.direction(), EventDirection::Increasing);
    }

    #[test]
    fn test_DLatitudeEvent_direction_decreasing() {
        setup_global_test_eop();
        let event = DLatitudeEvent::new(
            0.0,
            "Latitude",
            EventDirection::Decreasing,
            AngleFormat::Degrees,
        );
        assert_eq!(event.direction(), EventDirection::Decreasing);
    }

    #[test]
    fn test_DLatitudeEvent_direction_any() {
        setup_global_test_eop();
        let event = DLatitudeEvent::new(0.0, "Latitude", EventDirection::Any, AngleFormat::Degrees);
        assert_eq!(event.direction(), EventDirection::Any);
    }

    #[test]
    fn test_DLatitudeEvent_angle_format_radians() {
        setup_global_test_eop();
        let event = DLatitudeEvent::new(
            std::f64::consts::FRAC_PI_4,
            "Latitude",
            EventDirection::Any,
            AngleFormat::Radians,
        );
        assert_eq!(event.target_value(), std::f64::consts::FRAC_PI_4);
    }

    #[test]
    fn test_DLatitudeEvent_builder_chaining() {
        setup_global_test_eop();
        let callback: DEventCallback =
            Box::new(|_t, _state, _params| (None, None, EventAction::Stop));
        let event = DLatitudeEvent::new(
            35.6,
            "Tokyo Latitude",
            EventDirection::Decreasing,
            AngleFormat::Degrees,
        )
        .with_instance(3)
        .with_tolerances(1e-5, 1e-8)
        .with_step_reduction_factor(0.2)
        .with_callback(callback)
        .set_terminal();

        assert_eq!(event.name(), "Tokyo Latitude 3");
        assert_eq!(event.time_tolerance(), 1e-5);
        assert_eq!(event.step_reduction_factor(), 0.2);
        assert!(event.callback().is_some());
        assert_eq!(event.action(), EventAction::Stop);
    }

    // =========================================================================
    // Eclipse/Shadow Event Tests
    // =========================================================================

    // SUmbraEvent comprehensive tests

    #[test]
    fn test_SUmbraEvent_new_without_source() {
        let event = SUmbraEvent::<6, 0>::new("Enter Umbra", EdgeType::RisingEdge, None);
        assert_eq!(event.name(), "Enter Umbra");
    }

    #[test]
    fn test_SUmbraEvent_new_with_source() {
        let event = SUmbraEvent::<6, 0>::new(
            "Umbra DE440s",
            EdgeType::FallingEdge,
            Some(crate::propagators::EphemerisSource::DE440s),
        );
        assert_eq!(event.name(), "Umbra DE440s");
    }

    #[test]
    fn test_SUmbraEvent_with_ephemeris_sources() {
        use crate::propagators::EphemerisSource;

        let event1 = SUmbraEvent::<6, 0>::new(
            "Umbra LP",
            EdgeType::RisingEdge,
            Some(EphemerisSource::LowPrecision),
        );
        assert_eq!(event1.name(), "Umbra LP");

        let event2 = SUmbraEvent::<6, 0>::new(
            "Umbra DE440s",
            EdgeType::RisingEdge,
            Some(EphemerisSource::DE440s),
        );
        assert_eq!(event2.name(), "Umbra DE440s");

        let event3 = SUmbraEvent::<6, 0>::new(
            "Umbra DE440",
            EdgeType::RisingEdge,
            Some(EphemerisSource::DE440),
        );
        assert_eq!(event3.name(), "Umbra DE440");
    }

    #[test]
    fn test_SUmbraEvent_with_instance() {
        let event = SUmbraEvent::<6, 0>::new("Umbra", EdgeType::RisingEdge, None).with_instance(2);
        assert_eq!(event.name(), "Umbra 2");
    }

    #[test]
    fn test_SUmbraEvent_with_tolerances() {
        let event = SUmbraEvent::<6, 0>::new("Umbra", EdgeType::RisingEdge, None)
            .with_tolerances(1e-6, 1e-9);
        assert_eq!(event.time_tolerance(), 1e-6);
        assert_eq!(event.value_tolerance(), 1e-9);
    }

    #[test]
    fn test_SUmbraEvent_time_tolerance() {
        let event = SUmbraEvent::<6, 0>::new("Umbra", EdgeType::RisingEdge, None);
        assert!(event.time_tolerance() > 0.0);
    }

    #[test]
    fn test_SUmbraEvent_value_tolerance() {
        let event = SUmbraEvent::<6, 0>::new("Umbra", EdgeType::RisingEdge, None);
        assert!(event.value_tolerance() > 0.0);
    }

    #[test]
    fn test_SUmbraEvent_with_step_reduction_factor() {
        let event = SUmbraEvent::<6, 0>::new("Umbra", EdgeType::RisingEdge, None)
            .with_step_reduction_factor(0.2);
        assert_eq!(event.step_reduction_factor(), 0.2);
    }

    #[test]
    fn test_SUmbraEvent_step_reduction_factor() {
        let event = SUmbraEvent::<6, 0>::new("Umbra", EdgeType::RisingEdge, None);
        assert!(event.step_reduction_factor() > 0.0);
        assert!(event.step_reduction_factor() < 1.0);
    }

    #[test]
    fn test_SUmbraEvent_with_callback() {
        let callback: SEventCallback<6, 0> =
            Box::new(move |_t, _state, _params| (None, None, EventAction::Continue));
        let event =
            SUmbraEvent::<6, 0>::new("Umbra", EdgeType::RisingEdge, None).with_callback(callback);
        assert!(event.callback().is_some());
    }

    #[test]
    fn test_SUmbraEvent_callback_none() {
        let event = SUmbraEvent::<6, 0>::new("Umbra", EdgeType::RisingEdge, None);
        assert!(event.callback().is_none());
    }

    #[test]
    fn test_SUmbraEvent_set_terminal() {
        let event = SUmbraEvent::<6, 0>::new("Umbra", EdgeType::RisingEdge, None).set_terminal();
        assert_eq!(event.action(), EventAction::Stop);
    }

    #[test]
    fn test_SUmbraEvent_action_continue() {
        let event = SUmbraEvent::<6, 0>::new("Umbra", EdgeType::RisingEdge, None);
        assert_eq!(event.action(), EventAction::Continue);
    }

    #[test]
    fn test_SUmbraEvent_action_stop() {
        let event = SUmbraEvent::<6, 0>::new("Umbra", EdgeType::RisingEdge, None).set_terminal();
        assert_eq!(event.action(), EventAction::Stop);
    }

    #[test]
    fn test_SUmbraEvent_builder_chaining() {
        let callback: SEventCallback<6, 0> =
            Box::new(|_t, _state, _params| (None, None, EventAction::Stop));
        let event = SUmbraEvent::<6, 0>::new("Umbra", EdgeType::RisingEdge, None)
            .with_instance(1)
            .with_tolerances(1e-5, 1e-8)
            .with_step_reduction_factor(0.1)
            .with_callback(callback)
            .set_terminal();

        assert_eq!(event.name(), "Umbra 1");
        assert_eq!(event.time_tolerance(), 1e-5);
        assert_eq!(event.step_reduction_factor(), 0.1);
        assert!(event.callback().is_some());
        assert_eq!(event.action(), EventAction::Stop);
    }

    // DUmbraEvent comprehensive tests

    #[test]
    fn test_DUmbraEvent_new() {
        let event = DUmbraEvent::new("Umbra", EdgeType::AnyEdge, None);
        assert_eq!(event.name(), "Umbra");
    }

    #[test]
    fn test_DUmbraEvent_with_ephemeris_sources() {
        use crate::propagators::EphemerisSource;

        let event1 = DUmbraEvent::new(
            "Umbra LP",
            EdgeType::RisingEdge,
            Some(EphemerisSource::LowPrecision),
        );
        assert_eq!(event1.name(), "Umbra LP");

        let event2 = DUmbraEvent::new(
            "Umbra DE440s",
            EdgeType::RisingEdge,
            Some(EphemerisSource::DE440s),
        );
        assert_eq!(event2.name(), "Umbra DE440s");

        let event3 = DUmbraEvent::new(
            "Umbra DE440",
            EdgeType::RisingEdge,
            Some(EphemerisSource::DE440),
        );
        assert_eq!(event3.name(), "Umbra DE440");
    }

    #[test]
    fn test_DUmbraEvent_with_instance() {
        let event = DUmbraEvent::new("Umbra", EdgeType::RisingEdge, None).with_instance(3);
        assert_eq!(event.name(), "Umbra 3");
    }

    #[test]
    fn test_DUmbraEvent_with_tolerances() {
        let event =
            DUmbraEvent::new("Umbra", EdgeType::RisingEdge, None).with_tolerances(1e-6, 1e-9);
        assert_eq!(event.time_tolerance(), 1e-6);
        assert_eq!(event.value_tolerance(), 1e-9);
    }

    #[test]
    fn test_DUmbraEvent_time_tolerance() {
        let event = DUmbraEvent::new("Umbra", EdgeType::RisingEdge, None);
        assert!(event.time_tolerance() > 0.0);
    }

    #[test]
    fn test_DUmbraEvent_value_tolerance() {
        let event = DUmbraEvent::new("Umbra", EdgeType::RisingEdge, None);
        assert!(event.value_tolerance() > 0.0);
    }

    #[test]
    fn test_DUmbraEvent_with_step_reduction_factor() {
        let event =
            DUmbraEvent::new("Umbra", EdgeType::RisingEdge, None).with_step_reduction_factor(0.25);
        assert_eq!(event.step_reduction_factor(), 0.25);
    }

    #[test]
    fn test_DUmbraEvent_step_reduction_factor() {
        let event = DUmbraEvent::new("Umbra", EdgeType::RisingEdge, None);
        assert!(event.step_reduction_factor() > 0.0);
        assert!(event.step_reduction_factor() < 1.0);
    }

    #[test]
    fn test_DUmbraEvent_with_callback() {
        let callback: DEventCallback =
            Box::new(move |_t, _state, _params| (None, None, EventAction::Continue));
        let event = DUmbraEvent::new("Umbra", EdgeType::RisingEdge, None).with_callback(callback);
        assert!(event.callback().is_some());
    }

    #[test]
    fn test_DUmbraEvent_callback_none() {
        let event = DUmbraEvent::new("Umbra", EdgeType::RisingEdge, None);
        assert!(event.callback().is_none());
    }

    #[test]
    fn test_DUmbraEvent_set_terminal() {
        let event = DUmbraEvent::new("Umbra", EdgeType::RisingEdge, None).set_terminal();
        assert_eq!(event.action(), EventAction::Stop);
    }

    #[test]
    fn test_DUmbraEvent_action_continue() {
        let event = DUmbraEvent::new("Umbra", EdgeType::RisingEdge, None);
        assert_eq!(event.action(), EventAction::Continue);
    }

    #[test]
    fn test_DUmbraEvent_action_stop() {
        let event = DUmbraEvent::new("Umbra", EdgeType::RisingEdge, None).set_terminal();
        assert_eq!(event.action(), EventAction::Stop);
    }

    #[test]
    fn test_DUmbraEvent_builder_chaining() {
        let callback: DEventCallback =
            Box::new(|_t, _state, _params| (None, None, EventAction::Stop));
        let event = DUmbraEvent::new("Umbra", EdgeType::FallingEdge, None)
            .with_instance(2)
            .with_tolerances(1e-5, 1e-8)
            .with_step_reduction_factor(0.15)
            .with_callback(callback)
            .set_terminal();

        assert_eq!(event.name(), "Umbra 2");
        assert_eq!(event.time_tolerance(), 1e-5);
        assert_eq!(event.step_reduction_factor(), 0.15);
        assert!(event.callback().is_some());
        assert_eq!(event.action(), EventAction::Stop);
    }

    // SPenumbraEvent comprehensive tests

    #[test]
    fn test_SPenumbraEvent_new() {
        let event = SPenumbraEvent::<6, 0>::new("Penumbra", EdgeType::RisingEdge, None);
        assert_eq!(event.name(), "Penumbra");
    }

    #[test]
    fn test_SPenumbraEvent_with_ephemeris_sources() {
        use crate::propagators::EphemerisSource;

        let event1 = SPenumbraEvent::<6, 0>::new(
            "Penumbra LP",
            EdgeType::RisingEdge,
            Some(EphemerisSource::LowPrecision),
        );
        assert_eq!(event1.name(), "Penumbra LP");

        let event2 = SPenumbraEvent::<6, 0>::new(
            "Penumbra DE440s",
            EdgeType::RisingEdge,
            Some(EphemerisSource::DE440s),
        );
        assert_eq!(event2.name(), "Penumbra DE440s");

        let event3 = SPenumbraEvent::<6, 0>::new(
            "Penumbra DE440",
            EdgeType::RisingEdge,
            Some(EphemerisSource::DE440),
        );
        assert_eq!(event3.name(), "Penumbra DE440");
    }

    #[test]
    fn test_SPenumbraEvent_with_instance() {
        let event =
            SPenumbraEvent::<6, 0>::new("Penumbra", EdgeType::RisingEdge, None).with_instance(1);
        assert_eq!(event.name(), "Penumbra 1");
    }

    #[test]
    fn test_SPenumbraEvent_with_tolerances() {
        let event = SPenumbraEvent::<6, 0>::new("Penumbra", EdgeType::RisingEdge, None)
            .with_tolerances(1e-6, 1e-9);
        assert_eq!(event.time_tolerance(), 1e-6);
        assert_eq!(event.value_tolerance(), 1e-9);
    }

    #[test]
    fn test_SPenumbraEvent_time_tolerance() {
        let event = SPenumbraEvent::<6, 0>::new("Penumbra", EdgeType::RisingEdge, None);
        assert!(event.time_tolerance() > 0.0);
    }

    #[test]
    fn test_SPenumbraEvent_value_tolerance() {
        let event = SPenumbraEvent::<6, 0>::new("Penumbra", EdgeType::RisingEdge, None);
        assert!(event.value_tolerance() > 0.0);
    }

    #[test]
    fn test_SPenumbraEvent_with_step_reduction_factor() {
        let event = SPenumbraEvent::<6, 0>::new("Penumbra", EdgeType::RisingEdge, None)
            .with_step_reduction_factor(0.3);
        assert_eq!(event.step_reduction_factor(), 0.3);
    }

    #[test]
    fn test_SPenumbraEvent_step_reduction_factor() {
        let event = SPenumbraEvent::<6, 0>::new("Penumbra", EdgeType::RisingEdge, None);
        assert!(event.step_reduction_factor() > 0.0);
        assert!(event.step_reduction_factor() < 1.0);
    }

    #[test]
    fn test_SPenumbraEvent_with_callback() {
        let callback: SEventCallback<6, 0> =
            Box::new(move |_t, _state, _params| (None, None, EventAction::Continue));
        let event = SPenumbraEvent::<6, 0>::new("Penumbra", EdgeType::RisingEdge, None)
            .with_callback(callback);
        assert!(event.callback().is_some());
    }

    #[test]
    fn test_SPenumbraEvent_callback_none() {
        let event = SPenumbraEvent::<6, 0>::new("Penumbra", EdgeType::RisingEdge, None);
        assert!(event.callback().is_none());
    }

    #[test]
    fn test_SPenumbraEvent_set_terminal() {
        let event =
            SPenumbraEvent::<6, 0>::new("Penumbra", EdgeType::RisingEdge, None).set_terminal();
        assert_eq!(event.action(), EventAction::Stop);
    }

    #[test]
    fn test_SPenumbraEvent_action_continue() {
        let event = SPenumbraEvent::<6, 0>::new("Penumbra", EdgeType::RisingEdge, None);
        assert_eq!(event.action(), EventAction::Continue);
    }

    #[test]
    fn test_SPenumbraEvent_action_stop() {
        let event =
            SPenumbraEvent::<6, 0>::new("Penumbra", EdgeType::RisingEdge, None).set_terminal();
        assert_eq!(event.action(), EventAction::Stop);
    }

    #[test]
    fn test_SPenumbraEvent_builder_chaining() {
        let callback: SEventCallback<6, 0> =
            Box::new(|_t, _state, _params| (None, None, EventAction::Stop));
        let event = SPenumbraEvent::<6, 0>::new("Penumbra", EdgeType::AnyEdge, None)
            .with_instance(3)
            .with_tolerances(1e-5, 1e-8)
            .with_step_reduction_factor(0.1)
            .with_callback(callback)
            .set_terminal();

        assert_eq!(event.name(), "Penumbra 3");
        assert_eq!(event.time_tolerance(), 1e-5);
        assert_eq!(event.step_reduction_factor(), 0.1);
        assert!(event.callback().is_some());
        assert_eq!(event.action(), EventAction::Stop);
    }

    // DPenumbraEvent comprehensive tests

    #[test]
    fn test_DPenumbraEvent_new() {
        let event = DPenumbraEvent::new("Penumbra", EdgeType::FallingEdge, None);
        assert_eq!(event.name(), "Penumbra");
    }

    #[test]
    fn test_DPenumbraEvent_with_ephemeris_sources() {
        use crate::propagators::EphemerisSource;

        let event1 = DPenumbraEvent::new(
            "Penumbra LP",
            EdgeType::RisingEdge,
            Some(EphemerisSource::LowPrecision),
        );
        assert_eq!(event1.name(), "Penumbra LP");

        let event2 = DPenumbraEvent::new(
            "Penumbra DE440s",
            EdgeType::RisingEdge,
            Some(EphemerisSource::DE440s),
        );
        assert_eq!(event2.name(), "Penumbra DE440s");

        let event3 = DPenumbraEvent::new(
            "Penumbra DE440",
            EdgeType::RisingEdge,
            Some(EphemerisSource::DE440),
        );
        assert_eq!(event3.name(), "Penumbra DE440");
    }

    #[test]
    fn test_DPenumbraEvent_with_instance() {
        let event = DPenumbraEvent::new("Penumbra", EdgeType::RisingEdge, None).with_instance(4);
        assert_eq!(event.name(), "Penumbra 4");
    }

    #[test]
    fn test_DPenumbraEvent_with_tolerances() {
        let event =
            DPenumbraEvent::new("Penumbra", EdgeType::RisingEdge, None).with_tolerances(1e-6, 1e-9);
        assert_eq!(event.time_tolerance(), 1e-6);
        assert_eq!(event.value_tolerance(), 1e-9);
    }

    #[test]
    fn test_DPenumbraEvent_time_tolerance() {
        let event = DPenumbraEvent::new("Penumbra", EdgeType::RisingEdge, None);
        assert!(event.time_tolerance() > 0.0);
    }

    #[test]
    fn test_DPenumbraEvent_value_tolerance() {
        let event = DPenumbraEvent::new("Penumbra", EdgeType::RisingEdge, None);
        assert!(event.value_tolerance() > 0.0);
    }

    #[test]
    fn test_DPenumbraEvent_with_step_reduction_factor() {
        let event = DPenumbraEvent::new("Penumbra", EdgeType::RisingEdge, None)
            .with_step_reduction_factor(0.35);
        assert_eq!(event.step_reduction_factor(), 0.35);
    }

    #[test]
    fn test_DPenumbraEvent_step_reduction_factor() {
        let event = DPenumbraEvent::new("Penumbra", EdgeType::RisingEdge, None);
        assert!(event.step_reduction_factor() > 0.0);
        assert!(event.step_reduction_factor() < 1.0);
    }

    #[test]
    fn test_DPenumbraEvent_with_callback() {
        let callback: DEventCallback =
            Box::new(move |_t, _state, _params| (None, None, EventAction::Continue));
        let event =
            DPenumbraEvent::new("Penumbra", EdgeType::RisingEdge, None).with_callback(callback);
        assert!(event.callback().is_some());
    }

    #[test]
    fn test_DPenumbraEvent_callback_none() {
        let event = DPenumbraEvent::new("Penumbra", EdgeType::RisingEdge, None);
        assert!(event.callback().is_none());
    }

    #[test]
    fn test_DPenumbraEvent_set_terminal() {
        let event = DPenumbraEvent::new("Penumbra", EdgeType::RisingEdge, None).set_terminal();
        assert_eq!(event.action(), EventAction::Stop);
    }

    #[test]
    fn test_DPenumbraEvent_action_continue() {
        let event = DPenumbraEvent::new("Penumbra", EdgeType::RisingEdge, None);
        assert_eq!(event.action(), EventAction::Continue);
    }

    #[test]
    fn test_DPenumbraEvent_action_stop() {
        let event = DPenumbraEvent::new("Penumbra", EdgeType::RisingEdge, None).set_terminal();
        assert_eq!(event.action(), EventAction::Stop);
    }

    #[test]
    fn test_DPenumbraEvent_builder_chaining() {
        let callback: DEventCallback =
            Box::new(|_t, _state, _params| (None, None, EventAction::Stop));
        let event = DPenumbraEvent::new("Penumbra", EdgeType::FallingEdge, None)
            .with_instance(5)
            .with_tolerances(1e-5, 1e-8)
            .with_step_reduction_factor(0.2)
            .with_callback(callback)
            .set_terminal();

        assert_eq!(event.name(), "Penumbra 5");
        assert_eq!(event.time_tolerance(), 1e-5);
        assert_eq!(event.step_reduction_factor(), 0.2);
        assert!(event.callback().is_some());
        assert_eq!(event.action(), EventAction::Stop);
    }

    // SEclipseEvent comprehensive tests

    #[test]
    fn test_SEclipseEvent_new() {
        let event = SEclipseEvent::<6, 0>::new("Eclipse", EdgeType::AnyEdge, None);
        assert_eq!(event.name(), "Eclipse");
    }

    #[test]
    fn test_SEclipseEvent_with_ephemeris_sources() {
        use crate::propagators::EphemerisSource;

        let event1 = SEclipseEvent::<6, 0>::new(
            "Eclipse LP",
            EdgeType::RisingEdge,
            Some(EphemerisSource::LowPrecision),
        );
        assert_eq!(event1.name(), "Eclipse LP");

        let event2 = SEclipseEvent::<6, 0>::new(
            "Eclipse DE440s",
            EdgeType::RisingEdge,
            Some(EphemerisSource::DE440s),
        );
        assert_eq!(event2.name(), "Eclipse DE440s");

        let event3 = SEclipseEvent::<6, 0>::new(
            "Eclipse DE440",
            EdgeType::RisingEdge,
            Some(EphemerisSource::DE440),
        );
        assert_eq!(event3.name(), "Eclipse DE440");
    }

    #[test]
    fn test_SEclipseEvent_with_instance() {
        let event =
            SEclipseEvent::<6, 0>::new("Eclipse", EdgeType::RisingEdge, None).with_instance(1);
        assert_eq!(event.name(), "Eclipse 1");
    }

    #[test]
    fn test_SEclipseEvent_with_tolerances() {
        let event = SEclipseEvent::<6, 0>::new("Eclipse", EdgeType::RisingEdge, None)
            .with_tolerances(1e-6, 1e-9);
        assert_eq!(event.time_tolerance(), 1e-6);
        assert_eq!(event.value_tolerance(), 1e-9);
    }

    #[test]
    fn test_SEclipseEvent_time_tolerance() {
        let event = SEclipseEvent::<6, 0>::new("Eclipse", EdgeType::RisingEdge, None);
        assert!(event.time_tolerance() > 0.0);
    }

    #[test]
    fn test_SEclipseEvent_value_tolerance() {
        let event = SEclipseEvent::<6, 0>::new("Eclipse", EdgeType::RisingEdge, None);
        assert!(event.value_tolerance() > 0.0);
    }

    #[test]
    fn test_SEclipseEvent_with_step_reduction_factor() {
        let event = SEclipseEvent::<6, 0>::new("Eclipse", EdgeType::RisingEdge, None)
            .with_step_reduction_factor(0.4);
        assert_eq!(event.step_reduction_factor(), 0.4);
    }

    #[test]
    fn test_SEclipseEvent_step_reduction_factor() {
        let event = SEclipseEvent::<6, 0>::new("Eclipse", EdgeType::RisingEdge, None);
        assert!(event.step_reduction_factor() > 0.0);
        assert!(event.step_reduction_factor() < 1.0);
    }

    #[test]
    fn test_SEclipseEvent_with_callback() {
        let callback: SEventCallback<6, 0> =
            Box::new(move |_t, _state, _params| (None, None, EventAction::Continue));
        let event = SEclipseEvent::<6, 0>::new("Eclipse", EdgeType::RisingEdge, None)
            .with_callback(callback);
        assert!(event.callback().is_some());
    }

    #[test]
    fn test_SEclipseEvent_callback_none() {
        let event = SEclipseEvent::<6, 0>::new("Eclipse", EdgeType::RisingEdge, None);
        assert!(event.callback().is_none());
    }

    #[test]
    fn test_SEclipseEvent_set_terminal() {
        let event =
            SEclipseEvent::<6, 0>::new("Eclipse", EdgeType::RisingEdge, None).set_terminal();
        assert_eq!(event.action(), EventAction::Stop);
    }

    #[test]
    fn test_SEclipseEvent_action_continue() {
        let event = SEclipseEvent::<6, 0>::new("Eclipse", EdgeType::RisingEdge, None);
        assert_eq!(event.action(), EventAction::Continue);
    }

    #[test]
    fn test_SEclipseEvent_action_stop() {
        let event =
            SEclipseEvent::<6, 0>::new("Eclipse", EdgeType::RisingEdge, None).set_terminal();
        assert_eq!(event.action(), EventAction::Stop);
    }

    #[test]
    fn test_SEclipseEvent_builder_chaining() {
        let callback: SEventCallback<6, 0> =
            Box::new(|_t, _state, _params| (None, None, EventAction::Stop));
        let event = SEclipseEvent::<6, 0>::new("Eclipse", EdgeType::AnyEdge, None)
            .with_instance(2)
            .with_tolerances(1e-5, 1e-8)
            .with_step_reduction_factor(0.1)
            .with_callback(callback)
            .set_terminal();

        assert_eq!(event.name(), "Eclipse 2");
        assert_eq!(event.time_tolerance(), 1e-5);
        assert_eq!(event.step_reduction_factor(), 0.1);
        assert!(event.callback().is_some());
        assert_eq!(event.action(), EventAction::Stop);
    }

    // DEclipseEvent comprehensive tests

    #[test]
    fn test_DEclipseEvent_new() {
        let event = DEclipseEvent::new("Eclipse", EdgeType::FallingEdge, None);
        assert_eq!(event.name(), "Eclipse");
    }

    #[test]
    fn test_DEclipseEvent_with_ephemeris_sources() {
        use crate::propagators::EphemerisSource;

        let event1 = DEclipseEvent::new(
            "Eclipse LP",
            EdgeType::RisingEdge,
            Some(EphemerisSource::LowPrecision),
        );
        assert_eq!(event1.name(), "Eclipse LP");

        let event2 = DEclipseEvent::new(
            "Eclipse DE440s",
            EdgeType::RisingEdge,
            Some(EphemerisSource::DE440s),
        );
        assert_eq!(event2.name(), "Eclipse DE440s");

        let event3 = DEclipseEvent::new(
            "Eclipse DE440",
            EdgeType::RisingEdge,
            Some(EphemerisSource::DE440),
        );
        assert_eq!(event3.name(), "Eclipse DE440");
    }

    #[test]
    fn test_DEclipseEvent_with_instance() {
        let event = DEclipseEvent::new("Eclipse", EdgeType::RisingEdge, None).with_instance(6);
        assert_eq!(event.name(), "Eclipse 6");
    }

    #[test]
    fn test_DEclipseEvent_with_tolerances() {
        let event =
            DEclipseEvent::new("Eclipse", EdgeType::RisingEdge, None).with_tolerances(1e-6, 1e-9);
        assert_eq!(event.time_tolerance(), 1e-6);
        assert_eq!(event.value_tolerance(), 1e-9);
    }

    #[test]
    fn test_DEclipseEvent_time_tolerance() {
        let event = DEclipseEvent::new("Eclipse", EdgeType::RisingEdge, None);
        assert!(event.time_tolerance() > 0.0);
    }

    #[test]
    fn test_DEclipseEvent_value_tolerance() {
        let event = DEclipseEvent::new("Eclipse", EdgeType::RisingEdge, None);
        assert!(event.value_tolerance() > 0.0);
    }

    #[test]
    fn test_DEclipseEvent_with_step_reduction_factor() {
        let event = DEclipseEvent::new("Eclipse", EdgeType::RisingEdge, None)
            .with_step_reduction_factor(0.45);
        assert_eq!(event.step_reduction_factor(), 0.45);
    }

    #[test]
    fn test_DEclipseEvent_step_reduction_factor() {
        let event = DEclipseEvent::new("Eclipse", EdgeType::RisingEdge, None);
        assert!(event.step_reduction_factor() > 0.0);
        assert!(event.step_reduction_factor() < 1.0);
    }

    #[test]
    fn test_DEclipseEvent_with_callback() {
        let callback: DEventCallback =
            Box::new(move |_t, _state, _params| (None, None, EventAction::Continue));
        let event =
            DEclipseEvent::new("Eclipse", EdgeType::RisingEdge, None).with_callback(callback);
        assert!(event.callback().is_some());
    }

    #[test]
    fn test_DEclipseEvent_callback_none() {
        let event = DEclipseEvent::new("Eclipse", EdgeType::RisingEdge, None);
        assert!(event.callback().is_none());
    }

    #[test]
    fn test_DEclipseEvent_set_terminal() {
        let event = DEclipseEvent::new("Eclipse", EdgeType::RisingEdge, None).set_terminal();
        assert_eq!(event.action(), EventAction::Stop);
    }

    #[test]
    fn test_DEclipseEvent_action_continue() {
        let event = DEclipseEvent::new("Eclipse", EdgeType::RisingEdge, None);
        assert_eq!(event.action(), EventAction::Continue);
    }

    #[test]
    fn test_DEclipseEvent_action_stop() {
        let event = DEclipseEvent::new("Eclipse", EdgeType::RisingEdge, None).set_terminal();
        assert_eq!(event.action(), EventAction::Stop);
    }

    #[test]
    fn test_DEclipseEvent_builder_chaining() {
        let callback: DEventCallback =
            Box::new(|_t, _state, _params| (None, None, EventAction::Stop));
        let event = DEclipseEvent::new("Eclipse", EdgeType::FallingEdge, None)
            .with_instance(7)
            .with_tolerances(1e-5, 1e-8)
            .with_step_reduction_factor(0.25)
            .with_callback(callback)
            .set_terminal();

        assert_eq!(event.name(), "Eclipse 7");
        assert_eq!(event.time_tolerance(), 1e-5);
        assert_eq!(event.step_reduction_factor(), 0.25);
        assert!(event.callback().is_some());
        assert_eq!(event.action(), EventAction::Stop);
    }

    // SSunlitEvent comprehensive tests

    #[test]
    fn test_SSunlitEvent_new() {
        let event = SSunlitEvent::<6, 0>::new("Sunlit", EdgeType::RisingEdge, None);
        assert_eq!(event.name(), "Sunlit");
    }

    #[test]
    fn test_SSunlitEvent_with_ephemeris_sources() {
        use crate::propagators::EphemerisSource;

        let event1 = SSunlitEvent::<6, 0>::new(
            "Sunlit LP",
            EdgeType::RisingEdge,
            Some(EphemerisSource::LowPrecision),
        );
        assert_eq!(event1.name(), "Sunlit LP");

        let event2 = SSunlitEvent::<6, 0>::new(
            "Sunlit DE440s",
            EdgeType::RisingEdge,
            Some(EphemerisSource::DE440s),
        );
        assert_eq!(event2.name(), "Sunlit DE440s");

        let event3 = SSunlitEvent::<6, 0>::new(
            "Sunlit DE440",
            EdgeType::RisingEdge,
            Some(EphemerisSource::DE440),
        );
        assert_eq!(event3.name(), "Sunlit DE440");
    }

    #[test]
    fn test_SSunlitEvent_with_instance() {
        let event =
            SSunlitEvent::<6, 0>::new("Sunlit", EdgeType::RisingEdge, None).with_instance(1);
        assert_eq!(event.name(), "Sunlit 1");
    }

    #[test]
    fn test_SSunlitEvent_with_tolerances() {
        let event = SSunlitEvent::<6, 0>::new("Sunlit", EdgeType::RisingEdge, None)
            .with_tolerances(1e-6, 1e-9);
        assert_eq!(event.time_tolerance(), 1e-6);
        assert_eq!(event.value_tolerance(), 1e-9);
    }

    #[test]
    fn test_SSunlitEvent_time_tolerance() {
        let event = SSunlitEvent::<6, 0>::new("Sunlit", EdgeType::RisingEdge, None);
        assert!(event.time_tolerance() > 0.0);
    }

    #[test]
    fn test_SSunlitEvent_value_tolerance() {
        let event = SSunlitEvent::<6, 0>::new("Sunlit", EdgeType::RisingEdge, None);
        assert!(event.value_tolerance() > 0.0);
    }

    #[test]
    fn test_SSunlitEvent_with_step_reduction_factor() {
        let event = SSunlitEvent::<6, 0>::new("Sunlit", EdgeType::RisingEdge, None)
            .with_step_reduction_factor(0.5);
        assert_eq!(event.step_reduction_factor(), 0.5);
    }

    #[test]
    fn test_SSunlitEvent_step_reduction_factor() {
        let event = SSunlitEvent::<6, 0>::new("Sunlit", EdgeType::RisingEdge, None);
        assert!(event.step_reduction_factor() > 0.0);
        assert!(event.step_reduction_factor() < 1.0);
    }

    #[test]
    fn test_SSunlitEvent_with_callback() {
        let callback: SEventCallback<6, 0> =
            Box::new(move |_t, _state, _params| (None, None, EventAction::Continue));
        let event =
            SSunlitEvent::<6, 0>::new("Sunlit", EdgeType::RisingEdge, None).with_callback(callback);
        assert!(event.callback().is_some());
    }

    #[test]
    fn test_SSunlitEvent_callback_none() {
        let event = SSunlitEvent::<6, 0>::new("Sunlit", EdgeType::RisingEdge, None);
        assert!(event.callback().is_none());
    }

    #[test]
    fn test_SSunlitEvent_set_terminal() {
        let event = SSunlitEvent::<6, 0>::new("Sunlit", EdgeType::RisingEdge, None).set_terminal();
        assert_eq!(event.action(), EventAction::Stop);
    }

    #[test]
    fn test_SSunlitEvent_action_continue() {
        let event = SSunlitEvent::<6, 0>::new("Sunlit", EdgeType::RisingEdge, None);
        assert_eq!(event.action(), EventAction::Continue);
    }

    #[test]
    fn test_SSunlitEvent_action_stop() {
        let event = SSunlitEvent::<6, 0>::new("Sunlit", EdgeType::RisingEdge, None).set_terminal();
        assert_eq!(event.action(), EventAction::Stop);
    }

    #[test]
    fn test_SSunlitEvent_builder_chaining() {
        let callback: SEventCallback<6, 0> =
            Box::new(|_t, _state, _params| (None, None, EventAction::Stop));
        let event = SSunlitEvent::<6, 0>::new("Sunlit", EdgeType::AnyEdge, None)
            .with_instance(2)
            .with_tolerances(1e-5, 1e-8)
            .with_step_reduction_factor(0.1)
            .with_callback(callback)
            .set_terminal();

        assert_eq!(event.name(), "Sunlit 2");
        assert_eq!(event.time_tolerance(), 1e-5);
        assert_eq!(event.step_reduction_factor(), 0.1);
        assert!(event.callback().is_some());
        assert_eq!(event.action(), EventAction::Stop);
    }

    // DSunlitEvent comprehensive tests

    #[test]
    fn test_DSunlitEvent_new() {
        let event = DSunlitEvent::new("Sunlit", EdgeType::AnyEdge, None);
        assert_eq!(event.name(), "Sunlit");
    }

    #[test]
    fn test_DSunlitEvent_with_ephemeris_sources() {
        use crate::propagators::EphemerisSource;

        let event1 = DSunlitEvent::new(
            "Sunlit LP",
            EdgeType::RisingEdge,
            Some(EphemerisSource::LowPrecision),
        );
        assert_eq!(event1.name(), "Sunlit LP");

        let event2 = DSunlitEvent::new(
            "Sunlit DE440s",
            EdgeType::RisingEdge,
            Some(EphemerisSource::DE440s),
        );
        assert_eq!(event2.name(), "Sunlit DE440s");

        let event3 = DSunlitEvent::new(
            "Sunlit DE440",
            EdgeType::RisingEdge,
            Some(EphemerisSource::DE440),
        );
        assert_eq!(event3.name(), "Sunlit DE440");
    }

    #[test]
    fn test_DSunlitEvent_with_instance() {
        let event = DSunlitEvent::new("Sunlit", EdgeType::RisingEdge, None).with_instance(8);
        assert_eq!(event.name(), "Sunlit 8");
    }

    #[test]
    fn test_DSunlitEvent_with_tolerances() {
        let event =
            DSunlitEvent::new("Sunlit", EdgeType::RisingEdge, None).with_tolerances(1e-6, 1e-9);
        assert_eq!(event.time_tolerance(), 1e-6);
        assert_eq!(event.value_tolerance(), 1e-9);
    }

    #[test]
    fn test_DSunlitEvent_time_tolerance() {
        let event = DSunlitEvent::new("Sunlit", EdgeType::RisingEdge, None);
        assert!(event.time_tolerance() > 0.0);
    }

    #[test]
    fn test_DSunlitEvent_value_tolerance() {
        let event = DSunlitEvent::new("Sunlit", EdgeType::RisingEdge, None);
        assert!(event.value_tolerance() > 0.0);
    }

    #[test]
    fn test_DSunlitEvent_with_step_reduction_factor() {
        let event = DSunlitEvent::new("Sunlit", EdgeType::RisingEdge, None)
            .with_step_reduction_factor(0.55);
        assert_eq!(event.step_reduction_factor(), 0.55);
    }

    #[test]
    fn test_DSunlitEvent_step_reduction_factor() {
        let event = DSunlitEvent::new("Sunlit", EdgeType::RisingEdge, None);
        assert!(event.step_reduction_factor() > 0.0);
        assert!(event.step_reduction_factor() < 1.0);
    }

    #[test]
    fn test_DSunlitEvent_with_callback() {
        let callback: DEventCallback =
            Box::new(move |_t, _state, _params| (None, None, EventAction::Continue));
        let event = DSunlitEvent::new("Sunlit", EdgeType::RisingEdge, None).with_callback(callback);
        assert!(event.callback().is_some());
    }

    #[test]
    fn test_DSunlitEvent_callback_none() {
        let event = DSunlitEvent::new("Sunlit", EdgeType::RisingEdge, None);
        assert!(event.callback().is_none());
    }

    #[test]
    fn test_DSunlitEvent_set_terminal() {
        let event = DSunlitEvent::new("Sunlit", EdgeType::RisingEdge, None).set_terminal();
        assert_eq!(event.action(), EventAction::Stop);
    }

    #[test]
    fn test_DSunlitEvent_action_continue() {
        let event = DSunlitEvent::new("Sunlit", EdgeType::RisingEdge, None);
        assert_eq!(event.action(), EventAction::Continue);
    }

    #[test]
    fn test_DSunlitEvent_action_stop() {
        let event = DSunlitEvent::new("Sunlit", EdgeType::RisingEdge, None).set_terminal();
        assert_eq!(event.action(), EventAction::Stop);
    }

    #[test]
    fn test_DSunlitEvent_builder_chaining() {
        let callback: DEventCallback =
            Box::new(|_t, _state, _params| (None, None, EventAction::Stop));
        let event = DSunlitEvent::new("Sunlit", EdgeType::FallingEdge, None)
            .with_instance(9)
            .with_tolerances(1e-5, 1e-8)
            .with_step_reduction_factor(0.3)
            .with_callback(callback)
            .set_terminal();

        assert_eq!(event.name(), "Sunlit 9");
        assert_eq!(event.time_tolerance(), 1e-5);
        assert_eq!(event.step_reduction_factor(), 0.3);
        assert!(event.callback().is_some());
        assert_eq!(event.action(), EventAction::Stop);
    }

    // =========================================================================
    // SAOIEntryEvent Tests
    // =========================================================================

    #[test]
    fn test_SAOIEntryEvent_from_coordinates() {
        setup_global_test_eop();

        // Simple square polygon around (0, 0)
        let vertices = vec![
            (-10.0, -10.0),
            (10.0, -10.0),
            (10.0, 10.0),
            (-10.0, 10.0),
            (-10.0, -10.0), // closed
        ];

        let event =
            SAOIEntryEvent::<6, 0>::from_coordinates(&vertices, "AOI Entry", AngleFormat::Degrees);

        assert_eq!(event.name(), "AOI Entry");
        assert_eq!(event.action(), EventAction::Continue);
    }

    #[test]
    fn test_SAOIEntryEvent_from_polygon() {
        setup_global_test_eop();

        // Create polygon location - vertices are (lon, lat, alt) in degrees
        let polygon = PolygonLocation::new(vec![
            Vector3::new(-10.0, -10.0, 0.0),
            Vector3::new(10.0, -10.0, 0.0),
            Vector3::new(10.0, 10.0, 0.0),
            Vector3::new(-10.0, 10.0, 0.0),
            Vector3::new(-10.0, -10.0, 0.0),
        ])
        .unwrap();

        let event = SAOIEntryEvent::<6, 0>::from_polygon(&polygon, "AOI Entry");

        assert_eq!(event.name(), "AOI Entry");
    }

    #[test]
    #[serial]
    fn test_SAOIEntryEvent_evaluate_inside() {
        setup_global_test_eop();

        // At J2000 epoch, ECI x-axis position maps to approximately lon=79.54°, lat=0°
        // Create a polygon that covers this region (60-100° lon, ±30° lat)
        let vertices = vec![
            (60.0, -30.0),
            (100.0, -30.0),
            (100.0, 30.0),
            (60.0, 30.0),
            (60.0, -30.0),
        ];

        let event =
            SAOIEntryEvent::<6, 0>::from_coordinates(&vertices, "AOI Entry", AngleFormat::Degrees);

        let epoch = Epoch::from_jd(2451545.0, TimeSystem::UTC);

        // ECI position along x-axis at altitude 500 km
        // At this epoch, this maps to approximately lon=79.54°, lat=0°
        let r = R_EARTH + 500e3;
        let state_inside = Vector6::new(r, 0.0, 0.0, 0.0, 7.5e3, 0.0);
        let val = event.evaluate(epoch, &state_inside, None);
        // Inside polygon -> binary returns true -> should be positive
        assert!(
            val > 0.0,
            "Expected positive value for inside AOI, got {}",
            val
        );
    }

    #[test]
    #[serial]
    fn test_SAOIEntryEvent_evaluate_outside() {
        setup_global_test_eop();

        // Small polygon around (0, 0)
        let vertices = vec![
            (-5.0, -5.0),
            (5.0, -5.0),
            (5.0, 5.0),
            (-5.0, 5.0),
            (-5.0, -5.0),
        ];

        let event =
            SAOIEntryEvent::<6, 0>::from_coordinates(&vertices, "AOI Entry", AngleFormat::Degrees);

        let epoch = Epoch::from_jd(2451545.0, TimeSystem::UTC);

        // State at lon=90, lat=0 - directly over Indian Ocean
        // ECI position along y-axis at this epoch
        let r = R_EARTH + 500e3;
        let state_outside = Vector6::new(0.0, r, 0.0, 0.0, -7.5e3, 0.0);
        let val = event.evaluate(epoch, &state_outside, None);
        // Outside polygon -> binary returns false -> should be negative
        assert!(
            val < 0.0,
            "Expected negative value for outside AOI, got {}",
            val
        );
    }

    #[test]
    fn test_SAOIEntryEvent_with_instance() {
        setup_global_test_eop();

        let vertices = vec![
            (-10.0, -10.0),
            (10.0, -10.0),
            (10.0, 10.0),
            (-10.0, 10.0),
            (-10.0, -10.0),
        ];

        let event =
            SAOIEntryEvent::<6, 0>::from_coordinates(&vertices, "AOI Entry", AngleFormat::Degrees)
                .with_instance(3);

        assert_eq!(event.name(), "AOI Entry 3");
    }

    #[test]
    fn test_SAOIEntryEvent_set_terminal() {
        setup_global_test_eop();

        let vertices = vec![
            (-10.0, -10.0),
            (10.0, -10.0),
            (10.0, 10.0),
            (-10.0, 10.0),
            (-10.0, -10.0),
        ];

        let event =
            SAOIEntryEvent::<6, 0>::from_coordinates(&vertices, "AOI Entry", AngleFormat::Degrees)
                .set_terminal();

        assert_eq!(event.action(), EventAction::Stop);
    }

    #[test]
    fn test_SAOIEntryEvent_with_tolerances() {
        setup_global_test_eop();

        let vertices = vec![
            (-10.0, -10.0),
            (10.0, -10.0),
            (10.0, 10.0),
            (-10.0, 10.0),
            (-10.0, -10.0),
        ];

        let event =
            SAOIEntryEvent::<6, 0>::from_coordinates(&vertices, "AOI Entry", AngleFormat::Degrees)
                .with_tolerances(1e-4, 1e-6);

        assert_eq!(event.time_tolerance(), 1e-4);
    }

    #[test]
    fn test_SAOIEntryEvent_radians_input() {
        setup_global_test_eop();

        use std::f64::consts::PI;

        // Same polygon but in radians
        let vertices = vec![
            (-10.0 * PI / 180.0, -10.0 * PI / 180.0),
            (10.0 * PI / 180.0, -10.0 * PI / 180.0),
            (10.0 * PI / 180.0, 10.0 * PI / 180.0),
            (-10.0 * PI / 180.0, 10.0 * PI / 180.0),
            (-10.0 * PI / 180.0, -10.0 * PI / 180.0),
        ];

        let event = SAOIEntryEvent::<6, 0>::from_coordinates(
            &vertices,
            "AOI Entry Radians",
            AngleFormat::Radians,
        );

        assert_eq!(event.name(), "AOI Entry Radians");
    }

    // =========================================================================
    // SAOIExitEvent Tests
    // =========================================================================

    #[test]
    fn test_SAOIExitEvent_from_coordinates() {
        setup_global_test_eop();

        let vertices = vec![
            (-10.0, -10.0),
            (10.0, -10.0),
            (10.0, 10.0),
            (-10.0, 10.0),
            (-10.0, -10.0),
        ];

        let event =
            SAOIExitEvent::<6, 0>::from_coordinates(&vertices, "AOI Exit", AngleFormat::Degrees);

        assert_eq!(event.name(), "AOI Exit");
        assert_eq!(event.action(), EventAction::Continue);
    }

    #[test]
    fn test_SAOIExitEvent_from_polygon() {
        setup_global_test_eop();

        let polygon = PolygonLocation::new(vec![
            Vector3::new(-10.0, -10.0, 0.0),
            Vector3::new(10.0, -10.0, 0.0),
            Vector3::new(10.0, 10.0, 0.0),
            Vector3::new(-10.0, 10.0, 0.0),
            Vector3::new(-10.0, -10.0, 0.0),
        ])
        .unwrap();

        let event = SAOIExitEvent::<6, 0>::from_polygon(&polygon, "AOI Exit");

        assert_eq!(event.name(), "AOI Exit");
    }

    #[test]
    #[serial]
    fn test_SAOIExitEvent_evaluate_inside() {
        setup_global_test_eop();

        // At J2000 epoch, ECI x-axis position maps to approximately lon=79.54°, lat=0°
        // Create a polygon that covers this region (60-100° lon, ±30° lat)
        let vertices = vec![
            (60.0, -30.0),
            (100.0, -30.0),
            (100.0, 30.0),
            (60.0, 30.0),
            (60.0, -30.0),
        ];

        let event =
            SAOIExitEvent::<6, 0>::from_coordinates(&vertices, "AOI Exit", AngleFormat::Degrees);

        let epoch = Epoch::from_jd(2451545.0, TimeSystem::UTC);

        // State inside the AOI (equatorial position, maps to ~lon=79.54°)
        let r = R_EARTH + 500e3;
        let state_inside = Vector6::new(r, 0.0, 0.0, 0.0, 7.5e3, 0.0);
        let val = event.evaluate(epoch, &state_inside, None);
        // Inside polygon -> binary returns true -> should be positive
        assert!(
            val > 0.0,
            "Expected positive value for inside AOI, got {}",
            val
        );
    }

    #[test]
    #[serial]
    fn test_SAOIExitEvent_evaluate_outside() {
        setup_global_test_eop();

        let vertices = vec![
            (-5.0, -5.0),
            (5.0, -5.0),
            (5.0, 5.0),
            (-5.0, 5.0),
            (-5.0, -5.0),
        ];

        let event =
            SAOIExitEvent::<6, 0>::from_coordinates(&vertices, "AOI Exit", AngleFormat::Degrees);

        let epoch = Epoch::from_jd(2451545.0, TimeSystem::UTC);

        // State outside the AOI (at lon=90)
        let r = R_EARTH + 500e3;
        let state_outside = Vector6::new(0.0, r, 0.0, 0.0, -7.5e3, 0.0);
        let val = event.evaluate(epoch, &state_outside, None);
        // Outside polygon -> binary returns false -> should be negative
        assert!(
            val < 0.0,
            "Expected negative value for outside AOI, got {}",
            val
        );
    }

    #[test]
    fn test_SAOIExitEvent_with_instance() {
        setup_global_test_eop();

        let vertices = vec![
            (-10.0, -10.0),
            (10.0, -10.0),
            (10.0, 10.0),
            (-10.0, 10.0),
            (-10.0, -10.0),
        ];

        let event =
            SAOIExitEvent::<6, 0>::from_coordinates(&vertices, "AOI Exit", AngleFormat::Degrees)
                .with_instance(5);

        assert_eq!(event.name(), "AOI Exit 5");
    }

    #[test]
    fn test_SAOIExitEvent_set_terminal() {
        setup_global_test_eop();

        let vertices = vec![
            (-10.0, -10.0),
            (10.0, -10.0),
            (10.0, 10.0),
            (-10.0, 10.0),
            (-10.0, -10.0),
        ];

        let event =
            SAOIExitEvent::<6, 0>::from_coordinates(&vertices, "AOI Exit", AngleFormat::Degrees)
                .set_terminal();

        assert_eq!(event.action(), EventAction::Stop);
    }

    // =========================================================================
    // DAOIEntryEvent Tests
    // =========================================================================

    #[test]
    fn test_DAOIEntryEvent_from_coordinates() {
        setup_global_test_eop();

        let vertices = vec![
            (-10.0, -10.0),
            (10.0, -10.0),
            (10.0, 10.0),
            (-10.0, 10.0),
            (-10.0, -10.0),
        ];

        let event =
            DAOIEntryEvent::from_coordinates(&vertices, "Dynamic AOI Entry", AngleFormat::Degrees);

        assert_eq!(event.name(), "Dynamic AOI Entry");
        assert_eq!(event.action(), EventAction::Continue);
    }

    #[test]
    fn test_DAOIEntryEvent_from_polygon() {
        setup_global_test_eop();

        let polygon = PolygonLocation::new(vec![
            Vector3::new(-10.0, -10.0, 0.0),
            Vector3::new(10.0, -10.0, 0.0),
            Vector3::new(10.0, 10.0, 0.0),
            Vector3::new(-10.0, 10.0, 0.0),
            Vector3::new(-10.0, -10.0, 0.0),
        ])
        .unwrap();

        let event = DAOIEntryEvent::from_polygon(&polygon, "Dynamic AOI Entry");

        assert_eq!(event.name(), "Dynamic AOI Entry");
    }

    #[test]
    #[serial]
    fn test_DAOIEntryEvent_evaluate() {
        setup_global_test_eop();

        // At J2000 epoch, ECI x-axis position maps to approximately lon=79.54°, lat=0°
        // Create a polygon that covers this region (60-100° lon, ±30° lat)
        let vertices = vec![
            (60.0, -30.0),
            (100.0, -30.0),
            (100.0, 30.0),
            (60.0, 30.0),
            (60.0, -30.0),
        ];

        let event =
            DAOIEntryEvent::from_coordinates(&vertices, "Dynamic AOI Entry", AngleFormat::Degrees);

        let epoch = Epoch::from_jd(2451545.0, TimeSystem::UTC);

        // State inside the AOI (equatorial position, maps to ~lon=79.54°)
        let r = R_EARTH + 500e3;
        let state_inside = DVector::from_vec(vec![r, 0.0, 0.0, 0.0, 7.5e3, 0.0]);
        let val = event.evaluate(epoch, &state_inside, None);
        assert!(
            val > 0.0,
            "Expected positive value for inside AOI, got {}",
            val
        );

        // State outside the AOI (polar position at z-axis, maps to ~lat=90°)
        let state_outside = DVector::from_vec(vec![0.0, 0.0, r, 0.0, 7.5e3, 0.0]);
        let val = event.evaluate(epoch, &state_outside, None);
        assert!(
            val < 0.0,
            "Expected negative value for outside AOI, got {}",
            val
        );
    }

    #[test]
    fn test_DAOIEntryEvent_with_instance() {
        setup_global_test_eop();

        let vertices = vec![
            (-10.0, -10.0),
            (10.0, -10.0),
            (10.0, 10.0),
            (-10.0, 10.0),
            (-10.0, -10.0),
        ];

        let event =
            DAOIEntryEvent::from_coordinates(&vertices, "Dynamic AOI Entry", AngleFormat::Degrees)
                .with_instance(2);

        assert_eq!(event.name(), "Dynamic AOI Entry 2");
    }

    // =========================================================================
    // DAOIExitEvent Tests
    // =========================================================================

    #[test]
    fn test_DAOIExitEvent_from_coordinates() {
        setup_global_test_eop();

        let vertices = vec![
            (-10.0, -10.0),
            (10.0, -10.0),
            (10.0, 10.0),
            (-10.0, 10.0),
            (-10.0, -10.0),
        ];

        let event =
            DAOIExitEvent::from_coordinates(&vertices, "Dynamic AOI Exit", AngleFormat::Degrees);

        assert_eq!(event.name(), "Dynamic AOI Exit");
        assert_eq!(event.action(), EventAction::Continue);
    }

    #[test]
    fn test_DAOIExitEvent_from_polygon() {
        setup_global_test_eop();

        let polygon = PolygonLocation::new(vec![
            Vector3::new(-10.0, -10.0, 0.0),
            Vector3::new(10.0, -10.0, 0.0),
            Vector3::new(10.0, 10.0, 0.0),
            Vector3::new(-10.0, 10.0, 0.0),
            Vector3::new(-10.0, -10.0, 0.0),
        ])
        .unwrap();

        let event = DAOIExitEvent::from_polygon(&polygon, "Dynamic AOI Exit");

        assert_eq!(event.name(), "Dynamic AOI Exit");
    }

    #[test]
    #[serial]
    fn test_DAOIExitEvent_evaluate() {
        setup_global_test_eop();

        // At J2000 epoch, ECI x-axis position maps to approximately lon=79.54°, lat=0°
        // Create a polygon that covers this region (60-100° lon, ±30° lat)
        let vertices = vec![
            (60.0, -30.0),
            (100.0, -30.0),
            (100.0, 30.0),
            (60.0, 30.0),
            (60.0, -30.0),
        ];

        let event =
            DAOIExitEvent::from_coordinates(&vertices, "Dynamic AOI Exit", AngleFormat::Degrees);

        let epoch = Epoch::from_jd(2451545.0, TimeSystem::UTC);

        // State inside the AOI (equatorial position, maps to ~lon=79.54°)
        let r = R_EARTH + 500e3;
        let state_inside = DVector::from_vec(vec![r, 0.0, 0.0, 0.0, 7.5e3, 0.0]);
        let val = event.evaluate(epoch, &state_inside, None);
        assert!(
            val > 0.0,
            "Expected positive value for inside AOI, got {}",
            val
        );

        // State outside the AOI (polar position at z-axis, maps to ~lat=90°)
        let state_outside = DVector::from_vec(vec![0.0, 0.0, r, 0.0, 7.5e3, 0.0]);
        let val = event.evaluate(epoch, &state_outside, None);
        assert!(
            val < 0.0,
            "Expected negative value for outside AOI, got {}",
            val
        );
    }

    #[test]
    fn test_DAOIExitEvent_with_instance() {
        setup_global_test_eop();

        let vertices = vec![
            (-10.0, -10.0),
            (10.0, -10.0),
            (10.0, 10.0),
            (-10.0, 10.0),
            (-10.0, -10.0),
        ];

        let event =
            DAOIExitEvent::from_coordinates(&vertices, "Dynamic AOI Exit", AngleFormat::Degrees)
                .with_instance(4);

        assert_eq!(event.name(), "Dynamic AOI Exit 4");
    }

    // =========================================================================
    // AOI Event Anti-Meridian Tests
    // =========================================================================

    #[test]
    fn test_SAOIEntryEvent_antimeridian() {
        setup_global_test_eop();

        // Polygon crossing anti-meridian (e.g., Pacific region)
        let vertices = vec![
            (170.0, -10.0),
            (-170.0, -10.0), // crosses anti-meridian
            (-170.0, 10.0),
            (170.0, 10.0),
            (170.0, -10.0),
        ];

        let event = SAOIEntryEvent::<6, 0>::from_coordinates(
            &vertices,
            "Pacific AOI",
            AngleFormat::Degrees,
        );

        assert_eq!(event.name(), "Pacific AOI");
        // Event should be created successfully even for anti-meridian crossing polygons
    }

    #[test]
    fn test_SAOIExitEvent_antimeridian() {
        setup_global_test_eop();

        // Polygon crossing anti-meridian
        let vertices = vec![
            (170.0, -10.0),
            (-170.0, -10.0),
            (-170.0, 10.0),
            (170.0, 10.0),
            (170.0, -10.0),
        ];

        let event = SAOIExitEvent::<6, 0>::from_coordinates(
            &vertices,
            "Pacific AOI Exit",
            AngleFormat::Degrees,
        );

        assert_eq!(event.name(), "Pacific AOI Exit");
    }
}