Numpy 多项式构成的矩阵乘法

是否可以使用numpy/scipy对由多项式组成的矩阵进行乘法

具体地说，我希望将一个120乘120的稀疏矩阵相乘，它的条目本身看起来就像一个+7*b+c

老实说，我并没有尽力去做这件事。我看到numpy中有一个多项式模块，但我没有使用它的经验。我只是希望有人看到这一点，并说“很明显，这是可能的，这样做”

---

从我所看到的情况来看，有一个相关的问题曾经被问过：

我不知道稀疏，但是numpy对象数组工作得很好

In [1]: from numpy.polynomial import Polynomial as P

In [2]: a = np.array([[P([1,2]), P([3,4])]]*2)

In [3]: a
Out[3]: 
array([[Polynomial([ 1.,  2.], [-1,  1], [-1,  1]),
        Polynomial([ 3.,  4.], [-1,  1], [-1,  1])],
       [Polynomial([ 1.,  2.], [-1,  1], [-1,  1]),
        Polynomial([ 3.,  4.], [-1,  1], [-1,  1])]], dtype=object)

In [4]: np.dot(a, a)
Out[4]: 
array([[Polynomial([  4.,  14.,  12.], [-1.,  1.], [-1.,  1.]),
        Polynomial([ 12.,  34.,  24.], [-1.,  1.], [-1.,  1.])],
       [Polynomial([  4.,  14.,  12.], [-1.,  1.], [-1.,  1.]),
        Polynomial([ 12.,  34.,  24.], [-1.,  1.], [-1.,  1.])]], dtype=object)