Lines Matching refs:p34
258 int p34,
387 if (ind <= p34) {
390 if (ind + e4 < p34 + expmin) {
396 } else { // if (ind > p34)
398 // round R256 to p34 digits
399 x0 = ind - p34; // 1 <= x0 <= 34 as 35 <= ind <= 68
422 // the rounded result has p34 = 34 digits
439 // the number of digits in the significand is p34 = 34
444 ind = p34; // the number of decimal digits in the signifcand of res
453 // where q (significand) = ind <= p34
457 if (rnd_mode == ROUNDING_TO_NEAREST && (ind + e4) > (p34 + expmax)) {
487 } else if (x0 == ind) { // 1 <= x0 = ind <= p34 = 34
667 int p34 = P34; // used to modify the limit on the number of digits
1100 // where scale = min (p34-q3, (z_exp-p_exp) >> 49)
1101 scale = p34 - q3;
1293 if (q4 > p34) {
1295 // truncate C4 to p34 digits into res
1296 // x = q4-p34, 1 <= x <= 34 because 35 <= q4 <= 68
1297 x0 = q4 - p34;
1327 q4 = p34;
1330 } else { // if (q4 <= p34)
1334 if ((q4 + e4 <= p34 + expmax) && (e4 > expmax)) {
1362 if (q4 + e4 > p34 + expmax) {
1529 // where scale = min (p34-q4, (e4 - e3))
1530 scale = p34 - q4;
1614 // where scale = min (p34-q4, (p_exp-z_exp) >> 49)
1615 scale = p34 - q4;
1654 if (p34 <= delta - 1 || // Case (1')
1655 (p34 == delta && e3 + 6176 < p34 - q3)) { // Case (1''A)
1657 if ((q3 + e3) > (p34 + expmax) && p34 <= delta - 1) {
1658 // e3 > expmax implies p34 <= delta-1 and e3 > expmax is a necessary
1659 // condition for (q3 + e3) > (p34 + expmax)
1670 // q3 <= p34; if (q3 < p34) scale C3 up by 10^(p34-q3)
1671 scale = p34 - q3;
1689 // the coefficient in res has q3 + scale = p34 digits
1707 if (q3 < p34) { // the preferred exponent is z_exp - (p34 - q3)
1709 // where scale = min (p34-q3, z_exp-EMIN)
1710 scale = p34 - q3;
1734 if (scale + q3 < p34)
1855 // if (q3 + scale < p34)
1857 // else // if (q3 + scale = p34)
1880 if ((e3 == expmin && (q3 + scale) < p34) ||
1881 (e3 == expmin && (q3 + scale) == p34 &&
1902 } else if (p34 == delta) { // Case (1''B)
1904 // because Case (1''A) was treated above, e3 + 6176 >= p34 - q3
1905 // and C3 can be scaled up to p34 digits if needed
1907 // scale C3 to p34 digits if needed
1908 scale = p34 - q3; // 0 <= scale <= p34 - 1
1927 // now z_sign, z_exp, and res correspond to a z scaled to p34 = 34 digits
2284 } else if (((q3 <= delta && delta < p34 && p34 < delta + q4) || // Case (2)
2285 (q3 <= delta && delta + q4 <= p34) || // Case (3)
2286 (delta < q3 && p34 < delta + q4) || // Case (4)
2287 (delta < q3 && q3 <= delta + q4 && delta + q4 <= p34) || // Case (5)
2293 if ((q3 <= delta && delta < p34 && p34 < delta + q4) || // Case (2)
2294 (delta < q3 && p34 < delta + q4)) { // Case (4)
2296 // scale C3 up by scale = p34 - q3, 1 <= scale <= p34-1,
2299 scale = p34 - q3;
2300 x0 = delta + q4 - p34;
2329 // scale C3 up by scale = delta + q4 - q3, 1 <= scale <= p34-1,
2360 // round C4 to nearest to q4 - x0 digits, where x0 = delta + q4 - p34,
2361 // 1 <= x0 <= min (q4 - 1, 2 * p34 - 1) <=> 1 <= x0 <= min (q4 - 1, 67)
2362 // Also: 1 <= q4 - x0 <= p34 -1 => 1 <= q4 - x0 <= 33 (so the result of
2364 // x0 = delta + q4 - p34 (calculated before reaching case2_repeat)
2365 // because q3 + q4 - x0 <= P => x0 >= q3 + q4 - p34
2371 // 2 <= q4 <= 18, max(1, q3+q4-p34) <= x0 <= q4 - 1, 1 <= x0 <= 17
2382 // 19 <= q4 <= 38, max(1, q3+q4-p34) <= x0 <= q4 - 1, 1 <= x0 <= 37
2400 // 38 <= q4 <= 57, max(1, q3+q4-p34) <= x0 <= q4 - 1, 5 <= x0 <= 56
2424 // 58 <= q4 <= 68, max(1, q3+q4-p34) <= x0 <= q4 - 1, 25 <= x0 <= 67
2451 // the sum can result in [up to] p34 or p34 + 1 digits
2588 // the sum can result in [up to] p34 or p34 - 1 digits
2901 add_and_round (q3, q4, e4, delta, p34, z_sign, p_sign, C3, C4,
2918 if (p34 < q4 && q4 <= delta) { // Case (7)
2920 // truncate C4 to p34 digits into res
2921 // x = q4-p34, 1 <= x <= 34 because 35 <= q4 <= 68
2922 x0 = q4 - p34;
2954 // if C4 rounded to p34 digits is exact then the result is inexact,
2967 if (delta > p34 + 1) { // C3 < 1/2
2970 } else { // if (delta == p34 + 1)
3067 } else if ((q4 <= p34 && p34 <= delta) || // Case (8)
3068 (q4 <= delta && delta < p34 && p34 < delta + q3) || // Case (9)
3069 (q4 <= delta && delta + q3 <= p34) || // Case (10)
3070 (delta < q4 && q4 <= p34 && p34 < delta + q3) || // Case (13)
3071 (delta < q4 && q4 <= delta + q3 && delta + q3 <= p34) || // Case (14)
3072 (delta + q3 < q4 && q4 <= p34)) { // Case (18)
3104 } else if ((p34 <= delta && delta < q4 && q4 < delta + q3) || // Case (11)
3105 (delta < p34 && p34 < q4 && q4 < delta + q3)) { // Case (12)
3108 // 1 <= x0 <= q3 - 1 <= p34 - 1
3197 ind = nr_digits256 (R256); // ind >= p34
3198 // if R256 is sum, then ind > p34; if R256 is a difference, then
3199 // ind >= p34; this means that we can calculate the result rounded to
3204 if (ind < p34) {
3206 } else if (ind == p34) {
3212 } else { // if (ind > p34)
3214 // round R256 to p34 digits
3215 x0 = ind - p34; // 1 <= x0 <= 34 as 35 <= ind <= 68
3226 // round to p34 digits; the result fits in 113 bits
3252 // the rounded result has p34 = 34 digits
3332 // the number of digits in the significand is p34 = 34
3343 ind = p34; // the number of decimal digits in the signifcand of res
3348 // where q (significand) = ind = p34
3353 && (ind + e4) > (p34 + expmax)) {
3390 } else if (x0 == ind) { // 1 <= x0 = ind <= p34 = 34
3540 } else if ((p34 <= delta && delta + q3 <= q4) || // Case (15)
3541 (delta < p34 && p34 < delta + q3 && delta + q3 <= q4) || //Case (16)
3542 (delta + q3 <= p34 && p34 < q4)) { // Case (17)
3547 add_and_round (q3, q4, e4, delta, p34, z_sign, p_sign, C3, C4,