<html>
<head>
<style><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 10pt;
font-family:Tahoma
}
--></style></head>
<body class='hmmessage'><div dir='ltr'>
Hi,<br><br>I swapped two arguments to a function call. I have fixed it and it should appear soon<br>in the public repository. You can find the fix below.<br><br>Cheers,<br><br>Berk<br><br>diff --git a/src/kernel/topshake.c b/src/kernel/topshake.c<br>index c5f3957..78961c5 100644<br>--- a/src/kernel/topshake.c<br>+++ b/src/kernel/topshake.c<br>@@ -204,7 +204,7 @@ void make_shake (t_params plist[],t_atoms *atoms,gpp_atomtype_t at,int nshake)<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; p.C1 = pr-&gt;param[i].C2;<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; add_param_to_list (&amp;(plist[F_CONSTR]),&amp;p);<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; } else {<br>-&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; copy_bond(pr,i,j++); <br>+&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; copy_bond(pr,j++,i); <br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; }<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; }<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; pr-&gt;nr = j;<br><br><br><div><div id="SkyDrivePlaceholder"></div>&gt; Date: Mon, 21 May 2012 11:13:46 +0200<br>&gt; From: sebastian.waltz@physik.uni-freiburg.de<br>&gt; To: gmx-users@gromacs.org<br>&gt; Subject: [gmx-users] constraints hbonds with 4.6<br>&gt; <br>&gt; Hi all together,<br>&gt; <br>&gt; I am testing the 4.6 dev version of gromacs with heterogeneous<br>&gt; acceleration (force field: amber99-ildn) on a 95 residue protein (all<br>&gt; atom model) in SPCE and run into major troubles using in the .mdp file<br>&gt; option constraints = hbonds with the default lincs parameters. The<br>&gt; system is exploding with the standard lincs warning (relative constraint<br>&gt; deviation after LINCS ......). Everything is working fine with the<br>&gt; option constraints = all-bonds or when I constrain all bonds with<br>&gt; hydrogen atoms manually in the .top file (which is a lot of effort).<br>&gt; Using gmxdump to search the .tpr file for the constraints I can't find<br>&gt; the bond constraints using the constraints = hbonds. Therefore I am<br>&gt; wondering if the very useful option is not included in the 4.6 version<br>&gt; of grompp or if there is something equal I don't know anything about.<br>&gt; <br>&gt; Thank a lot for any response<br>&gt; <br>&gt; Sebastian<br>&gt;  <br>&gt; <br>&gt; -- <br>&gt; gmx-users mailing list    gmx-users@gromacs.org<br>&gt; http://lists.gromacs.org/mailman/listinfo/gmx-users<br>&gt; Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/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/Support/Mailing_Lists<br></div>                                               </div></body>
</html>