is_hpd function

Check if a matrix is Hermitian positive definite

is_hpd(A, tol = 1e-15)