X-Git-Url: https://pintos-os.org/cgi-bin/gitweb.cgi?a=blobdiff_plain;ds=sidebyside;f=lib%2Ftukey%2Fqtukey.c;h=54c174be2074be0378ffaed41c6903106d47e92f;hb=refs%2Fheads%2Fpivot-titles-5;hp=ec146001420eca7d4876e6af12a813efe3f206da;hpb=cb00cc9f9dc738f2850cd4b5a9776031f87de8a2;p=pspp diff --git a/lib/tukey/qtukey.c b/lib/tukey/qtukey.c index ec14600142..54c174be20 100644 --- a/lib/tukey/qtukey.c +++ b/lib/tukey/qtukey.c @@ -143,8 +143,8 @@ static double qinv(double p, double c, double v) ps = 0.5 - 0.5 * p; yi = sqrt (log (1.0 / (ps * ps))); - t = yi + (((( yi * p4 + p3) * yi + p2) * yi + p1) * yi + p0) - / (((( yi * q4 + q3) * yi + q2) * yi + q1) * yi + q0); + t = yi + ((((yi * p4 + p3) * yi + p2) * yi + p1) * yi + p0) + / ((((yi * q4 + q3) * yi + q2) * yi + q1) * yi + q0); if (v < vmax) t += (t * t * t + t) / v / 4.0; q = c1 - c2 * t; if (v < vmax) q += -c3 / v + c4 * t / v; @@ -192,14 +192,14 @@ double qtukey(double p, double rr, double cc, double df, } /* df must be > 1 ; there must be at least two values */ - /* ^^ + /* ^^ JMD: The comment says 1 but the code says 2. Which is correct? */ assert (df >= 2); assert (rr >= 1); assert (cc >= 2); - + R_Q_P01_boundaries (p, 0, ML_POSINF);