X-Git-Url: https://pintos-os.org/cgi-bin/gitweb.cgi?a=blobdiff_plain;ds=sidebyside;f=lib%2Ftukey%2Fptukey.c;h=07923d98bd7b5a363d7969847a333950243d5a6e;hb=7a99952fb25df826b584a3652cac9fd8b74c307a;hp=9a89104cf90e78f206a002e611d850b5ad922eb6;hpb=cb00cc9f9dc738f2850cd4b5a9776031f87de8a2;p=pspp diff --git a/lib/tukey/ptukey.c b/lib/tukey/ptukey.c index 9a89104cf9..07923d98bd 100644 --- a/lib/tukey/ptukey.c +++ b/lib/tukey/ptukey.c @@ -59,7 +59,7 @@ * Vol.30, pp.1-15, 1988. */ - + #include @@ -71,7 +71,7 @@ #include #include -#define R_D__0 (log_p ? ML_NEGINF : 0.) /* 0 */ +#define R_D__0 (log_p ? -INFINITY : 0.) /* 0 */ #define R_D__1 (log_p ? 0. : 1.) /* 1 */ #define R_DT_0 (lower_tail ? R_D__0 : R_D__1) /* 0 */ #define R_DT_1 (lower_tail ? R_D__1 : R_D__0) /* 1 */ @@ -84,13 +84,13 @@ #define ME_PRECISION 8 -static inline double +static inline double pnorm(double x, double mu, double sigma, int lower_tail, int log_p) { assert (lower_tail == 1); assert (log_p == 0); assert (sigma == 1.0); - + return gsl_cdf_gaussian_P (x - mu, sigma); } @@ -272,7 +272,6 @@ wprob (double w, double rr, double cc) double ptukey (double q, double rr, double cc, double df, int lower_tail, int log_p) { - const double ML_NEGINF = -1.0 / 0.0; /* function ptukey() [was qprob() ]: q = value of studentized range @@ -327,7 +326,7 @@ ptukey (double q, double rr, double cc, double df, int lower_tail, int log_p) f(x) = .5 + erf(x / sqrt(2)) / 2 for x > 0 - f(x) = erfc( -x / sqrt(2)) / 2 for x < 0 + f(x) = erfc(-x / sqrt(2)) / 2 for x < 0 where f(x) is standard normal c. d. f.