<html>
<head>
<style>
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 10pt;
font-family:Verdana
}
</style>
</head>
<body class='hmmessage'>
Hi,<br><br>I fixed the bug for version 4.0.3.<br>There was also a bug with writing the forces at the end<br>with nstfout=0.<br><br>But if you are minimizing for starting a md simulation,<br>you do not need cg, steepest descents is enough<br>(with em_tol of 1000 or 100).<br><br>Note that in 4.0 you can now minimize with constraints,<br>box for water and other molecules.<br><br>If you really need to run parallel cg now, you can use<br>the -pd switch of mdrun.<br><br>Berk<br><br><br><br>&gt; Date: Sat, 22 Nov 2008 00:37:18 +0530<br>&gt; From: suman@sscu.iisc.ernet.in<br>&gt; To: gmx-users@gromacs.org<br>&gt; Subject: [gmx-users] Gromacs-4.0.2: energy minimization error in parallel        version<br>&gt; <br>&gt; Dear all,<br>&gt; <br>&gt; I am trying to use mdrun_mpi for energy minimization over 2/4/8 nodes.<br>&gt; But in every case the run gets aborted with the following error:<br>&gt; <br>&gt; -------------------------------------------------------<br>&gt; Program mdrun_mpi, VERSION 4.0.2<br>&gt; Source code file: minimize.c, line: 404<br>&gt; <br>&gt; Software inconsistency error:<br>&gt; state mismatch in do_em_step<br>&gt; -------------------------------------------------------<br>&gt; <br>&gt; What could be the reason behind this? As far as I remember I have not<br>&gt; faced this problem in 3.3.3 version. Not that I regularly use parallel<br>&gt; version for energy minimization, but it would be nice to be able to<br>&gt; sort out this problem. The PR or MD run goes fine in the same<br>&gt; settings.<br>&gt; <br>&gt; The mdp file is:<br>&gt; <br>&gt; title           = cg+steep energymin spce watbox<br>&gt; cpp             = /lib/cpp<br>&gt; define          = -DFLEXIBLE<br>&gt; <br>&gt; integrator      = cg<br>&gt; nsteps          = 5000<br>&gt; nstcgsteep      = 1000<br>&gt; comm_mode       = linear<br>&gt; nstcomm         = 1<br>&gt; <br>&gt; emtol           = 1.0<br>&gt; emstep          = 0.01<br>&gt; <br>&gt; nstxout         = 100<br>&gt; nstvout         = 100<br>&gt; nstfout         = 100<br>&gt; nstlog          = 100<br>&gt; nstenergy       = 100<br>&gt; energygrps      = protein sol<br>&gt; <br>&gt; nstlist         = 5<br>&gt; ns_type         = grid<br>&gt; pbc             = xyz<br>&gt; rlist           = 0.9<br>&gt; <br>&gt; coulombtype     = PME<br>&gt; rcoulomb        = 0.9<br>&gt; <br>&gt; vdwtype         = cut-off<br>&gt; rvdw            = 1.2<br>&gt; <br>&gt; fourier_spacing = 0.12<br>&gt; pme_order       = 6<br>&gt; ewald_rtol      = 1e-6<br>&gt; optimize_fft    = yes<br>&gt; <br>&gt; constraints     = none<br>&gt; <br>&gt; <br>&gt; <br>&gt; Thanks,<br>&gt; Suman.<br>&gt; _______________________________________________<br>&gt; gmx-users mailing list    gmx-users@gromacs.org<br>&gt; http://www.gromacs.org/mailman/listinfo/gmx-users<br>&gt; Please search the archive at http://www.gromacs.org/search before posting!<br>&gt; Please don't post (un)subscribe requests to the list. Use the <br>&gt; www interface or send it to gmx-users-request@gromacs.org.<br>&gt; Can't post? Read http://www.gromacs.org/mailing_lists/users.php<br><br /><hr />Express yourself instantly with MSN Messenger! <a href='http://clk.atdmt.com/AVE/go/onm00200471ave/direct/01/' target='_new'>MSN Messenger</a></body>
</html>