演算加速器を有する計算機クラスタ向け固有値解法の並列実装に関する研究

具有计算加速器的计算机集群特征值求解方法的并行实现研究

• 批准号: 16J02020
• 负责人: 矢野 貴大
• 金额: $ 0.83万
• 依托单位: University of Tsukuba
• 依托单位国家: 日本
• 项目类别: Grant-in-Aid for JSPS Fellows
• 财政年份: 2016
• 资助国家: 日本
• 起止时间: 2016-04-22 至 2018-03-31
• 项目状态: 已结题
• 来源: https://kaken.nii.ac.jp/grant/KAKENHI-PROJECT-16J02020/
• 关键词: 固有値問題; 固有値計算; Sakurai-Sugiura法; 並列計算; Krylov部分空間反復法; 大規模計算

本年度は研究課題3.を達成するために，ナノシミュレーションにおいて現れる大規模な密エルミート一般化固有値問題に対し，大規模Xeon Phi(Knights Landing,KNL)クラスタOakforest-PACSにおける周回積分型固有値解法の適用，および性能評価を行うことを目的に研究を行った．研究課題1.について，前年度の成果であるGPUクラスタ向け分散メモリ並列な行列-複数右辺ベクトル積のKNLクラスタ向け移植を行った．また，複数右辺ベクトルのデータ保持形式を使用して，直交化手法であるCholesky QR法の実装を行った．これらを利用してKrylov部分空間反復法の一種であるBlock BiCGrQ法の実装を行った．分散メモリ並列数値計算パッケージScaLAPACKの直接法による連立一次方程式ソルバであるPZGETRFおよびPZGETRSと提案法を，周回積分型固有値解法で現れる連立一次方程式に対して適用した．Oakforest-PACS上で性能評価を行い，特に多数ノードでの実験において提案法が高速であることを確認し，学会で口頭発表を行った．研究課題3.について，Sakurai-Sugiura法のパッケージz-Paresに前述のBlock BiCGrQ法を組み込むことで，KNLクラスタにおける分散メモリ並列な一般化固有値問題解法を実現した．ナノシミュレーションにおいて実際に現れる行列に対して，ScaLAPACKのエルミート一般化固有値問題向けソルバであるPZHEGVXと提案法を適用し，性能評価を行った．PZHEGVXの多数ノードでのスケーラビリティを確認し，少数ノードでの提案法の実行結果から多数ノードでの計算時間の見積もり，多数ノード時にはPZHEGVXよりも高速になる可能性が示された．この成果に関して，国際会議SIAM PP18にて口頭発表を行った．

は research topic this year 3. を reached す る た め に, ナ ノ シ ミ ュ レ ー シ ョ ン に お い て now れ る large-scale な dense エ ル ミ ー ト generalized numerical problems inherent に し polices, Large-scale Xeon Phi (Knights Landing, KNL) ク ラ ス タ Oakforest - PACS に お け る weeks back to integral type の applicable, inherent numerical method お よ び performance evaluation 価 を line う こ と を purpose に を line っ た. Research topic 1. に つ い て, former annual の results で あ る GPU ク ラ ス タ メ け decentralized モ リ parallel な ranks - the plural right 辺 ベ ク ト ル product の KNL ク ラ ス タ line to け transplantation を っ た. ま た, plural right 辺 ベ ク ト ル の デ ー タ keep form を use し て, straight minimization algorithm technique で あ る Cholesky QR method の line be loaded を っ た. <s:1> れらを by using the <s:1> てKrylov partial space repetition method <s:1> a であるBlock BiCGrQ method を row った. Scattered メ モ リ tied for the numerical calculation パ ッ ケ ー ジ ScaLAPACK の direct method に よ る even set a formula ソ ル バ で あ る PZGETRF お よ び PZGETRS と を proposal method, inherent numerical solution で weeks back to integral type now れ る even set a formula に し seaborne て applicable し た. Oakforest - で performance evaluation on PACS 価 を い, most especially に ノ ー ド で の be 験 に お い て proposal method が high-speed で あ る こ と を confirm し, learn で oral 発 table line を っ た. Research topic 3. に つ い て, Sakurai - Sugiura method の パ ッ ケ ー ジ z - Pares に aforesaid の Block BiCGrQ を group み 込 む こ と で, KNL ク ラ ス タ に お け る scattered メ モ を リ tied for numerical problems inherent な generalization method be presently し た. ナ ノ シ ミ ュ レ ー シ ョ ン に お い て be interstate に now れ る ranks に し seaborne て, ScaLAPACK の エ ル ミ ー ト generalized numerical problems inherent to the け ソ ル バ で あ る PZHEGVX と を applicable し proposal method, performance evaluation of 価 を line っ た. PZHEGVX の most ノ ー ド で の ス ケ ー ラ ビ リ テ ィ を confirm し, minority ノ ー ド で の の proposal method be line results か ら most ノ ー ド で の computing time の see product も り, most ノ ー ド when に は PZHEGVX よ り も high-speed に な が る possibility in さ れ た. <s:1> results に related to て, international conference SIAM PP18にて oral presentation を line った.

项目成果

密行列一般化シフト線形方程式向けKrylov部分空間反復法の分散並列実装と性能評価

稠密矩阵广义平移线性方程Krylov子空间迭代方法的分布式并行实现及性能评估

• 发表时间: 2017
• 作者: Md. Imran Hossain;Shinya Hanashima;Takuto Nomura;Sebastien Lethu;Hiroshi Tsuchikawa;Michio Murata;Hiroki Kusaka;Shunsuke Kita;Katsumi Maenaka;矢野貴大，二村保徳，櫻井鉄也;Md. Imran Hossain;矢野貴大，二村保徳，櫻井鉄也
• 通讯作者: 矢野貴大，二村保徳，櫻井鉄也

A real-valued block conjugate gradient type method for solving complex symmetric linear systems with multiple right-hand sides

求解多右边复对称线性系统的实值分块共轭梯度型方法

• DOI: 10.21136/am.2017.0023-17
• 发表时间: 2017
• 期刊: Applications of Mathematics
• 影响因子: 0.7
• 作者: Yasunori Futamura;Takahiro Yano;Akira Imakura;Tetsuya Sakurai
• 通讯作者: Tetsuya Sakurai

Distributed Parallel Implementation of the Sakurai-Sugiura Method for Large Dense Generalized Eigenvalue Problems

大稠密广义特征值问题 Sakurai-Sugiura 方法的分布式并行实现

• 发表时间: 2018
• 作者: Takahiro Yano;Yasunori Futamura;Tetsuya Sakurai
• 通讯作者: Tetsuya Sakurai

密行列に対する周回積分型固有値解法のGPUクラスタにおける性能評価

稠密矩阵循环积分型特征值求解方法GPU集群性能评估

• 发表时间: 2016
• 作者: Md. Imran Hossain;Shinya Hanashima;Takuto Nomura;Sebastien Lethu;Hiroshi Tsuchikawa;Michio Murata;Hiroki Kusaka;Shunsuke Kita;Katsumi Maenaka;矢野貴大，二村保徳，櫻井鉄也;Md. Imran Hossain;矢野貴大，二村保徳，櫻井鉄也;矢野貴大，二村保徳，櫻井鉄也
• 通讯作者: 矢野貴大，二村保徳，櫻井鉄也

密一般化固有値問題向け周回積分型固有値解法における連立一次方程式解法の分散並列実装手法

稠密广义特征值问题圆积分型特征值求解方法中求解联立线性方程组的分布式并行实现方法

• 发表时间: 2017
• 作者: Md. Imran Hossain;Shinya Hanashima;Takuto Nomura;Sebastien Lethu;Hiroshi Tsuchikawa;Michio Murata;Hiroki Kusaka;Shunsuke Kita;Katsumi Maenaka;矢野貴大，二村保徳，櫻井鉄也
• 通讯作者: 矢野貴大，二村保徳，櫻井鉄也