矩阵求解泛谈

         矩阵求解AX = b 在工程上使用相当广泛，尤其是在simulation领域。矩阵求解的速度相当程度地决定一款模拟软件的命运。矩阵求解的方法很多，如一般的迭代法，投影法，AMG等。

         1）迭代法

               不实用，尤其是求解大型矩阵。但可以作为投影法的prediction算子存在。如CGJacobi, CGSSOR，这是把Jacobi, SSOR方法当作CG方法的Prediction.

         2) 投影法

              主要是CG(Conjugate gradient method) 和BiCG. CG适用的前提条件是对称正定，而BiCG只需要正定。BiCG这一类型方法中一般有CGS，BiCGStable。一般而言，BiCGStable比CGS还是要稳定些的。

        3）AMG

              个人认为，对于矩阵求解，AMG是一个超级武器。关键点是如何寻找AMG中的Coarse DOF in each coarser level，关于这方面的文章相当多。