[論文レビュー] Accelerating the Hybrid Monte Carlo algorithm with ILU preconditioning

Mike Peardon|ArXiv.org|Nov 24, 2000

Nuclear reactor physics and engineering被引用数 4

ひとこと要約

本稿では、格子QCDシミュレーションにおけるハイブリッドモンテカルロ（HMC）アルゴリズムの高速化を目的として、二段階のILU前処理法—eo-ILU—を提案する。偶数奇数サブラティスにまず不完全LU分解を適用し、その後各サブラティス内のサイトに対して再びILU分解を施すことで、確率的ノイズを低減し、分子動力学ステップサイズを2倍に拡大可能となる。これにより、2フレーバーのシュヴィンガー模型において計算コストを2倍に削減できる。

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シミュレーションの計算コストを低減するため、フェルミオン行列の反転収束を改善すること。
• 格子QCDシミュレーションにおいて計算負荷の大部分を占める行列反転処理の高コスト問題に対処すること。
• ILU前処理がフェルミオン行列に限らず、擬スピンフェルミオン作用に直接適用可能かどうかを検討し、アルゴリズム性能の向上を図ること。
• サイト順序の選択がHMC効率に与える影響を評価し、並列計算環境に最適な前処理スキームを同定すること。
• ILU前処理の適用範囲を、改良されたフェルミオン作用や他のダイナミカルフェルミオンアルゴリズムへと拡張すること。

提案手法

• 作用項における逆フェルミオン行列を変更することで、擬スピンフェルミオン作用に直接ILU前処理を組み込む。
• 二段階の前処理を実施：まず偶数奇数（赤黒）分解を行い、その後各サブラティス内で局所的なサイト順序を用いてILU因子分解を実行する。
• 行列反転における最適収束を実現するため、サブラティス上でグローバルな辞書式順序を用いる一方で、HMCの受容率向上の観点から局所的な「4フレーバー」分解を併用する。
• 前処理系の反復解法としてBiCGStabを採用し、前処理子を (I−L)⁻¹M(I−U)⁻¹ の形で構築する。
• 16×16格子を用いた2フレーバーのシュヴィンガー模型において、標準的な偶数奇数HMCおよび修正なしHMCと比較して、本手法のHMC性能を評価する。
• アルゴリズム性能の指標として、エルミートフェルミオン行列Qの条件数を用い、ステップサイズおよび受容率と関連付ける。

実験結果

リサーチクエスチョン

• RQ1 ILU前処理を擬スピンフェルミオン作用に直接効果的に適用することで、確率的ノイズを低減し、HMC効率を向上させることができるか？
• RQ2 サイト順序の選択が、ILU前処理を用いたHMCアルゴリズムの性能に与える影響は何か？
• RQ3 二段階の前処理スキーム（eo-ILU）が、HMCシミュレーションにおける分子動力学ステップサイズおよび受容率に与える影響は何か？
• RQ4 eo-ILU前処理による性能向上効果は、4次元QCDや改良されたフェルミオン作用へ一般化可能か？
• RQ5 フェルミオン行列の条件数とHMCアルゴリズムの受容率およびステップサイズの間に相関関係があるか？

主な発見

• 2段階のeo-ILU前処理スキームにより、シュヴィンガー模型における分子動力学ステップサイズが、標準的な偶数奇数前処理と比較して2倍に拡大された。
• HMC性能は、標準的な偶数奇数擬スピンフェルミオンHMCと比較して2倍向上し、修正なしHMCと比較しては3倍向上した。
• メトロポリステストの受容率はサイト順序に弱い依存性を示すが、フェルミオン行列Qの条件数と強く相関していた。
• サブラティス上で局所的なサイト順序（例：「4フレーバー」分解）を用いることで、グローバルな辞書式順序よりも優れたHMC性能が得られた。ただし、行列反転の最適収束には辞書式順序が優れていた。
• 逆前処理行列のヘッジパラメータ展開が、元のO(κ)から2段階のeo-ILUスキームではO(κ⁴)にシフトしており、ノイズ低減と安定性向上を示している。
• 本手法は4次元ゲージ理論（例：QCD）にも直接適用可能であり、シェイクホルスラミ＝ウォーライトやD234作用など、改良されたフェルミオン作用へも拡張可能である。