phononの分散関係の計算 －カイラルナノチューブ(18,3)－

Presentation on theme: "phononの分散関係の計算 －カイラルナノチューブ(18,3)－"— Presentation transcript:

1 phononの分散関係の計算 －カイラルナノチューブ(18,3)－
　－カイラルナノチューブ(18,3)－ 齋藤グループ 小林 直樹

2 プログラムの実行 (1) 座標計算 プログラム tube-xyz1.f (2) 最近接情報データ プログラム saikin1.f
この計算に 非常に時間が 掛かる (2) 最近接情報データ 　　プログラム saikin1.f (3) tube の phonon の分散関係 　　プログラム tu-phonon1.f tube　 　Pentium4　 4時間28分 Xeon 時間38分 wire　 　Alpha 時間24分 計算時間 tube GNU Fortran 77 wire DIGITAL Fortran 77 コンパイラ

3 MPIへの対応① －新しく追加－ include 'mpif.h' integer myrank, tnode, error
include ‘nk-size’ の下に追加 integer myrank, tnode, error include ‘mpif.h’ の下に追加 call mpi_init(error) call mpi_comm_rank(MPI_COMM_WORLD, myrank, error) call mpi_comm_size(MPI_COMM_WORLD, tnode, error) 変数宣言の下に追加 call mpi_finalize(error) メインプログラムの stop の上に追加

4 MPIへの対応② －編集－ do 20 i=1,nj lb = nj/tnode lc = mod(nj,tnode)
if(myrank+1.le.lc) then is = myrank*lb+myrank+1 ie = is + lb else is = myrank*lb+lc+1 ie = is + lb - 1 end if do 20 i=is,ie 計算回数をプロセス数に分割 doループで計算開始 Fortranソースファイル

5 MPI実行時のエラー　－２プロセス－ T = do 20 i,nj = do 20 i,nj = back deigch nnv= e(j)<0 ,i , j = D end of do 20 do 20 i,nj = do 20 i,nj = e(j)<0 ,i , j = D do 20 i,nj =

6 MPI実行時のエラー　－４プロセス－ T = do 20 i,nj = do 20 i,nj = do 20 i,nj = do 20 i,nj = back deigch nnv= e(j)<0 ,i , j = D end of do 20 do 20 i,nj =

7 エラーの比較 FORTRAN77 mpif 77 での 計算結果
e(j)<0 ,i , j = E e(j)<0 ,i , j = E e(j)<0 ,i , j = E e(j)<0 ,i , j = e(j)<0 ,i , j = e(j)<0 ,i , j = e(j)<0 ,i , j = e(j)<0 ,i , j = e(j)<0 ,i , j = e(j)<0 ,i , j = e(j)<0 ,i , j = e(j)<0 ,i , j = FORTRAN77 での計算結果 mpif 77 での 計算結果 e(j)<0 ,i , j = D e(j)<0 ,i , j = D e(j)<0 ,i , j = D e(j)<0 ,i , j = D e(j)<0 ,i , j = D e(j)<0 ,i , j = D e(j)<0 ,i , j = D e(j)<0 ,i , j = D e(j)<0 ,i , j = D

8 プロセス数と計算時間の関係 np h m s 1 2 11 43 7 15 3 45 4 33 19 5 26 42 6 23 8 18 56 9 22

9 問題点 解決法及び今後の課題 各プロセスで計算 出力ファイルが各プロセスに分割 正しい出力ファイルが得られない
サブルーチン　『MPI_GATHER』を用いてデータの結合 (サブルーチン機能 ： 全プロセスから1プロセスにメッセージ送信)