The R Project SVN R

Rev

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

Rev 68947 Rev 69670
Line 72... Line 72...
72 
    bi = (double *) R_alloc((size_t) nb, sizeof(double));
 72 
    bi = (double *) R_alloc((size_t) nb, sizeof(double));
73 
#endif
 73 
#endif
74 
    I_bessel(&x, &alpha, &nb, &ize, bi, &ncalc);
 74 
    I_bessel(&x, &alpha, &nb, &ize, bi, &ncalc);
75 
    if(ncalc != nb) {/* error input */
 75 
    if(ncalc != nb) {/* error input */
76 
	if(ncalc < 0)
 76 
	if(ncalc < 0)
77 
	    MATHLIB_WARNING4(_("bessel_i(%g): ncalc (=%ld) != nb (=%ld); alpha=%g. Arg. out of range?\n"),
 77 
	    MATHLIB_WARNING4(_("bessel_i(%g): ncalc (=%d) != nb (=%d); alpha=%g. Arg. out of range?\n"),
78 
			     x, ncalc, nb, alpha);
 78 
			     x, ncalc, nb, alpha);
79 
	else
 79 
	else
80 
	    MATHLIB_WARNING2(_("bessel_i(%g,nu=%g): precision lost in result\n"),
 80 
	    MATHLIB_WARNING2(_("bessel_i(%g,nu=%g): precision lost in result\n"),
81 
			     x, alpha+(double)nb-1);
 81 
			     x, alpha+(double)nb-1);
82 
    }
 82 
    }
Line 117... Line 117...
117 
    nb = 1 + (int)na;/* nb-1 <= alpha < nb */
 117 
    nb = 1 + (int)na;/* nb-1 <= alpha < nb */
118 
    alpha -= (double)(nb-1);
 118 
    alpha -= (double)(nb-1);
119 
    I_bessel(&x, &alpha, &nb, &ize, bi, &ncalc);
 119 
    I_bessel(&x, &alpha, &nb, &ize, bi, &ncalc);
120 
    if(ncalc != nb) {/* error input */
 120 
    if(ncalc != nb) {/* error input */
121 
	if(ncalc < 0)
 121 
	if(ncalc < 0)
122 
	    MATHLIB_WARNING4(_("bessel_i(%g): ncalc (=%ld) != nb (=%ld); alpha=%g. Arg. out of range?\n"),
 122 
	    MATHLIB_WARNING4(_("bessel_i(%g): ncalc (=%d) != nb (=%d); alpha=%g. Arg. out of range?\n"),
123 
			     x, ncalc, nb, alpha);
 123 
			     x, ncalc, nb, alpha);
124 
	else
 124 
	else
125 
	    MATHLIB_WARNING2(_("bessel_i(%g,nu=%g): precision lost in result\n"),
 125 
	    MATHLIB_WARNING2(_("bessel_i(%g,nu=%g): precision lost in result\n"),
126 
			     x, alpha+(double)nb-1);
 126 
			     x, alpha+(double)nb-1);
127 
    }
 127 
    }