gonzui

t2ex

bsd_source/lib/libc/src_bsd/math/e_lgamma_r.c - 12.2KB - 304 lines

218:                 *signgamp = -1;
219:                 return - log(-x);
220:             } else return - log(x);
221:         }
226:             if(t==zero) return one/zero; /* -integer */
227:             nadj = log(pi/fabs(t*x));
228:             if(t<zero) *signgamp = -1;
236:             if(ix<=0x3feccccc) {       /* lgamma(x) = lgamma(x+1)-log(x) */
237:                 r = - log(x);
238:                 if(ix>=0x3FE76944) {y = one-x; i= 0;}
281:             case 3: z *= (y+2.0);      /* FALLTHRU */
282:                     r += log(z); break;
283:             }
285:         } else if (ix < 0x43900000) {
286:             t = log(x);
287:             z = one/x;
292:     /* 2**58 <= x <= inf */
293:             r =  x*(log(x)-one);
294:         if(hx<0) r = nadj - r;