fn log_det_psd(m: &Array2<f64>) -> f64
Cholesky-decomposition log-determinant of a symmetric PSD matrix.
Returns 2 * sum(log diag(L)) where M = L L^T. Returns NaN if the matrix is not positive definite.