Message: 1
Date: Wed, 21 Nov 2007 21:54:57 +0100
From: David van der Spoel <[EMAIL PROTECTED]>
Subject: Re: [gmx-users] gromacs preformance versus amber
To: Discussion list for GROMACS users <gmx-users@gromacs.org>
Message-ID: <[EMAIL PROTECTED]>
Content-Type: text/plain; charset=ISO-8859-1; format=flowed

servaas michielssens wrote:
I am considering to use gromacs for MD-simulation for speed reasons, so
I tried to compare a amber simulation to a gromacs simulation. I ran a
20ps simulation with amber. I converted my topology and coordinate file
from amber to gromacs using and tried the same run in gromacs.
amb2gmx.pl(http://chemistry.csulb.edu/ffamber/). This is my amber input:


&cntrl
  imin=0,nstlim=10000,nmropt=0
  irest=1,ntx=5,
  iwrap=1,
  nscm=1000,
  dt=0.002,
  ntb=2,ntp=1,
  ntt=1,tautp=1.0,temp0=300.0,
  ntc=2,ntf=2,
  cut=8.0,
  ntwx=250,ntpr=250,
/

This is my gromacs input:

title           = FWS
cpp             = /usr/bin/cpp
constraints     = all-bonds
integrator      = md
dt              = 0.002 ; ps !
nsteps          = 10000 ; total 20 ps.
nstcomm         = 1000
nstxout         = 250 ; collect data every 1 ps
nstlog          = 0
nstenergy       = 250
nstvout         = 0
nstfout         = 0
nstlist         = 10
ns_type         = grid
rlist           = 1.0
coulombtype     = PME
rcoulomb        = 1.0
vdwtype         = cutoff
rvdw            = 1.4
fourierspacing  = 0.12
fourier_nx      = 0
fourier_ny      = 0
fourier_nz      = 0
pme_order       = 4
ewald_rtol      = 1e-5
optimize_fft    = yes
; Berendsen temperature coupling is on
Tcoupl          = berendsen
tau_t           = 1.0 1.0
tcgrps          = protein non-protein
ref_t           = 300 300
; Pressure coupling is on
Pcoupl          = parrinellorahman
tau_p           = 0.5
compressibility = 4.5e-5
ref_p           = 1.0
; Generate velocites is on at 300 K.
gen_vel         = no
gen_temp        = 300.0
gen_seed        = 937475


I ran it on 4 AMD Opteron(tm) Processor 250. It turned out that amber
could simulate 1.5ns/day and gromacs only 1ns/day. This was a surprise
to me, I expected that gromacs would be faster.  Could it be that using
gromacs with the amberFF is not faster than just using amber? I am new
to gromacs and perhaps I didn't compare them correct.

Did you run gromacs in parallel at all?
Did you get four md.log files, one for each processor?


Thanks for you answer,

Yes this worked fine, I have a log file for every processor, their must be some other problem. I also tried to run it on 8 or 16 processors, here the run failed. Amber runs worked fine, so I don't think their are problems with the processors

kind regards,

servaas michielssens.

_______________________________________________
gmx-users mailing list    gmx-users@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 [EMAIL PROTECTED]
Can't post? Read http://www.gromacs.org/mailing_lists/users.php

Reply via email to