fn percentile_95(mag: &[f64], mask: &[u8]) -> f64
Compute 95th percentile of magnitude within mask (matching ROMEO.jl).