The R Project SVN R

Rev

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

Rev 13183 Rev 13435
Line 27... Line 27...
27
#include "nmath.h"
27
#include "nmath.h"
28
#include "dpq.h"
28
#include "dpq.h"
29
 
29
 
30
double dnchisq(double x, double df, double lambda, int give_log)
30
double dnchisq(double x, double df, double lambda, int give_log)
31
{
31
{
32
    const static double eps = 5e-15;
32
    const double eps = 5e-15;
33
 
33
 
34
    double i, lambda2, term, sum, q, mid, dfmid, imax, errorbound;
34
    double i, lambda2, term, sum, q, mid, dfmid, imax, errorbound;
35
 
35
 
36
#ifdef IEEE_754
36
#ifdef IEEE_754
37
    if (ISNAN(x) || ISNAN(df) || ISNAN(lambda))
37
    if (ISNAN(x) || ISNAN(df) || ISNAN(lambda))