現在位置: ホーム / Forums / PHASE Forums / Install, Compile, Usage / PHASE 2019.02, total energyがNaNになる件

PHASE 2019.02, total energyがNaNになる件

Up to Install, Compile, Usage

PHASE 2019.02, total energyがNaNになる件

Posted by H.I. at December 31. 2020

お世話になります。

PHASE 2019.02の2次元版にて、
install.shのほぼデフォルト設定の通りmakeしてテスト実行したところ、
total energyがNaNになる結果が出てきました。
テスト環境、入力ファイル、結果などは下記の通りです。
どの様に対処したら良いかご教示頂けないでしょうか。

●install.sh設定
Supported platforms: 1) GNU Linux (EM64T/AMD64)
Supported compilers:0) GNU compiler collection (gfortran)
Supported programming-models 1) MPI parallel
Supported MPI libraries:0) Open MPI
Supported BLAS/LAPACK:0) Netlib BLAS/LAPACK
Supported FFT libraries:0) Built-in FFT subroutnes

●入力ファイル→phase0_2019.02/test/scf中のファイル一式
●実行コマンド→time mpirun -np 1 ../../bin/bin/phase ne=1 nk=1
   binの下にbinをもう一つ作り、その中に実行ファイルを保存しています。
   また途中で計算を打ち切っています(ctrl+C)。

●出力ファイル(output)→2019.02での実行結果を添付します。

●実行環境
windows10, 64bit, cygwin 64bit
gcc ver.9.3.0
gfortran ver.9.3.0
mpirun (Open MPI) 4.0.5

何卒よろしくお願いいたします。

Attachments

Re: PHASE 2019.02, total energyがNaNになる件

Posted by H.I. at December 31. 2020

補足します。

添付したoutput000中に

DEBUG @ m_ES_occup.F90 [AB002]

といった行が出力されていますが、私の計算環境由来かどうか調べる時に一時的に埋め込んだwrite文の出力です(特に何もわからなかったのですが…)。

普段の計算では埋め込まれていないコードです。

Re: PHASE 2019.02, total energyがNaNになる件

Posted by 宇佐見護 at January 07. 2021

直接的な回答はできませんが、現在の主要開発メンバーはPHASE/0実行(テスト)環境としてcygwinを用いておりません。
過去にも、cygwinでMPI並列まで行った例を、私は知りません。

cygwinをお使いであれば、まずは「MPI並列なし」をお試しください。

それよりもむしろ、Windows10でPHASE/0ご利用に際しては、VirtualBoxやWSLなどの仮想Linux環境のご利用をお勧めします。
以下、最新の情報ではありませんが、ご参考になれば幸いです。
https://qiita.com/USAMI_Mamoru/items/385e4eaf8fc30fbb4db2
https://qiita.com/USAMI_Mamoru/items/d7cde2391ab69a1c61ac

Re: PHASE 2019.02, total energyがNaNになる件

Posted by H.I. at January 10. 2021

宇佐見様

お世話になっております。
ご回答ありがとうございます。

ご回答いただいた
・cygwinでMPI並列なし
・VirtualBoxの利用
で試したところ、どちらも付属のtest/Run.sh中に書かれたエネルギー値
eref="-7.897015064593"
と一致する結果を得ました。

・cygwinでMPI並列なし
-7.897015064593

・VirtualBoxの利用:OpenMPI並列, -np 1
-7.897015064593

・VirtualBoxの利用:OpenMPI並列, -np 2
-7.897015064593

ひとまずtestのscfまで確認した段階ですが、
cygwinでMPI並列なしかVirtualBoxの利用で計算してみたいと思います。

 

Powered by Ploneboard