Xavier,<br><br><div><span class="gmail_quote"><br>
</span><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">I basically rerun a simulation, both in xtc and trr<br>format,<br>using slightly different tpr file, varying the gro file
<br>used to construct it, and different gmx version, 3.2 and<br>3.3. From those rerun I extracted and compared the bonded<br>and non-bonded energies.<br><br>The most stricking is that I do not see differences<br>between<br>
the two tpr files, only small differences between the two<br>gmx versions, and more significant changes when reruning<br>on trr and xtc files.<br><br>Those results fit with the already identified bug in<br>gmx3.2<br>with rerun fixed for
gmx3.3. The only problem is the lack<br>of<br>differences between the two gmx versions, and the fact<br>that<br>xtc and trr file give such differences from the "correct"<br>energies extracted from the original energy file.
</blockquote><div><br>
<br>
I, too, am slightly confused about what exactly the differences are
you're seeing. Could you perhaps make a little table or something?
Partly, I just don't get what you mean about the difference between the
two versions... Are you taking trajectories you originally ran with 3.2
and reprocessing them with 3.3? I think I had tried this at one point
and concluded that it wasn't a good idea since there are some things
the two versions do differently, so I shouldn't necessarily get the
same energies even using the same (full precision) trajectory and tpr
files. <br>
<br>
The other issue is using mdrun -rerun with xtc files. Since these are
reduced precision, you should expect to get different energies even
reprocessing a 3.3 trajectory with 3.3 using the same tpr file, as some
of the energies are quite sensitive to distances. How much difference
you will see depends on the precision of your xtc files. I *never* use
xtc files for reprocessing for this reason.<br>
<br>
Also, make sure you set the list and energies to update every step in
mdrun -rerun, otherwise you will see different energies for other
reasons.<br>
<br>
If you do all that, and see *significantly* different energies, do let
me know, as I'd be concerned in that case. Minor differences due to
finite precision are to be expected.<br>
<br>
David<br>
</div><br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">Best<br>XAvier<br><br>>><br>>> The energies indicate that:<br>>> 0) I see a systematic difference with the original
<br>>>energies (not<br>>> significant), but using<br>>> the trr file helps getting better.<br>>> 1) I do not see ANY difference between the rerun using<br>>>the two tpr file<br>>> (Unexpected ?!?!)
<br>>> 2) I do not see significant differences between gmx3.2<br>>>and gmx3.3 (New<br>>> ?!?!) in all cases.<br>>> 3) I see significant differences in energies when<br>>>generated with xtc and
<br>>> trr files (New ?!?!).<br>>> trr files help.<br>>><br>>> All this seems in agreement with the past remarks of<br>>>David and Mickael<br>>> until I compare<br>>> the energies obtained with the twp tpr files ?! Am I
<br>>>missing something ?<br>>> Note that my system is quite big (in cubic (!) box about<br>>>12.0, 12.0,<br>>> 8.0, in nm).<br>>><br>>> Does this make sense to any of you ?<br>><br>
> I'm having difficulty understanding what you've done,<br>>but some things that might help you understand what<br>>you've done are<br>><br>> 1) The .trr and .xtc versions of the trajectories will<br>>correspond to slightly different positions and so you
<br>>should expect slightly different energies to result<br>> 2) You need to use .edr files with grompp/tpbconv for<br>>exact restarts to be exact<br>> 3) Different versions of the mdrun code should be<br>>expected to be slightly different. A default value of a
<br>>parameter may have changed, of course, but in general the<br>>execution order of the two versions need not be the same.<br>>If you add up a few thousand floating point numbers in a<br>>slightly different order you will get different results.
<br>>Evaluation of LJ and Coulomb terms requires this.<br>> 4) Different versions may produce identical .tpr files,<br>>or might not. You can check this with gmxdump.<br>><br>> Mark<br>> _______________________________________________
<br>> gmx-users mailing list <a href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a><br>> <a href="http://www.gromacs.org/mailman/listinfo/gmx-users">http://www.gromacs.org/mailman/listinfo/gmx-users</a>
<br>> Please don't post (un)subscribe requests to the list.<br>>Use the www interface or send it to<br>><a href="mailto:gmx-users-request@gromacs.org">gmx-users-request@gromacs.org</a>.<br>> Can't post? Read<br>
><a href="http://www.gromacs.org/mailing_lists/users.php">http://www.gromacs.org/mailing_lists/users.php</a><br><br>_______________________________________________<br>gmx-users mailing list <a href="mailto:gmx-users@gromacs.org">
gmx-users@gromacs.org</a><br><a href="http://www.gromacs.org/mailman/listinfo/gmx-users">http://www.gromacs.org/mailman/listinfo/gmx-users</a><br>Please don't post (un)subscribe requests to the list. Use the<br>www interface or send it to
<a href="mailto:gmx-users-request@gromacs.org">gmx-users-request@gromacs.org</a>.<br>Can't post? Read <a href="http://www.gromacs.org/mailing_lists/users.php">http://www.gromacs.org/mailing_lists/users.php</a><br></blockquote>
</div><br>