[gmx-developers] velocity: gmxdump vs. interface, precision?
Mehmet Suezen
suzen at theochem.tu-muenchen.de
Thu Aug 29 15:16:58 CEST 2002
Hi,
I've just compared results of interface and gmxdump. *But they are not
quite the same.
In case of positions, they are very close.
Any idea, what could be the reason?
cheers,
Mehmet
gmxdump gives;
v (1536x3):
v[ 0]={ 7.36503e-01, -8.57067e-01, 7.37718e-01}
v[ 1]={ 6.71643e-01, 3.62874e-01, 6.77159e-01}
v[ 2]={ 1.50657e-01, -8.60839e-01, -9.45431e-02}
v[ 3]={ 1.38135e+00, -7.58055e-01, -8.14016e-01}
v[ 4]={ 4.33980e-01, -2.16724e+00, 1.24819e+00}
v[ 5]={ 1.77014e-01, 3.58265e-01, -2.77815e+00}
v[ 6]={-3.01506e-02, 7.85357e-01, -5.22182e-01}
v[ 7]={ 1.20010e+00, -1.38198e-01, -4.33350e-01}
v[ 8]={-1.24630e+00, 9.64708e-01, 5.24720e-01}
v[ 9]={ 3.95388e-01, -4.95744e-01, 6.13505e-01}
and the corresponding output of the interface;
-----------
if(trn.v_size) {
vxdyn[myindex]=fr.v[i][XX];
vydyn[myindex]=fr.v[i][YY];
vzdyn[myindex]=fr.v[i][ZZ];
printf("nframe %d Velocity x %f Velocity y %f Velocity
z %f \n",
nframe,vxdyn[myindex],vydyn[myindex],vzdyn[myindex]);
}
----------
nframe 0 Velocity x 0.804340 Velocity y -0.719031 Velocity z 0.780158
nframe 0 Velocity x 0.519643 Velocity y 0.370223 Velocity z 0.704696
nframe 0 Velocity x -0.090481 Velocity y -0.715459 Velocity z -0.181907
nframe 0 Velocity x 1.041339 Velocity y -0.620405 Velocity z -0.673146
nframe 0 Velocity x 0.514215 Velocity y -2.010526 Velocity z 1.178636
nframe 0 Velocity x 0.326832 Velocity y 0.273683 Velocity z -2.329099
nframe 0 Velocity x -0.004167 Velocity y 0.498045 Velocity z -0.372272
nframe 0 Velocity x 0.780044 Velocity y 0.008229 Velocity z -0.480785
nframe 0 Velocity x -1.115026 Velocity y 0.908866 Velocity z 0.331761
More information about the gromacs.org_gmx-developers
mailing list