Sympy: Matrix.doit()不适用于矩阵元素

Question

在下面的示例中，我将一个函数替换为一个矩阵表达式。我希望在.doit()调用之后进行一些简化：

x = symbols('x', real=True)
f = symbols('f', real=True)(x)
v = Matrix([f * sin(x), f * cos(x)])
v1 = v.diff(x)

print v1.subs(f, x)
# prints: Matrix([[x*cos(x) + sin(x)*Derivative(x, x)], [-x*sin(x) + cos(x)*Derivative(x, x)]])
print v1.subs(f, x).doit()
# prints: Matrix([[x*cos(x) + sin(x)*Derivative(x, x)], [-x*sin(x) + cos(x)*Derivative(x, x)]])
print Matrix([e.doit() for e in v1.subs(f, x)])
# prints: Matrix([[x*cos(x) + sin(x)], [-x*sin(x) + cos(x)]])

但这种事不会发生。如果我对矩阵的每个元素调用.doit()，那么一切都是正确的。

问题：什么是不正确的，我的代码或.doit()函数的行为？怎么修呢？

Answer 1

在SymPy的开发版本中，这个问题似乎已经解决了。一个新的版本(1.0)应该会很快发布。