pub fn correlation_matrix(data: &Array2<f64>) -> Result<Array2<f64>>
Compute the pairwise correlation matrix (from original tools/utils.rs)