search for: r_d_clog

Displaying 1 result from an estimated 1 matches for "r_d_clog".

Did you mean: r_dt_clog
2004 Apr 11
3
pcauchy precision (PR#6756)
...Submission from: (NULL) (65.213.85.129) Two things are wrong. 1. There is nan test outside IEEE_754. 2. The meat part of the function should really be something like... if (!lower_tail) x = -x; if (fabs (x) > 1) { double temp = atan (1 / x) / M_PI; return (x > 0) ? R_D_Clog (temp) : R_D_val (-temp); } else return R_D_val (0.5 + atan (x) / M_PI); ...instead of the current heavily truncated series expansion. The above is much simpler and more precise. (Thanks to Ian Smith for pointing the in-hindsight- obvious 1/x trick out.)