fn gaussian_mi(x: &[Vec<f64>], y: &[Vec<f64>]) -> f64
Gaussian mutual information: MI(X;Y) = 0.5 * ln(det(Cov_X) * det(Cov_Y) / det(Cov_XY))