<!DOCTYPE html PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html>
<head>
  <meta content="text/html;charset=ISO-8859-1" http-equiv="Content-Type">
</head>
<body bgcolor="#ffffff" text="#000000">
Hi Berk,<br>
<br>
Thanks for pointing out my mistake!<br>
What I worry about is how big the deviation is after several
continuations? If it is unpredictable, it is a sad news.<br>
<br>
regards,<br>
Baofu Qiao<br>
<br>
Berk Hess wrote:
<blockquote cite="mid:COL113-W22D2942A621F4B6F7EB31A8E130@phx.gbl"
 type="cite">
  <pre wrap="">Hi,

Your command line is incorrect.
-reprod should be used without yes.
-reprod is yes
-noreprod is no

mdrun -cpi has no effect on your results at all.
But always remember that MD is chaotic, there is no point in trying
to reproduce trajectories exactly, unless for a very special purpose.

Berk

  </pre>
  <blockquote type="cite">
    <pre wrap="">Date: Thu, 30 Jul 2009 16:53:13 +0200
From: <a class="moz-txt-link-abbreviated" href="mailto:qiaobf@gmail.com">qiaobf@gmail.com</a>
To: <a class="moz-txt-link-abbreviated" href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a>
Subject: Re: [gmx-users] question about continuation using tpbconv and mdrun        -cpi

Hi Mark,

Sorry to trouble you again!

I made two tests, by using one-processor and one 32-processors on the
same cluster. I used 4000SPC/E waters (OPLS-AA ff.) The former
(one-processor) gives the exactly the same potential. However, the
latter still shown some deviation of the potential. When I use "-reprod
yes", the potential is -190098+-457.199 for "Comparison", and
-190088+-467.225 for "Reference". When I use "-dlb no", -190116+-476.512
for "Comparison", and -190114+-483.749 for "Reference".

The following is the job script lines when using "-reprod yes":
********************************************************************************
# Reference
grompp -f full.mdp -c initial.gro -p system.top -o full
mpirun -np 32 mdrun -deffnm full -reprod yes

# Comparison
grompp -f full.mdp -c initial.gro -p system.top -o full2
mpirun -np 32 mdrun -deffnm full2 -reprod yes -cpt 2 -maxh 0.2
mpirun -np 32 mdrun -deffnm full2 -reprod yes -cpi full2.cpt -append
********************************************************************************

Do you have some further ideas?

best wishes,
Baofu Qiao


Mark Abraham wrote:
    </pre>
    <blockquote type="cite">
      <pre wrap="">Baofu Qiao wrote:
      </pre>
      <blockquote type="cite">
        <pre wrap="">Hi Mark,

Thanks!
Because the maximum time for one single job is set to be 24hours on the
cluster I'm using, I want to make sure which is the best way to continue
the gmx jobs. I wonder how strong effect "mdrun -cpi" has?  From the
introduction of mdrun, it seems that there are some EXTRA energy frames,
but for the trajectory file (.xtc), there is no extra frames? Am I
right?

"mdrun -h
--&gt; The contents will be binary identical (unless you use dynamic load
balancing), but for technical reasons there might be some extra energy
frames when using checkpointing (necessary for restarts without
appending)."
        </pre>
      </blockquote>
      <pre wrap="">The intent with GROMACS 4.x is for a user to be able to construct a
.tpr with a very long simulation time, and perhaps constrain mdrun
with -maxh (or rely on the cluster killing the job), and to use the
information in the checkpoint file to restart correctly and perhaps to
then use mdrun -append so that when the simulation is running
smoothly, only one set of files needs to exist. Thus one doesn't need
to trouble with using tpbconv correctly, crashes can restart
transparently, etc. The old-style approach still works, however.

Obviously you should (be able to) verify with mdrun -reprod that
whatever approach you use when you construct your job scripts leads to
simulations that are in principle reproducible. For production, don't
use -reprod because you will want the better speed from dynamic load
balancing, etc.

Mark
_______________________________________________
gmx-users mailing list    <a class="moz-txt-link-abbreviated" href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a>
<a class="moz-txt-link-freetext" href="http://lists.gromacs.org/mailman/listinfo/gmx-users">http://lists.gromacs.org/mailman/listinfo/gmx-users</a>
Please search the archive at <a class="moz-txt-link-freetext" href="http://www.gromacs.org/search">http://www.gromacs.org/search</a> before
posting!
Please don't post (un)subscribe requests to the list. Use the www
interface or send it to <a class="moz-txt-link-abbreviated" href="mailto:gmx-users-request@gromacs.org">gmx-users-request@gromacs.org</a>.
Can't post? Read <a class="moz-txt-link-freetext" href="http://www.gromacs.org/mailing_lists/users.php">http://www.gromacs.org/mailing_lists/users.php</a>

      </pre>
    </blockquote>
    <pre wrap="">
_______________________________________________
gmx-users mailing list    <a class="moz-txt-link-abbreviated" href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a>
<a class="moz-txt-link-freetext" href="http://lists.gromacs.org/mailman/listinfo/gmx-users">http://lists.gromacs.org/mailman/listinfo/gmx-users</a>
Please search the archive at <a class="moz-txt-link-freetext" href="http://www.gromacs.org/search">http://www.gromacs.org/search</a> before posting!
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to <a class="moz-txt-link-abbreviated" href="mailto:gmx-users-request@gromacs.org">gmx-users-request@gromacs.org</a>.
Can't post? Read <a class="moz-txt-link-freetext" href="http://www.gromacs.org/mailing_lists/users.php">http://www.gromacs.org/mailing_lists/users.php</a>
    </pre>
  </blockquote>
  <pre wrap=""><!---->
_________________________________________________________________
Express yourself instantly with MSN Messenger! Download today it's FREE!
<a class="moz-txt-link-freetext" href="http://messenger.msn.click-url.com/go/onm00200471ave/direct/01/">http://messenger.msn.click-url.com/go/onm00200471ave/direct/01/</a>
  </pre>
  <pre wrap="">
<hr size="4" width="90%">
_______________________________________________
gmx-users mailing list    <a class="moz-txt-link-abbreviated" href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a>
<a class="moz-txt-link-freetext" href="http://lists.gromacs.org/mailman/listinfo/gmx-users">http://lists.gromacs.org/mailman/listinfo/gmx-users</a>
Please search the archive at <a class="moz-txt-link-freetext" href="http://www.gromacs.org/search">http://www.gromacs.org/search</a> before posting!
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to <a class="moz-txt-link-abbreviated" href="mailto:gmx-users-request@gromacs.org">gmx-users-request@gromacs.org</a>.
Can't post? Read <a class="moz-txt-link-freetext" href="http://www.gromacs.org/mailing_lists/users.php">http://www.gromacs.org/mailing_lists/users.php</a></pre>
</blockquote>
<br>
</body>
</html>