第一原理分子動力学プログラム STATE Senri Wiki
開始行:
* xc_gga.f90 [#o190ebca]
** subroutine xc_gga_chg [#od1e2fca]
- if文でxctypeに関する比較を順に行う際には最も使われやす...
現状ではggapbeを使う際のif文の回数がxctypeの数に比例して...
- xctypeが間違っている場合にelseでggapbeが呼ばれ計算が実...
xctypeが間違っているならエラーかxctypeのリストを出力してa...
** subroutine exch*_wb [#dd31e70e]
- LDAの交換エネルギーの計算でrho^{1/3}が必要になるので
rho = max(fft_grad_l(ig,0),eps_chg)
によりrhoが負になるのを防いでいるが、rhoをゼロでなくeps_c...
しかしrhoが振動し負の値が現れる真空領域は本来はほぼ密度一...
その副作用としてrhoが負の位置ではsの分母に微小量eps_chgが...
一方vdW-DF(の派生版)ではLieb-Oxford boundが効かない交換...
* 全般 [#mc7f00d4]
多くの関数が大量に引数を持つのは、定数を引数にしているた...
* Makefile [#oe5a3941]
- 本来-O1でコンパイルするファイルを
NOVECT = opgr.o filter.o inputf.o util.o
で区別していたはずが、5.5.*辺りから
NOPARA = opgr.o filter.o inputf.o util.o inopts...
は-O3でコンパイルするように変更されている。~
%%5.5.3で例えばnum_space_group=221とすると計算が止まって...
→num_space_groupで失敗するのはブリルアンゾーンの対称性が...
終了行:
* xc_gga.f90 [#o190ebca]
** subroutine xc_gga_chg [#od1e2fca]
- if文でxctypeに関する比較を順に行う際には最も使われやす...
現状ではggapbeを使う際のif文の回数がxctypeの数に比例して...
- xctypeが間違っている場合にelseでggapbeが呼ばれ計算が実...
xctypeが間違っているならエラーかxctypeのリストを出力してa...
** subroutine exch*_wb [#dd31e70e]
- LDAの交換エネルギーの計算でrho^{1/3}が必要になるので
rho = max(fft_grad_l(ig,0),eps_chg)
によりrhoが負になるのを防いでいるが、rhoをゼロでなくeps_c...
しかしrhoが振動し負の値が現れる真空領域は本来はほぼ密度一...
その副作用としてrhoが負の位置ではsの分母に微小量eps_chgが...
一方vdW-DF(の派生版)ではLieb-Oxford boundが効かない交換...
* 全般 [#mc7f00d4]
多くの関数が大量に引数を持つのは、定数を引数にしているた...
* Makefile [#oe5a3941]
- 本来-O1でコンパイルするファイルを
NOVECT = opgr.o filter.o inputf.o util.o
で区別していたはずが、5.5.*辺りから
NOPARA = opgr.o filter.o inputf.o util.o inopts...
は-O3でコンパイルするように変更されている。~
%%5.5.3で例えばnum_space_group=221とすると計算が止まって...
→num_space_groupで失敗するのはブリルアンゾーンの対称性が...
ページ名: