The R Project SVN R

Rev

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

Rev 68947 Rev 69670
Line 66... Line 66...
66 
    bk = (double *) R_alloc((size_t) nb, sizeof(double));
 66 
    bk = (double *) R_alloc((size_t) nb, sizeof(double));
67 
#endif
 67 
#endif
68 
    K_bessel(&x, &alpha, &nb, &ize, bk, &ncalc);
 68 
    K_bessel(&x, &alpha, &nb, &ize, bk, &ncalc);
69 
    if(ncalc != nb) {/* error input */
 69 
    if(ncalc != nb) {/* error input */
70 
      if(ncalc < 0)
 70 
      if(ncalc < 0)
71 
	MATHLIB_WARNING4(_("bessel_k(%g): ncalc (=%ld) != nb (=%ld); alpha=%g. Arg. out of range?\n"),
 71 
	MATHLIB_WARNING4(_("bessel_k(%g): ncalc (=%d) != nb (=%d); alpha=%g. Arg. out of range?\n"),
72 
			 x, ncalc, nb, alpha);
 72 
			 x, ncalc, nb, alpha);
73 
      else
 73 
      else
74 
	MATHLIB_WARNING2(_("bessel_k(%g,nu=%g): precision lost in result\n"),
 74 
	MATHLIB_WARNING2(_("bessel_k(%g,nu=%g): precision lost in result\n"),
75 
			 x, alpha+(double)nb-1);
 75 
			 x, alpha+(double)nb-1);
76 
    }
 76 
    }
Line 103... Line 103...
103 
    nb = 1+ (int)floor(alpha);/* nb-1 <= |alpha| < nb */
 103 
    nb = 1+ (int)floor(alpha);/* nb-1 <= |alpha| < nb */
104 
    alpha -= (double)(nb-1);
 104 
    alpha -= (double)(nb-1);
105 
    K_bessel(&x, &alpha, &nb, &ize, bk, &ncalc);
 105 
    K_bessel(&x, &alpha, &nb, &ize, bk, &ncalc);
106 
    if(ncalc != nb) {/* error input */
 106 
    if(ncalc != nb) {/* error input */
107 
      if(ncalc < 0)
 107 
      if(ncalc < 0)
108 
	MATHLIB_WARNING4(_("bessel_k(%g): ncalc (=%ld) != nb (=%ld); alpha=%g. Arg. out of range?\n"),
 108 
	MATHLIB_WARNING4(_("bessel_k(%g): ncalc (=%d) != nb (=%d); alpha=%g. Arg. out of range?\n"),
109 
			 x, ncalc, nb, alpha);
 109 
			 x, ncalc, nb, alpha);
110 
      else
 110 
      else
111 
	MATHLIB_WARNING2(_("bessel_k(%g,nu=%g): precision lost in result\n"),
 111 
	MATHLIB_WARNING2(_("bessel_k(%g,nu=%g): precision lost in result\n"),
112 
			 x, alpha+(double)nb-1);
 112 
			 x, alpha+(double)nb-1);
113 
    }
 113 
    }