[論文レビュー] FEAST as Subspace Iteration Accelerated by Approximate Spectral Projection

Ping Tang, Eric Polizzi|arXiv (Cornell University)|Feb 2, 2013

Matrix Theory and Algorithms被引用数 4

ひとこと要約

本論文は、FEAST固有値アルゴリズムの最初の厳密な収束解析を提供し、それがスペクトル射影子への有理関数近似によって加速された部分空間反復法として機能することを示している。この解析により、近似的なスペクトル射影子に依存しているにもかかわらずFEASTが収束することを証明し、エルミート固有値問題を解く際のその頑健性と並列効率の理論的基盤を確立している。

ABSTRACT

The calculation of a segment of eigenvalues and their corresponding eigenvectors of a Hermitian matrix or matrix pencil has many applications. A new density-matrix-based algorithm has been proposed recently and a software package FEAST has been developed. The density-matrix approach allows FEAST's implementation to exploit a key strength of modern computer architectures, namely, multiple levels of parallelism. Consequently, the software package has been well received, especially in the electronic structure community. Nevertheless, theoretical analysis of FEAST has lagged. For instance, the FEAST algorithm has not been proven to converge. This paper offers a detailed numerical analysis of FEAST. In particular, we show that the FEAST algorithm can be understood as an accelerated subspace iteration algorithm in conjunction with the Rayleigh-Ritz procedure. The novelty of FEAST lies in its accelerator which is a rational matrix function that approximates the spectral projector onto the eigenspace in question. Analysis of the numerical nature of this approximate spectral projector and the resulting subspaces generated in the FEAST algorithm establishes the algorithm's convergence. This paper shows that FEAST is resilient against rounding errors and establishes properties that can be leveraged to enhance the algorithm's robustness. Finally, we propose an extension of FEAST to handle non-Hermitian problems and suggest some future research directions.

研究の動機と目的

• FEAST固有値アルゴリズムに厳密な理論的基盤を提供すること。これは、これまで形式的な収束解析が欠落していたためである。
• FEASTのメカニズムを、目的固有空間へのスペクトル射影子の近似による部分空間反復の加速として説明し、その数値的挙動を明確化すること。
• 丸め誤差に対するFEASTの耐性を分析し、アルゴリズムの頑健性を高める要因を同定すること。
• FEASTの理論枠組みを非エルミート固有値問題へ拡張することにより、今後の研究方向性を示唆すること。

提案手法

• FEASTを、各反復で目的固有空間へのスペクトル射影子を近似する有理行列関数によって加速する部分空間反復法として再定式化する。
• 各ステップで生成された部分空間からRitz対を抽出するためにレイリー・リツの手続きを用い、最適なRitzベクトルと固有値を保証する。
• スペクトル射影子への有理関数近似の数値的性質を分析し、それが望ましい固有空間を効果的に捉えていることを示す。
• 標準的な仮定の下で、FEASTが生成する部分空間が目的不変部分空間に近づくことを証明することで収束を確立する。
• 後退誤差と摂動の影響を分析することで、有限精度算術下でもアルゴリズムが安定していることを示す。
• 複素曲線と有理関数へのスペクトル射影子近似の適応により、非エルミート問題へのFEASTの拡張を提案する。

実験結果

リサーチクエスチョン

• RQ1 FEASTアルゴリズムは、近似的なスペクトル射影子に依存しているにもかかわらず、エルミート固有値問題に対して収束するのか？
• RQ2 スペクトル射影子への有理関数近似が、部分空間反復の収束速度と安定性にどのように影響するか？
• RQ3 丸め誤差に対する耐性を高めるために、近似的なスペクトル射影子が有する性質は何か？
• RQ4 FEASTの理論枠組みを非エルミート固有値問題へ拡張することは可能か？
• RQ5 部分空間反復に有理関数加速を施したという解釈から、どのような理論的知見が得られるか？

主な発見

• FEASTがスペクトル射影子への有理関数近似によって加速された部分空間反復法として収束することが証明され、長年の理論的空白が解消された。
• 近似的なスペクトル射影子の良好な条件数とレイリー・リツ手続きのおかげで、丸め誤差に対して頑健であることが示された。
• 収束速度は、スペクトル射影子への有理関数近似の品質に依存し、より良い近似がより速い収束をもたらす。
• 部分空間反復の枠組みが、実装上観察されるFEASTの頑健性と並列スケーラビリティを説明する。
• 非エルミート問題への拡張は理論的に妥当であり、有理ケイリーロフ法の分野における新たな研究方向性を示唆する。
• 本分析により、スペクトル射影子のより良い近似と誤差制御の強化を通じたFEASTの頑健性向上の基盤が提供された。