class:`cvxopt.base.matrix` class).
"""
+from copy import copy
from math import sqrt
from cvxopt import matrix
from cvxopt.lapack import gees, gesdd, syevr
domain_dim = symmat.size[0]
eigs = matrix(0, (domain_dim, 1), tc='d')
- syevr(symmat, eigs)
+
+ # Create a copy of ``symmat`` here because ``syevr`` clobbers it.
+ dummy = copy(symmat)
+ syevr(dummy, eigs)
return list(eigs)