Print this page

        

*** 879,889 **** /* * gamma(x+i) for 0 <= x < 1 */ static struct LDouble gam_n(int i, long double x) { ! struct LDouble rr, yy; long double r1, r2, t2, z, xh, xl, yh, yl, zh, z1, z2, zl, x5, wh, wl; /* compute yy = gamma(x+1) */ if (x > 0.2845L) { if (x > 0.6374L) { --- 879,889 ---- /* * gamma(x+i) for 0 <= x < 1 */ static struct LDouble gam_n(int i, long double x) { ! struct LDouble rr = {0.0L, 0.0L}, yy; long double r1, r2, t2, z, xh, xl, yh, yl, zh, z1, z2, zl, x5, wh, wl; /* compute yy = gamma(x+1) */ if (x > 0.2845L) { if (x > 0.6374L) {
*** 1030,1040 **** ww = large_gam(x, &m); w = ww.h + ww.l; return (scalbnl(w, m)); } ! if (hx > 0) { /* x from 0 to 8 */ i = (int) x; ww = gam_n(i, x - (long double) i); return (ww.h + ww.l); } /* INDENT OFF */ --- 1030,1040 ---- ww = large_gam(x, &m); w = ww.h + ww.l; return (scalbnl(w, m)); } ! if (hx > 0) { /* 0 < x < 8 */ i = (int) x; ww = gam_n(i, x - (long double) i); return (ww.h + ww.l); } /* INDENT OFF */