计算化学公社
标题: 关于CP2k-8.2使用OT方法时SCF卡住的解决方法 [打印本页]
作者Author: 丁越 时间: 2021-10-3 11:21
标题: 关于CP2k-8.2使用OT方法时SCF卡住的解决方法
本帖最后由 丁越 于 2021-12-6 22:31 编辑
关于CP2k-8.2使用OT方法时SCF卡住的解决方法
最近在论坛上好多小伙伴反映在编译openmpi的CP2K-8.2时明明一切编译的很顺利，没有任何报错，但是使用OT算法时SCF却莫名其妙卡死。于是去google group 上咨询了一下，Hutter教授是这样回复的：
there is a good chance that this behavior is due to an interface problem with the ELPA library. You can find information about it
in this mailing list or on the github repository.
Easy fix is to use the SCALAPACK library by adding
&GLOBAL
...
PREFERRED_DIAG_LIBRARY SL
&END GLOBAL
to your input.
regards
Juerg Hutter
由于CP2K-8.2中toolchain编译默认安装的ELPA库是elpa-2020.11.001.tar.gz，目前这个库太新了导致的bug还真不少
所以当下的解决办法要么是在&GLOBAL 中替换掉默认的对角化库ELPA改用SCALAPACK库，要么在编译的时候直接写--with-elpa=no。
PS：openmpi需要自行编译4.1.1版本，结合--with-elpa=no可以顺利解决OT卡住问题。
作者Author: ChemG 时间: 2021-10-7 14:46
本帖最后由 ChemG 于 2021-10-8 12:39 编辑
经过楼主的提示，编译的时候--with-elpa=no确实可以消除使用export OMP_NUM_THREADS=1 cp2k.psmp第一个离子步后卡在（OPTIMIZATION STEP: 1）的问题，简单的测试了一下速度变化不多。我的系统是centos 8.4 gcc版本8.4.1 编译方法按照社长的博文 http://sobereva.com/586。
作者Author: 丁越 时间: 2021-10-8 20:09
这也是没得办法的办法
另外使用popt版，不要用psmp
作者Author: ChemG 时间: 2021-10-9 09:39
能解决问题就是好办法哈哈，对对忘记用popt了，重新编译测试完就赶紧来评论了
作者Author: djjj148 时间: 2021-10-9 10:15
添加PREFERRED_DIAG_LIBRARY SL亲测有效，运行速度会有一点损失。按照楼主另个帖子的方法编译8.2
使用intel全家桶编译cp2k8.2
http://bbs.keinsci.com/forum.php ... 25162&fromuid=15641
(出处: 计算化学公社)
速度比intel编译的7.1还快10%左右，但是第一个离子步会卡住，添加PREFERRED_DIAG_LIBRARY SL后不会卡住，但是速度会变得和7.1一样。
希望官方快点修复这个bug吧
作者Author: zx4562430 时间: 2022-1-7 17:26
丁老师，不知道您有没有试过刚发布的cp2k 9.1，这个bug修复了吗
作者Author: ChemG 时间: 2022-1-7 21:54
修复了，不需要添加PREFERRED_DIAG_LIBRARY SL了
作者Author: 丁越 时间: 2022-1-8 09:04
我还没装呢，目前还用的8.2。等出9.2了我再装吧
作者Author: zx4562430 时间: 2022-1-8 16:22
谢谢老师，我去试试
| 欢迎光临 计算化学公社 (http://bbs.keinsci.com/) |
Powered by Discuz! X3.3 |