[論文レビュー] A fast minimal residual solver for overlap fermions

Artan Boriçi, Alban Allkoci|ArXiv.org|Feb 9, 2006

Particle physics theoretical and experimental studies参考文献 12被引用数 4

ひとこと要約

本稿では、格子QCDにおけるシフトされたユニタリ系に特化した、Arnoldi法の新規3項再帰を導入し、2つの新しいクリロフ部分空間ソルバー（SUOMおよびSHUMR）を構築する。SHUMRアルゴリズムは、全テストされた格子サイズおよびクォーク質量において、以前で最も速かった方法であるSUMRよりも10–15%高速に収束し、大規模なオーバーラップフェルミオン線形系を解く上で優れた効率性を示している。

ABSTRACT

Computing quark propagators with overlap fermions requires the solution of a shifted unitary linear system. Jagels and Reichel have shown that for such systems it is possible to construct a minimal residual algorithm by short recurrences. The Jülich-Wuppertal group have found this algorithm to be the fastest among overlap solvers. In this paper we present a three-term recurrence for the Arnoldi unitary process. Using the new recurrence we construct a minimal residual solver which is the fastest among all Krylov subspace algorithms considered so far for the overlap inversion.

研究の動機と目的

• 格子QCDにおけるオーバーラップフェルミオンのクォーク伝播関数計算のための、より効率的な反復的ソルバーの開発。
• 標準的なGMRESがシフトされたユニタリ系に適用される際の長い再帰による計算ボトルネックの解消。
• オーバーラップフェルミオン系において、残差の最小化の最適性を保ちつつ、短い再帰を実現するクリロフ部分空間法の構築。
• 既存のSUMRアルゴリズムを、3項再帰を用いた新しいユニタリArnoldiプロセスを導入することで改善。
• CGNE、CR、SUMRといった既存の手法と比較して、現実的な格子構成における新規ソルバーの妥当性の検証。

提案手法

• ユニタリArnoldiプロセスのための新規3項再帰を提案し、標準的なArnoldi法と比較して計算コストとメモリ使用量を削減。
• 残差ベクトルをクリロフ部分空間に直交させる、シフトされたユニタリ直交法（SUOM）を導出。
• 新規再帰を用いて、クリロフ部分空間上で残差ノルムを最小化する、シフトされたユニタリ最小残差法（SHUMR）を導入。
• 過去のベクトルを完全に蓄積することを避けるために、短い再帰により基底ベクトルの正規直交性を維持する修正Arnoldiプロセスを採用。
• Neubergerのオーバーラップディラック作用素（シフトされたユニタリ行列）の構造を活用し、ユニタリ性を活かして短い再帰を可能に。
• MATLABでアルゴリズムを実装し、高性能計算用にC++に移植可能なuBLASライブラリを用いたコードを提供。

実験結果

リサーチクエスチョン

• RQ1 ユニタリArnoldiプロセスに対して3項再帰を導出可能か。これにより、シフトされたユニタリ系のための短い再帰クリロフソルバーが実現可能か。
• RQ2 新規SHUMRアルゴリズムは、オーバーラップフェルミオン系において、以前で最も最適であったSUMRソルバーを上回る収束速度を達成するか。
• RQ3 新規ソルバー（SUOMおよびSHUMR）は、さまざまなクォーク質量および格子結合定数において、収束速度およびロバスト性の観点でどのように比較されるか。
• RQ4 クリロフ部分空間法におけるオーバーラップフェルミオンの収束挙動に、幾何的最適性と代数的最適性の影響は何か。
• RQ5 新規再帰法は、大規模な格子QCDシミュレーションにおける数値的安定性および効率性を向上させられるか。

主な発見

• 新規3項再帰により、ユニタリArnoldiプロセスに対して初めて短い再帰クリロフ部分空間法が実現され、計算およびメモリのオーバーヘッドが削減された。
• SHUMRアルゴリズムは、全テストされた構成（軽いクォーク質量 m=0.005、β=5.7およびβ=6.0）において、SUMRよりも10–15%高速に収束した。
• SHUMRは、軽いクォーク質量において超線形収束を示し、残差ノルムが10⁻³（β=6.0）および10⁻⁵（β=5.7）の周辺で収束速度が増加した。これは、CRの線形収束とは対照的である。
• SUOMとSUMRはほぼ同等の収束を示したが、SHUMRは一貫して両者を上回り、クリロフ部分空間の構築方法が収束効率に影響を与えることを示している。
• 新規ソルバー（SHUMRおよびSUOM）は、CGNEおよびCG-CHIよりも効率的であり、重いクォーク質量でもそれらを上回った。
• SHUMRは幾何的に最適であり、両者とも残差ノルムを最小化しているにもかかわらず、SUMRよりも高速に収束する。これは、SHUMRがクリロフ部分空間の探索がより効率的であることを示唆している。