[論文レビュー] Real dqds for the nonsymmetric tridiagonal eigenvalue problem

Carla Ferreira, Beresford Ν. Parlett|arXiv (Cornell University)|Jan 24, 2012

Matrix Theory and Algorithms参考文献 18被引用数 3

ひとこと要約

この論文は、実数演算のみを用いて、実非対称三重対角行列の固有値を計算するための実数 dqds アルゴリズムを導入し、複素数シフトを三重ステップ変換により効率的に実装している。性能は複素数 dqds コードと同等であり、Ehrlich-Aberth 法に比べて最大100倍高速であり、要素の成長を制御し、動的変換拒否機構を用いることで安定性と高速性を両立している。

ABSTRACT

We present a new transform, triple dqds, to help to compute the eigenvalues of a real tridiagonal matrix C using real arithmetic. The algorithm uses the real dqds transform to shift by a real number and triple dqds to shift by a complex conjugate pair. We present what seems to be a new criteria for splitting the current pair L,U. The algorithm rejects any transform which suffers from excessive element growth and then tries a new transform. Our numerical tests show that the algorithm is about 100 times faster than the Ehrlich-Aberth method of D. A. Bini, L. Gemignani and F. Tisseur. Our code is comparable in performance to a complex dqds code and is sometimes 3 times faster.

研究の動機と目的

• 実数演算のみを用いて、実非対称三重対角行列の固有値を高速かつ高精度に計算するソルバーを開発すること。
• 複素数シフト（共役対を用いて）を実数演算で実装し、数値的安定性と効率性を保ちながら、複素数演算を回避すること。
• 要素の成長が著しく増加するステップを検出し、破棄する機構を設計し、アルゴリズムの堅牢性を高めること。
• 複素数 dqds アルゴリズムと同等の性能を達成しながら、実数演算の利点を維持すること。
• Ehrlich-Aberth や複素数 QR/dqds と同様の既存手法に対する実用的で高性能な代替手法を提供すること。

提案手法

• アルゴリズムは、実数 dqds と複素数シフト（共役対を用いて）を組み合わせた三重ステップ変換である tridqds を用い、実数演算で二重シフトを模倣する。
• 行列 $ Z_i $, $ ilde{L}_i^{-1} $, $ ilde{L}_i $, $ Y_i^{-1} $, $ Y_i $ を用いた類似変換の系列を適用し、状態変数 $ x_r, y_r, z_r, x_l, y_l $ を再帰的に更新する。
• 要素の成長は比 $ \hat{u}_i / \text{max}(u_i, l_i) $ を用いて監視され、成長がしきい値を超える場合は変換が拒否される。
• 現在の $ L,U $ 組に基づく動的シフト戦略を用い、収束を最大化しつつ不安定化を避けるようにシフトを選択する。
• MATLAB で実装され、収束した固有値を分離するためのデフラスティケーションおよびスプリッティング基準が含まれる。
• コア計算では、$ x, y, z $ 状態ベクトルの再帰的更新を用い、各ステップで正規化とスケーリングを実施して数値的安定性を維持する。

実験結果

リサーチクエスチョン

• RQ1 実数演算のみを用いたアルゴリズムは、非対称三重対角固有値問題において、複素数 dqds と同等の性能を達成できるか？
• RQ2 複素数演算を避けながら、実数演算で効果的に複素数シフトを適用する方法は何か？ その際に数値的不安定性が生じないか？
• RQ3 要素の成長が著しく増加する変換を信頼性高く検出し、拒否するための基準は何か？
• RQ4 三重ステップ tridqds 変換は、Ehrlich-Aberth や複素数 dqds といった既存手法と比較して、速度と精度の面でどのように差がつくか？
• RQ5 複素数演算の計算コストを回避しつつ、高い相対的精度を維持できるか？

主な発見

• 提案された実数 tridqds アルゴリズムは、Bini, Gemignani, Tisseur が開発した Ehrlich-Aberth 法に比べて約100倍高速である。
• アルゴリズムの性能は複素数 dqds コードと同等であり、一部のケースでは最大3倍も速い。
• 要素の成長が著しく増加する変換の動的拒否は、数値的安定性と収束の信頼性を顕著に向上させた。
• アルゴリズムは、固有値が複素共役対であっても、実数演算のみで実非対称三重対角行列のすべての固有値を正しく計算できた。
• スプリッティングおよびデフラスティケーション戦略は、収束した固有値を効果的に分離でき、全スペクトルの効率的計算を可能にした。
• 本手法は高い相対的精度を維持しており、特異値計算における dqds の既知の利点を非対称固有値問題へと拡張した。