[论文解读] Krylov space solvers for shifted linear systems

Beat Jegerlehner|ArXiv.org|Dec 15, 1996

Matrix Theory and Algorithms参考文献 8被引用 100

一句话总结

该论文提出了一种统一框架，用于构建移位Krylov子空间求解器——特别针对CG、CR、BiCGstab和CG-M方法——通过利用移位多项式，仅使用单个系统的矩阵-向量运算即可同时求解多个移位线性系统 $(A + \sigma)x = b$。关键贡献在于实现了格点QCD计算中的最优性能，尤其是在阶梯费米子情况下，当需求解多个夸克质量时可获得显著加速。

ABSTRACT

We investigate the application of Krylov space methods to the solution of shifted linear systems of the form (A+σ) x - b = 0 for several values of σsimultaneously, using only as many matrix-vector operations as the solution of a single system requires. We find a suitable description of the problem, allowing us to understand known algorithms in a common framework and developing shifted methods basing on short recurrence methods, most notably the CG and the BiCGstab solvers. The convergence properties of these shifted solvers are well understood and the derivation of other shifted solvers is easily possible. The application of these methods to quark propagator calculations in quenched QCD using Wilson and Clover fermions is discussed and numerical examples in this framework are presented. With the shifted CG method an optimal algorithm for staggered fermions is available.

研究动机与目标

• 开发一种统一框架，用于构建移位Krylov子空间求解器，以最小计算成本同时求解多个移位系统 $(A + \sigma)x = b$。
• 解决在冻结QCD中夸克传播子计算的计算瓶颈问题，其中必须高效地对多个 $\sigma$ 值（夸克质量）进行求逆。
• 将短递推Krylov方法（如CG、BiCGstab和CR）扩展为移位变体，每个质量值仅需额外两个向量。
• 确保在舍入误差和32位浮点数运算下，移位系统的数值稳定性和收敛性。
• 在使用威尔逊、Clover和阶梯费米子的真实格点QCD模拟中，展示方法的实际效率和可扩展性。

提出的方法

• 定义移位多项式 $P_n^\sigma(z + \sigma) = c_n^\sigma P_n(z)$，使得 $P_n^\sigma(A)$ 生成的向量是 $P_n(A)$ 生成向量的缩放版本，从而实现在不同 $\sigma$ 值间复用。
• 利用Krylov子空间的不变性 $\mathcal{K}_n(A, v_0) = \mathcal{K}_n(A + \sigma, v_0)$，确保移位多项式能为每个 $\sigma$ 生成有效迭代向量，而无需额外的矩阵-向量乘积。
• 通过推导传播移位多项式结构的递推关系，构建CG、CR、BiCGstab和CG-M的移位变体。
• 实施线性预条件（如公式4.63中的 $a=0$），以稳定收敛性，尤其在小夸克质量时，并确保条件(4.62)得以保持，避免停滞。
• 将该方法应用于使用威尔逊、Clover和阶梯费米子的夸克传播子计算，采用tadpole改进参数，并使用非局部源和点源进行测试。
• 通过监控移位系统的残差，验证收敛性和稳定性，特别是在32位算术中舍入效应可能影响移位残差的情况下。

实验结果

研究问题

• RQ1 能否通过统一的多项式框架，系统地从标准短递推求解器（如CG和BiCGstab）推导出移位Krylov子空间方法？
• RQ2 如何仅使用求解最慢单个系统所需的矩阵-向量运算，求解多个移位系统 $(A + \sigma)x = b$？
• RQ3 舍入误差对移位残差收敛性的影响是什么？在实际实现中如何缓解？
• RQ4 在具有非局部源的格点QCD中，移位求解器的性能如何随质量数量及其间距变化而扩展？
• RQ5 高阶多项式预条件是否可有效应用于极小夸克质量下的移位求解器？

主要发现

• 对于未提供高效预条件的阶梯费米子，移位CG-M和CR-M方法为最优，可实现显著加速。
• 当内存允许时，BiCGstab-M方法是威尔逊和Clover费米子的首选，提供稳定收敛和高效率。
• 对于非局部源，当夸克质量彼此接近时，性能提升因子可能非常显著，其中 $\sum_{i=1}^{n}N^{\rm cont}_{i} - 2N^{\rm zero\text{-}guess}_{n} \gg 0$ 表明收益巨大。
• 数值测试表明，由于舍入误差，某些情况下移位系统的残差可能停滞在约 $10^{-2}$ 水平，尤其是在大尺寸格点中，因此必须进行残差检查。
• 线性预条件可稳定小质量下的收敛性，且在正确实现下，该方法在32位算术中仍具可行性。
• 该框架可轻松扩展至其他Krylov方法，并支持对极小质量的高阶预条件。