config/libbid/bid128_fma.c

258                UINT64 p_sign,
336   if (p_sign == z_sign) { // R256 = C4 + R256
340   } else { // if (p_sign != z_sign) { // R256 = C4 - R256
353       // flip p_sign too, because the result has the sign of z
354       p_sign = z_sign;
365         p_sign = 0x0000000000000000ull;
367         p_sign = 0x8000000000000000ull;
372       res.w[1] = p_sign | ((UINT64) (e4 + 6176) << 49);
382   // the exact result is (-1)^p_sign * R256 * 10^e4 where q (R256) = ind;
391     res.w[1] = p_sign | ((UINT64) (e4 + 6176) << 49) | R256.w[1];
429       // (-1)^p_sign * significand * 10^0
430       P128.w[1] = p_sign | 0x3040000000000000ull | R128.w[1];
443     res.w[1] = p_sign | ((UINT64) (e4 + 6176) << 49) | R128.w[1]; // RN
450   // res = (-1)^p_sign * significand * 10^e4,
456     res.w[1] = p_sign | 0x7800000000000000ull;
467     //   (-1)^p_sign * significand_res * 10^e4
482       res.w[1] = p_sign | 0x0000000000000000ull;
527       res.w[1] = p_sign | res.w[1];
556         p_sign | ((UINT64) (e4 + 6176) << 49) | (res.w[1] & MASK_COEFF);
657   UINT64 x_sign, y_sign, z_sign, p_sign, tmp_sign;
855   p_sign = x_sign ^ y_sign; // sign of the product
862         if (p_sign == z_sign) {
873         res.w[1] = p_sign | MASK_INF;
878         if (p_sign == z_sign) {
889         res.w[1] = p_sign | MASK_INF;
908       if ((p_sign != z_sign)
928       res.w[1] = p_sign | MASK_INF;
961     if (p_sign == z_sign) {
1362         res.w[1] = p_sign | 0x7800000000000000ull; // +/-inf
1366         res.w[1] = p_sign | res.w[1];
1558       res.w[1] = p_sign | ((UINT64) (e4 + 6176) << 49) | res.w[1];
1575     res.w[1] = p_sign | ((UINT64) (e4 + 6176) << 49) | res.w[1];
1629           res.w[1] = p_sign | (p_exp & MASK_EXP) | res.w[1];
1633           res.w[1] = p_sign | (p_exp & MASK_EXP) | res.w[1];
1663           if (p_sign == z_sign) {
1743       if ((p_sign != z_sign) && (delta == (q3 + scale + 1))) {
1839         if (p_sign == z_sign) {
1842         } else { // if (p_sign != z_sign)
1882           z_sign != p_sign && ((!z_sign && rnd_mode != ROUNDING_UP) ||
1985       if (p_sign == z_sign) {
2040       } else { // if (p_sign != z_sign)
2239             // p_sign != z_sign so swap gt_half_ulp and lt_half_ulp
2273       } // end if (p_sign != z_sign)
2287         !(delta <= 1 && p_sign != z_sign)) { // Case (2), (3), (4), (5) or (6)
2447       if (z_sign == p_sign) {
2583       } else { // if (z_sign != p_sign)
2585         // used to swap rounding indicators if p_sign != z_sign
2875         // (17) if we swap (C3, C4), (q3, q4), (e3, e4), (z_sign, p_sign)
2891         z_sign = p_sign;
2892         p_sign = tmp_sign;
2899       add_and_round (q3, q4, e4, delta, p34, z_sign, p_sign, C3, C4,
2954         if (p_sign == z_sign) {
2956         } else { // if (p_sign != z_sign)
3003         if (z_sign != p_sign) {
3018         } else { // if (z_sign == p_sign)
3023         if (z_sign == p_sign) {
3030         } else { // if (z_sign != p_sign)
3039         res.w[1] = p_sign | 0x7800000000000000ull;
3044         res.w[1] = p_sign | (p_exp & MASK_EXP) | res.w[1];
3079       // swap (C3, C4), (q3, q4), (e3, 34), (z_sign, p_sign), (z_exp, p_exp)
3095       z_sign = p_sign;
3096       p_sign = tmp_sign;
3138       if (p_sign == z_sign) { // R256 = C4 + R256
3140       } else { // if (p_sign != z_sign) { // R256 = C4 - R256
3207         // is (-1)^p_sign * R256 * 10^e4
3321         // (-1)^p_sign * res * 10^0
3322         P128.w[1] = p_sign | 0x3040000000000000ull | res.w[1];
3337       // is (-1)^p_sign * res * 10^e4
3338       res.w[1] = p_sign | ((UINT64) (e4 + 6176) << 49) | res.w[1]; // RN
3345       // res = (-1)^p_sign * significand * 10^e4,
3352         res.w[1] = p_sign | 0x7800000000000000ull;
3370         //   (-1)^p_sign * significand_res * 10^e4
3385           res.w[1] = p_sign | 0x0000000000000000ull;
3430           res.w[1] = p_sign | res.w[1];
3460             p_sign | ((UINT64) (e4 + 6176) << 49) | (res.
3545       add_and_round (q3, q4, e4, delta, p34, z_sign, p_sign, C3, C4,