pub fn cholesky(a: &Matrix) -> Option<Matrix>
Cholesky decomposition for symmetric positive-definite matrices. Returns lower triangular L such that A = L * L^T. Returns None if not SPD.