gonzui

t2ex

bsd_source/lib/libc/src_bsd/math/b_tgamma.c - 10.9KB - 346 lines

179:         v.a = (x -= .5);
180:         TRUNC(v.a);
181:         v.b = x - v.a;
208:         r.a = y;
209:         TRUNC(r.a);
210:         yy.a = r.a - one;
217:                 r.a = t;
218:                 TRUNC(r.a);
219:                 r.b += (t - r.a);
238:                 t = x;
239:                 TRUNC(t);
240:                 d = (t+x)*(x-t);
242:                 xx.a = (t + x);
243:                 TRUNC(xx.a);
244:                 xx.b = x - xx.a; xx.b += t; xx.b += d;
249:                 xx.a =  x;
250:                 TRUNC(xx.a);
251:                 xx.b = x - xx.a;
256:         d = r.a/x;
257:         TRUNC(d);
258:         r.a -= d*xx.a; r.a -= d*xx.b; r.a += r.b;
277:         t.a = z;
278:         TRUNC(t.a);                    /* t ~= z + c */
279:         t.b = (z - t.a) + c;
281:         q = (t.a *= t.a);              /* t = (z+c)^2 */
282:         TRUNC(t.a);
283:         t.b += (q - t.a);
284:         r.a = p;
285:         TRUNC(r.a);                    /* r = P/Q */
286:         r.b = p - r.a;
289:         r.a = t.a + a0_hi;
290:         TRUNC(r.a);
291:         r.b = ((a0_hi-r.a) + t.a) + t.b;