is_hpd function

Check if a matrix is Hermitian positive definite