前处理 — CAE用语解说

这是个好问题。虽然名字容易混淆，但两者完全不同。网格生成等前处理（preprocessing）是为求解器准备输入数据的工程。而前处理（preconditioner）是在求解线性方程组 $\mathbf{Ax}=\mathbf{b}$ 时，加速迭代法收敛的数学变换技巧。

如果矩阵A的条件数很大，CG法或GMRES等迭代法的收敛会非常缓慢。使用前处理矩阵M将方程变换为 $\mathbf{M}^{-1}\mathbf{A}\mathbf{x} = \mathbf{M}^{-1}\mathbf{b}$ ，条件数会大幅改善，迭代次数急剧减少。在实际应用中，有无前处理的计算时间差异有时超过10倍。

理想情况是M ≈ A且M^{-1}的计算成本低廉。代表性的选择有：雅可比（只用对角元素·超简单但效果有限）、ILU（不完全LU分解·通用性强）、AMG（代数多重网格·大规模问题最有力）等，这些是实际应用的标准做法。

取决于问题规模。对于数十万自由度的问题，ILU通常就足够了。但当超过1000万自由度时，ILU的迭代次数会过多，AMG会显著更快。AMG相对于问题规模的扩展性几乎是线性的，这是它的优势。

雅可比完全可并行，但效果较弱。ILU本质上是串行处理，并行性不高。AMG具有中等并行性。基于区域分割法（DDM）的前处理在并行中最优，在数千核以上的HPC环境中成为主流。

这种认识很重要。很多人收敛缓慢时都以为是"网格不好"，但只改变前处理就能大幅改善的例子不在少数。养成经常检查计算日志中迭代次数的习惯很有帮助。