Python 以更快的方式乘以numpy矩阵

我有3个矩阵，我想乘以A*B*C

A是维数为k×d的非常稀疏的矩阵

B是d×d的阿达玛矩阵

C是d×d的对角矩阵

d的数量级为100000，而k的数量级为1000

---

是否有人可以建议是否有更好的方法来表示稀疏矩阵，并更快地执行此操作

只是为了澄清一下，您是否正在寻找改进scipy.sparse实现性能的方法？我认为对于你的维度，scipy.sparse.csr_矩阵应该相当快。我不知道scipy.sparse，我正在研究它。谢谢