From cacf9dee745e70d4e8af8b2fcef53c8f39c1ef2d Mon Sep 17 00:00:00 2001
From: "P.M. Sutter" <sutter@iap.fr>
Date: Sun, 2 Jun 2013 20:48:07 -0500
Subject: [PATCH] corrected ellipticity calculation

---
 c_tools/stacking/pruneVoids.cpp | 6 ++++--
 1 file changed, 4 insertions(+), 2 deletions(-)

diff --git a/c_tools/stacking/pruneVoids.cpp b/c_tools/stacking/pruneVoids.cpp
index 11398d8..49c2ec1 100644
--- a/c_tools/stacking/pruneVoids.cpp
+++ b/c_tools/stacking/pruneVoids.cpp
@@ -616,9 +616,11 @@ int main(int argc, char **argv) {
     float c = sqrt(2.5*(gsl_vector_get(voids[iVoid].eval,0) +
                         gsl_vector_get(voids[iVoid].eval,1) - 
                         gsl_vector_get(voids[iVoid].eval,2)));
-    float ca = c/a;
+    float ca;
     float cb = c/b;
-    voids[iVoid].ellip = 1.0 - c/a;
+    if (a < c)  ca = a/c;
+    if (a >= c) ca = c/a;
+    voids[iVoid].ellip = fabs(1.0 - ca);
 
   } // iVoid