Lines Matching defs:ww
739 * t = rr.h*ww.h+hln2pi_h ... exact
740 * w = (hln2pi_l + ((x-0.5)*ww.l+rr.l*ww.h)) + p
1275 * compute ww = (x-.5)*(log(x)-1) + .5*(log(2pi)-1) + 1/x*(P(1/x^2)))
1533 struct Double ss, ww;
1554 ww = large_gam(x, &m);
1555 w = ww.h + ww.l;
1561 ww = gam_n(i, x - (double) i);
1562 return (ww.h + ww.l);
1632 /* Then compute ww = gamma(1+y), note that result scale to 2**m */
1635 ww = gam_n(j + 1, z);
1639 ww = large_gam(w, &m);
1643 ww = large_gam(w, &m);
1646 ww = gam_n(j, z);
1648 ww = large_gam(y, &m);
1649 t4 = ww.h + ww.l;
1653 ww.l = y * (ww.l - (t2 - ww.h)) + (y - t1) * t2;
1654 ww.h = t1 * t2;
1659 /* compute 1/(ss*ww) */
1661 t4 = ww.h + ww.l;
1665 z2 = ww.l - (t2 - ww.h); /* (t2,z2) = ww */
1666 t3 = t3 * t4; /* t3 = ss*ww */
1667 z3 = one / t3; /* z3 = 1/(ss*ww) */
1669 z5 = z1 * t4 + t1 * z2; /* (t5,z5) = ss*ww */
1670 t1 = (double) ((float) t3); /* (t1,z1) = ss*ww */
1672 t2 = (double) ((float) z3); /* leading 1/(ss*ww) */