"""
if not A.is_symmetric():
+ # This function is more or less internal, so blow up if passed
+ # something unexpected.
raise ValueError('The matrix ``A`` must be symmetric.')
r = rank(A)
- n = A.nrows() # Columns would work, too, since ``A`` is symmetric.
+ n = ZZ(A.nrows()) # Columns would work, too, since ``A`` is symmetric.
if r == 0:
# Zero is in the doubly-nonnegative cone.