プログラム構造(サブルーチン呼び出しリスト、各論編用)


以下に、公開が可能と考えられるサブルーチン部分を一挙に公開したいと思 います(10/10、1997版)。(全てのコードが筆者の完全なオリジナ ルではないところ〔他の共同研究者による部分〕もありますが、一応許可を得 られたもののみ公開しています。もし問題があれば連絡して下さい。)

公開されていないものには、元々公開が不可能なものもあります、ベッセル 関数や、対角化、時間表示などのサブルーチンは、非公開の汎用ライブラリー 等のものを利用しているので、これらに関しては各自で確保してもらうしかな いと思います。



  文字色が青色っぽいのが既に公開されているサブルーチンです。

 +-main-----TIME------CLOCKM  時間計測用サブルーチン(なくてもよい)
         |--SIMP              シンプソン積分用係数計算
         |--INPUT1            座標データ入力用
         |--INFOUT            計算情報出力(なくてもよい)
         |--LATTIC----HPSORT  基本逆格子ベクトル計算(LATTIC)
         |--GSTEP1            逆格子座標計算(2Gmax以下)
         |--HPSORT            逆格子座標を小さい順に並べる
         |--SYMM------OPGR------OPMTRX    対称性計算(SYMM)
         |--GSTEPF            FFT用マッピング
         |--KSTEP-----KPMSF   対称性なし(表面用)サンプリングk点計算
         |         |--KPMWBZ  対称性なし(バルク用)
         |         |--SCCM    単純立方用
         |         |--BCCM    体心立方用
         |         |--FCCM    面心立方用
         |         |--HEXM    六方最密充填用
         |         |--TETRAH  正方晶(ルチル構造)用
         |         |--APBO2   斜方晶(αPbO2構造)用
         |--BASNUM            平面波用逆格子計算(Gmax以下)
         |--LATSCA            LATTICストレス用(定圧)
         |--GSTSCA            GSTEP1ストレス用(定圧)
         |--SYMSCA            SYMMストレス用(定圧)
         |--GSFSCA            GSTEPFストレス用(定圧)
         |--FORM              フォームファクター計算用
         |--PCC-------DSJNV   部分内殻補正用(PCC)
         |--PSEUDO----DSJNV   局所擬ポテンシャル用(入力、フーリエ変換)(PSEUDO)
         |--PSELMD            PSEUDO(非MDループ用)
         |--EWVEC             エバルト計算用
         |--EWVMD             EWVEC(非MDループ用)
         |--INTCHG            初期電荷密度設定用
         |--EVIN              パラメータ入力(なくてもよい)
         |--XCFFT-----C3FFT   <---- MFFT  交換相関項計算用(XCFTT)
         |--KBMAT-----DSJNV   非局所擬ポテンシャル部分(入力、フーリエ変換)
         |--DIAGON----CHOBSD  初期波動関数(固有ベクトル)設定〔対角化〕(DIAGON)
         |--KBINT-----DSJNV   非局所擬ポテンシャル(非初期用、フーリエ変換)
         |--MSD-------KBINT-----DSJNV  カー・パリネロ本体(MSD)
         |         |--FORZFB    非局所部分計算用(非MDループ用)       
         |         |--C3FFT   <---- MFFT  高速フーリエ変換
         |         |--TIME------CLOCKM    時間計測(なくてもよい)
         |--FERMI-----WIDTH2  フェルミエネルギー計算用(FERMI)
         |--FORCE             力の計算(非局所部分)用
         |--FORZFB            非局所部分計算用(非MDループ用)
         |--CHAVER----C3FFT   <---- MFFT  電荷密度計算用(CHAVER)
         |         |--CHGAVR              電荷密度対称性計算用
         |--STRESS----XCFFT-----C3FFT <---- MFFT  ストレス計算用(STRESS)
         |         |--XSTPC-----C3FFT <---- MFFT  (局所部分用)
         |         |--STRNL   ストレス計算(非局所部分用)
         |--FORLOC            力の計算用(局所部分)
         |--FORCES            力の対称性計算用
         |--MD                古典的分子動力学部分用
         |--ENERGY----XCFFT-----C3FFT <---- MFFT  全エネルギー計算用(ENERGY)
         |--CONV2             収束判定用
         |--EVOUT             計算結果、パラメーター出力用(その1)
         |--EVOU2             計算結果、パラメーター出力用(その2)



[公開プログラム](←やや重めのページ、およそ5000行)
[サブルーチン群][作り方][各論編][著作権][Copyright][To Guide][Top]