Gamess 做能量分解的时候遇到了 floating point error

elvisng · 发表于 Post on 2020-5-28 21:42:59

本帖最后由 elvisng 于 2020-5-28 23:25 编辑

第一次做能量分解，一直发生各种错误。经过各种试错之后，现在遇到了floating point error （节录在下面，这些错误内容只显示了在命令窗口，没有在输出文件里显示），实在不知道如何入手。

「DDIProcess 9: trapped a floating point error (SIGFPE).

DDI Process12: trapped a floating point error (SIGFPE).

DDI Process15: trapped a floating point error (SIGFPE).

DDI Process5: trapped a floating point error (SIGFPE).

DDI Process13: trapped a floating point error (SIGFPE).

DDI Process11: trapped a floating point error (SIGFPE).

DDI Process10: trapped a floating point error (SIGFPE).

DDI Process14: trapped a floating point error (SIGFPE).

DDI Process4: trapped a floating point error (SIGFPE).

DDI Process3: trapped a floating point error (SIGFPE).

DDI Process6: trapped a floating point error (SIGFPE).

DDI Process17: trapped a floating point error (SIGFPE).

DDI Process8: trapped a floating point error (SIGFPE).

DDI Process7: trapped a floating point error (SIGFPE).

DDI Process2: trapped a floating point error (SIGFPE).

DDI Process16: trapped a floating point error (SIGFPE).

DDI Process0: trapped a floating point error (SIGFPE).

DDI Process1: trapped a floating point error (SIGFPE).

ddikick.x:Fatal error detected.

The error ismost likely to be in the application, so check for

inputerrors, disk space, memory needs, application bugs, etc.

ddikick.xwill now clean up all processes, and exit...

unset echo

ls: No match.

ls: No match.」

EDA_simple.log (110.41 KB, 下载次数 Times of downloads: 3)

EDA_simple.inp (2.5 KB, 下载次数 Times of downloads: 8)

--更新--

为了方便大家测试，把输入文件进一步简化了，但依然同样的报错。

补充进一步信息：我曾经试过把核心数调成1，也试过调整 Mwords 和 Memddi，但依然报错（我的内存有64GB）。

也试了运行测试文件的EDA例子（exam28），那是成功没报错的。

$CONTRLSCFTYP=RHF RUNTYP=eda

DFTTYP=B3LYPV1R ISPHER=1 nosym=1 $END

$BASIS GBASIS=N31 NGAUSS=6 NDFUNC=1 NPFUNC=0DIFFSP=.F. DIFFS=.F. $END

$SYSTEM MWORDS=125 MEMDDI=4000 $END

$lmoeda matom(1)=20,6 mcharg(1)=0,0mmult(1)=1,1 $end

$DATA

syn-1c

C1

C 6.0 2.17846900 0.76043200 -0.35904400

C 6.0 1.71866500 -0.67787500 -0.43591500

C 6.0 1.14688200 1.60696000 0.40540400

Ni 28.0 -0.65316300 0.53172100 0.06140000

C 6.0 0.43521900 -0.96837300 -0.16637100

C 6.0 -2.25885200 -0.44569800 -0.01044700

N 7.0 -3.28242100 -0.26354700 -0.88964900

C 6.0 -4.30954500 -1.16772100 -0.66288200

C 6.0 -3.92978200 -1.92913000 0.39496300

N 7.0 -2.67551000 -1.47523700 0.77914900

H 1.0 2.25916900 1.14757400 -1.38919200

H 1.0 1.19842600 1.35866100 1.47582400

H 1.0 -5.20430800 -1.18530000 -1.26524000

H 1.0 -4.43433500 -2.73196000 0.90930800

H 1.0 0.10112800 -2.00967200 -0.23360900

H 1.0 1.44883000 2.66750400 0.34774900

H 1.0 2.45910624 -1.38580161 -0.74492421

H 1.0 3.15453431 0.80994804 0.07655291

H 1.0 -3.26541175 0.40685257 -1.63145430

H 1.0 -2.16844938 -1.82556068 1.56665319

O 8.0 -1.33257700 2.40040500 -0.05458900

C 6.0 -1.87858200 2.92626800 1.14147600

H 1.0 -2.82131600 2.41729300 1.38922800

H 1.0 -1.19334500 2.80279100 1.99393100

H 1.0 -2.08431500 3.99856000 1.01909600

H 1.0 -0.13162200 2.17944600 0.05021700

$END

zjxitcc · 发表于 Post on 2020-5-29 00:07:34

本帖最后由 zjxitcc 于 2020-5-29 00:09 编辑

我用GAMESS 2017算就好好的，没问题。不知道是不是版本问题。为节省空间压缩了文件，需用7-zip解压

EDA_simple.7z (2.45 MB, 下载次数 Times of downloads: 10)

paramecium86 · 发表于 Post on 2020-5-29 00:39:43

试着用2016和2019版都算了算，可以正常结束。是不是有可能是编译的问题。以前编译的时候gcc gfortran 选的不好编译出来看起来没问题算的时候会出类似这样的bug

elvisng · 发表于 Post on 2020-5-29 14:10:58

paramecium86 发表于 2020-5-29 00:39
试着用2016和2019版都算了算，可以正常结束。是不是有可能是编译的问题。以前编译的时候gcc gfortran 选的 ...

我把 gamess 重新編譯了一下, 這次用 gfortran (之前是 ifort），而且這次不用 mkl, 結果還是一樣的錯誤。不過也謝謝你的幫忙, 我要再慢慢找問題在哪裡。

zjxitcc · 发表于 Post on 2020-5-29 14:16:21

elvisng 发表于 2020-5-29 14:10
我把 gamess 重新編譯了一下, 這次用 gfortran (之前是 ifort），而且這次不用 mkl, 結果還是一樣的錯誤 ...

我就是用ifort + MKL（都包含在INTEL全家桶parallel studio里）编译GAMESS的，没有问题。

paramecium86 · 发表于 Post on 2020-5-29 19:35:14

elvisng 发表于 2020-5-29 14:10
我把 gamess 重新編譯了一下, 這次用 gfortran (之前是 ifort），而且這次不用 mkl, 結果還是一樣的錯誤 ...

嗯我用的是gfortran+mkl的组合我之前出这类bug是因为比如我的 gfortran比它要求的差了比如0.1的版本比如2019版推荐4.9-6.3 我用4.8 编译会出这种可以算但是中途某一步报这类错的信息。

elvisng · 发表于 Post on 2020-5-29 22:38:15

paramecium86 发表于 2020-5-29 19:35
嗯我用的是gfortran+mkl的组合我之前出这类bug是因为比如我的 gfortran比它要求的差了比如0.1的版本 ...

我嘗試用新版 gfortran 7.3 再編譯一次, 仍然是同樣的錯誤。我留意到你們的是 2019 June 版本, 我的是 2019 Sept, 不知道是不是新版反而有 bug。不知道能否分享一下舊版的 gamess?

zjxitcc · 发表于 Post on 2020-5-29 23:18:47

elvisng 发表于 2020-5-29 22:38
我嘗試用新版 gfortran 7.3 再編譯一次, 仍然是同樣的錯誤。我留意到你們的是 2019 June 版本, 我的是 20 ...

你的gfortran可能版本太高了。另外，旧版的GAMESS在版上搜搜就有。

elvisng · 发表于 Post on 2020-5-29 23:23:54

zjxitcc 发表于 2020-5-29 23:18
你的gfortran可能版本太高了。另外，旧版的GAMESS在版上搜搜就有。

是的, 謝謝你!

elvisng · 发表于 Post on 2020-5-30 00:03:43

謝謝各位的幫忙! 問題已經解決, 用2019 R1 版本就沒有這個問題了! 這個問題在 Gamess 官方討論區也有人問到, 不過最終還是沒法解決。這次有賴大家幫忙, 解決了一個官方討論區也沒答案的問題!

paramecium86 · 发表于 Post on 2020-5-30 00:24:58

elvisng 发表于 2020-5-30 00:03
謝謝各位的幫忙! 問題已經解決, 用2019 R1 版本就沒有這個問題了! 這個問題在 Gamess 官方討論區也有人問到 ...

哦哦原来是因为GAMESS 19R2这个地方有bug呀。我之前用19R2 算MCQDPT也出了bug。怎么编译都不行。用回到R1 就好了。

喵星大佬 · 发表于 Post on 2021-3-12 16:01:52

本帖最后由喵星大佬于 2021-3-12 16:07 编辑

paramecium86 发表于 2020-5-30 00:24
哦哦原来是因为GAMESS 19R2这个地方有bug呀。我之前用19R2 算MCQDPT也出了bug。怎么编译都不行。用回到 ...

2020R2版本里面这个已经修复了

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[GAMESS-US] Gamess 做能量分解的时候遇到了 floating point error

评分 Rate

评分 Rate

浏览过的版块