From 64adea4a9bf1d0260daf4750d2afe07c52edd1af Mon Sep 17 00:00:00 2001
From: nhamaus <hamaus@illinois.edu>
Date: Tue, 18 Dec 2012 16:28:20 +0100
Subject: [PATCH] Enforced pbc for void-barycenters

---
 c_tools/stacking/pruneVoids.cpp | 3 +++
 1 file changed, 3 insertions(+)

diff --git a/c_tools/stacking/pruneVoids.cpp b/c_tools/stacking/pruneVoids.cpp
index f418c5d..eceab99 100644
--- a/c_tools/stacking/pruneVoids.cpp
+++ b/c_tools/stacking/pruneVoids.cpp
@@ -330,6 +330,9 @@ int main(int argc, char **argv) {
     voids[iVoid].barycenter[0] += voids[iVoid].center[0];
     voids[iVoid].barycenter[1] += voids[iVoid].center[1];
     voids[iVoid].barycenter[2] += voids[iVoid].center[2];
+    voids[iVoid].barycenter[0] = fmod(voids[iVoid].barycenter[0], boxLen[0]);
+    voids[iVoid].barycenter[1] = fmod(voids[iVoid].barycenter[1], boxLen[1]);
+    voids[iVoid].barycenter[2] = fmod(voids[iVoid].barycenter[2], boxLen[2]);
 
     // compute central density
     centralRad = voids[iVoid].radius/args_info.centralRadFrac_arg;