研究課題
平成26年度に繰越して2分法と逆反復法による大規模スパース行列の部分特異対計算法について研究した。前年度までに再直交化つきGolub-Kahan-Lanczos(GKL)の前処理アルゴリズムによる部分特異対を与える近似行列生成法を開発しているが、平成27年度は、まず、GKLアルゴリズムを共有メモリ並列計算機(マルチコアプロセッサ)で実行するにあたって、ボトルネックであった再直交化プロセスの高速化を目指し、キャッシュメモリの効果的な利用によってデータ再利用性の高めた。大きい方または小さい方からいくつかの特異値と対応する特異ベクトルのみが必要な場合には、特異値分解を正定値な大規模スパース行列の固有値分解に変換した上で米国のライブラリARPACK(ARnoldi PACKage)を利用することができるが、部分固有対を与える近似行列計算のGKLアルゴリズムでは大規模行列の乗算が頻出し、キャッシュメモリのオーバーフローが起きやすい。石上裕之氏、高田雅美氏と共同提案する実装法では、行列の行ベクトルのデータ再利用によってキャッシュヒット率を高めることで、ARPACKの計算時間を80%に削減することができた。さらに、GKLアルゴリズムと同様に前処理法となるべき古典Gram-Schmidt法、compact WY法等を比較検討し、より大規模なスパース行列の部分特異対計算法では、高速な前処理にはthick restart GKLアルゴリズムを採用すべしという今後の研究の方針が得られるに至った。
27年度が最終年度であるため、記入しない。
Our laboratory continually develops numerical computation package LAPIS (Linear Algebra Package by Integrable Systems) which is based on integrable systems.
すべて 2016 2015 その他
すべて 雑誌論文 (2件) (うち査読あり 2件、 オープンアクセス 2件) 学会発表 (7件) (うち国際学会 2件、 招待講演 1件) 備考 (1件)
Proceedings of the 10th International Conference on P2P, Parallel, Grid, Cloud and Internet Computing (3PGCIC-2015)
巻: 1 ページ: 216, 223
Proceedings of The 2015 International Conference on Parallel and Distributed Processing Techniques and Applications (PDPTA2015)
巻: 2 ページ: 243, 249
http://www-is.amp.i.kyoto-u.ac.jp/lab/en/isvd/