fn cholesky(m: &Array2<f64>) -> Array2<f64>
Cholesky decomposition: returns L lower-triangular such that L L^T = M for symmetric PSD M. No checks on PSD-ness; if M is not PSD, the diagonal of L will contain a NaN.