Lines Matching defs:ww
689 * compute ww = (x-.5)*(log(x)-1) + .5*(log(2pi)-1) + 1/x*(P(1/x^2)))
1011 struct LDouble ss, ww;
1028 ww = large_gam(x, &m);
1029 w = ww.h + ww.l;
1035 ww = gam_n(i, x - (long double) i);
1036 return (ww.h + ww.l);
1119 /* Then compute ww = gamma(1+y), note that result scale to 2**m */
1122 ww = gam_n(j + 1, z);
1126 ww = large_gam(w, &m);
1130 ww = large_gam(w, &m);
1133 ww = gam_n(j, z);
1135 ww = large_gam(y, &m);
1136 t4 = ww.h + ww.l;
1140 ww.l = y * (ww.l - (t2 - ww.h)) + (y - t1) * t2;
1141 ww.h = t1 * t2;
1146 /* compute 1/(ss*ww) */
1148 t4 = ww.h + ww.l;
1152 z2 = ww.l - (t2 - ww.h); /* (t2,z2) = ww */
1153 t3 = t3 * t4; /* t3 = ss*ww */
1154 z3 = one / t3; /* z3 = 1/(ss*ww) */
1156 z5 = z1 * t4 + t1 * z2; /* (t5,z5) = ss*ww */
1157 t1 = CHOPPED((t3)); /* (t1,z1) = ss*ww */
1159 t2 = CHOPPED((z3)); /* leading 1/(ss*ww) */