PunctureTracker: add READS to workaround poison for position arrrays

Author: lwJi

PunctureTracker/schedule.ccl

@@ -23,6 +23,7 @@ SCHEDULE PunctureTracker_Setup IN PunctureTracker_SetupGroup
   OPTIONS: GLOBAL
   READS: pt_loc
   READS: pt_vel
+  READS: BoxInBox::positions
   WRITES: BoxInBox::positions
 } "Calculate initial location of punctures"
 
@@ -31,6 +32,7 @@ SCHEDULE PunctureTracker_Track AT evol AFTER ODESolvers_Solve
   LANG: C
   OPTIONS: GLOBAL
   READS: ADMBaseX::shift(everywhere)
+  READS: BoxInBox::positions
   WRITES: pt_loc
   WRITES: pt_vel
   WRITES: BoxInBox::positions

PunctureTracker/src/puncture_tracker.cxx

@@ -83,17 +83,11 @@ extern "C" void PunctureTracker_Setup(CCTK_ARGUMENTS) {
   if (track_boxes) {
     const std::array<std::vector<CCTK_REAL>, Loop::dim> &location =
         g_punctures->getLocation();
-    for (int n = 0; n < 3; ++n) { // BoxInBox currently hardcodes position[3]
-      if (n < nPunctures) {
-        CCTK_VINFO("Writing punc coords to box %d.", n);
-        position_x[n] = location[0][n];
-        position_y[n] = location[1][n];
-        position_z[n] = location[2][n];
-      } else {
-        position_x[n] = 0.0;
-        position_y[n] = 0.0;
-        position_z[n] = 0.0;
-      }
+    for (int n = 0; n < nPunctures; ++n) {
+      CCTK_VINFO("Writing punc coords to box %d.", n);
+      position_x[n] = location[0][n];
+      position_y[n] = location[1][n];
+      position_z[n] = location[2][n];
     }
   }
 }
@@ -193,16 +187,10 @@ extern "C" void PunctureTracker_Track(CCTK_ARGUMENTS) {
   }
 
   if (track_boxes) {
-    for (int i = 0; i < 3; ++i) { // BoxInBox currently hardcodes position[3]
-      if (i < nPunctures) {
-        position_x[i] = location[0][i];
-        position_y[i] = location[1][i];
-        position_z[i] = location[2][i];
-      } else {
-        position_x[i] = 0.0;
-        position_y[i] = 0.0;
-        position_z[i] = 0.0;
-      }
+    for (int i = 0; i < nPunctures; ++i) {
+      position_x[i] = location[0][i];
+      position_y[i] = location[1][i];
+      position_z[i] = location[2][i];
     }
   }
 }