From f9fc9e899055db4fa1d6c2abf331e6b0df082e30 Mon Sep 17 00:00:00 2001 From: "Paul M. Sutter" Date: Tue, 25 Mar 2025 14:49:41 -0400 Subject: [PATCH] cleaned up procedure for checking for 'obviously bad' voids; no longer filters on density contrast since this only affects edge voids anyway --- c_source/pruning/pruneVoids.cpp | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/c_source/pruning/pruneVoids.cpp b/c_source/pruning/pruneVoids.cpp index bbf3878..d84cef6 100644 --- a/c_source/pruning/pruneVoids.cpp +++ b/c_source/pruning/pruneVoids.cpp @@ -277,7 +277,7 @@ int main(int argc, char **argv) { volNorm = numPartTot/(boxLen[0]*boxLen[1]*boxLen[2]); printf(" VOL NORM = %f\n", volNorm); - printf(" CENTRAL DEN = %f\n", args.maxCentralDen_arg); + printf(" MAX CENTRAL DEN = %e\n", args.maxCentralDen_arg); fread(&dummy, 1, 4, fp); fread(temp, numPartTot, 4, fp); @@ -823,8 +823,8 @@ int main(int argc, char **argv) { // toss out voids that are obviously wrong int iGood = 0; for (iVoid = 0; iVoid < voids.size(); iVoid++) { - if (voids[iVoid].densCon > 1.e4 || isnan(voids[iVoid].vol) || - isinf(voids[iVoid].vol)) { + + if (isnan(voids[iVoid].vol) || isinf(voids[iVoid].vol)) { numWrong++; } else { voids[iGood++] = voids[iVoid];