from os.path import abspath
from site import addsitedir
addsitedir(abspath('../../'))
-from mjo.symbolic import matrix_simplify_full
def is_symmetric_psd(A):
sage: A = matrix(SR, [[0, 2, 3], [2, 0, 0], [3, 0, 0]])
sage: X = factor_psd(A)
- sage: A2 = matrix_simplify_full(X*X.transpose())
+ sage: A2 = (X*X.transpose()).simplify_full()
sage: A == A2
True