sympy.factor() doesn't have an effect on matrix?

[mgeier]

When I call sympy.expand() on a matrix, all matrix elements get expanded, as expected.

However, when I call sympy.factor(), nothing happens?

>>> import sympy as sp
>>> a, b = sp.symbols('a b')
>>> one = sp.Matrix([(a - b)**2])
>>> one
Matrix([[(a - b)**2]])
>>> two = sp.expand(one)
>>> two
Matrix([[a**2 - 2*a*b + b**2]])
>>> sp.factor(two)
Matrix([[a**2 - 2*a*b + b**2]])

Expected result:

>>> sp.factor(two)
Matrix([[(a - b)**2]])

Am I doing something wrong?

I'm using the master branch of SymPy.

[jksuom]

The method M.applyfunc(factor) could probably be used.

[oscarbenjamin]

See #18181

[mgeier]

Thanks for the quick answers!

It's good to know that this is a known issue. I will try to use applyfunc() in the meantime.

[sidhu1012]

@oscarbenjamin I would like to add .factor() to matrix if needed because after reading #18181 I'm not sure if it's needed or not . Shall I?