Lines Matching defs:u1
381 volatile double u1;
419 /* case 1: u1 = u rounded to 2^-43 absolute. Since u < 2^-8,
420 * u1 has at most 35 bits, and F*u1 is exact, as F has < 8 bits.
424 u1 = u + 513, u1 -= 513;
427 * u1 = u to 24 bits.
430 u1 = u, TRUNC(u1);
431 u2 = (2.0*(f - F*u1) - u1*f) * g;
432 /* u1 + u2 = 2f/(2F+f) to extra precision. */
435 /* (m*log2_hi+logF_head[j]+u1) + (m*log2_lo+logF_tail[j]+q); */
438 u1 += m*logF_head[N] + logF_head[j]; /* exact */
441 return (u1 + u2);
453 volatile double u1;
474 u1 = u + 513, u1 -= 513;
476 u1 = u, TRUNC(u1);
477 u2 = (2.0*(f - F*u1) - u1*f) * g;
479 u1 += m*logF_head[N] + logF_head[j];
483 r.a = u1 + u2; /* Only difference is here */
485 r.b = (u1 - r.a) + u2;