Matrix expressions generally follow the noncommutative algebra and I'm not sure that sympy have much framework built for noncommutative polynomials yet. Someone could try to implement such ...