fn cholesky_solve(a: &[Vec<f64>], b: &[f64]) -> Vec<f64>
Solve Ax = b where A is symmetric positive definite via Cholesky.