计算化学公社

标题: Gamess 做能量分解的时候遇到了 floating point error [打印本页]
作者
Author: elvisng    时间: 2020-5-28 21:42
标题: Gamess 做能量分解的时候遇到了 floating point error
本帖最后由 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.
ls: No match.」



(, 下载次数 Times of downloads: 3)


(, 下载次数 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



作者
Author: zjxitcc    时间: 2020-5-29 00:07
本帖最后由 zjxitcc 于 2020-5-29 00:09 编辑

我用GAMESS 2017算就好好的，没问题。不知道是不是版本问题。为节省空间压缩了文件，需用7-zip解压
(, 下载次数 Times of downloads: 10)


作者
Author: paramecium86    时间: 2020-5-29 00:39
试着用2016和2019版都算了算，可以正常结束。是不是有可能是编译的问题。 以前编译的时候gcc gfortran 选的不好编译出来看起来没问题 算的时候会出类似这样的bug
作者
Author: elvisng    时间: 2020-5-29 14:10
paramecium86 发表于 2020-5-29 00:39
试着用2016和2019版都算了算，可以正常结束。是不是有可能是编译的问题。 以前编译的时候gcc gfortran 选的 ...

我把 gamess 重新編譯了一下, 這次用 gfortran (之前是 ifort），而且這次不用 mkl, 結果還是一樣的錯誤。不過也謝謝你的幫忙, 我要再慢慢找問題在哪裡。
作者
Author: zjxitcc    时间: 2020-5-29 14:16
elvisng 发表于 2020-5-29 14:10
我把 gamess 重新編譯了一下, 這次用 gfortran (之前是 ifort），而且這次不用 mkl, 結果還是一樣的錯誤 ...

我就是用ifort + MKL（都包含在INTEL全家桶parallel studio里）编译GAMESS的，没有问题。
作者
Author: paramecium86    时间: 2020-5-29 19:35
elvisng 发表于 2020-5-29 14:10
我把 gamess 重新編譯了一下, 這次用 gfortran (之前是 ifort），而且這次不用 mkl, 結果還是一樣的錯誤 ...

嗯 我用的是gfortran+mkl的组合  我之前出这类bug是因为 比如我的 gfortran比它要求的差了比如0.1的版本 比如2019版推荐4.9-6.3 我用4.8 编译  会出这种可以算但是中途某一步报这类错的信息。
作者
Author: elvisng    时间: 2020-5-29 22:38
paramecium86 发表于 2020-5-29 19:35
嗯 我用的是gfortran+mkl的组合  我之前出这类bug是因为 比如我的 gfortran比它要求的差了比如0.1的版本  ...

我嘗試用新版 gfortran 7.3 再編譯一次, 仍然是同樣的錯誤。我留意到你們的是 2019 June 版本, 我的是 2019 Sept, 不知道是不是新版反而有 bug。不知道能否分享一下舊版的 gamess?
作者
Author: zjxitcc    时间: 2020-5-29 23:18
elvisng 发表于 2020-5-29 22:38
我嘗試用新版 gfortran 7.3 再編譯一次, 仍然是同樣的錯誤。我留意到你們的是 2019 June 版本, 我的是 20 ...

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

是的, 謝謝你!
作者
Author: elvisng    时间: 2020-5-30 00:03
謝謝各位的幫忙! 問題已經解決, 用2019 R1 版本就沒有這個問題了! 這個問題在 Gamess 官方討論區也有人問到, 不過最終還是沒法解決。這次有賴大家幫忙, 解決了一個官方討論區也沒答案的問題!

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

哦哦 原来是因为GAMESS 19R2这个地方有bug呀。 我之前用19R2 算MCQDPT也出了bug。怎么编译都不行。用回到R1 就好了。
作者
Author: 喵星大佬    时间: 2021-3-12 16:01
本帖最后由 喵星大佬 于 2021-3-12 16:07 编辑
paramecium86 发表于 2020-5-30 00:24
哦哦 原来是因为GAMESS 19R2这个地方有bug呀。 我之前用19R2 算MCQDPT也出了bug。怎么编译都不行。用回到 ...

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




欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3