<html>
<head>
<style>
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 10pt;
font-family:Verdana
}
</style>
</head>
<body class='hmmessage'>
Thinking about this, I will probably not fix it in 4.0,<br>but only in git master for 4.1.<br><br>There are several issues with pulling "infinitely" long/far.<br>One of them is that you can no longer determine the virial contribution.<br>Also, grompp should check that you do not pressure coupl in the direction<br>that you pull.<br><br>Berk<br><br>&gt; Date: Fri, 4 Sep 2009 12:45:04 +0200<br>&gt; From: alexander.herz@mytum.de<br>&gt; To: gmx-users@gromacs.org<br>&gt; Subject: Re: [gmx-users] t_trxframe speed units<br>&gt; <br>&gt; Great!<br>&gt; <br>&gt; Berk Hess schrieb:<br>&gt; &gt; Ah, now I see what the problem is.<br>&gt; &gt; The pull code works correctly when the distance is shorter or longer<br>&gt; &gt; than half the box length. But when the distance crosses half the box<br>&gt; &gt; length,<br>&gt; &gt; the code might take the "wrong" periodic distance.<br>&gt; &gt; I'll see if I can fix this.<br>&gt; &gt;<br>&gt; &gt; Berk<br>&gt; &gt;<br>&gt; &gt; ------------------------------------------------------------------------<br>&gt; &gt; From: gmx3@hotmail.com<br>&gt; &gt; To: gmx-users@gromacs.org<br>&gt; &gt; Subject: RE: [gmx-users] t_trxframe speed units<br>&gt; &gt; Date: Fri, 4 Sep 2009 10:25:09 +0200<br>&gt; &gt;<br>&gt; &gt; Hi,<br>&gt; &gt;<br>&gt; &gt; Could you mail me your input files?<br>&gt; &gt;<br>&gt; &gt; Berk<br>&gt; &gt;<br>&gt; &gt; &gt; Date: Fri, 4 Sep 2009 10:05:27 +0200<br>&gt; &gt; &gt; From: alexander.herz@mytum.de<br>&gt; &gt; &gt; To: gmx-users@gromacs.org<br>&gt; &gt; &gt; Subject: Re: [gmx-users] t_trxframe speed units<br>&gt; &gt; &gt;<br>&gt; &gt; &gt; Hi<br>&gt; &gt; &gt;<br>&gt; &gt; &gt; yes, I tried pull_start=yes (you also suggested trying constraint<br>&gt; &gt; &gt; instead of umbrella, which I tried as well).<br>&gt; &gt; &gt; my box is rectangular.<br>&gt; &gt; &gt;<br>&gt; &gt; &gt; Alex<br>&gt; &gt; &gt;<br>&gt; &gt; &gt; Berk Hess schrieb:<br>&gt; &gt; &gt; &gt; Hi,<br>&gt; &gt; &gt; &gt;<br>&gt; &gt; &gt; &gt; The only thing I suggested was using pull_start=yes,<br>&gt; &gt; &gt; &gt; otherwise you will start at a "random" distance and you will have<br>&gt; &gt; enormous forces.<br>&gt; &gt; &gt; &gt; Did you try that?<br>&gt; &gt; &gt; &gt;<br>&gt; &gt; &gt; &gt; Do you have a rectangular or a triclinic box?<br>&gt; &gt; &gt; &gt;<br>&gt; &gt; &gt; &gt; Berk<br>&gt; &gt; &gt; &gt;<br>&gt; &gt; &gt; &gt;<br>&gt; &gt; &gt; &gt;&gt; Date: Thu, 3 Sep 2009 18:00:44 +0200<br>&gt; &gt; &gt; &gt;&gt; From: alexander.herz@mytum.de<br>&gt; &gt; &gt; &gt;&gt; To: gmx-users@gromacs.org<br>&gt; &gt; &gt; &gt;&gt; Subject: Re: [gmx-users] t_trxframe speed units<br>&gt; &gt; &gt; &gt;&gt;<br>&gt; &gt; &gt; &gt;&gt; Hey,<br>&gt; &gt; &gt; &gt;&gt;<br>&gt; &gt; &gt; &gt;&gt; I tried everything you suggested and a couple more things and<br>&gt; &gt; still it<br>&gt; &gt; &gt; &gt;&gt; is not possible to get a reasonable speed<br>&gt; &gt; &gt; &gt;&gt; for the pulled groups (always way beyond the 0.01nm/ps set in the mdp<br>&gt; &gt; &gt; &gt;&gt; file and mostly containing some high frequence patttern probably<br>&gt; &gt; related<br>&gt; &gt; &gt; &gt;&gt; to pbc). (I tried: constraint, pbcatoms, direction,position and<br>&gt; &gt; using no<br>&gt; &gt; &gt; &gt;&gt; ref group).<br>&gt; &gt; &gt; &gt;&gt;<br>&gt; &gt; &gt; &gt;&gt; Can I send you a small sample system with the problem or is there<br>&gt; &gt; &gt; &gt;&gt; something else I can do?<br>&gt; &gt; &gt; &gt;&gt;<br>&gt; &gt; &gt; &gt;&gt; Thx,<br>&gt; &gt; &gt; &gt;&gt; Alex<br>&gt; &gt; &gt; &gt;&gt;<br>&gt; &gt; &gt; &gt;&gt; Berk Hess schrieb:<br>&gt; &gt; &gt; &gt;&gt;<br>&gt; &gt; &gt; &gt;&gt;&gt; But you can not talk about popping back, since there is pbc.<br>&gt; &gt; &gt; &gt;&gt;&gt; What matters is how the pull code determines the distance.<br>&gt; &gt; &gt; &gt;&gt;&gt; At least in 4.0.4 and 4.0.5 it add the rate*t before doing pbc,<br>&gt; &gt; &gt; &gt;&gt;&gt; such that it pulls smoothly even if the distance is larger than<br>&gt; &gt; the box<br>&gt; &gt; &gt; &gt;&gt;&gt;<br>&gt; &gt; &gt; &gt;&gt;&gt;<br>&gt; &gt; ------------------------------------------------------------------------<br>&gt; &gt; &gt; &gt;&gt;&gt;<br>&gt; &gt; &gt; &gt;&gt;&gt; _______________________________________________<br>&gt; &gt; &gt; &gt;&gt;&gt; gmx-users mailing list gmx-users@gromacs.org<br>&gt; &gt; &gt; &gt;&gt;&gt; http://lists.gromacs.org/mailman/listinfo/gmx-users<br>&gt; &gt; &gt; &gt;&gt;&gt; Please search the archive at http://www.gromacs.org/search<br>&gt; &gt; before posting!<br>&gt; &gt; &gt; &gt;&gt;&gt; Please don't post (un)subscribe requests to the list. Use the<br>&gt; &gt; &gt; &gt;&gt;&gt; www interface or send it to gmx-users-request@gromacs.org.<br>&gt; &gt; &gt; &gt;&gt;&gt; Can't post? Read http://www.gromacs.org/mailing_lists/users.php<br>&gt; &gt; &gt;<br>&gt; &gt; &gt; _______________________________________________<br>&gt; &gt; &gt; gmx-users mailing list gmx-users@gromacs.org<br>&gt; &gt; &gt; http://lists.gromacs.org/mailman/listinfo/gmx-users<br>&gt; &gt; &gt; Please search the archive at http://www.gromacs.org/search before<br>&gt; &gt; posting!<br>&gt; &gt; &gt; Please don't post (un)subscribe requests to the list. Use the<br>&gt; &gt; &gt; www interface or send it to gmx-users-request@gromacs.org.<br>&gt; &gt; &gt; Can't post? Read http://www.gromacs.org/mailing_lists/users.php<br>&gt; &gt;<br>&gt; &gt; ------------------------------------------------------------------------<br>&gt; &gt; Express yourself instantly with MSN Messenger! MSN Messenger<br>&gt; &gt; &lt;http://clk.atdmt.com/AVE/go/onm00200471ave/direct/01/&gt;<br>&gt; &gt; ------------------------------------------------------------------------<br>&gt; &gt; What can you do with the new Windows Live? Find out<br>&gt; &gt; &lt;http://www.microsoft.com/windows/windowslive/default.aspx&gt;<br>&gt; &gt; ------------------------------------------------------------------------<br>&gt; &gt;<br>&gt; &gt; _______________________________________________<br>&gt; &gt; gmx-users mailing list    gmx-users@gromacs.org<br>&gt; &gt; http://lists.gromacs.org/mailman/listinfo/gmx-users<br>&gt; &gt; Please search the archive at http://www.gromacs.org/search before posting!<br>&gt; &gt; Please don't post (un)subscribe requests to the list. Use the <br>&gt; &gt; www interface or send it to gmx-users-request@gromacs.org.<br>&gt; &gt; Can't post? Read http://www.gromacs.org/mailing_lists/users.php<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/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>