[論文レビュー] A FEAST Algorithm for Generalized non-Hermitian Eigenvalue Problems

Guojian Yin, Raymond H. Chan+1 more|arXiv (Cornell University)|Apr 7, 2014

Matrix Theory and Algorithms被引用数 12

ひとこと要約

本稿は、一般化非エルミート固有値問題にFEASTアルゴリズムを拡張し、複素平面内のユーザーが定義した領域内で、信頼性のある固有値および固有ベクトルの計算を可能にする。数値実験を通じて、実用性と効率性・安定性を示し、領域内の固有値数を推定するための数学的枠組みを導入する。

ABSTRACT

Recently, a contour-integral based solver for generalized Hermitian eigenvalue problems, called FEAST, was proposed by Polizzi. It calculates the eigenvalues inside a given interval and their associated eigenvectors. This novel algorithm is reliable and accurate, however its applicability is limited since in practical implementation, it has to know in advance the number of eigenvalues inside the interval. In this work, we generalize the FEAST algorithm to non-Hermitian problems to compute eigenvalues inside a region in the complex plane along with their eigenvectors. We establish the mathematical framework and present a way to find the number of eigenvalues inside the given region. Our algorithm is applicable in practical implementation and numerical experiments illustrate that our algorithm is efficient and stable.

研究の動機と目的

• 元々一般化エルミート問題を対象として設計されたFEASTアルゴリズムを、非エルミート固有値問題に対応させる。
• 複素平面内の指定された領域内での固有値および固有ベクトルの計算を可能にする数学的枠組みを確立する。
• 与えられた領域内の固有値数を推定する手法を開発し、元のFEASTアルゴリズムの主な限界を克服する。
• 計算中に固有値数の決定が現実的に行えるようにすることで、実装の実用性を確保する。

提案手法

• アルゴリズムは、複素平面における線積分を用いて、一般化固有値問題を指定された領域内の固有値を含む部分空間に射影する。
• コーシーの積分公式を用いた有理関数フィルタリングを適用し、線積分に沿った領域内の固有ペアを抽出することで、エルミート版FEASTのアプローチを非エルミートケースに一般化する。
• 領域内の固有値数を推定するための新規技術を導入し、線積分されたリソルベント行列のトレースを分析することで実現する。
• 線積分の数値的評価を確実にするために、線積分に基づく求積則を採用し、ロバスト性と正確性を確保する。
• 個々の固有ペアを抽出しながら数値的安定性を維持するためのデフラッティング戦略を統合する。
• 実装は、大規模な非エルミート問題に対する実用的応用に適した効率的かつスケーラブルな設計である。

実験結果

リサーチクエスチョン

• RQ1 FEASTアルゴリズムは、非エルミート一般化固有値問題を解くためにどのように一般化できるか？
• RQ2 ユーザーが定義した複素平面内の領域内での固有値数の正確な推定を可能にする数学的枠組みは何か？
• RQ3 拡張されたFEASTアルゴリズムは、非エルミート問題に適用された際にも高い正確性と安定性を維持できるか？
• RQ4 既存手法と比較して、一般化FEASTアルゴリズムの効率性およびロバスト性はどのように評価できるか？

主な発見

• 一般化FEASTアルゴリズムは、非エルミート問題において、指定された領域内の固有値およびそれに関連する固有ベクトルを効果的に計算できる。
• 領域内の固有値数を信頼性を持って推定する方法を提供し、元のFEASTアルゴリズムの主要な限界を解消する。
• 数値実験により、アルゴリズムがさまざまなテスト問題において効率的かつ数値的に安定であることが確認された。
• 線積分と有理関数フィルタリングの使用により、非エルミート行列に対しても高精度な固有ペア抽出が保証される。
• 収束性とスケーラビリティが優れており、大規模応用に適している。