The following sympy code outputs False instead of True
from sympy import symbols, Identity, MatrixSymbol
d = symbols("d")
I = Identity(d)
A = MatrixSymbol("A", d, d)
print(I - A@(A + I).inv() == (A + I).inv())
Question. How to get sympy to identity that basic matricial equivalence and output True ?