1 / 15

高性能並列固有値計算アルゴリズムの開発

高性能並列固有値計算アルゴリズムの開発. コンピューティクス による物質デザイン：複合相関と非平衡 ダイナミクス. 櫻井鉄也 1 ), 2) , 二村保 徳 1) ， 前田恭 行 1) ， 矢野 貴大 1) ， Du Lei 1 ), 2) ， 今倉 暁 1) ， 多田野 寛人 1) 1) 筑波大学大学院システム情報工学研究科 2) 科学技術振興機構 CREST. 研究課題. 研究課題名 計算量子科学のニーズに対応した高性能固有値計算アルゴリズムの開発 体制 代表者　　　　櫻井鉄也（筑波大学） 連携研究者　多田野寛人，今倉暁（筑波大学）

teranika-fullerton
Télécharger la présentation

高性能並列固有値計算アルゴリズムの開発

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. 高性能並列固有値計算アルゴリズムの開発 コンピューティクスによる物質デザイン：複合相関と非平衡ダイナミクス 櫻井鉄也1), 2), 二村保徳1)，前田恭行1)，矢野貴大1)，Du Lei1), 2)，今倉暁1)，多田野寛人1)1)筑波大学大学院システム情報工学研究科2)科学技術振興機構CREST

  2. 研究課題 • 研究課題名 • 計算量子科学のニーズに対応した高性能固有値計算アルゴリズムの開発 • 体制 • 代表者　　　　櫻井鉄也（筑波大学） • 連携研究者　多田野寛人，今倉暁（筑波大学） • 研究協力者　Du Lei，二村保徳，前田恭行，矢野貴大（筑波大学） • 主な実施課題 • 階層型固有値解法によるバンド構造計算法の高性能化 • 確率的推定法を利用した状態密度計算法の開発 • 計算量子科学分野で現れるさまざまな固有値計算への展開 • 実問題での性能評価とボトルネックの改善

  3. 発表内容 • はじめに • 研究背景と目的 • 超並列固有値解析エンジン • CRESTプロジェクトによるソフトウェア開発 • 周回積分による疎行列向け固有値解法 • 解法の概要と実行例 • 提案法が有効な固有値問題 • おわりに

  4. はじめに アプリ B アプリ C アプリ A アプリ X ・・・ 固有値解析エンジン アプリケーションを支える基盤的技術 • 研究背景と目的 • 固有値問題が現れるアプリケーションの高性能化 • 本課題では特に計算量子科学のニーズに対応 • ペタスケール超規模のハードウェア性能を引き出す並列固有値計算アルゴリズムの研究とその高性能実装技術を開発する

  5. はじめに Cray XT4 (LBNL)における計算時間 時間(秒) 提案法 (SS法*) 2500 従来法 (Shift-invert Lanczos法) 1000 500 並列度が低いとき 並列度が高いとき 100 コア数 *Sakurai, et al. 2003 8 32 128 512 2048 • 大規模並列で性能を発揮する解法 • 数値例：高エネルギー加速器設計で現れる固有値問題(110万次元) • SS法*，Lanczos法ともに，内部で同じ線形ソルバーを使用

  6. CRESTプロジェクト超並列固有値解析エンジンの開発CRESTプロジェクト超並列固有値解析エンジンの開発 • ポストペタスケール計算環境でのハードウェア性能を十分に引き出す超並列固有値解析エンジンを開発 (H23年度〜H27年度) • 主要メンバー • 櫻井(筑波大)，張(名古屋大)，今村(理研)，山本(神戸大)，藏増(筑波大)，星(鳥取大) • アーキテクチャの超並列性を活用する階層型並列アルゴリズム • 密行列・疎行列に対応 • 共通基盤技術としての線形計算手法・実装技術開発 • 高性能な線形計算手法，性能予測モデルによる評価，GPU／メニーコア向け実装 • 実問題での利用技術の開発 • 基礎科学・ナノ物質分野での実問題に適用し，高度利用技術を開発

  7. 周回積分を用いた固有値解法 Im Re 基本アイディア：有理式に対する周回積分 行列Aのレゾルベント: 固有値, : 固有ベクトル(ここでは簡単のためAは実対称とし，固有値は相異なるとする) 周回積分によるスペクトル射影

  8. 周回積分を用いた固有値解法 • 疎行列向け超並列固有値解法の概略 • 複素平面上の特定領域内部の固有ベクトル成分を抽出 • 周回積分は数値積分で近似（各積分点では線形方程式を解く）

  9. 領域内固有値数の確率的推定法 GEP: Futamura, et al. (2010), NEP: Maeda et al. (2011) • パラメータの最適化のために確率的固有値数推定を用いる • G 内の固有値数： • 逆行列のトレースを以下のように近似するここでベクトルの要素は1 または-1 を等確率にとる • この推定値を利用して必要な部分空間サイズを決める

  10. 計算例 Futamura, et al. (2013) • 実空間密度汎関数法(RSDFT)によるバンド図計算 • シリコンナノワイヤー9,924原子 • 行列サイズ：8,719,488次元，利用リソース：「京」 27,648コア • バンドギャップ付近の固有値を計算

  11. 開発ソフトウェア Eigen-Supercomputing Sparse Engine Application Interface Eigensolver Block Linear Solver Base Modules: PETSc, Trilinos, MUMPS BLAS, SpBLAS, LAPACK, ScaLAPACK MatVec, Data Converter • 疎行列向けモジュール • 固有値ソルバ，線形ソルバ，基本線形計算ライブラリ，アプリケーションインタフェイス • 開発ソフトウェア • MATLAB版：sseig.m • Webからダウンロード可 • 多項式，非線型版も利用可能 • SLEPc版 (C言語)：eps_ciss • SLEPc ver. 3.4 に実装 • Fortran90版： • 京向けを近日公開予定RSDFTのバンド図計算で利用 • Trilinos版 (C++)：開発中

  12. GPU/メニーコアへの対応 GPUクラスタ HA-PACS (Tsukuba)での実験 Yano, et al. accepted • 密実対称行列の一般化固有値問題 • 10,000次元，最小固有値から10％の固有対を計算 • BLAS/LAPACKルーチンとしてMAGMAを利用 • miscの大部分(SVD, Rayleigh-Ritz）は，現状では単一GPU or CPU実装

  13. 提案法が有効な問題 Im Re 指定した範囲 複素平面上の指定領域 Im 標準固有値問題 一般化固有値問題 Re 多項式固有値問題 非線形固有値問題 リング状領域 （ex. 平方根，有理式，etc.） 固有値を求める領域求めたい範囲やその範囲内にある固有値数が　　ある程度予想できる問題に対して効果が高い 固有値問題の種類

  14. おわりに • 大規模並列で性能を発揮する解法の開発 • 高性能固有値解析ソフトウェアを開発 • 京版は近日中に公開予定，今後はポストペタへ • アルゴリズムレベルでの階層構造 • 周回積分により固有成分をフィルタリング • 独立な複数の線形方程式に帰着 • 実アプリケーション • RSDFTのバンド図計算 • 今後の課題 • 実用性を考慮したアプリケーション向けインターフェイスの開発 • 計算量子科学分野で現れるさまざまな固有値計算への展開

  15. References Y. Futamura, T. Sakurai, S. Furuya and J.-I. Iwata, Proc. VECPAR 2012, LNCS 7851, 226–235 (2013). T. Mizusaki, K. Kaneko, M. Honma, and T. Sakurai, Phys. Rev. C82, 024310 (2010). Y. Nagai, Y. Shinohara, Y. Futamura, Y. Ota and T. Sakurai, arXiv:1303.3683 (2013). H. Ohno, Y. Kuramashi, H. Tadano and T. Sakurai, JSIAM Letters2, 115-118 (2010). T. Sakurai and H. Sugiura, J. Comput. Appl. Math.159, 119-128 (2003). T. Sakurai, H. Tadano, T. Ikegami and U. Nagashima, Taiwanese J. Math.14, 855-867 (2010). T. Sakurai, H. Tadano and Y. Futamura, J. Alg. Comput. Tech.7 (in print). T. Yano, Y. Futamura and T. Sakurai, Proc. 3PGCIC2013 (accepted).

More Related