pub fn column_means(array: &Array2<f64>, config: &SimdConfig) -> Array1<f64>
SIMD-optimized column-wise mean calculation