[論文レビュー] On a weighted quasi-residual minimization strategy of the QMR method for solving complex symmetric shifted linear systems

Tomohiro Sogabe, Takeo Hoshi+2 more|ArXiv.org|Feb 16, 2009

Matrix Theory and Algorithms被引用数 4

ひとこと要約

本稿では、電子構造シミュレーションに生じる大規模な複素対称シフト付き線形系を解くために、QMR_SYM法に重み付き準最小残差（WQMR）戦略を提案する。最適化された重み行列を導入することで、近似解の更新にかかる計算コストを顕著に低減し、シード系の選択を必要としないものの、シフトCOCG法と同等の性能を達成する。

ABSTRACT

We consider the solution of complex symmetric shifted linear systems. Such systems arise in large-scale electronic structure simulations and there is a strong need for the fast solution of the systems. With the aim of solving the systems efficiently, we consider a special case of the QMR method for non-Hermitian shifted linear systems and propose its weighted quasi-minimal residual approach. A numerical algorithm, referred to as shifted QMR\_SYM($B$), is given by the choice of a particularly cost-effective weight. Numerical examples are presented to show the performance of the shifted QMR\_SYM($B$) method.

研究の動機と目的

• 大規模な数のシフト付き線形系を解く際、シフトQMR_SYM法における近似解の更新処理にかかる計算コストの高さに対処すること。
• 収束失敗を避けるために注意深くシード系を選択する必要がある、シフトCOCG法の制限を克服すること。
• Krylov部分空間のシフト不変性特性を活用し、複数のシフト系間で基底ベクトルを再利用すること。
• 残差の更新コストを最小限に抑えつつ収束特性を損なわない、費用対効果の高い重み戦略を開発すること。
• 電子構造計算における複素対称シフト付き系に対して、既存のKrylov法の代替として、強固で効率的な手法を提供すること。

提案手法

• 非エルミートなシフト付き系に特化したQMR法を、複素対称系に限定することで、シフトQMR_SYM法を導出する。
• 数が多数のシフト付き線形系を解く際、近似解の更新処理が、シフトQMR_SYM法において最も計算コストがかかる部分であることを特定する。
• 特定の重み行列を導入することで、近似解の更新コストを低減する重み付き準最小残差（WQMR）アプローチを提案する。
• 実数演算を可能にする費用対効果の高い重み行列Bを選択し、浮動小数点演算数を削減する。
• WQMR戦略をQMR_SYMフレームワークに統合し、収束性を維持する新しいアルゴリズム、すなわちシフトQMR_SYM(B)を構築する。
• Krylov部分空間のシフト不変性を利用することで、シフトCOCG法とは異なり、シード系の選択を必要としない。

実験結果

リサーチクエスチョン

• RQ1 大規模な数のシフト付き線形系を解く際、シフトQMR_SYM法における主要な計算ボトルネックは何か？
• RQ2 重み付き準最小残差戦略は、収束性を損なわせることなく、近似解の更新コストを低減できるか？
• RQ3 計算時間と頑健性の観点から、提案されたシフトQMR_SYM(B)法はシフトCOCG法と比較してどの程度の性能を示すか？
• RQ4 新しい手法は、シフトCOCG法の主要な制限要因であるシード系選択を不要にするか？
• RQ5 提案手法は、シフトCOCG法と同等の効率性を達成できるか、また、より広範な複素対称シフト付き系に適用可能か？

主な発見

• シフトQMR_SYM法は、近似解の更新処理に複素数演算が使用されるため、シフトCOCG法に比べて反復ごとの計算コストが顕著に高い。
• シフト系の数が閾値（例：m > 200）を超えると、近似解の更新処理がシフトQMR_SYM法の最も計算コストがかかる部分となる。
• 実数重み戦略を用いることで、提案されたシフトQMR_SYM(B)法は、計算コストをシフトCOCG法と同等の水準にまで低減する。
• m = 1501の場合、シフトQMR_SYM(B)はシフトCOCG法と比較して計算時間比が約1.0に達し、同等の性能を示す。
• シード系の選択を必要とせず、シフトCOCG法とは異なり頑健性を維持する。
• 大規模な電子構造問題に対する数値実験により、シフトQMR_SYM(B)が多くのシフトを持つ系に対して実用的で効率的な代替手段であることが確認された。