fn chol_logdet(chol: &Cholesky<f64, Dyn>, n: usize) -> f64
log-determinant of an SPD matrix from its Cholesky factor (2 Σ log L_ii).
2 Σ log L_ii