config/libbid/bid128_fma.c

663   int e1, e2, e3, e4;
1134   e3 = (z_exp >> 49) - 6176; // unbiased exponent of z
1524         if (e3 < e4) {
1525           // if (e3 < e4) the preferred exponent is e3
1527           // where scale = min (p34-q4, (e4 - e3))
1529           ind = e4 - e3;
1648   delta = q3 + e3 - q4 - e4;
1653         (p34 == delta && e3 + 6176 < p34 - q3)) { // Case (1''A)
1655       if ((q3 + e3) > (p34 + expmax) && p34 <= delta - 1) {
1656         // e3 > expmax implies p34 <= delta-1 and e3 > expmax is a necessary
1657         // condition for (q3 + e3) > (p34 + expmax)
1689           e3 = e3 - scale;
1695         		       e3, &res, pfpsf);
1709         ind = e3 + 6176;
1730         e3 = e3 - scale;
1736         res.w[1] = z_sign | ((UINT64) (e3 + 6176) << 49) | C3.w[1];
1755           // ok from above e3 = (z_exp >> 49) - 6176;
1803           } else if ((e3 == expmin) ||
1825             e3 = e3 - 1;
1826             res.w[1] = z_sign | ((UINT64) (e3 + 6176) << 49) | res.w[1];
1828           if (e3 == expmin) {
1878       if ((e3 == expmin && (q3 + scale) < p34) ||
1879           (e3 == expmin && (q3 + scale) == p34 &&
1891         		     e3, &res, pfpsf);
1902       // because Case (1''A) was treated above, e3 + 6176 >= p34 - q3
1924       e3 = e3 - scale;
1999             e3 = e3 + 1;
2001             z_exp = ((UINT64) (e3 + 6176) << 49) & MASK_EXP;
2021         if (e3 > expmax && rnd_mode == ROUNDING_TO_NEAREST) {
2037         		       e3, &res, pfpsf);
2067           if (e3 > expmax) {
2077         			   is_midpoint_gt_even, e3, &res,
2094         			 is_midpoint_gt_even, e3, &res, pfpsf);
2098           e3 = (z_exp >> 49) - 6176;
2099           if (e3 > expmin) {
2108               e3 = e3 - 1;
2151         	e3 = e3 - 1;
2166         	// 100...0 * 10^e3 - 0.41 * 10^e3 =
2167         	// 0999...9.59 * 10^e3 -> rounds to 99...96*10^(e3-1)
2168         	// note that (e3 > expmin}
2177         	e3 = e3 - 1;
2195         	if (e3 > expmax) {
2205         				 is_midpoint_gt_even, e3, &res,
2218         	  z_sign | ((UINT64) (e3 + 6176) << 49) | res.w[1];
2224         			       is_midpoint_gt_even, e3, &res,
2230           } else { // if (e3 = emin)
2231             // if e3 = expmin the result is also tiny (the condition for
2263                   is_midpoint_gt_even, e3, &res,
2267           } // end e3 = emin
2352       // e3 is already calculated
2353       e3 = e3 - scale;
2354       // now res = C3 * 10^scale and e3 = e3 - scale
2526           e3 = e3 + 1;
2552         	e3++;
2568         	// the exponent is max (e3, expmin)
2594         if (e3 > expmin && ((res.w[1] < 0x0000314dc6448d93ull ||
2602           // first restore e3, otherwise it will be too small
2603           e3 = e3 + scale;
2643               e3++;
2657               // the exponent is max (e3, expmin)
2675       if (e3 == expmin) { // and if significand < 10^33 => result is tiny
2681       } else if (e3 < expmin) {
2684         x0 = expmin - e3;
2721         // digit (q3 = 1) equal to 1 (C3 = 1), e3 is expmin (e3 = expmin),
2814         e3 = e3 + x0;
2839         e3 = e3 + 1;
2841       res.w[1] = z_sign | ((UINT64) (e3 + 6176) << 49) | res.w[1];
2843       if (rnd_mode == ROUNDING_TO_NEAREST && e3 > expmax) {
2853         		     e3, &res, pfpsf);
2875         // (17) if we swap (C3, C4), (q3, q4), (e3, e4), (z_sign, p_sign)
2887         ind = e3;
2888         e3 = e4;
3079       // swap (C3, C4), (q3, q4), (e3, 34), (z_sign, p_sign), (z_exp, p_exp)
3091       ind = e3;
3092       e3 = e4;
3105       // round C3 to nearest to q3 - x0 digits, where x0 = e4 - e3,
3107       x0 = e4 - e3; // or x0 = delta + q3 - q4
3131       e3 = e3 + x0; // this is e4