The R Project SVN R

Rev

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

Rev 68947 Rev 69670
Line 75... Line 75...
75 
    bj = (double *) R_alloc((size_t) nb, sizeof(double));
 75 
    bj = (double *) R_alloc((size_t) nb, sizeof(double));
76 
#endif
 76 
#endif
77 
    J_bessel(&x, &alpha, &nb, bj, &ncalc);
 77 
    J_bessel(&x, &alpha, &nb, bj, &ncalc);
78 
    if(ncalc != nb) {/* error input */
 78 
    if(ncalc != nb) {/* error input */
79 
      if(ncalc < 0)
 79 
      if(ncalc < 0)
80 
	MATHLIB_WARNING4(_("bessel_j(%g): ncalc (=%ld) != nb (=%ld); alpha=%g. Arg. out of range?\n"),
 80 
	MATHLIB_WARNING4(_("bessel_j(%g): ncalc (=%d) != nb (=%d); alpha=%g. Arg. out of range?\n"),
81 
			 x, ncalc, nb, alpha);
 81 
			 x, ncalc, nb, alpha);
82 
      else
 82 
      else
83 
	MATHLIB_WARNING2(_("bessel_j(%g,nu=%g): precision lost in result\n"),
 83 
	MATHLIB_WARNING2(_("bessel_j(%g,nu=%g): precision lost in result\n"),
84 
			 x, alpha+(double)nb-1);
 84 
			 x, alpha+(double)nb-1);
85 
    }
 85 
    }
Line 122... Line 122...
122 
    nb = 1 + (int)na; /* nb-1 <= alpha < nb */
 122 
    nb = 1 + (int)na; /* nb-1 <= alpha < nb */
123 
    alpha -= (double)(nb-1); // ==> alpha' in [0, 1)
 123 
    alpha -= (double)(nb-1); // ==> alpha' in [0, 1)
124 
    J_bessel(&x, &alpha, &nb, bj, &ncalc);
 124 
    J_bessel(&x, &alpha, &nb, bj, &ncalc);
125 
    if(ncalc != nb) {/* error input */
 125 
    if(ncalc != nb) {/* error input */
126 
      if(ncalc < 0)
 126 
      if(ncalc < 0)
127 
	MATHLIB_WARNING4(_("bessel_j(%g): ncalc (=%ld) != nb (=%ld); alpha=%g. Arg. out of range?\n"),
 127 
	MATHLIB_WARNING4(_("bessel_j(%g): ncalc (=%d) != nb (=%d); alpha=%g. Arg. out of range?\n"),
128 
			 x, ncalc, nb, alpha);
 128 
			 x, ncalc, nb, alpha);
129 
      else
 129 
      else
130 
	MATHLIB_WARNING2(_("bessel_j(%g,nu=%g): precision lost in result\n"),
 130 
	MATHLIB_WARNING2(_("bessel_j(%g,nu=%g): precision lost in result\n"),
131 
			 x, alpha+(double)nb-1);
 131 
			 x, alpha+(double)nb-1);
132 
    }
 132 
    }