Added debugging code to locate where velocities get trashed.

1 files changed, 16 insertions, 0 deletions

diff --git a/src/tests/md_openmp_util.c b/src/tests/md_openmp_util.c
index 99ab60e..41ccff5 100644
--- a/src/tests/md_openmp_util.c
+++ b/src/tests/md_openmp_util.c
@@ -269,6 +269,21 @@ int main ( int argc, char *argv[] )
                 printf("Error adding data. %s: %d\n", __FILE__, __LINE__);
                 exit(1);
             }
+#if 1
+	    {
+	      int j,i;
+    for ( j = 0; j < np; j++ )
+    {
+      printf("vel in md for %d %d:",step,j);
+        for ( i = 0; i < nd; i++ )
+        {
+	  printf (" %g",vel[i+j*nd]);
+	}
+	printf("\n");
+    }
+	    }
+#endif
+
             if(tng_util_vel_write(traj, step, vel) != TNG_SUCCESS)
             {
                 printf("Error adding data. %s: %d\n", __FILE__, __LINE__);
@@ -742,5 +757,6 @@ void update ( int np, int nd, float pos[], float vel[], float f[],
         }
     }
 
+
     return;
 }