Hi Mark,<br><br>Yes I have compared the runtime for single as well as 8 processor . Could you please tell me the right steps to configure GROMACS-MPI version?.<br><br>Thanks.<br><br>thamu<br><br><div class="gmail_quote">2009/6/11  <span dir="ltr">&lt;<a href="mailto:gmx-users-request@gromacs.org">gmx-users-request@gromacs.org</a>&gt;</span><br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">Send gmx-users mailing list submissions to<br>
        <a href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a><br>
<br>
To subscribe or unsubscribe via the World Wide Web, visit<br>
        <a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/mailman/listinfo/gmx-users</a><br>
or, via email, send a message with subject or body &#39;help&#39; to<br>
        <a href="mailto:gmx-users-request@gromacs.org">gmx-users-request@gromacs.org</a><br>
<br>
You can reach the person managing the list at<br>
        <a href="mailto:gmx-users-owner@gromacs.org">gmx-users-owner@gromacs.org</a><br>
<br>
When replying, please edit your Subject line so it is more specific<br>
than &quot;Re: Contents of gmx-users digest...&quot;<br>
<br>
<br>
Today&#39;s Topics:<br>
<br>
   1. Re: gromacs-4.0.5 parallel run in 8 cpu: slow speed (Mark Abraham)<br>
<br>
<br>
----------------------------------------------------------------------<br>
<br>
Message: 1<br>
Date: Thu, 11 Jun 2009 13:56:15 +0000 (GMT)<br>
From: Mark Abraham &lt;<a href="mailto:mark.abraham@anu.edu.au">mark.abraham@anu.edu.au</a>&gt;<br>
Subject: Re: [gmx-users] gromacs-4.0.5 parallel run in 8 cpu: slow<br>
        speed<br>
To: Discussion list for GROMACS users &lt;<a href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a>&gt;<br>
Message-ID: &lt;<a href="mailto:fbe6b07c8ea.4a310cff@anu.edu.au">fbe6b07c8ea.4a310cff@anu.edu.au</a>&gt;<br>
Content-Type: text/plain; charset=&quot;iso-8859-1&quot;<br>
<br>
On 06/11/09, Thamu  &lt;<a href="mailto:asthamu@gmail.com">asthamu@gmail.com</a>&gt; wrote:<br>
&gt;<br>
&gt; Hi Mark,<br>
&gt;<br>
&gt; The top md.log is below. The mdrun command was &quot;mpirun -np 8 ~/software/bin/mdrun_mpi -deffnm md&quot;<br>
In my experience, correctly-configured MPI gromacs running in parallel reports information about the number of nodes and the identity of the node writing the .log file. This is missing, so something is wrong with your setup.<br>

<br>
I&#39;ve assumed that you&#39;ve compared this &quot;8-processor&quot; runtime with a single-processor runtime and found them comparable...<br>
<br>
Mark<br>
<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt;                           :-)   G  R  O  M  A  C  S  (-:<br>
&gt;<br>
&gt;                       GROup of MAchos and Cynical Suckers<br>
&gt;<br>
&gt;                              :-)   VERSION 4.0.5  (-:<br>
&gt;<br>
&gt;<br>
&gt;       Written by David van der Spoel, Erik Lindahl, Berk Hess, and others.<br>
&gt;        Copyright (c) 1991-2000, University of Groningen, The Netherlands.<br>
&gt;              Copyright (c) 2001-2008, The GROMACS development team,<br>
&gt;             check out <a href="http://www.gromacs.org" target="_blank">http://www.gromacs.org</a> (<a href="http://www.gromacs.org" target="_blank">http://www.gromacs.org</a>) for more information.<br>
&gt;<br>
&gt;          This program is free software; you can redistribute it and/or<br>
&gt;           modify it under the terms of the GNU General Public License<br>
&gt;          as published by the Free Software Foundation; either version 2<br>
&gt;              of the License, or (at your option) any later version.<br>
&gt;<br>
&gt;                    :-)   /home/thamu/software/bin/mdrun_mpi  (-:<br>
&gt;<br>
&gt;<br>
&gt; ++++ PLEASE READ AND CITE THE FOLLOWING REFERENCE ++++<br>
&gt; B. Hess and C. Kutzner and D. van der Spoel and E. Lindahl<br>
&gt; GROMACS 4: Algorithms for highly efficient, load-balanced, and scalable<br>
&gt; molecular simulation<br>
&gt; J. Chem. Theory Comput. 4 (2008) pp. 435-447<br>
&gt; -------- -------- --- Thank You --- -------- --------<br>
&gt;<br>
&gt;<br>
&gt; ++++ PLEASE READ AND CITE THE FOLLOWING REFERENCE ++++<br>
&gt; D. van der Spoel, E. Lindahl, B. Hess, G. Groenhof, A. E. Mark and H. J. C.<br>
&gt; Berendsen<br>
&gt; GROMACS: Fast, Flexible and Free<br>
&gt; J. Comp. Chem. 26 (2005) pp. 1701-1719<br>
&gt; -------- -------- --- Thank You --- -------- --------<br>
&gt;<br>
&gt;<br>
&gt; ++++ PLEASE READ AND CITE THE FOLLOWING REFERENCE ++++<br>
&gt; E. Lindahl and B. Hess and D. van der Spoel<br>
&gt; GROMACS 3.0: A package for molecular simulation and trajectory analysis<br>
&gt; J. Mol. Mod. 7 (2001) pp. 306-317<br>
&gt; -------- -------- --- Thank You --- -------- --------<br>
&gt;<br>
&gt;<br>
&gt; ++++ PLEASE READ AND CITE THE FOLLOWING REFERENCE ++++<br>
&gt; H. J. C. Berendsen, D. van der Spoel and R. van Drunen<br>
&gt; GROMACS: A message-passing parallel molecular dynamics implementation<br>
&gt; Comp. Phys. Comm. 91 (1995) pp. 43-56<br>
&gt; -------- -------- --- Thank You --- -------- --------<br>
&gt;<br>
&gt; Input Parameters:<br>
&gt;    integrator           = md<br>
&gt;    nsteps               = 10000000<br>
&gt;    init_step            = 0<br>
&gt;    ns_type              = Grid<br>
&gt;    nstlist              = 10<br>
&gt;    ndelta               = 2<br>
&gt;    nstcomm              = 1<br>
&gt;    comm_mode            = Linear<br>
&gt;    nstlog               = 100<br>
&gt;    nstxout              = 1000<br>
&gt;    nstvout              = 0<br>
&gt;    nstfout              = 0<br>
&gt;    nstenergy            = 100<br>
&gt;    nstxtcout            = 0<br>
&gt;    init_t               = 0<br>
&gt;    delta_t              = 0.002<br>
&gt;    xtcprec              = 1000<br>
&gt;    nkx                  = 70<br>
&gt;    nky                  = 70<br>
&gt;    nkz                  = 70<br>
&gt;    pme_order            = 4<br>
&gt;    ewald_rtol           = 1e-05<br>
&gt;    ewald_geometry       = 0<br>
&gt;    epsilon_surface      = 0<br>
&gt;    optimize_fft         = TRUE<br>
&gt;    ePBC                 = xyz<br>
&gt;    bPeriodicMols        = FALSE<br>
&gt;    bContinuation        = FALSE<br>
&gt;    bShakeSOR            = FALSE<br>
&gt;    etc                  = V-rescale<br>
&gt;    epc                  = Parrinello-Rahman<br>
&gt;    epctype              = Isotropic<br>
&gt;    tau_p                = 0.5<br>
&gt;    ref_p (3x3):<br>
&gt;       ref_p[    0]={ 1.00000e+00,  0.00000e+00,  0.00000e+00}<br>
&gt;       ref_p[    1]={ 0.00000e+00,  1.00000e+00,  0.00000e+00}<br>
&gt;       ref_p[    2]={ 0.00000e+00,  0.00000e+00,  1.00000e+00}<br>
&gt;    compress (3x3):<br>
&gt;       compress[    0]={ 4.50000e-05,  0.00000e+00,  0.00000e+00}<br>
&gt;       compress[    1]={ 0.00000e+00,  4.50000e-05,  0.00000e+00}<br>
&gt;       compress[    2]={ 0.00000e+00,  0.00000e+00,  4.50000e-05}<br>
&gt;    refcoord_scaling     = No<br>
&gt;    posres_com (3):<br>
&gt;       posres_com[0]= 0.00000e+00<br>
&gt;       posres_com[1]= 0.00000e+00<br>
&gt;       posres_com[2]= 0.00000e+00<br>
&gt;    posres_comB (3):<br>
&gt;       posres_comB[0]= 0.00000e+00<br>
&gt;       posres_comB[1]= 0.00000e+00<br>
&gt;       posres_comB[2]= 0.00000e+00<br>
&gt;    andersen_seed        = 815131<br>
&gt;    rlist                = 1<br>
&gt;    rtpi                 = 0.05<br>
&gt;    coulombtype          = PME<br>
&gt;    rcoulomb_switch      = 0<br>
&gt;    rcoulomb             = 1<br>
&gt;    vdwtype              = Cut-off<br>
&gt;    rvdw_switch          = 0<br>
&gt;    rvdw                 = 1.4<br>
&gt;    epsilon_r            = 1<br>
&gt;    epsilon_rf           = 1<br>
&gt;    tabext               = 1<br>
&gt;    implicit_solvent     = No<br>
&gt;    gb_algorithm         = Still<br>
&gt;    gb_epsilon_solvent   = 80<br>
&gt;    nstgbradii           = 1<br>
&gt;    rgbradii             = 2<br>
&gt;    gb_saltconc          = 0<br>
&gt;    gb_obc_alpha         = 1<br>
&gt;    gb_obc_beta          = 0.8<br>
&gt;    gb_obc_gamma         = 4.85<br>
&gt;    sa_surface_tension   = 2.092<br>
&gt;    DispCorr             = No<br>
&gt;    free_energy          = no<br>
&gt;    init_lambda          = 0<br>
&gt;    sc_alpha             = 0<br>
&gt;    sc_power             = 0<br>
&gt;    sc_sigma             = 0.3<br>
&gt;    delta_lambda         = 0<br>
&gt;    nwall                = 0<br>
&gt;    wall_type            = 9-3<br>
&gt;    wall_atomtype[0]     = -1<br>
&gt;    wall_atomtype[1]     = -1<br>
&gt;    wall_density[0]      = 0<br>
&gt;    wall_density[1]      = 0<br>
&gt;    wall_ewald_zfac      = 3<br>
&gt;    pull                 = no<br>
&gt;    disre                = No<br>
&gt;    disre_weighting      = Conservative<br>
&gt;    disre_mixed          = FALSE<br>
&gt;    dr_fc                = 1000<br>
&gt;    dr_tau               = 0<br>
&gt;    nstdisreout          = 100<br>
&gt;    orires_fc            = 0<br>
&gt;    orires_tau           = 0<br>
&gt;    nstorireout          = 100<br>
&gt;    dihre-fc             = 1000<br>
&gt;    em_stepsize          = 0.01<br>
&gt;    em_tol               = 10<br>
&gt;    niter                = 20<br>
&gt;    fc_stepsize          = 0<br>
&gt;    nstcgsteep           = 1000<br>
&gt;    nbfgscorr            = 10<br>
&gt;    ConstAlg             = Lincs<br>
&gt;    shake_tol            = 0.0001<br>
&gt;    lincs_order          = 4<br>
&gt;    lincs_warnangle      = 30<br>
&gt;    lincs_iter           = 1<br>
&gt;    bd_fric              = 0<br>
&gt;    ld_seed              = 1993<br>
&gt;    cos_accel            = 0<br>
&gt;    deform (3x3):<br>
&gt;       deform[    0]={ 0.00000e+00,  0.00000e+00,  0.00000e+00}<br>
&gt;       deform[    1]={ 0.00000e+00,  0.00000e+00,  0.00000e+00}<br>
&gt;       deform[    2]={ 0.00000e+00,  0.00000e+00,  0.00000e+00}<br>
&gt;    userint1             = 0<br>
&gt;    userint2             = 0<br>
&gt;    userint3             = 0<br>
&gt;    userint4             = 0<br>
&gt;    userreal1            = 0<br>
&gt;    userreal2            = 0<br>
&gt;    userreal3            = 0<br>
&gt;    userreal4            = 0<br>
&gt; grpopts:<br>
&gt;    nrdf:     6706.82      106800<br>
&gt;    ref_t:         300         300<br>
&gt;    tau_t:         0.1         0.1<br>
&gt; anneal:          No          No<br>
&gt; ann_npoints:           0           0<br>
&gt;    acc:               0           0           0<br>
&gt;    nfreeze:           N           N           N<br>
&gt;    energygrp_flags[  0]: 0 0 0<br>
&gt;    energygrp_flags[  1]: 0 0 0<br>
&gt;    energygrp_flags[  2]: 0 0 0<br>
&gt;    efield-x:<br>
&gt;       n = 0<br>
&gt;    efield-xt:<br>
&gt;       n = 0<br>
&gt;    efield-y:<br>
&gt;       n = 0<br>
&gt;    efield-yt:<br>
&gt;       n = 0<br>
&gt;    efield-z:<br>
&gt;       n = 0<br>
&gt;    efield-zt:<br>
&gt;       n = 0<br>
&gt;    bQMMM                = FALSE<br>
&gt;    QMconstraints        = 0<br>
&gt;    QMMMscheme           = 0<br>
&gt;    scalefactor          = 1<br>
&gt; qm_opts:<br>
&gt;    ngQM                 = 0<br>
&gt; Table routines are used for coulomb: TRUE<br>
&gt; Table routines are used for vdw:     FALSE<br>
&gt; Will do PME sum in reciprocal space.<br>
&gt;<br>
&gt; ++++ PLEASE READ AND CITE THE FOLLOWING REFERENCE ++++<br>
&gt; U. Essman, L. Perela, M. L. Berkowitz, T. Darden, H. Lee and L. G. Pedersen<br>
&gt; A smooth particle mesh Ewald method<br>
&gt; J. Chem. Phys. 103 (1995) pp. 8577-8592<br>
&gt; -------- -------- --- Thank You --- -------- --------<br>
&gt;<br>
&gt; Using a Gaussian width (1/beta) of 0.320163 nm for Ewald<br>
&gt; Cut-off&#39;s:   NS: 1   Coulomb: 1   LJ: 1.4<br>
&gt; System total charge: -0.000<br>
&gt; Generated table with 1200 data points for Ewald.<br>
&gt; Tabscale = 500 points/nm<br>
&gt; Generated table with 1200 data points for LJ6.<br>
&gt; Tabscale = 500 points/nm<br>
&gt; Generated table with 1200 data points for LJ12.<br>
&gt; Tabscale = 500 points/nm<br>
&gt; Generated table with 1200 data points for 1-4 COUL.<br>
&gt; Tabscale = 500 points/nm<br>
&gt; Generated table with 1200 data points for 1-4 LJ6.<br>
&gt; Tabscale = 500 points/nm<br>
&gt; Generated table with 1200 data points for 1-4 LJ12.<br>
&gt; Tabscale = 500 points/nm<br>
&gt;<br>
&gt; Enabling TIP4p water optimization for 17798 molecules.<br>
&gt;<br>
&gt; Configuring nonbonded kernels...<br>
&gt; Testing x86_64 SSE support... present.<br>
&gt;<br>
&gt;<br>
&gt; Removing pbc first time<br>
&gt;<br>
&gt; Initializing LINear Constraint Solver<br>
&gt;<br>
&gt; ++++ PLEASE READ AND CITE THE FOLLOWING REFERENCE ++++<br>
&gt; B. Hess and H. Bekker and H. J. C. Berendsen and J. G. E. M. Fraaije<br>
&gt; LINCS: A Linear Constraint Solver for molecular simulations<br>
&gt; J. Comp. Chem. 18 (1997) pp. 1463-1472<br>
&gt; -------- -------- --- Thank You --- -------- --------<br>
&gt;<br>
&gt; The number of constraints is 3439<br>
&gt;<br>
&gt; ++++ PLEASE READ AND CITE THE FOLLOWING REFERENCE ++++<br>
&gt; S. Miyamoto and P. A. Kollman<br>
&gt; SETTLE: An Analytical Version of the SHAKE and RATTLE Algorithms for Rigid<br>
&gt; Water Models<br>
&gt; J. Comp. Chem. 13 (1992) pp. 952-962<br>
&gt; -------- -------- --- Thank You --- -------- --------<br>
&gt;<br>
&gt; Center of mass motion removal mode is Linear<br>
&gt; We have the following groups for center of mass motion removal:<br>
&gt;   0:  rest<br>
&gt;<br>
&gt; ++++ PLEASE READ AND CITE THE FOLLOWING REFERENCE ++++<br>
&gt; G. Bussi, D. Donadio and M. Parrinello<br>
&gt; Canonical sampling through velocity rescaling<br>
&gt; J. Chem. Phys. 126 (2007) pp. 014101<br>
&gt; -------- -------- --- Thank You --- -------- --------<br>
&gt;<br>
&gt; There are: 56781 Atoms<br>
&gt; There are: 17798 VSites<br>
&gt; Max number of connections per atom is 59<br>
&gt; Total number of connections is 216528<br>
&gt; Max number of graph edges per atom is 4<br>
&gt; Total number of graph edges is 113666<br>
&gt;<br>
&gt; Constraining the starting coordinates (step 0)<br>
&gt;<br>
&gt; Constraining the coordinates at t0-dt (step 0)<br>
&gt; RMS relative constraint deviation after constraining: 3.77e-05<br>
&gt; Initial temperature: 299.838 K<br>
&gt;<br>
&gt;<br>
&gt; &gt; &gt; Recently I successfully installed the gromacs-4.0.5 mpi version.<br>
&gt; &gt; &gt; I could run in 8 cpu. but the speed is very slow.<br>
&gt; &gt; &gt; Total number of atoms in the system is 78424.<br>
&gt; &gt; &gt; while running all 8 cpu showing 95-100% CPU.<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; How to speed up the calculation.<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; Thanks<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;<br>
&gt; &gt; That&#39;s normal for a system that atoms/cpu ratio.<br>
&gt; &gt; What&#39;s your system and what mdp file are you using?<br>
&gt; &gt; --<br>
&gt; &gt; ------------------------------------------------------<br>
&gt; &gt; You haven&#39;t given us any diagnostic information. The problem could be<br>
&gt; &gt; that you&#39;re not running an MPI GROMACS (show us your configure line,<br>
&gt; &gt; your mdrun command line and the top 50 lines of your .log file).<br>
&gt; &gt;<br>
&gt; &gt; Mark<br>
&gt; &gt;<br>
&gt; &gt;<br>
&gt; &gt;<br>
&gt;<br>
&gt;<br>
<br>
<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <a href="http://lists.gromacs.org/pipermail/gmx-users/attachments/20090611/b60097b8/attachment.html" target="_blank">http://lists.gromacs.org/pipermail/gmx-users/attachments/20090611/b60097b8/attachment.html</a><br>

<br>
------------------------------<br>
<br>
_______________________________________________<br>
gmx-users mailing list<br>
<a href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a><br>
<a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/mailman/listinfo/gmx-users</a><br>
Please search the archive at <a href="http://www.gromacs.org/search" target="_blank">http://www.gromacs.org/search</a> before posting!<br>
<br>
End of gmx-users Digest, Vol 62, Issue 57<br>
*****************************************<br>
</blockquote></div><br>