[论文解读] Deation and Flexible SAP-Preconditioning of GMRES in Lattice QCD Simulation

Andreas Frommer, Andrea Nobile+1 more|arXiv (Cornell University)|Apr 24, 2012

Matrix Theory and Algorithms参考文献 30被引用 7

一句话总结

本文提出一种灵活的GMRES-DR方法，结合降维和乘法Schwarz预条件化（SAP），以求解格点QCD模拟中出现的病态线性系统。通过在GMRES重启过程中保留信息，该方法在标准重启GMRES失效时仍能实现收敛，显著提升了在并行架构上对威尔逊-狄拉克算子求逆的鲁棒性。

ABSTRACT

The simulation of lattice QCD on massively parallel computers stimulated the development of scalable algorithms for the solution of sparse linear systems. We tackle the problem of the Wilson-Dirac operator inversion by combining a Schwarz alternating procedure (SAP) in multiplicative form with a exible variant of the GMRES-DR algorithm. We show that restarted GMRES is not able to converge when the system is poorly conditioned. By adding deation in the form of the FGMRES-DR algorithm, an important fraction of the information produced by the iterates is kept between successive restarts leading to convergence in cases in which FGMRES stagnates.

研究动机与目标

• 解决在格点QCD模拟中求解大规模、稀疏、病态线性系统时，重启GMRES方法的收敛失败问题。
• 提升在大规模并行计算机上对威尔逊-狄拉克算子求逆的迭代求解器的鲁棒性。
• 将乘法Schwarz预条件化（SAP）与灵活GMRES-DR结合，以在重启过程中保留有用信息。
• 证明降维可防止停滞，并在标准FGMRES失效的情况下实现收敛。

提出的方法

• 使用乘法Schwarz交替过程（SAP）作为预条件器，以改善系统条件。
• 应用GMRES-DR的灵活变体（FGMRES-DR），以允许在迭代过程中使用变化的预条件器。
• 实施降维，以保留来自先前Krylov子空间的近似不变子空间信息。
• 维持一种灵活的预条件器结构，可随Krylov子空间基的变化而自适应调整。
• 将SAP与FGMRES-DR结合，以增强在条件较差系统上的收敛行为。
• 采用带降维的重启策略，在控制内存和计算成本的同时保持收敛稳定性。

实验结果

研究问题

• RQ1 在格点QCD的病态线性系统求解中，带降维的FGMRES-DR能否克服停滞问题？
• RQ2 SAP与灵活GMRES-DR的结合相比标准重启GMRES，如何提升收敛性能？
• RQ3 通过降维在重启之间保留信息，对收敛鲁棒性有何影响？
• RQ4 灵活预条件化策略在大规模并行架构上是否保持可扩展性？
• RQ5 该方法能否在标准方法失效的情况下，可靠地实现格点QCD模拟中威尔逊-狄拉克算子的求逆？

主要发现

• 带降维的FGMRES-DR在标准重启GMRES因条件差而停滞的系统上成功实现收敛。
• 降维保留了先前Krylov子空间中的关键信息，防止了重启过程中收敛进度的丢失。
• SAP与灵活GMRES-DR的结合为威尔逊-狄拉克算子求逆提供了鲁棒且可扩展的解决方案。
• 通过仅在重启之间保留一个小型相关子空间，该方法保持了计算效率。
• 该方法在标准迭代求解器失效的困难格点QCD场景中实现了收敛。