fn percentile(sorted_scratch: &mut [f32], q: f64) -> f64
np.percentile(sorted, q) with linear interpolation between order statistics.
np.percentile(sorted, q)