pub fn phi_star(data: &[Vec<f64>], tau: usize) -> f64
Compute Phi* for a multi-channel time series. data: [n_channels][n_timesteps], row-major.