xTAPP開発版170619

[yoshizk07]

Posted by 吉本 芳英 at June 19. 2017

xTAPP開発版170619を公開しました。

• 任意の形状の計算セルを最局在wannier関数の計算で利用できるようになりました。
• 最局在ワニエ関数の計算においてinner windowを利用できるようになりました。
• wannier関数で補間したバンド図を計算するhmatr2bndが利用する中間ファイル.hmatrのフォーマットが変更されました。hmatr2bndは補間に用いる実格子点位置をWigner Seitzセルの中にとるようになりました。

[yoshizk07]

Posted by 吉本 芳英 at June 19. 2017

hmatr2bndの変更によって、これまでの平行６面体内に実格子をとる場合に比べて、補間されたバンドのなめらかさが向上するはずです。k点の濃度が大きい極限で両方の方法は一致した結果になります。

[yoshizk07]

Posted by 藤堂 眞治 at October 26. 2017

xTAPP開発版170619をgfortran 4.7.2でコンパイルしようとすると

mpif90 -ffree-line-length-none -O3 -fopenmp -msse3 -fno-inline-functions -O3 -DPWB_SORTED_BY_INDEX -DESM -DINCONSISTENT_THREAD_COMP -I/usr/include -c vbchrg.f90

vbchrg.f90:837.11:

if (srg(irdkpt%irrg(irg))) then

1

Error: IF clause at (1) requires a scalar LOGICAL expression

というエラーで怒られます。487行目
if (srg(irdkpt%irrg(irg)).eq.1) then

と同じに修正すればよいでしょうか？

[yoshizk07]

Posted by 藤堂 眞治 at October 26. 2017

また、wannier.F90 についても

mpif90 -ffree-line-length-none -O3 -fopenmp -msse3 -fno-inline-functions -O3 -DPWB_SORTED_BY_INDEX -DESM -DINCONSISTENT_THREAD_COMP -I/usr/include -c wannier.F90

wannier.F90:232.14:

mod(nk,mmm(3,2,1).ne.0) ) then

1

Error: 'a' and 'p' arguments of 'mod' intrinsic at (1) must have the same type

と怒られますが、230行目から232行目は

if ( mod(nk,mmm(1,2,1)).ne.0 .or. &

mod(nk,mmm(2,2,1)).ne.0 .or. &

mod(nk,mmm(3,2,1)).ne.0 ) then

と修正したので良いでしょうか？

[yoshizk07]

Posted by 吉本 芳英 at October 27. 2017

藤堂先生、どうもありがとうございます。この修正で正しいです。バグです。

[yoshizk07]

Posted by 吉本 芳英 at October 27. 2017

藤堂先生、修正どうもありがとうございます。これで正しいです。バグになります。

[yoshizk07]

Posted by 吉本 芳英 at October 27. 2017

Previously 吉本 芳英 wrote:

藤堂先生、修正どうもありがとうございます。これで正しいです。バグになります。

---

これは最初のコメントに対する物です。

[yoshizk07]

Posted by 藤堂 眞治 at November 01. 2017

gfortran 6.3.0 では以下のようなエラーがでます。

makeqg.f90:419:26:

stop'GETQJ'

1

Error: Blank required in STOP statement near (1)

makeqg.f90以外にも、bsplpot.F90、esm_util.F90、esm_2d.F90、stropt_cnstr.F90、subrpot.f90、bunpu.F90、esm_manage.F90、doss.f90、sdoss.f90 の修正が必要です。上のコメントの修正も含めたパッチは https://github.com/cmsi/ma-xtapp/blob/develop/xtapp/debian/patches/fortran.patch にあります。

[yoshizk07]

Posted by 吉本 芳英 at November 01. 2017

藤堂先生

どうもありがとうございます。

gfortran 6.3.0で文法が厳格になったのでしょうか。

新しい版を作る際に取り込みます。

[yoshizk07]

Posted by 藤堂 眞治 at November 02. 2017

xtapp-util も一箇所だけ修正が必要なようです。

https://github.com/cmsi/ma-xtapp/blob/develop/xtapp-util/debian/patches/fortran.patch

[yoshizk07]

Posted by 吉本 芳英 at November 02. 2017

藤堂先生

どうもありがとうございます

こちらも直します