* Vol.30, pp.1-15, 1988.
*/
-
+
#include <config.h>
#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);
}
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;