关于CP2k-8.2使用OT方法时SCF卡住的解决方法

丁越

关于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还真不少

https://github.com/cp2k/cp2k/issues/1444

https://github.com/cp2k/cp2k/pull/1515

所以当下的解决办法要么是在&GLOBAL 中替换掉默认的对角化库ELPA改用SCALAPACK库，要么在编译的时候直接写--with-elpa=no。

PS：openmpi需要自行编译4.1.1版本，结合--with-elpa=no可以顺利解决OT卡住问题。

ChemG

经过楼主的提示，编译的时候--with-elpa=no确实可以消除使用export OMP_NUM_THREADS=1 cp2k.psmp第一个离子步后卡在（OPTIMIZATION STEP:      1）的问题，简单的测试了一下速度变化不多。我的系统是centos 8.4 gcc版本8.4.1 编译方法按照社长的博文 http://sobereva.com/586。

丁越

ChemG 发表于 2021-10-7 14:46
经过楼主的提示，编译的时候--with-elpa=no确实可以消除使用export OMP_NUM_THREADS=1 cp2k.psmp第一个离子 ...

这也是没得办法的办法另外使用popt版，不要用psmp

ChemG

丁越 发表于 2021-10-8 20:09
这也是没得办法的办法另外使用popt版，不要用psmp

能解决问题就是好办法哈哈，对对忘记用popt了，重新编译测试完就赶紧来评论了

djjj148

添加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吧

zx4562430

丁老师，不知道您有没有试过刚发布的cp2k 9.1，这个bug修复了吗

ChemG

zx4562430 发表于 2022-1-7 17:26
丁老师，不知道您有没有试过刚发布的cp2k 9.1，这个bug修复了吗

修复了，不需要添加PREFERRED_DIAG_LIBRARY SL了

丁越

zx4562430 发表于 2022-1-7 17:26
丁老师，不知道您有没有试过刚发布的cp2k 9.1，这个bug修复了吗

我还没装呢，目前还用的8.2。等出9.2了我再装吧

zx4562430

ChemG 发表于 2022-1-7 21:54
修复了，不需要添加PREFERRED_DIAG_LIBRARY SL了

谢谢老师，我去试试