pub fn bitonic_argsort(data: &[f64]) -> Vec<usize>
Bitonic sort that returns the original indices (argsort variant).
Pads with (f64::INFINITY, usize::MAX) pairs and trims back.
(f64::INFINITY, usize::MAX)