[论文解读] Accelerating the Hybrid Monte Carlo algorithm with ILU preconditioning

Mike Peardon|ArXiv.org|Nov 24, 2000

Nuclear reactor physics and engineering被引用 4

一句话总结

该论文提出了一种两步ILU预条件化方法——称为eo-ILU——用于加速格点QCD模拟中的混合蒙特卡罗（HMC）算法。通过先对偶-奇子晶格进行不完全LU分解，再对每个子晶格内的格点进行分解，该方法减少了随机噪声，并使分子动力学步长增加两倍，从而在两味 Schwinger 模型中将计算成本降低了两倍。

ABSTRACT

The pseudofermion action of the Hybrid Monte Carlo (HMC) algorithm for dynamical fermions is modified to directly incorporate Incomplete LU (ILU) factorisation. This reduces the stochastic noise and allows a larger molecular dynamics step-size to be taken, cutting the computational cost. Numerical tests using the two-flavour Schwinger model are presented, where a two-step ILU preconditioning of the even-odd fermion matrix allows the step-size to be increased by a factor of two over the standard even-odd formulation.

研究动机与目标

• 通过改善费米子矩阵求逆的收敛性，降低含动态费米子的HMC模拟的计算成本。
• 解决HMC中矩阵求逆的高成本问题，该问题在格点QCD模拟中主导了计算工作量。
• 探究是否可将ILU预条件化直接应用于伪费米子作用量，而不仅限于费米子矩阵，以提升算法性能。
• 评估不同格点排序对HMC效率的影响，并在并行计算环境中识别最优预条件化方案。
• 将ILU预条件化的适用性扩展至改进的费米子作用量及其他含动态费米子的算法。

提出的方法

• 通过修改作用量项中的逆费米子矩阵，将ILU预条件化直接引入伪费米子作用量。
• 采用两级预条件化：首先进行偶-奇（红-黑）分解，然后在每个子晶格上使用局部格点排序进行ILU分解。
• 在子晶格上使用全局字典序以实现矩阵求逆的最优收敛性，而局部的“四费米子”分解则提升了HMC的接受率。
• 采用BiCGStab作为预条件化系统的迭代求解器，预条件子形式为 (I−L)⁻¹M(I−U)⁻¹。
• 在16×16晶格的两味Schwinger模型中测试该方法，将HMC性能与标准偶-奇HMC和未修改的HMC进行比较。
• 分析厄米特费米子矩阵Q的条件数作为算法性能的代理指标，将其与步长和接受率关联。

实验结果

研究问题

• RQ1 能否有效将ILU预条件化直接应用于伪费米子作用量，以减少随机噪声并提升HMC效率？
• RQ2 在使用ILU预条件化时，格点排序的选择如何影响HMC算法的性能？
• RQ3 两步预条件化方案（eo-ILU）对HMC模拟中分子动力学步长和接受率有何影响？
• RQ4 eo-ILU预条件化带来的性能提升是否可推广至四维QCD和改进的费米子作用量？
• RQ5 费米子矩阵的条件数与HMC算法的接受率和步长之间是否存在相关性？

主要发现

• 在Schwinger模型中，与标准偶-奇预条件化相比，两步eo-ILU预条件化方案使分子动力学步长提高了两倍。
• HMC性能相比标准偶-奇伪费米子HMC提升了两倍，相比未修改的HMC提升了三倍。
• Metropolis检验的接受率对格点排序的依赖性较弱，但与厄米特费米子矩阵Q的条件数强相关。
• 尽管全局字典序在矩阵求逆中为最优，但在子晶格上采用局部格点排序（如“四费米子”分解）可获得比其更好的HMC性能。
• 逆预条件化矩阵的跳跃参数展开从原始情况的O(κ)变为两步eo-ILU方案中的O(κ⁴)，表明噪声减少且稳定性提高。
• 该方法可直接应用于四维规范理论（如QCD），并可扩展至如Sheikholeslami-Wohlert和D234等改进的费米子作用量。