The R Project SVN R

Rev

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

Rev 19500 Rev 22929
Line 35... Line 35...
35
 
35
 
36
double bessel_i(double x, double alpha, double expo)
36
double bessel_i(double x, double alpha, double expo)
37
{
37
{
38
    long nb, ncalc, ize;
38
    long nb, ncalc, ize;
39
    double *bi;
39
    double *bi;
-
 
40
#ifndef MATHLIB_STANDALONE
40
    char *vmax;
41
    char *vmax;
-
 
42
#endif
41
 
43
 
42
#ifdef IEEE_754
44
#ifdef IEEE_754
43
    /* NaNs propagated correctly */
45
    /* NaNs propagated correctly */
44
    if (ISNAN(x) || ISNAN(alpha)) return x + alpha;
46
    if (ISNAN(x) || ISNAN(alpha)) return x + alpha;
45
#endif
47
#endif