pub fn cholesky_determinant(n: usize, l: &[f64]) -> f64
Determinant of an SPD matrix from its Cholesky factor: det(A) = Π L[i][i]².
det(A) = Π L[i][i]²