pub fn psi(reference: &[f64], current: &[f64], n_bins: usize) -> Result<f64>
Compute PSI between a reference and current 1D distribution.