[论文解读] Fast Krylov Space Methods for Calculation of Quark Propogator

Artan Boriçi, Philippe de Forcrand|ArXiv.org|May 2, 1994

Particle physics theoretical and experimental studies被引用 7

一句话总结

本文评估了用于求解格点QCD狄拉克方程的Krylov子空间方法，证明稳定双共轭梯度法（BiCGstab2）是计算威尔逊夸克传播子的鲁棒且高效的求解器。该方法通过减少迭代次数并提高在大格点上的收敛稳定性，显著加速了数值模拟。

ABSTRACT

Different recently developed Krylov space methods for solving linear systems are studied and compared for the solution of the Dirac equation on the lattice. Stabilized Biconjugate Gradient (BiCGstab2) is shown to be a robust and efficient solver for the calculation of the Wilson quark propagator in lattice QCD.

研究动机与目标

• 识别并评估格点QCD模拟中狄拉克方程的高效迭代求解器。
• 根据收敛速度和数值稳定性，比较各种Krylov子空间方法的性能。
• 确定在实际格点QCD应用中计算威尔逊夸克传播子的最合适方法。
• 提供一种计算成本低且保持精度的鲁棒数值解。

提出的方法

• 本研究将Krylov子空间方法（包括BiCGstab2）应用于威尔逊费米子离散化产生的大型稀疏线性系统。
• 选择BiCGstab2是因其稳定化特性，可减少标准BiCG方法中常见的振荡行为。
• 在使用标准QCD参数和典型规范场组态的格点配置上测试了这些方法。
• 通过残差范数和迭代次数监控收敛性，以评估效率和鲁棒性。
• 实现使用了针对稀疏矩阵-向量运算优化的标准线性代数原语。
• 在多个测试案例上进行了对比分析，以评估稳定性和可扩展性。

实验结果

研究问题

• RQ1 哪种Krylov子空间方法在威尔逊夸克传播子计算中展现出最稳定且最快的收敛性？
• RQ2 BiCGstab2在迭代次数和残差减少方面相较于其他迭代求解器表现如何？
• RQ3 BiCGstab2能否在不同格点尺寸和规范配置下保持数值稳定性？
• RQ4 在真实的QCD模拟中，BiCGstab2的计算成本相较于其他方法如何？

主要发现

• 在测试配置中，BiCGstab2展现出优于标准BiCG及其他Krylov方法的收敛行为。
• 该方法在较大格点上实现了显著减少的迭代次数，并保持了稳定的收敛性。
• BiCGstab2对由物理规范场引起的病态矩阵表现出强健性。
• 该算法在多个测试案例中均保持了较低的残差范数，表明具有高数值精度。