[論文レビュー] Subspace Iteration with Approximate Spectral Projection

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

Matrix Theory and Algorithms被引用数 10

ひとこと要約

この論文は、FEAST固有値アルゴリズムの最初の厳密な収束解析を提供し、スペクトル射影子への有理近似を用いた加速された部分空間反復として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の理論的枠組みを拡張し、新たな応用分野を開く。

提案手法

• FEASTを、目的固有空間へのスペクトル射影子を近似する有理行列関数によって加速された部分空間反復に再定式化する。
• 各反復で生成された部分空間からレイリー・リッツ法を用いてリッツ対を抽出する。
• 近似スペクトル射影子のスペクトル特性を分析することで、部分空間反復の収束を研究する。
• 真のスペクトル射影子への有理近似の精度に基づいて、残差および収束速度の上限を確立する。
• 有理行列関数の構造を活用して、有限精度算術における浮動小数点誤差に対する頑健性を証明する。
• 複素積分路を用いた有理射影への適応により、非エルミート行列へのFEASTフレームワークの一般化を提案する。

実験結果

リサーチクエスチョン

• RQ1 FEASTアルゴリズムは有限精度算術下でも収束するのか。その収束の理論的根拠は何か。
• RQ2 スペクトル射影子への有理近似は、FEASTの収束速度および安定性にどのように影響するか。
• RQ3 近似スペクトル射影子のどのような性質が、丸め誤差に対する数値的頑健性を保証するか。
• RQ4 FEASTフレームワークは、収束性および効率性を保ったまま非エルミート固有値問題に拡張可能か。
• RQ5 FEASTを有理射影を伴う加速された部分空間反復として解釈することで、どのような理論的知見が得られるか。

主な発見

• FEASTが、スペクトル射影子への有理近似によって加速された部分空間反復として形式的に収束することを証明した。
• FEASTの収束速度は、スペクトル射影子への有理近似の精度に依存し、近似が良いほど収束が速くなる。
• 有理行列関数の構造とレイリー・リッツ手順のおかげで、FEASTは丸め誤差に対して耐性を示し、数値的安定性を確保する。
• 理想状態では二次収束の挙動を示すが、有限精度算術下でも収束が保たれる。
• 理論的枠組みにより、複素積分路と有理射影を用いた非エルミート問題へのFEASTの拡張が可能である。
• 解析により、より良い有理近似戦略およびプリコンディショニングを用いたFEASTの頑健性向上の基盤が得られた。