Displaying 1 result from an estimated 1 matches for "me_rang".
Did you mean:
me_range
2002 Oct 17
1
underflow handling in besselK (PR#2179)
...re is an
*underflow*). The answer may be as simple as returning 0 rather than +Inf
in this case.
If this is appropriate, the following patch should do the job:
*** bessel_k.c Thu Oct 17 13:55:55 2002
--- bessel_k.c.dist Thu Oct 17 13:48:53 2002
***************
*** 216,224 ****
ML_ERROR(ME_RANGE);
*ncalc = *nb;
for(i=0; i < *nb; i++)
! if (ex<=0)
! bk[i] = ML_POSINF;
! else bk[i]=0.0;
return;
}
k = 0;
--- 216,222 ----
ML_ERROR(ME_RANGE);
*ncalc = *nb;
for(i=0; i < *nb; i++)
! bk[i] = ML_POSINF;
return;...