docs

[pspp] / lib / tukey / ptukey.c

diff --git a/lib/tukey/ptukey.c b/lib/tukey/ptukey.c
index 9a89104cf90e78f206a002e611d850b5ad922eb6..07923d98bd7b5a363d7969847a333950243d5a6e 100644 (file)
--- a/lib/tukey/ptukey.c
+++ b/lib/tukey/ptukey.c
@@ -59,7 +59,7 @@
  *    Vol.30, pp.1-15, 1988.
  */
 
- 
+
 
 #include <config.h>
 
@@ -71,7 +71,7 @@
 #include <assert.h>
 #include <math.h>
 
-#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.