pub fn compute_spin_drift_with_decay( time_s: f64, stability: f64, twist_rate: f64, is_twist_right: bool, decay_factor: Option<f64>, ) -> f64
Calculate spin drift with optional spin decay modeling