Print this page
*** 1404,1414 ****
/* INDENT ON */
/* gamma(x+i) for 0 <= x < 1 */
static struct Double
gam_n(int i, double x) {
! struct Double rr, yy;
double r1, r2, t2, z, xh, xl, yh, yl, zh, z1, z2, zl, x5, wh, wl;
/* compute yy = gamma(x+1) */
if (x > 0.2845) {
if (x > 0.6374) {
--- 1404,1414 ----
/* INDENT ON */
/* gamma(x+i) for 0 <= x < 1 */
static struct Double
gam_n(int i, double x) {
! struct Double rr = {0.0L, 0.0L}, yy;
double r1, r2, t2, z, xh, xl, yh, yl, zh, z1, z2, zl, x5, wh, wl;
/* compute yy = gamma(x+1) */
if (x > 0.2845) {
if (x > 0.6374) {
*** 1556,1566 ****
ww = large_gam(x, &m);
w = ww.h + ww.l;
__HI(w) += m << 20;
return (w);
}
! if (hx > 0) { /* x from 0 to 8 */
i = (int) x;
ww = gam_n(i, x - (double) i);
return (ww.h + ww.l);
}
--- 1556,1566 ----
ww = large_gam(x, &m);
w = ww.h + ww.l;
__HI(w) += m << 20;
return (w);
}
! if (hx > 0) { /* 0 < x < 8 */
i = (int) x;
ww = gam_n(i, x - (double) i);
return (ww.h + ww.l);
}