The file 'particle_dump.dat' now also contains the particle masses.

examples/plot.py

@@ -39,17 +39,17 @@ print "Plotting..."
 # Read Quickshed accelerations
 data=loadtxt("particle_dump.dat")
 id = data[:,0]
-accx_e=data[:,4]
-accy_e=data[:,5]
-accz_e=data[:,6]
+accx_e=data[:,5]
+accy_e=data[:,6]
+accz_e=data[:,7]
 
-accx_bh=data[:,7]
-accy_bh=data[:,8]
-accz_bh=data[:,9]
+accx_bh=data[:,8]
+accy_bh=data[:,9]
+accz_bh=data[:,10]
 
-accx_new=data[:,10]
-accy_new=data[:,11]
-accz_new=data[:,12]
+accx_new=data[:,11]
+accy_new=data[:,12]
+accz_new=data[:,13]
 
 # Sort accelerations
 rank = argsort(id)

examples/test_bh_sorted.c

@@ -1900,10 +1900,10 @@ void test_bh(int N, int nr_threads, int runs, char *fileName) {
 
   /* Dump the particles to a file */
   file = fopen("particle_dump.dat", "w");
-  fprintf(file, "# x y z a_exact.x   a_exact.y    a_exact.z    a_legacy.x    "
+  fprintf(file, "# m x y z a_exact.x   a_exact.y    a_exact.z    a_legacy.x    "
                 "a_legacy.y    a_legacy.z    a_new.x     a_new.y    a_new.z\n");
   for (k = 0; k < N; ++k)
-    fprintf(file, "%d %e %e %e %e %e %e %e %e %e %e %e %e\n", parts[k].id,
+    fprintf(file, "%d %e %e %e %e %e %e %e %e %e %e %e %e %e\n", parts[k].id, parts[k].mass,
             parts[k].x[0], parts[k].x[1], parts[k].x[2], parts[k].a_exact[0],
             parts[k].a_exact[1], parts[k].a_exact[2], parts[k].a_legacy[0],
             parts[k].a_legacy[1], parts[k].a_legacy[2], parts[k].a[0],