[gmx-users] Gromacs slow for 23000 atom DPPC bilayer on a (1 x 4) node: 50 ps in 10 hours

maria goranovic mariagoranovic at gmail.com
Mon Mar 24 15:08:51 CET 2008


Dear All,

My apologies. I had too big a simulation cell, and too few atoms, hence the
problem.

No particular reason to choose order 5. I will try with pme_order 4 and see
if it improves performance anyway.

thanks !

-maria

On Mon, Mar 24, 2008 at 2:19 PM, Carsten Kutzner <ckutzne at gwdg.de> wrote:

> Am 24.03.2008 um 10:17 schrieb maria goranovic:
>
> > Hi Folks,
> >
> > My simulation is running too slow. It took 10 wall clock hours (40
> > cpu hours) for a short 50 ps simulation of a ~ 23000 atom DPPC
> > bilayer. The hardware is a 4-cpu core. The installation is gromacs
> > 3.3.1. I have run much larger systems (~ 160000 atoms) using the
> > same gromacs installation on the same hardware, and they run much
> > faster than this (200 ps per 40 cpu hours).
> >
> > Can anybody suggest why this is happening ? Is it because of latency
> > in the cpu communication? If so, what is the workaround ?
> Is there a special reason for using pme_order=5? I would use the
> default, 4 instead, or at least an even number.
>
> Carsten
>
> > My .mdp script is below.
> > These are the run commands.
> >
> > grompp -np 4 -v -f heat.mdp -c minim4.gro -p dppc.top -n dppc.ndx -o
> > heat.tpr
> > mpirun -np 4 ~/bin/mdrun_mpi -np 4 -v -s heat.tpr -o heat.trr -c
> > heat.gro -e heat -g heat.log > & heat.out
> >
> > ;###########################################
> > ; heat.mdp
> > ;
> > title               =  heating
> > cpp                 =  /usr/bin/cpp
> > constraints         =  hbonds
> > constraint_algorithm = lincs
> > unconstrained_start =  yes
> > integrator          =  md
> > nsteps              =  25000
> > dt                  =  0.002
> > comm_mode           =  linear
> > nstxout             =  5000
> > nstvout             =  5000
> > nstlog              =  5000
> > nstenergy           =  5000
> > nstlist             =  10
> > ns_type             =  grid
> > pbc                 =  xyz
> > ; ------------------
> > coulombtype         = PME
> > rcoulomb            = 1.0
> > vdwtype             = cut-off
> > rlist               = 1.0
> > rvdw                = 1.0
> > fourierspacing      = 0.1
> > pme_order           = 5
> > ewald_rtol          = 1e-5
> > ; ---------------------------
> > ; Berendsen temperature and preasure coupling
> > Tcoupl              =  berendsen
> > tc-grps             =  DPPC  SOL
> > tau_t               =  0.6      0.6
> > ; i have also tried a tau value of 0.1, but no speed up
> > ref_t               =  323.0     323.0
> > Pcoupl              =  berendsen
> > Pcoupltype          =  semiisotropic
> > tau_p               =  1.0        1.0
> > compressibility     =  4.5e-5  4.5e-5
> > ref_p               =  1.0        1.0
> > ; ---------------------------
> > gen_vel             =  yes
> > gen_temp            =  323.0
> > gen_seed            =  194040
> > ;###########################################
> >
> > --
> > Maria G.
> > Technical University of Denmark
> > Copenhagen _______________________________________________
> > gmx-users mailing list    gmx-users at gromacs.org
> > http://www.gromacs.org/mailman/listinfo/gmx-users
> > Please search the archive at http://www.gromacs.org/search before
> > posting!
> > Please don't post (un)subscribe requests to the list. Use the
> > www interface or send it to gmx-users-request at gromacs.org.
> > Can't post? Read http://www.gromacs.org/mailing_lists/users.php
>
> _______________________________________________
> gmx-users mailing list    gmx-users at gromacs.org
> http://www.gromacs.org/mailman/listinfo/gmx-users
> Please search the archive at http://www.gromacs.org/search before posting!
> Please don't post (un)subscribe requests to the list. Use the
> www interface or send it to gmx-users-request at gromacs.org.
> Can't post? Read http://www.gromacs.org/mailing_lists/users.php
>



-- 
Maria G.
Technical University of Denmark
Copenhagen
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://maillist.sys.kth.se/pipermail/gromacs.org_gmx-users/attachments/20080324/53824fd4/attachment.html>


More information about the gromacs.org_gmx-users mailing list