Lattice Tool Kit in f90 Project在F90的项目格工具包.pptVIP

Lattice Tool Kit in f90 Project 野中千穂、中村純 ２０００年１０月１２日 航空宇宙技術研究所 動機と目的 格子は場の量子論を研究するための道具。 得にハドロン物理の研究に有用 格子QCDの基本的な部分の分かりやすいコードを書く それを基に、数値計算の非専門科でも、自分の目的に合ったプログラムが書けるように 「モジュール」機能を持ったFortran90はC＋＋と同様のデータ構造記述力 プロジェクトメンバー S.Choe (崔勝皓) T.Hatou (羽藤隆夫) S.Muroya (室谷心) C.Nonaka (野中千穂) A.Nakamura (中村純) I.Pushkina (イリーナ?プシキーナ) 用意する予定のプログラム Metropolis and Pseudo-Heatbath gauge update Standard Wilson Action (1x1+1x2) Imroved Action Anisotropic Lattice Hybrid Monte Carlo for dynamical fermions (Wilson Fermions) Quark Propagators Wilson Fermion + Clover MPIで並列化 コードは white box モジュール:データ構造＋演算規則 c-----------------------------c program ex1 c-----------------------------c use vec_def TYPE(vector) p, q, r p = vector(0.5, -1.5); q = vector(0.5, 0.5) r = p + q write(*,*) r end 格子QCD(1) プラケット 格子QCD（２）クォークプロパゲータ　 格子QCD（３）改良された作用 並列化　のりしろを付ける * * MODULE v e c_def TYPE v e c t o r r e a l x, y END TYPE INTERFACE OPERATOR(+) MODULE PROCEDURE vadd END INTERFACE CONTAINS FUNCTION vadd(a,b) RESULT(c) TYPE(vector), INTENT(IN):: a, b TYPE(vector) c c %x = a%x + b%x c %y = a%y + b%y END FUNCTION END MODULE : SU(3)行列 t emp1 = u(nu) t emp2 = nu.gshift.u(mu) t emp3 = temp1 * temp2 t emp1 = mu.gshift.u(nu) staple = staple + (temp3.prod1.temp1) TYPE(g_field1) staple, temp1, temp2, temp3 TYPE g_field1 SEQUENCE COMPLEX*16, DIMENSION(NC,NC,NV/NBUSH) :: g INTEGER parity, direction END TYPE parameter( NC=3 ) TYPE g_field0 SEQUENCE COMPLEX*16, DIMENSION (NC,NC,-NDW+1:NX+NDW,-NDW+1:NY+NDW, -NDW+1:NZ+NDW,-NDW+1:NT+NDW) :: g INTEGER parity, direction END TYPE c---------------------------------------------------------------