The R Project SVN R

(root)/trunk/src/nmath/bessel_i.c @ 28362 – Rev 88401 → 89291

Rev 88401 | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 88401		Rev 89291
Line 243...		Line 243...
243	`------------------------------------------------------------------- */`	243	`------------------------------------------------------------------- */`
244	`if (nb > 0 && x >= 0. && (0. <= nu && nu < 1.) &&`	244	`if (nb > 0 && x >= 0. && (0. <= nu && nu < 1.) &&`
245	`(1 <= ize && ize <= 2) ) {`	245	`(1 <= ize && ize <= 2) ) {`
246		246
247	`ncalc = nb;`	247	`ncalc = nb;`
248	`if(ize == 1 && x > exparg_BESS) {`	248	`if(ize == 1 && x > exparg_BESS) { // x > 709`
249	`for(k=1; k <= *nb; k++)`	249	`for(k=1; k <= *nb; k++)`
250	`bi[k]=ML_POSINF; /* the limit is = Inf */`	250	`bi[k]=ML_POSINF; /* the limit is = Inf */`
251	`return;`	251	`return;`
252	`}`	252	`}`
253	`if(ize == 2 && x > xlrg_BESS_IJ) {`	253	`if(ize == 2 && x > xlrg_BESS_IJ) {`