Lines Matching refs:a
156 return(__exp__D(u.a, u.b));
190 u.a -= one;
191 v.a = (x -= .5);
192 TRUNC(v.a);
193 v.b = x - v.a;
194 t.a = v.a*u.a; /* t = (x-.5)*(log(x)-1) */
195 t.b = v.b*u.a + x*u.b;
198 u.a = lns2pi_hi + t.b; u.a += t.a;
199 u.b = t.a - u.a;
217 return (yy.a + yy.b);
219 r.a = y;
220 TRUNC(r.a);
221 yy.a = r.a - one;
223 yy.b = r.b = y - yy.a;
225 for (ym1 = y-one; ym1 > LEFT + x0; y = ym1--, yy.a--) {
226 t = r.a*yy.a;
227 r.b = r.a*yy.b + y*r.b;
228 r.a = t;
229 TRUNC(r.a);
230 r.b += (t - r.a);
234 y = r.b*(yy.a + yy.b) + r.a*yy.b;
235 y += yy.a*r.a;
251 xx.a = (t + x), TRUNC(xx.a);
252 xx.b = x - xx.a; xx.b += t; xx.b += d;
255 x = xx.a + xx.b;
257 xx.a = x, TRUNC(xx.a);
258 xx.b = x - xx.a;
263 d = r.a/x, TRUNC(d);
264 r.a -= d*xx.a; r.a -= d*xx.b; r.a += r.b;
265 return (d + r.a/x);
282 t.a = z, TRUNC(t.a); /* t ~= z + c */
283 t.b = (z - t.a) + c;
284 t.b *= (t.a + z);
285 q = (t.a *= t.a); /* t = (z+c)^2 */
286 TRUNC(t.a);
287 t.b += (q - t.a);
288 r.a = p, TRUNC(r.a); /* r = P/Q */
289 r.b = p - r.a;
290 t.b = t.b*p + t.a*r.b + a0_lo;
291 t.a *= r.a; /* t = (z+c)^2*(P/Q) */
292 r.a = t.a + a0_hi, TRUNC(r.a);
293 r.b = ((a0_hi-r.a) + t.a) + t.b;
326 lg.a -= lsine.a; /* exact (opposite signs) */
328 y = -(lg.a + lg.b);
329 z = (y + lg.a) + lg.b;