Function positive_semidefinite 

fn positive_semidefinite<const N: usize>(m: &[[f64; N]; N]) -> bool

Whether a symmetric matrix is positive semi-definite, by an LDLᵀ walk that allows a zero pivot only on a zero remainder.