"""
domain_dim = real_matrix.size[0] # Assume ``real_matrix`` is square.
- W = matrix(0, (domain_dim, 1), tc='d')
- syev(real_matrix, W)
- return list(W)
+ eigs = matrix(0, (domain_dim, 1), tc='d')
+ syev(real_matrix, eigs)
+ return list(eigs)
def identity(domain_dim):