Hi GMX people!
I'm trying to check that the energies I can calculate using the trajectory
output are the same as those given by g_energy. I'm running a diatomic for only
a couple of steps and already there are inconsistencies. I suspect that one
possible factor could be that the leap-frog algorithm uses trajectory data
(i.e. positions and velocities) different from those printed to the trajectory
file.
Could you please direct me to the file and line in the source code where the
kinetic and potential energies are calculated, just before they are printed to
the edr file?
Thanks in advance,
--
Inon Sharony
ינון שרוני
+972(3)6407634
atto.TAU.ac.IL/~inonshar
Please consider your environmental responsibility before printing this e-mail.
_______________________________________________
gmx-users mailing list [email protected]
http://lists.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the
www interface or send it to [email protected].
Can't post? Read http://www.gromacs.org/mailing_lists/users.php