pub fn drift_persistence(norms: &[f64], k: usize, window: usize) -> f64
Compute rolling drift persistence: fraction of last W windows with drift > 0