[论文解读] An adaptive aggregation based domain decomposition multilevel method for the lattice wilson dirac operator: multilevel results

Andreas Frommer, Karsten Kahl+3 more|arXiv (Cornell University)|Jul 23, 2013

Matrix Theory and Algorithms参考文献 14被引用 12

一句话总结

本论文提出了一种基于自适应聚合的多水平区域分解代数多重网格方法（DD-α AMG），用于格点QCD中的格点威尔逊狄拉克算符。通过将两水平DD-α AMG求解器扩展为使用自适应插值和乘法Schwarz平滑的真正多水平框架，该方法在保持强可扩展性和γ₅对称性的同时，相较于标准Krylov方法实现了700倍的加速，相较于两水平版本提升了1.7倍。

ABSTRACT

In lattice QCD computations a substantial amount of work is spent in solving linear systems arising in Wilson's discretization of the Dirac equations. We show first numerical results of the extension of the two-level DD-αAMG method to a true multilevel method based on our parallel MPI-C implementation. Using additional levels pays off, allowing to cut down the core minutes spent on one system solve by a factor of approximately 700 compared to standard Krylov subspace methods and yielding another speed-up of a factor of 1.7 over the two-level approach.

研究动机与目标

• 将基于自适应聚合的两水平区域分解代数多重网格方法（DD-α AMG）扩展为真正多水平框架，以求解格点威尔逊狄拉克算符。
• 通过减少威尔逊费米子离散化产生的大型线性系统求解时间，提升格点QCD模拟的计算效率。
• 通过保持对称性的伽辽金粗网格转移方法，在所有粗网格层级上维持狄拉克算符的γ₅对称性。
• 通过使用MPI和自适应设置例程，实现在大规模并行架构上的可扩展性能。
• 探究在多重网格层次中增加额外层级是否能带来相较于两水平方法的可测量性能提升。

提出的方法

• 该方法采用V循环多重网格策略，所有层级均使用乘法Schwarz（SAP）平滑，确保良好的数据局部性和强可扩展性。
• 采用自适应聚合构造插值算子 $ P_\nu $，分别对自旋-1,2和自旋-3,4自由度进行分组，以保持γ₅对称性。
• 粗网格算符通过伽辽金投影构造：$ D_{\nu+1} = P_\nu^\top D_\nu P_\nu $，确保对称性保持和误差传播的一致性。
• 提出一种新颖的设置例程，结合反幂迭代与自举技术，高效生成用于粗网格校正子空间构造的测试向量。
• 该方法作为预条件子应用于灵活Krylov子空间求解器（如FGMRES）中，粗网格系统采用近似求解以降低设置开销。
• 多级层次结构递归应用，直至最粗网格系统小到可直接求解，从而实现完整的多重网格递归。

实验结果

研究问题

• RQ1 将两水平DD-α AMG方法扩展为真正多水平框架，是否能在求解格点威尔逊狄拉克算符时带来可测量的性能提升？
• RQ2 自适应聚合、SAP平滑与对称性保持的粗网格转移相结合，是否能在多级之间保持鲁棒的收敛性？
• RQ3 相较于两水平版本，增加额外层级对求解时间与可扩展性有何影响？
• RQ4 γ₅对称性保持对多水平方法的有效性与稳定性有何影响？
• RQ5 多级设置能否足够高效，以证明额外层级带来的计算开销是合理的？

主要发现

• 三水平DD-α AMG方法在128个进程配置下，单次求解相较于标准Krylov子空间方法实现了约700倍的加速。
• 与两水平DD-α AMG方法相比，该多级方法性能提升了1.7倍，表明额外层级确实带来了可测量的性能增益。
• 该方法在所有层级上均保持了γ₅对称性，这对格点QCD模拟中的物理一致性至关重要。
• 三水平方法的设置时间略高于两水平变体，但更短的求解时间足以弥补，即使仅针对单个右端向量也是如此。
• 数值实验表明，该多级方法在增大格点尺寸时表现出有利的可扩展性，并在临界夸克质量附近展现出更强的鲁棒性。
• 尽管由于粗网格复杂度降低，不精确降维结合不精确投影在设置阶段更快，但该方法在速度和对称性保持方面仍优于后者。