分散メモリ型スパースソルバの开发と评似 - 计算基盘研究室.pptVIP

応用数理工学特論線形計算とハイパフォーマンスコンピューティング 第8回 名古屋大学　計算理工学専攻 山本有作 １.　はじめに 本研究で対象とする問題 標準固有値問題　Ax = lx A ： n×n 非対称密行列 応用分野 MHD 化学工学 量子化学 流体力学 Cf. Bai and Demmel: A test matrix collection for non-Hermitian eigenvalue problems. 非対称行列の固有値計算の流れ 各部分の実行時間 全固有値を求める場合の演算量 ヘッセンベルグ 化：　(10/3) n3 QR法： 　10n3 （経験値） 演算量（時間）の大部分をQR法が占める。 QR法の高速化が必要 QR法の特徴と高性能計算 QR法の基礎 アルゴリズム（Francis, 1961 Kublanovskaya, 1961） 行列 A0 から始めて次のように QR 分解と相似変換を繰り返す。 A0 = Q1R1 A1 = R1Q1 (= Q1–1A0 Q1 ) A1 = Q2R2 A2 = R2Q2 (= Q2–1A1 Q2 = Q2–1Q1–1A0 Q1Q2 ) 収束定理 適当な条件の下で，Ak は（ブロック）上三角行列に収束 A0 の固有値を絶対値の大きい順に l1, l2, … , ln とすると，対角要素 aij