mmasub_mma_

移动渐进法对稀疏矩阵结构的操作则比较复杂，和排序、消去数等图论技巧关系密切，比如稀疏矩阵的求和与三角分解会产生填充（fill-in），两个稀疏矩阵的乘积的高效实现至今仍是一个具有挑战性的问题。另一方面，稀疏矩阵运算，比如三角分解，还要考虑数值稳定性问题，稀疏性和数值稳定性往往是一对矛盾。MMA