做一个稀疏-稀疏矩阵乘积几乎总是低效的$C=AB$. 简单地考虑几乎总是更有效$C$作为可以应用于向量的运算符，即计算$y=Cx$. 每当需要这样的操作时，您可以分两步进行计算：$y=A(Bx)$，即，您只需连续执行两个矩阵向量乘积即可得到结果。为此，不必实际形成乘积矩阵$AB$.

如果您绝对必须形成产品，请查看https://svn.dealii.org/trunk/deal.II/include/deal.II/lac/sparse_matrix.templates.h并搜索 mmult() 函数。