[论文解读] Deflation Methods in Fermion Inverters

Walter Wilcox|ArXiv.org|Oct 9, 2007

Quantum Chromodynamics and Particle Interactions被引用 3

一句话总结

本论文提出并比较了多种先进的降压方法——Morgan/Wilcox的GMRES-DR与BiCGStab-DR、Stathopoulos/Orginos的基于CG的降压方法，以及Lüscher的区域分解GCR方法——以加速格点QCD中的费米子求逆。通过投影出狄拉克算子的小特征值模，这些方法减少了Krylov子空间迭代次数，在多个夸克质量下实现了高达6.6倍的加速，并显著改善了随格点体积增长的缩放特性，尤其在小夸克质量下，标准求解器会遭遇临界减速时表现更优。

ABSTRACT

I will review recent developments in matrix deflation methods, by Ronald Morgan/Walter Wilcox, Andreas Stathopoulos/Konstantinos Orginos, and Martin Lüscher, with application to lattice QCD fermion inversion. I will begin with a short review of deflation-related work in the field. The Morgan/Wilcox algorithms using GMRES and BiCGStab for deflation will be described. Typical results for quenched (Wilson and twisted mass) and dynamical configurations (CP-PACS and ETMC) will be displayed and discussed. I will outline how inclusion of multiple right-hand sides and multiple shifts can be accommodated within versions of these methods. I will also review deflation methods based on Conjugate Gradient introduced by Stathopoulos and Orginos for hermitian systems with multiple right-hand sides as well as Lüscher's non-hermitian domain decomposition algorithm using Generalized Conjugate Residual.}

研究动机与目标

• 解决因小特征值导致的病态性而引起的格点QCD费米子求逆在小夸克质量下的临界减速问题。
• 开发并比较投影出低能级特征模的降压技术，以加速Krylov子空间求解器。
• 通过重用降压后的特征向量，实现对多个右端向量和多个夸克质量的高效求解。
• 评估不同格点尺寸和费米子实现下，降压方法的可扩展性与性能表现。
• 评估降压对真空与动力学配置中迭代次数与计算成本的影响。

提出的方法

• 使用GMRES-DR与BiCGStab-DR方法，通过在Krylov子空间中引入谐波Ritz向量，同时计算近似特征向量并求解线性系统。
• 采用基于CG的降压方法处理多个右端向量，在每个周期开始和重启时投影特征向量以改善收敛性。
• 实现Lüscher的非厄米区域分解方法，结合广义共轭残差（GCR）与代数多重网格预处理（SAP），以提升可扩展性。
• 通过从Krylov子空间中投影出特征模，有效消除其对迭代求解器收敛性的影响。
• 通过在不同线性系统间重用降压后的特征向量，支持多个移位量和多个右端向量的求解。
• 在降压子空间中使用谐波Ritz向量作为近似特征向量，确保其为矩阵与子空间的Ritz对。

实验结果

研究问题

• RQ1 小特征模的降压如何改善格点QCD费米子矩阵在Krylov子空间求解器中的收敛性？
• RQ2 与标准求解器（如BiCGStab）相比，GMRES-DR与基于CG的降压方法在多个夸克质量下的性能提升如何？
• RQ3 在动力学费米子模拟中，降压方法随格点体积增大时的缩放特性如何，尤其在大体积下？
• RQ4 降压方法能否有效应用于非厄米系统，如格点QCD中的威尔逊与扭曲质量狄拉克算子？
• RQ5 在求解多个右端向量时，降压的计算成本与所获得加速的相对关系如何，尤其在考虑开销的情况下？

主要发现

• 在大格点（32³×48）上，包含降压开销时，12个夸克质量求解的综合加速因子达到6.6，相比标准BiCGStab方法。
• Morgan/Wilcox的GMRES-DR方法在体积增大3.2倍时，外层GCR迭代次数仅增加13%，而标准BiCGStab方法则增加了30%。
• Lüscher的区域分解方法在相同体积变化下，外层迭代次数仅增加13%，显示出更优的缩放性能。
• 对于5个质量求解，降压开销为184秒，但包含此成本后加速比仍保持在3.1的显著水平。
• 在相同体积变化下，Morgan/Wilcox与Lüscher方法的迭代次数增长比分别约为0.45和0.43，表明其算法加速效果相近。
• 降压方法显著减少了所需Krylov迭代次数，尤其在小夸克质量下，标准求解器会遭遇临界减速。