プログラム構造(サブルーチン呼び出しリスト、各論編用)
以下に、これまでのものを含めて、公開が可能と考えられるサブルーチン部
分を一挙に公開したいと思います。(全てのコードが筆者の完全なオリジナル
ではないところ〔他の共同研究者による部分〕もありますが、それらはごく一
部で尚かつ非常に一般的なルーチンで、既にどこかで公開されていると考えら
れるコードです。もし問題があれば連絡して下さい。)
点滅している(ネットスケープでない人、ごめんなさい^^;)または
文字色が金色(黄)っぽいのが既に(2/28現在)公開されているサブルーチンです。
+-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)
|------DSJNV 局所擬ポテンシャル用(入力、フーリエ変換)(PSEUDO)
|--PSELMD PSEUDO(非MDループ用)
|--EWVEC エバルト計算用
|--EWVMD EWVEC(非MDループ用)
|--INTCHG 初期電荷密度設定用
|--EVIN パラメータ入力(なくてもよい)
|--XCFFT-----C3FFT <---- MFFT 交換相関項計算用(XCFTT)
|-------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 古典的分子動力学部分用
|------XCFFT-----C3FFT <---- MFFT 全エネルギー計算用(ENERGY)
|--CONV2 収束判定用
|--EVOUT 計算結果、パラメーター出力用(その1)
|--EVOU2 計算結果、パラメーター出力用(その2)
[To Guide][Top]