When using the sympy
library, what is the most efficient way of operating with matrices whose entries are polynomials?
I tried both storing symbols and storing Poly
instances as entries of a Matrix
object. In each case, the performance that I get for operations like matrix multiplication seems to be an order of magnitude slower than what I get with sage
.