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); }