powered by
this function checks if the input is a symmetric positive definite matrix
is_sympd_matrix(x, n)
is the input
is the input matrix dimension (assuming a square matrix)