HPCI戦略プログラム分野5 » ユーザー支援事例報告 http://www.jicfus.jp/field5/jp Tue, 01 Mar 2016 05:48:22 +0000 ja hourly 1 http://wordpress.org/?v=4.2.7 原子核の見せかけの励起状態を求める ../../../160216houkoku/ ../../../160216houkoku/#comments Tue, 16 Feb 2016 06:01:44 +0000 ../../../?p=5491 続きを読む ]]> HPCI戦略プログラム分野5では素粒子・原子核・宇宙分野の研究者、計算機科学の研究者らからなるユーザー支援チームを組織し、分野内の情報共有、計算の並列化や最適化に関するユーザーへのアドバイス、ユーザー支援を行なっています。最新のユーザー支援事例報告をご紹介いたします。

計算の概略 原子核の見せかけの励起状態を求める
使用言語 Fortran90, MPI
問題 上記計算のため、十万次元くらいのエルミート行列を用いた連立線形方程式を解くプログラムがほしい。
解決法、参考情報 詳しくはこちらをご覧ください(PDF)

コメントや質問がありましたら、 jicfus-support-at-ccs.tsukuba.ac.jp (-at-を@にして下さい)までお送りください。

]]>
../../../160216houkoku/feed/ 0
磁気流体シミュレーションコードの空間解像度を部分的に挙げるために非一様格子を導入する ../../../150715houkoku/ ../../../150715houkoku/#comments Wed, 15 Jul 2015 06:39:20 +0000 ../../../?p=5263 続きを読む ]]> HPCI戦略プログラム分野5では素粒子・原子核・宇宙分野の研究者、計算機科学の研究者らからなるユーザー支援チームを組織し、分野内の情報共有、計算の並列化や最適化に関するユーザーへのアドバイス、ユーザー支援を行なっています。最新のユーザー支援事例報告をご紹介いたします。

計算の概略 磁気流体シミュレーションコードの空間解像度を部分的に挙げるために非一様格子を導入する
使用言語 FORTRAN
問題 2013年の千葉大学の磁気流体シミュレーションサマースクールで作成したHLLD法によるMHDコードへの、非一様メッシュの導入に伴う諸問題。
解決法、参考情報 詳しくはこちらをご覧ください(PDF)

コメントや質問がありましたら、 jicfus-support-at-ccs.tsukuba.ac.jp (-at-を@にして下さい)までお送りください。

]]>
../../../150715houkoku/feed/ 0
流束制限拡散近似を用いた輻射輸送方程式(線形拡散方程式+非線形方程式)を陰的に解く ../../../150519houkoku/ ../../../150519houkoku/#comments Tue, 19 May 2015 06:14:16 +0000 ../../../?p=5194 続きを読む ]]> HPCI戦略プログラム分野5では素粒子・原子核・宇宙分野の研究者、計算機科学の研究者らからなるユーザー支援チームを組織し、分野内の情報共有、計算の並列化や最適化に関するユーザーへのアドバイス、ユーザー支援を行なっています。最新のユーザー支援事例報告をご紹介いたします。

計算の概略 流束制限拡散近似を用いた輻射輸送方程式(線形拡散方程式+非線形方程式)を陰的に解く。
ニュートン・ラプソン法とマルチグリッド法(スムーザ=SOR 法)を用いている。
使用言語 FORTRAN
問題 MPI を用いて領域分割で並列化しているが、並列効率が悪いので、問題点を特定の上、改良したい。
解決法、参考情報 詳しくはこちらをご覧ください(PDF)

コメントや質問がありましたら、 jicfus-support-at-ccs.tsukuba.ac.jp (-at-を@にして下さい)までお送りください。

]]>
../../../150519houkoku/feed/ 0
スパコン上でC言語を用いた複素指数関数計算について(原子核の励起状態の計算) ../../../150511houkoku/ ../../../150511houkoku/#comments Mon, 11 May 2015 09:20:20 +0000 ../../../?p=5120 続きを読む ]]> HPCI戦略プログラム分野5では素粒子・原子核・宇宙分野の研究者、計算機科学の研究者らからなるユーザー支援チームを組織し、分野内の情報共有、計算の並列化や最適化に関するユーザーへのアドバイス、ユーザー支援を行なっています。最新のユーザー支援事例報告をご紹介いたします。

計算の概略 スパコン上でC言語を用いた複素指数関数計算について(原子核の励起状態の計算)
使用言語 c
問題 基研のSR16000 上のIBM 製xlc で複素数の計算を行いたいと考えています。しかし同等のFortran のコードと比較して、C のコードは4倍以上も遅くなると言う問題が生じています。原因を調べた結果、特に標準ライブラリlibm.a 中の複素数cexp がFortran と比べて異常なぐらい遅い事がわかりました。おそらくcexp を自分で置き換えたりするのが一般的な方法でしょうが、その場合は複素数cexp が提供されているその他のライブラリをご存知ないでしょうか?またこのような場合、他に何か解決策があるのでしょうか?同じくC をお使いの方が、このような問題が生じた時に、一般的にどのように対処されているのか、お知恵を拝借させてください。
解決法、参考情報 詳しくはこちらをご覧ください(PDF)

コメントや質問がありましたら、 jicfus-support-at-ccs.tsukuba.ac.jp (-at-を@にして下さい)までお送りください。

]]>
../../../150511houkoku/feed/ 0
メインプログラムからの MPI サブルーチンの呼び出し、行列積に関する MPI 化について ../../../150217houkoku_2/ ../../../150217houkoku_2/#comments Tue, 17 Mar 2015 12:17:03 +0000 ../../../?p=5014 続きを読む ]]> HPCI戦略プログラム分野5では素粒子・原子核・宇宙分野の研究者、計算機科学の研究者らからなるユーザー支援チームを組織し、分野内の情報共有、計算の並列化や最適化に関するユーザーへのアドバイス、ユーザー支援を行なっています。最新のユーザー支援事例報告をご紹介いたします。

計算の概略 メインプログラムからの MPI サブルーチンの呼び出し、行列積に関する MPI 化について
使用言語 FORTRAN
問題 bosonic BFSS 行列模型のコードで、MPI を用いないシリアルなコードからMPI を用いて並列化することを目指している。この際、行列の積A*B=C 演算部分が計算時間の大半を占めるため、この部分を並列化することを目指す。しかしながらそのMPI 化した行列積演算部分をサブルーチンとしてメインプログラムから呼び出そうと試行錯誤したがうまくいかない。MPI 版を正しく動作するようにしたい。
解決法、参考情報 詳しくはこちらをご覧ください(PDF)

コメントや質問がありましたら、 jicfus-support-at-ccs.tsukuba.ac.jp (-at-を@にして下さい)までお送りください。

]]>
../../../150217houkoku_2/feed/ 0
Lattice QCD における gauge fixing ../../../150217houkoku_1/ ../../../150217houkoku_1/#comments Tue, 17 Mar 2015 12:12:34 +0000 ../../../?p=5010 続きを読む ]]> HPCI戦略プログラム分野5では素粒子・原子核・宇宙分野の研究者、計算機科学の研究者らからなるユーザー支援チームを組織し、分野内の情報共有、計算の並列化や最適化に関するユーザーへのアドバイス、ユーザー支援を行なっています。最新のユーザー支援事例報告をご紹介いたします。

計算の概略 Lattice QCD における gauge fixing
使用言語
問題 Lattice QCD の計算を行っているのですが、現在、(Coulomb) gauge fixing を、lattice QCD
library の一つ、CPS (The Columbia Physics System) を用いて行っています。しかし、コードの実効
効率が非常に低いため (例えば 32^3×48 lattice を KEK SR16000 1node で計算する時の実効
効率は 2-3%程度)、もっと効率の良いコードを使いたいと思っています。お薦めのコードなどあれば紹介
していただけるでしょうか。
特に、GPGPU (HA-PACS)で高速な gauge fixing ができれば非常に助かるので、その辺りも調査し
ていただければと思っています。
よろしくお願いします。
解決法、参考情報 詳しくはこちらをご覧ください(PDF)

コメントや質問がありましたら、 jicfus-support-at-ccs.tsukuba.ac.jp (-at-を@にして下さい)までお送りください。

]]>
../../../150217houkoku_1/feed/ 0
コードのバージョン管理、メンバーでの情報共有について ../../../150217houkoku/ ../../../150217houkoku/#comments Tue, 17 Feb 2015 07:20:05 +0000 ../../../?p=4919 続きを読む ]]> HPCI戦略プログラム分野5では素粒子・原子核・宇宙分野の研究者、計算機科学の研究者らからなるユーザー支援チームを組織し、分野内の情報共有、計算の並列化や最適化に関するユーザーへのアドバイス、ユーザー支援を行なっています。最新のユーザー支援事例報告をご紹介いたします。

計算の概略 コードのバージョン管理、メンバーでの情報共有について
使用言語
問題 依頼内容 コード開発に関わるメンバーが多くなり、コードのバージョン管理を効率的に行い、研究や開発上の情報を後で参照/集約する方法を探している。またコードのモジュール化を進めて、複数のメンバーで開発やチューニングを行うための体制を作りたい。これまでの格子QCDコード開発では共同作業体制がうまくいっていると聞いているため、そこで得られたノウハウを教えて欲しい。
解決法、参考情報 詳しくはこちらをご覧ください(PDF)

コメントや質問がありましたら、 jicfus-support-at-ccs.tsukuba.ac.jp (-at-を@にして下さい)までお送りください。

]]>
../../../150217houkoku/feed/ 0
ユニタリ行列のオイラー角分解 ../../../150105houkoku/ ../../../150105houkoku/#comments Mon, 05 Jan 2015 02:55:15 +0000 ../../../?p=4836 続きを読む ]]> HPCI戦略プログラム分野5では素粒子・原子核・宇宙分野の研究者、計算機科学の研究者らからなるユーザー支援チームを組織し、分野内の情報共有、計算の並列化や最適化に関するユーザーへのアドバイス、ユーザー支援を行なっています。最新のユーザー支援事例報告をご紹介いたします。

計算の概略 ユニタリ行列のオイラー角分解
応用としては格子QCDコードのメモリ量/帯域ないし演算負荷の軽減を考えていますが、割と基礎的な内容の数学の質問となります
使用言語 c, c++, N=2,3ではMaximaを使ってEuler角から合成はしてみましたが、分解まではできていません
問題 UをU(N)ないしSU(N)のNxN複素行列とします。
このとき、次のような分解が可能であるように思えます。
U = ADB,ただしA, BはSO(N)の実回転行列、Dはexp(i diag(φ_1, φ_2,…,φ_N))となるような複素対角行列で、UがSU(N)のときはφ_iの総和は0とします。
A, D, B ∈ U(N) or SU(N)なのは明らかなのでこれらの積もU(N) or SU(N)になります。
自由度を数えるとA, BについてはN(N-1)/2、DについてはN or N-1で、合わせてN^2 or N^2 -1で、数合わせの上ではこれでよいように見えます。肝心の質問ですが
・このような分解は常に可能か、可能ならABの符号を除いて一意か?
・このような分解をするには、実際にどんな計算をすればよいのか?
です。
使用する計算機(予定も含む):「京」、その他B/Fの低い計算機全般。
解決法、参考情報 詳しくはこちらをご覧ください(PDF)

コメントや質問がありましたら、 jicfus-support-at-ccs.tsukuba.ac.jp (-at-を@にして下さい)までお送りください。

]]>
../../../150105houkoku/feed/ 0
小サイズの行列-行列積の高速化on SR16K ../../../141107houkoku/ ../../../141107houkoku/#comments Thu, 06 Nov 2014 07:50:03 +0000 ../../../?p=4746 続きを読む ]]> HPCI戦略プログラム分野5では素粒子・原子核・宇宙分野の研究者、計算機科学の研究者らからなるユーザー支援チームを組織し、分野内の情報共有、計算の並列化や最適化に関するユーザーへのアドバイス、ユーザー支援を行なっています。最新のユーザー支援事例報告をご紹介いたします。

計算の概略 随伴表現フェルミオンを含むラージNゲージ理論のシミュレーション
使用言語 FORTRAN, Fortran90
問題 随伴表現フェルミオンを含むラージNゲージ理論のシミュレーションをしています。そこでは、N×Nサイズの行列の行列-行列積を何度も計算します。
N=64の場合に、行列-行列積を高速に計算する方法が知りたいです。
あるいは、1000より小さいNで、計算効率が良いN
の値があるならば知りたいです。
解決法、参考情報 詳しくはこちらをご覧ください(PDF)

コメントや質問がありましたら、 jicfus-support-at-ccs.tsukuba.ac.jp (-at-を@にして下さい)までお送りください。

]]>
../../../141107houkoku/feed/ 0
QCDのシミュレーションにおいて、BlueGeneで大きなゲージ配位を書き出す ../../../141017houkoku/ ../../../141017houkoku/#comments Fri, 17 Oct 2014 09:08:48 +0000 ../../../?p=4699 続きを読む ]]> HPCI戦略プログラム分野5では素粒子・原子核・宇宙分野の研究者、計算機科学の研究者らからなるユーザー支援チームを組織し、分野内の情報共有、計算の並列化や最適化に関するユーザーへのアドバイス、ユーザー支援を行なっています。最新のユーザー支援事例報告をご紹介いたします。

計算の概略 QCDのシミュレーションにおいて、BlueGeneで大きなゲージ配位を書き出す
使用言語 Fortran
問題 高速化したい。配位が大きい(32^^3 *64)の時、IOに非常に時間がかかって困っている。
他の格子サイズのものと同じ型の出力ファイルにしたいので、中身を変えることは避けたい。
多分、MPIをうまく使ってやる方法があると思うのですがわからないので、具体的に書き方を教えて欲しい。
解決法、参考情報 詳しくはこちらをご覧ください(PDF)

コメントや質問がありましたら、 jicfus-support-at-ccs.tsukuba.ac.jp (-at-を@にして下さい)までお送りください。

]]>
../../../141017houkoku/feed/ 0