The R Project SVN R

Rev

Rev 59170 | Rev 69128 | Go to most recent revision | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 59170 Rev 60268
Line 30... Line 30...
30 
double dnchisq(double x, double df, double ncp, int give_log)
 30 
double dnchisq(double x, double df, double ncp, int give_log)
31 
{
 31 
{
32 
    const static double eps = 5e-15;
 32 
    const static double eps = 5e-15;
33 
 
 33 
 
34 
    double i, ncp2, q, mid, dfmid, imax;
 34 
    double i, ncp2, q, mid, dfmid, imax;
35 
    long double sum, term;
 35 
    LDOUBLE sum, term;
36 
 
 36 
 
37 
#ifdef IEEE_754
 37 
#ifdef IEEE_754
38 
    if (ISNAN(x) || ISNAN(df) || ISNAN(ncp))
 38 
    if (ISNAN(x) || ISNAN(df) || ISNAN(ncp))
39 
	return x + df + ncp;
 39 
	return x + df + ncp;
40 
#endif
 40 
#endif