Lines Matching defs:u1
370 volatile double u1;
400 /* case 1: u1 = u rounded to 2^-43 absolute. Since u < 2^-8,
401 * u1 has at most 35 bits, and F*u1 is exact, as F has < 8 bits.
405 u1 = u + 513, u1 -= 513;
408 * u1 = u to 24 bits.
411 u1 = u, TRUNC(u1);
412 u2 = (2.0*(f - F*u1) - u1*f) * g;
413 /* u1 + u2 = 2f/(2F+f) to extra precision. */
416 /* (m*log2_hi+logF_head[j]+u1) + (m*log2_lo+logF_tail[j]+q); */
419 u1 += m*logF_head[N] + logF_head[j]; /* exact */
422 return (u1 + u2);
439 volatile double u1;
460 u1 = u + 513, u1 -= 513;
462 u1 = u, TRUNC(u1);
463 u2 = (2.0*(f - F*u1) - u1*f) * g;
465 u1 += m*logF_head[N] + logF_head[j];
469 r.a = u1 + u2; /* Only difference is here */
471 r.b = (u1 - r.a) + u2;
