gonzui

t2ex

bsd_source/lib/libc/src_bsd/math/e_lgammaf_r.c - 8.1KB - 232 lines

154:                 *signgamp = -1;
155:                 return - logf(-x);
156:             } else return - logf(x);
157:         }
162:             if(t==zero) return one/zero; /* -integer */
163:             nadj = logf(pi/fabsf(t*x));
164:             if(t<zero) *signgamp = -1;
172:             if(ix<=0x3f666666) {       /* lgamma(x) = lgamma(x+1)-log(x) */
173:                 r = - logf(x);
174:                 if(ix>=0x3f3b4a20) {y = one-x; i= 0;}
217:             case 3: z *= (y+(float)2.0);       /* FALLTHRU */
218:                     r += logf(z); break;
219:             }
221:         } else if (ix < 0x5c800000) {
222:             t = logf(x);
223:             z = one/x;
228:     /* 2**58 <= x <= inf */
229:             r =  x*(logf(x)-one);
230:         if(hx<0) r = nadj - r;