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 */