Displaying 1 result from an estimated 1 matches for "l_rt".
Did you mean:
_rt
2004 Apr 19
2
pgeom accuracy (PR#6792)
Full_Name: Morten Welinder
Version: snapshot
OS:
Submission from: (NULL) (65.213.85.218)
This should fix the remaining two 1-p cancellation issues.
double l_rt = log1p (-p) * (x + 1);
if (log_p)
return R_DT_Clog (l_rt);
else
return lower_tail ? -expm1 (l_rt) : exp (l_rt);