from math import sqrt
from cvxopt import matrix
-from cvxopt.lapack import syev
+from cvxopt.lapack import syevr
import options
domain_dim = symmat.size[0]
eigs = matrix(0, (domain_dim, 1), tc='d')
- syev(symmat, eigs)
+ syevr(symmat, eigs)
return list(eigs)