[X,fX,eig_val,eig_vec]=NCeigOpt(f,S,d) computes minimizers and the minimum of the polynomial f on D_S.

f is an NCpoly representing a polynomial
S is a set of nc polynomails defining D_S
d is a starting degree for the hierarchy (even number)

X: from GNS - a matrix where each of its rows represents a square matrix
fX: f(X) where X is from GNS
eig_val: eigenvalues of fX
eig_vec: corresponding eigenvectors

NCeigOpt(f,S,d)