计算化学公社

 找回密码 Forget password
 注册 Register
Views: 18918|回复 Reply: 17
打印 Print 上一主题 Last thread 下一主题 Next thread

[GROMACS] GPU加速出错问题

[复制链接 Copy URL]

43

帖子

0

威望

142

eV
积分
185

Level 3 能力者

老师您好,我之前在做蛋白的模拟的时候,正式MD过程使用了GPU加速,但是现在做蛋白-配体的模拟,发现使用GPU会报错,不使用GPU则没有问题Command line:
  gmx_mpi mdrun -deffnm md01 -nb gpu

Reading file md01.tpr, VERSION 2016.4 (single precision)
Note: file tpx version 110, software tpx version 112
Multiple energy groups is not implemented for GPUs, falling back to the CPU. For better performance, run on the GPU without energy groups and then do gmx mdrun -rerun option on the trajectory with an energy group .tpr file.

-------------------------------------------------------
Program:     gmx mdrun, version 2018.3
Source file: src/gromacs/taskassignment/decidegpuusage.cpp (line 292)
Function:    bool gmx::decideWhetherToUseGpusForNonbonded(const gmx::TaskTarget, const std::vector<int> &, const EmulateGpuNonbonded, const bool, const bool, const bool)

Inconsistency in user input:
Nonbonded interactions on the GPU were required, but not supported for these
simulation settings. Change your settings, or do not require using GPUs.

这是我的mdp文件参数
integrator               = md
dt                       = 0.002
nsteps                   = 250000

nstxout                 = 500
nstvout                 = 500
nstfout                 = 500
nstenergy               = 500
nstlog                  = 500
energygrps              = Protein_CAL NA_SOL

nstlist                  = 1000
ns-type                  = Grid
pbc                      = xyz
rlist                    = 300.0
cutoff-scheme = verlet
coulombtype              = PME
pme_order                = 4   
fourierspacing           = 3
rcoulomb                 = 300.0
vdw-type                 = Cut-off
rvdw                     = 300.0

Tcoupl                   = V-rescale
tc-grps                  = Protein_CAL NA_SOL
tau_t                    = 3        3
ref_t                    = 370      370

DispCorr                 = EnerPres

Pcoupl                   = no

gen_vel                  = no

constraints              = all-bonds
continuation             = yes
constraint_algorithm     = lincs
lincs_iter               = 1
lincs_order              = 4
E-z                      = 1 0 0

请问老师这是说不能设置energygrps吗?
还有老师,如果使用nose-hoover热浴,会出现这个错误,是要扩大tau—t的值吗?

QQ图片20190713153704.jpg (17.04 KB, 下载次数 Times of downloads: 56)

QQ图片20190713153704.jpg

43

帖子

0

威望

142

eV
积分
185

Level 3 能力者

2#
 楼主 Author| 发表于 Post on 2019-7-13 20:46:29 | 只看该作者 Only view this author
老师,我发现如果不分开设置能量组,就可以使用GPU加速,请问可以这样设置吗?

5万

帖子

99

威望

5万

eV
积分
112384

管理员

公社社长

3#
发表于 Post on 2019-7-14 01:15:07 | 只看该作者 Only view this author
用GPU做MD的时候去掉能量组设置,之后再用纯CPU版把轨迹rerun一遍,这时候再设能量组

热浴问题我看你在群里已经问过了,我也已经回了。同一个问题不要24小时内同时在两个地方问,这点在群文件里群规以及本论坛置顶的新人必读里都已经说了。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

43

帖子

0

威望

142

eV
积分
185

Level 3 能力者

4#
 楼主 Author| 发表于 Post on 2019-7-14 12:06:36 | 只看该作者 Only view this author
sobereva 发表于 2019-7-14 01:15
用GPU做MD的时候去掉能量组设置,之后再用纯CPU版把轨迹rerun一遍,这时候再设能量组

热浴问题我看你在 ...

不好意思老师,我下次一定注意!
老师,我还有两个问题
1.看了您在群里关于热浴的回复,由于我看文献中做的这个体系是用nose-hoover,所以如果我一定要用它的话,是不是得按照报错提示,扩大tau-t的值啊?
2.老师,像跑这种蛋白-配体的体系,请问能量组一定要分开设定吗?可以不rerun吗?谢谢老师

5万

帖子

99

威望

5万

eV
积分
112384

管理员

公社社长

5#
发表于 Post on 2019-7-15 01:00:49 | 只看该作者 Only view this author
mjluan 发表于 2019-7-14 12:06
不好意思老师,我下次一定注意!
老师,我还有两个问题
1.看了您在群里关于热浴的回复,由于我看文献中 ...

1 NH热浴一般用tau-t=0.5就行了。如果你没有设一些额外东西,应该不会有那种提示
2 你要想考察蛋白质与配体的相互作用,肯定得定义两个energy group。不rerun的话你只能MD的过程用纯CPU来跑
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

43

帖子

0

威望

142

eV
积分
185

Level 3 能力者

6#
 楼主 Author| 发表于 Post on 2019-7-15 10:18:08 | 只看该作者 Only view this author
sobereva 发表于 2019-7-15 01:00
1 NH热浴一般用tau-t=0.5就行了。如果你没有设一些额外东西,应该不会有那种提示
2 你要想考察蛋白质与 ...

感谢老师!

8

帖子

0

威望

165

eV
积分
173

Level 3 能力者

7#
发表于 Post on 2020-1-17 13:48:04 | 只看该作者 Only view this author

楼主您好,我刚学GMX,请教一下:我的电脑是2核E5-2670+1块1080Ti,在跑GROMACS教程的例子“漏斗网蜘蛛毒素肽的溶剂化研究" 时,提示错误“Multiple energy groups is not implemented for GPUs, falling back to the CPU. For better performance, run on the GPU without energy groups and then do gmx mdrun -rerun option on the trajectory with an energy group .tpr file.”
我把mdp中的“energygrps = Protein Non-Protein”删掉也不管用,请问你是怎么屏蔽的能量组啊?谢谢!

5万

帖子

99

威望

5万

eV
积分
112384

管理员

公社社长

8#
发表于 Post on 2020-1-18 16:50:38 | 只看该作者 Only view this author
aken5555 发表于 2020-1-17 13:48
楼主您好,我刚学GMX,请教一下:我的电脑是2核E5-2670+1块1080Ti,在跑GROMACS教程的例子“漏斗网蜘蛛毒 ...

跑轨迹时候照常用GPU跑,不设能量组
之后rerun时候再用CPU跑,并且设能量组
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

8

帖子

0

威望

165

eV
积分
173

Level 3 能力者

9#
发表于 Post on 2020-1-27 15:49:02 | 只看该作者 Only view this author
sobereva 发表于 2020-1-18 16:50
跑轨迹时候照常用GPU跑,不设能量组
之后rerun时候再用CPU跑,并且设能量组

明白了,谢谢社长,新年快乐!

1

帖子

0

威望

9

eV
积分
10

Level 1 能力者

10#
发表于 Post on 2021-1-8 11:11:32 | 只看该作者 Only view this author
aken5555 发表于 2020-1-27 15:49
明白了,谢谢社长,新年快乐!

你好,请问你最后是怎么把能量组去掉的呀?

3

帖子

0

威望

25

eV
积分
28

Level 2 能力者

11#
发表于 Post on 2021-5-6 23:04:35 | 只看该作者 Only view this author
DS00HY 发表于 2021-1-8 11:11
**** 作者被禁止或删除 内容自动屏蔽 ****

您好,想问一下您最后是怎么删除能量组的牙

14

帖子

0

威望

109

eV
积分
123

Level 2 能力者

12#
发表于 Post on 2021-5-29 12:12:30 | 只看该作者 Only view this author
tomato0301 发表于 2021-5-6 23:04
您好,想问一下您最后是怎么删除能量组的牙

去除能量组就把.mdp文件里这一行删除“energygrps  = Protein  XX  ; Which energy group(s) to write to disk”

17

帖子

0

威望

215

eV
积分
232

Level 3 能力者

13#
发表于 Post on 2022-11-17 10:30:12 | 只看该作者 Only view this author
tomato0301 发表于 2021-5-6 23:04
**** 作者被禁止或删除 内容自动屏蔽 ****

请问怎么rerun,是在gmx mdrun后面吗?但是他老提示我“使用mdrun-rerun时,输入轨迹文件md_0_1.xtc的名称不能与输出文件的名称相同”这个就需要改一下文件名吗

13

帖子

0

威望

55

eV
积分
68

Level 2 能力者

14#
发表于 Post on 2024-6-23 10:23:57 | 只看该作者 Only view this author
Alpaca 发表于 2022-11-17 10:30
请问怎么rerun,是在gmx mdrun后面吗?但是他老提示我“使用mdrun-rerun时,输入轨迹文件md_0_1.xtc的名 ...

你好,请问你知道怎么rerun了吗,命令应该怎么写呀

251

帖子

0

威望

426

eV
积分
677

Level 4 (黑子)

A Student

15#
发表于 Post on 2024-6-23 11:04:25 | 只看该作者 Only view this author
本帖最后由 student0618 于 2024-6-23 13:21 编辑

可參考教程 http://www.mdtutorials.com/gmx/complex/09_analysis.html

  1. gmx mdrun -deffnm ie -rerun md_0_10.xtc -nb cpu
复制代码

以上指令rerun的軌跡是md_0_10.xtc
tpr文件是ie.tpr
輸出文件會是ie.xxx

Yet to be strong in theory, yet to have enough practical skills.
Still I am having fun with MD simulation.

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2024-11-25 11:31 , Processed in 0.274255 second(s), 29 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list