[論文レビュー] Recursive Integral Method with Cayley Transformation

Ruihao Huang, Jiguang Sun+1 more|arXiv (Cornell University)|May 3, 2017

Matrix Theory and Algorithms参考文献 6被引用数 7

ひとこと要約

本稿では、複素平面内の指定された領域に含まれる大規模なスパース非エルミート行列ペンティッドのすべての固有値を計算するための改善版固有値ソルバー、RIM-Cを提案する。Cayley変換とアーノルド法を再帰的積分法（RIM）に統合することで、計算コストを著しく低減し、新たな指標を用いることで安定性が向上し、特にスペクトル情報が入手できない状況下でもMATLABの'eigs'よりも安定性と信頼性に優れる。

ABSTRACT

Recently, a non-classical eigenvalue solver, called RIM, was proposed to compute (all) eigenvalues in a region on the complex plane. Without solving any eigenvalue problem, it tests if a region contains eigenvalues using an approximate spectral projection. Regions that contain eigenvalues are subdivided and tested recursively until eigenvalues are isolated with a specified precision. This makes RIM an eigensolver distinct from all existing methods. Furthermore, it requires no a priori spectral information. In this paper, we propose an improved version of {\bf RIM} for non-Hermitian eigenvalue problems. Using Cayley transformation and Arnoldi's method, the computation cost is reduced significantly. Effectiveness and efficiency of the new method are demonstrated by numerical examples and compared with 'eigs' in Matlab.

研究の動機と目的

• 事前のスペクトル情報が得られない大規模非エルミート固有値問題に対して信頼性の高い固有値ソルバーを開発すること。
• 元のRIM手法が密度行列ソルバーを繰り返し用いることに起因する計算コストを低減すること。
• スペクトル射影に基づく新しい指標を用いることで、与えられた領域における固有値検出の堅牢性と効率性を向上させること。
• 特に従来の方法（例：'eigs'）がユーザーが指定するシフトに強く依存する状況下でも、多様な問題に対して安定的かつ一貫した性能を発揮できるようにすること。

提案手法

• RIM-Cはケーリー変換を用いて、関心領域を単位円に写像し、一般化固有値問題をKrylov部分空間法に適した形に変換する。
• スペクトル射影近似に必要なシフト付き線形系 $(A - z_j B)^{-1}f$ を効率的に解くために、アーノルド法を用いる。
• 再帰的領域分割戦略を適用：固有値を含む領域は、指定された許容誤差 $\epsilon$ 未満のサイズになるまで再帰的に分割する。
• 固有値検出の感度を高めるために、新しい堅牢な指標 $\delta_S = \left| P\left( \frac{Pf}{\|Pf\|} \right) \right|$ を用いる。
• 領域を破棄するか再分割するかを判断するための固定しきい値 $\delta_0 = 0.1$ を使用し、最終的な領域の中心を固有値の近似値とする。
• 全固有値問題を解くのではなく、閉曲線積分と数値積分を用いてスペクトル射影を近似する。

実験結果

リサーチクエスチョン

• RQ1 Krylov部分空間法を用いることで、大規模非エルミート固有値問題に対する再帰的積分法を著しく高速化できるか？
• RQ2 ケーリー変換を用いることで、RIMフレームワーク内でのパrameter化された線形系の条件数と効率性がどのように向上するか？
• RQ3 新しいスペクトル射影指標は、弱いスペクトル成分を有する領域においても、固有値検出の堅牢性を向上させられるか？
• RQ4 適切なシフトが得られない状況下で、RIM-Cは'eigs'と比較して性能と安定性に優れているか？
• RQ5 RIM-Cは、事前のスペクトル情報がなくても大規模問題において信頼性高く固有値を計算できるか？

主な発見

• RIM-Cは、すべてのテストケースで少なくとも10桁の精度で固有値に収束することを示し、与えられた領域内の固有値を高い精度で計算可能である。
• 例1では、RIM-Cが0.185秒で10個の固有値を計算したのに対し、'eigs'は0.201秒を要し、同等の効率性を示した。
• 例2では、RIM-Cが0.558687秒で6個の固有値を計算したが、'eigs'は0.333513秒を要し、やや高いコストであったが一貫性のある結果を得た。
• 例3では、RIM-Cはシフトに依存せず2.57秒で計算を完了したのに対し、'eigs'の性能は0.39〜7.18秒と大きく変動し、RIM-Cの堅牢性が顕著に現れた。
• 例4では、RIM-Cが104.23秒で処理を完了したが、'eigs'は悪くないシフトを用いた場合でも1696.7秒を要し、RIM-Cの劣悪な条件下での安定性が明確に示された。
• 新しい指標は、4つの数値積分点しか使用しなくても、固有値を含む領域（例：$\delta_S \approx 0.997$）と含まない領域（例：$\delta_S \approx 0.0377$）を一貫して区別できた。