corrected ellipticity calculation

This commit is contained in:
P.M. Sutter 2013-06-02 20:48:07 -05:00
parent ef5384c54b
commit cacf9dee74

View file

@ -616,9 +616,11 @@ int main(int argc, char **argv) {
float c = sqrt(2.5*(gsl_vector_get(voids[iVoid].eval,0) + 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,1) -
gsl_vector_get(voids[iVoid].eval,2))); gsl_vector_get(voids[iVoid].eval,2)));
float ca = c/a; float ca;
float cb = c/b; 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 } // iVoid