projects
/
pspp
/ blobdiff
commit
grep
author
committer
pickaxe
?
search:
re
summary
|
shortlog
|
log
|
commit
|
commitdiff
|
tree
raw
|
inline
| side by side
Merge remote-tracking branch 'origin/master' into sheet
[pspp]
/
lib
/
tukey
/
ptukey.c
diff --git
a/lib/tukey/ptukey.c
b/lib/tukey/ptukey.c
index c86ef23613c0766e90681dfdabc25c026257c58f..2af46cf89378fa6840502d3e986c6e38211420b2 100644
(file)
--- a/
lib/tukey/ptukey.c
+++ b/
lib/tukey/ptukey.c
@@
-59,7
+59,7
@@
* Vol.30, pp.1-15, 1988.
*/
* Vol.30, pp.1-15, 1988.
*/
-
+
#include <config.h>
#include <config.h>
@@
-84,13
+84,13
@@
#define ME_PRECISION 8
#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);
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);
}
return gsl_cdf_gaussian_P (x - mu, sigma);
}
@@
-371,9
+371,7
@@
ptukey (double q, double rr, double cc, double df, int lower_tail, int log_p)
double ans, f2, f21, f2lf, ff4, otsum, qsqz, rotsum, t1, twa1, ulen, wprb;
int i, j, jj;
double ans, f2, f21, f2lf, ff4, otsum, qsqz, rotsum, t1, twa1, ulen, wprb;
int i, j, jj;
-#ifdef IEEE_754
- abort (! (ISNAN (q) || ISNAN (rr) || ISNAN (cc) || ISNAN (df)));
-#endif
+ assert (! (isnan (q) || isnan (rr) || isnan (cc) || isnan (df)));
if (q <= 0)
return R_DT_0;
if (q <= 0)
return R_DT_0;