[gmx-users] "nstlist=-1" and performance problem
LuLanyuan
lulanyuan at msn.com
Thu Jan 29 01:40:55 CET 2009
Hi,
No warning was found. Is it because the use of "User" vdw potential?
The mdp file is the following:
title = A
cpp = cpp
include = -I../top
define =
integrator = md
dt = 0.001
nsteps = 10000
nstxout = 1000
nstvout = 1000
nstlog = 1000
nstenergy = 1000
nstxtcout = 0
energygrps = CGW CBB HAL NBB OBB
coulombtype = User
rcoulomb = 1.1
nstlist = -1
ns_type = grid
pbc = xyz
rlist = 1.1
vdwtype = User
rvdw = 1.1
table-extension = 0.
energygrp_table = CGW CGW CGW HAL CGW NBB CGW OBB CGW CBB CBB CBB CBB NBB CBB OBB CBB HAL HAL HAL HAL NBB HAL OBB NBB NBB NBB OBB OBB OBB
tcoupl = nose-hoover
tc_grps = System
tau_t = 0.1
ref_t = 310
pcoupl = no
gen_vel = yes
gen_temp = 310
gen_seed = 173529
Thanks.
Lanyuan
From: gmx3 at hotmail.com
To: gmx-users at gromacs.org
Subject: RE: [gmx-users] "nstlist=-1" and performance problem
Date: Wed, 28 Jan 2009 23:22:30 +0100
Hi,
But grompp should have given a warning about this.
You should have only obtained a tpr file when using the -maxwarn option of grompp.
Or did grompp not give a warning about the buffer being too small?
Berk
From: lulanyuan at msn.com
To: gmx-users at gromacs.org
Subject: RE: [gmx-users] "nstlist=-1" and performance problem
Date: Wed, 28 Jan 2009 17:12:00 -0500
Hi Berk,
Thanks very much. I've found the problem. As you said I forgot to increase rvdw. So the buffer is zero and the system needs to update the nblist every step, which is obviously wrong. After I corrected the mistake, the speed for "nstlist=-1" is roughly three times faster as I expected.
Lanyuan
From: gmx3 at hotmail.com
To: gmx-users at gromacs.org
Subject: RE: [gmx-users] "nstlist=-1" and performance problem
Date: Wed, 28 Jan 2009 10:42:59 +0100
Hi,
I assume you did not only change nstlist, but also other mdp settings, or not?
Changing only nstlist from 1 to -1 should result in nearly equal or much faster speeds.
For efficient simulations with nslist=-1 you need to set an appropriate
buffer region (rlist > rcoulomb=rvdw). In 4.0 you have to optimize this by hand.
It would be nice to have this automated.
At the end of the log file (a run of a few minutes is enough) you can see what
the neighborlist lifetime is. Optimal is usually something between 10 and 20.
Play with rlist to optimize the performance (also printed at the end of the log file).
Berk
From: lulanyuan at msn.com
To: gmx-users at gromacs.org
Date: Mon, 26 Jan 2009 17:31:43 -0500
Subject: [gmx-users] "nstlist=-1" and performance problem
Hello,
Could Berk or someone else answer my question regarding the new "nstlist=-1" option in gmx4? In my understanding, the algorithm for this option is the one used in DLPOLY/LAMMPS. However, When I was trying two gromacs runs with "nstlist=-1" and "nstlist=1", I found the speeds were basically the same. Isn't this automatic nblist checking supposed to be faster than updating the nblist every step? I also did similar tests on LAMMPS for the same system and found the speed for automatical updating was several times faster. As a result, although Gromacs is about 60% faster if we choose nstlist=1, it's much slower than LAMMPS if we choose the automatic nblist updating.
Thanks in advance for any ideas about it.
Lanyuan Lu
不加好友也能聊?试试MSN在线通! 现在就下载!
Express yourself instantly with MSN Messenger! MSN Messenger
八卦娱乐包打听,MSN资讯速递帮你忙! 了解详细!
See all the ways you can stay connected to friends and family
_________________________________________________________________
MSN热搜榜,每天最In的信息资讯和热点排行让您一览无余!
http://top.msn.com.cn
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://maillist.sys.kth.se/pipermail/gromacs.org_gmx-users/attachments/20090128/e994a58a/attachment.html>
More information about the gromacs.org_gmx-users
mailing list