is_spd function

Check if a matrix is symmetric positive definite

is_spd(A, tol = 1e-05)