ユーザー支援報告
寺崎順
筑波大計算科学研究センター

１．何をやるのか
２．どういう人達がやるのか
３．依頼者にとって何がいいのか
４．どういう依頼が来たか、どう回答したか
５．ウェブページ
６．まとめ

2013年3月6日
東京

1. 何をやるのか？
• すでにお持ちの研究用数値計算プログラムの最適化、• これからつくろうとしているプログラムを将来の発展まで考えてどのような構造にすればよいか、またどのような手法を用いればよいか、• シリアルプログラムの並列化、さらにその大規模化等の問題あるいは作業過程に関して、ユーザー支援チームメンバーが個々の支援要請に応じて検討の上、助言・提案を行います。

2. どういう人達がやるのか？
ユーザー支援チームメンバーリスト １/２
橋本省二
入江敦子
松古栄夫
野秋淳一
花田政範
伊藤悦子
上田 悟
Guido Cossu
濱口信行
清水裕也
朴 泰祐
櫻井鉄也

高エネルギー加速器研究機構
高エネルギー加速器研究機構

教授、チームリーダー
研究支援員

高エネルギー加速器研究機構計算科学センター 助教高エネルギー加速器研究機構

特任助教

高エネルギー加速器研究機構
高エネルギー加速器研究機構
高エネルギー加速器研究機構
高エネルギー加速器研究機構

特任助教
特任助教
研究員
研究員

高エネルギー加速器研究機構
研究支援員
高エネルギー加速器研究機構
研究員
筑波大学計算科学研究センター 副センター長／教授
筑波大学大学院システム情報工学研究科教授

ユーザー支援チームメンバーリスト ２/２
吉江友照
根村英克
寺崎 順
多田野寛人
佐々木健志
今倉 暁
船木靖郎
松元亮治
松本洋介
武田隆顕
馬場淳一
市川隆敏
平松尚志
鎌野寛之

筑波大学計算科学研究センター 准教授
筑波大学計算科学研究センター 准教授
筑波大学計算科学研究センター 准教授、連絡係
筑波大学計算科学研究センター 助教
筑波大学計算科学研究センター 研究員
筑波大学計算科学研究センター 研究員
理化学研究所仁科加速器研究センター
協力研究員
千葉大学大学院理学研究科教授
千葉大学大学院理学研究科
特任助教
国立天文台
特任助教
東京工業大学大学院理学研究科 理学研究流動機構 特任助教
京都大学基礎物理学研究所
特任助教／特定研究員
京都大学基礎物理学研究所
特定研究員
大阪大学核物理研究センター特任助教

総勢26名、うち戦略プログラムによる被雇用人員19名。
自分の研究を行いながら、ユーザー支援を行っている。
専門分野
Lattice QCD
計算システム

人数11名
5名

宇宙
原子核

5名
4名

技術支援

1名

計算システムはアルゴリズムとプログラミングを含む

3. 依頼者にとって何がいいのか？
いろいろな専門の豊富な経験をもつ支援員が結集している― Fortran, C, C++, MPI, OpenMP, Cuda, IDL などをカバー、― 物理はもとより応用数学、可視化もカバー

ひとりまたは限られたグループでの議論では解けないプログラ
ミング上の問題が解決できることがある。解けなくても何らかの情報が得られる。
･ 用いるプラットフォームはPCから超並列計算機まで何でも可、･ 大規模計算や並列化に限らず、どんなに些細な問題でも扱う、・ 最先端のプログラムを開発中の研究者の方々はもとより院生の皆さんからの依頼も歓迎。

4. どんな依頼が来たか、どう回答したか？
抜粋を以下に示す。
問題1解決法・参考 上記サブルーチンの作成を外注し納入された。それをもとに完成形を情報の要約 作成中。
問題

25次元行列を用いる格子QCD計算のため、様々な線形ソルバの収束の可否やスピードなどを知りたい。

解決法・参考 上記データの収集を外注し、レポートが納入された。それをもとに公表
情報の要約 できる最終レポートを作成中。
問題

3複素対称行列の一般化固有値問題の並列サブルーチンが必要だがどこにもない。つくるにはどうすればよいか。

格子QCDの線形ソルバへの応用のためmultigrid法一般の考え方と具体的なやり方の解説を求む。

解決法・参考 この方法の一般的説明と格子QCDへの応用の二つの解説レポートを情報の要約 作成し依頼者に提出した。

問題
4

原因の可能性としてメモリアクセスなどを指摘し、プロファイリングの注
解決法・参考意点を述べ、改善方法として、メモリアクセスをできるだけ速くする方法、情報の要約ロードバランスの向上、リダクション処理の試みを提案した。
問題

5

宇宙ダストの低速衝突での合体、破壊、変形を調べるため, 付着力を入れた弾性球のN体計算コードをOpenMPを用いて並列化しているがシリアルコードに比べ実時間が短縮されない。

アスキーデータを高速に数値化する方法を知りたい。

ファイル読み込みをブロックごとに分け、次のブロックを読み込んでいる
解決法・参考間に以前に読み込んだデータの数値化を別スレッドで行うことで高速化情報の要約が行える。その詳細をレポートで示した。。
問題

３次元Magnetohydrodynamics シミュレーションの結果から、Interactive
data language (IDL)を用いて磁力線構造の可視化を行いたい。また、可能であれば、磁力線構造が変化する様子を示すムービーを作成したい。

6
IDLを用いてサポート担当者が開発し公開している3次元可視化ツール
解決法・参考及び磁力線描画ツールを、サンプルプログラムと共に渡し、使い方を説情報の要約明した。

問題
7

粒子と核子の相互作用が中心力ポテンシャルで与えられている場合の、ΛN散乱断面積の計算の手順を教えてほしい。

解決法・参考 重心系における Schrödinger 方程式の動径部分を差分化した式を示し、情報の要約 その解き方の概略と位相のずれ、ならびに断面積の計算式を示した。

問題
8低次元の超対称ゲージ理論の並列計算用コードで、質量シフトを伴う疎行列の逆行列を解くアルゴリズム”multi-mass solver”を要する。このためのMultimass shift BICGStab法より高速なアルゴリズムがあれば知りたい

解決法・参考 Multimass shift BICGStab法の特性を説明し、この方法の前に2,3の方法情報の要約 を試すことを推奨して、参考文献を紹介した。
問題
9

2Gbを超えるデータをCで一度にファイルに出力するにはどうすればよいか。

解決法・参考 fopen, fread, fwrite, fcloseを使えばよい。そうでないと割り込みの排除が情報の要約 保証されない。

2012年度 計１９案件、（2011年度 計１５案件）

依頼内容
応用数学
プログラミング技術

件数(2012年度)
７５件数(2011年度)
５４

物理の方法
可視化

４
２４０

ソフトウェアの情報

１２

依頼は、内容、分野、難易度において極めてさまざま。

依頼者の専門
宇宙素粒子
原子核
その他

件数(2012年度)
５９
５０件数(2011年度)
３
３５１

各案件が終了した後、回答をレポートにして、ウェブサイトにアップロード。
JICFuSレポート一覧• 今までに行われたユーザー支援事例の分類
I. 分野別１．Lattice-QCD
２．宇宙
III. 数学の手法
３．原子核
１．線型方程式
II. プログラミング技術による分類
２．微分方程式
１．MPI
３．他
２．Open MP
IV. 理論物理の方法
３．Fortran 90
V. ライブラリ、ソフトウェア
４．C, C++

ユーザー支援チームメンバーの計算経験レポート
プログラム名テンソルくりこみ群の数値シミュレーションプログラム

プログラムの目的

テンソルくりこみ群を使って分配関数を計算する使っている言語

Fortran

基本的なアルゴリズム使っている、または使える行列の掛け算，行列の特異値分解
Linux デスクトップ

プラットフォーム、
詳しくはこちらアーキテクチャプログラム名プログラムの目的

Bridge++(NBSWF)
ハイペロン相互作用を研究するための格子 QCD による４点相関関数を計算する使っている言語

C++, MPI

基本的なアルゴリズム

CG, FFT使っている、または使えるプラットフォーム、アーキテクチャ

MPI 環境で C++ コンパイラの使える Linux 計算機(x86_64)
を利用して作業を行った。
詳しくはこちら６. まとめ
ユーザー支援活動の定義、実行メンバー、今年度の依頼
内容と回答の抜粋、昨年度との統計的比較、ウェブサイトについて述べた。
今回強調したいこと:発足以来、分野、難易度、内容において極めてさまざまな３２案件を受け付け、お手上げだったことは一度もない。

ぜひ、引き続きユーザー支援をご利用下さい。