Symbolic matrix calculation
Is it possible to do symbolic matrix calculation with sage? What I like to do is somehow declare that U_n are unitary matrices (without giving their dimension or explicit representation). Once declared I'd like sage to exploit the unitarity when it is evaluating expressions involving U_n, like:
A=U_1.H * U_1 * U_2
should yield
A=U_2
Is this possible?