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
Hitachi SR16000での高速化 ../../../../141211hpcdoor-2/ ../../../../141211hpcdoor-2/#comments Thu, 11 Dec 2014 10:11:46 +0000 ../../../../?p=4820 続きを読む ]]> 京コンピュータで計算を行ったり、京コンピュータにコードを移植した際にユーザの方が気づいたことを参考情報としてまとめました。

Hitachi SR16000を利用する際の注意点(PDF)
土井琢身氏(理化学研究所・HPCI戦略プログラム分野5課題1)
行った計算: 格子QCDによる核力の計算
使用言語、コンパイラ: C/C++, IBM XL C/C++ compiler for AIX

HPCI戦略プログラム分野5ではユーザー支援の一環として、最適なアルゴリズムや手法を探すための指針となるような情報をJICFuS WIKIに蓄積しています。項目については高性能計算の扉をご覧ください。

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

]]>
../../../../141211hpcdoor-2/feed/ 0
京コンピュータ上での高速化 ../../../../141211hpcdoor-1/ ../../../../141211hpcdoor-1/#comments Thu, 11 Dec 2014 10:11:39 +0000 ../../../../?p=4819 続きを読む ]]> 京コンピュータで計算を行ったり、京コンピュータにコードを移植した際にユーザの方が気づいたことを参考情報としてまとめました。

京/FX10を利用する際の注意点(PDF)
土井琢身氏 (理化学研究所・HPCI戦略プログラム分野5課題1)
行った計算: 格子QCDによる核力の計算
使用言語、コンパイラ: C/C++, Fujitsu C/C++ compiler

HPCI戦略プログラム分野5ではユーザー支援の一環として、最適なアルゴリズムや手法を探すための指針となるような情報をJICFuS WIKIに蓄積しています。項目については高性能計算の扉をご覧ください。

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

]]>
../../../../141211hpcdoor-1/feed/ 0