config/libbid/bid64_div.c

523 UINT128 CX, CY, T128, CQ, CQ2, CR, CA, TP128, Qh, Ql, Tmp;
656   CQ.w[0] = CQ.w[1] = 0;
660   // get CQ = CX/CY
661   __div_128_by_128 (&CQ, &CR, CX, CY);
663   // get number of decimal digits in CQ
666   fx.d = (float) CQ.w[1] * f64.d + (float) CQ.w[0];
667   // binary expon. of CQ
673   if (__unsigned_compare_ge_128 (CQ, TP128))
679 		       CQ.w[0], rnd_mode, pfpsf);
690     CQ.w[0] *= T128.w[0];
695     __mul_128x128_to_256 (P256, CQ, T128);
697     CQ.w[0] = (P256.w[2] >> amount) | (P256.w[3] << (64 - amount));
698     CQ.w[1] = 0;
700     __mul_64x64_to_128 (CQ2, CQ.w[0], (power10_table_128[ed2].w[0]));
717   __div_256_by_128 (&CQ, &CA4, CY);
749       __mul_128x128_full (Qh, Ql, CQ, reciprocals10_128[nzeros]);
753       __shr_128_long (CQ, Qh, amount);
758       Q_low = CQ.w[0];
809 	  __mul_128x128_full (Qh, Ql, CQ, reciprocals10_128[nzeros]);
813 	  __shr_128 (CQ, Qh, amount);
821       fast_get_BID64_check_OF (sign_x ^ sign_y, diff_expon, CQ.w[0],
841   carry64 = (1 + (((SINT64) CA4r.w[1]) >> 63)) & ((CQ.w[0]) | D);
843   CQ.w[0] += carry64;
856   CQ.w[0] += carry64;
857   if (CQ.w[0] < carry64)
858     CQ.w[1]++;
872     carry64 = (1 + (((SINT64) CA4r.w[1]) >> 63)) & ((CQ.w[0]) | D);
873     CQ.w[0] += carry64;
874     if (CQ.w[0] < carry64)
875       CQ.w[1]++;
886     CQ.w[0] += carry64;
887     if (CQ.w[0] < carry64)
888       CQ.w[1]++;
894     CQ.w[0]++;
895     if (!CQ.w[0])
896       CQ.w[1]++;
903       fast_get_BID64_check_OF (sign_x ^ sign_y, diff_expon, CQ.w[0], rnd_mode,
920       get_BID64_UF (sign_x ^ sign_y, diff_expon, CQ.w[0], CA4.w[1] | CA4.w[0], rmode, pfpsf);
937 UINT128 CX, CY, T128, CQ, CQ2, CR, CA, TP128, Qh, Ql, Tmp;
1079   CQ.w[0] = CQ.w[1] = 0;
1083   // get CQ = CX/CY
1084   __div_128_by_128 (&CQ, &CR, CX, CY);
1086   // get number of decimal digits in CQ
1089   fx.d = (float) CQ.w[1] * f64.d + (float) CQ.w[0];
1090   // binary expon. of CQ
1096   if (__unsigned_compare_ge_128 (CQ, TP128))
1102 		       CQ.w[0], rnd_mode, pfpsf);
1113     CQ.w[0] *= T128.w[0];
1118     __mul_128x128_to_256 (P256, CQ, T128);
1120     CQ.w[0] = (P256.w[2] >> amount) | (P256.w[3] << (64 - amount));
1121     CQ.w[1] = 0;
1123     __mul_64x64_to_128 (CQ2, CQ.w[0], (power10_table_128[ed2].w[0]));
1144   __div_256_by_128 (&CQ, &CA4, CY);
1175       __mul_128x128_full (Qh, Ql, CQ, reciprocals10_128[nzeros]);
1176       //__mul_128x128_to_256(P256, CQ, reciprocals10_128[nzeros]);Qh.w[1]=P256.w[3];Qh.w[0]=P256.w[2];
1180       __shr_128_long (CQ, Qh, amount);
1186       Q_low = CQ.w[0];
1237 	  __mul_128x128_full (Qh, Ql, CQ, reciprocals10_128[nzeros]);
1241 	  __shr_128 (CQ, Qh, amount);
1250       fast_get_BID64_check_OF (sign_x ^ sign_y, diff_expon, CQ.w[0],
1270   carry64 = (1 + (((SINT64) CA4r.w[1]) >> 63)) & ((CQ.w[0]) | D);
1272   CQ.w[0] += carry64;
1273   //if(CQ.w[0]<carry64)
1274   //CQ.w[1] ++;
1287   CQ.w[0] += carry64;
1288   if (CQ.w[0] < carry64)
1289     CQ.w[1]++;
1303     carry64 = (1 + (((SINT64) CA4r.w[1]) >> 63)) & ((CQ.w[0]) | D);
1304     CQ.w[0] += carry64;
1305     if (CQ.w[0] < carry64)
1306       CQ.w[1]++;
1317     CQ.w[0] += carry64;
1318     if (CQ.w[0] < carry64)
1319       CQ.w[1]++;
1325     CQ.w[0]++;
1326     if (!CQ.w[0])
1327       CQ.w[1]++;
1335       fast_get_BID64_check_OF (sign_x ^ sign_y, diff_expon, CQ.w[0], rnd_mode,
1352       get_BID64_UF (sign_x ^ sign_y, diff_expon, CQ.w[0], CA4.w[1] | CA4.w[0], rmode, pfpsf);
1374 UINT128 CX, CY, T128, CQ, CQ2, CR, CA, TP128, Qh, Ql, Tmp;
1512   CQ.w[0] = CQ.w[1] = 0;
1516   // get CQ = CX/CY
1517   __div_128_by_128 (&CQ, &CR, CX, CY);
1519   // get number of decimal digits in CQ
1522   fx.d = (float) CQ.w[1] * f64.d + (float) CQ.w[0];
1523   // binary expon. of CQ
1529   if (__unsigned_compare_ge_128 (CQ, TP128))
1535 		       CQ.w[0], rnd_mode, pfpsf);
1546     CQ.w[0] *= T128.w[0];
1551     __mul_128x128_to_256 (P256, CQ, T128);
1553     CQ.w[0] = (P256.w[2] >> amount) | (P256.w[3] << (64 - amount));
1554     CQ.w[1] = 0;
1556     __mul_64x64_to_128 (CQ2, CQ.w[0], (power10_table_128[ed2].w[0]));
1576   __div_256_by_128 (&CQ, &CA4, CY);
1609       __mul_128x128_full (Qh, Ql, CQ, reciprocals10_128[nzeros]);
1610       //__mul_128x128_to_256(P256, CQ, reciprocals10_128[nzeros]);Qh.w[1]=P256.w[3];Qh.w[0]=P256.w[2];
1614       __shr_128_long (CQ, Qh, amount);
1620       Q_low = CQ.w[0];
1671 	  __mul_128x128_full (Qh, Ql, CQ, reciprocals10_128[nzeros]);
1675 	  __shr_128 (CQ, Qh, amount);
1684       fast_get_BID64_check_OF (sign_x ^ sign_y, diff_expon, CQ.w[0],
1705   carry64 = (1 + (((SINT64) CA4r.w[1]) >> 63)) & ((CQ.w[0]) | D);
1707   CQ.w[0] += carry64;
1708   //if(CQ.w[0]<carry64)
1709   //CQ.w[1] ++;
1722   CQ.w[0] += carry64;
1723   if (CQ.w[0] < carry64)
1724     CQ.w[1]++;
1738     carry64 = (1 + (((SINT64) CA4r.w[1]) >> 63)) & ((CQ.w[0]) | D);
1739     CQ.w[0] += carry64;
1740     if (CQ.w[0] < carry64)
1741       CQ.w[1]++;
1752     CQ.w[0] += carry64;
1753     if (CQ.w[0] < carry64)
1754       CQ.w[1]++;
1760     CQ.w[0]++;
1761     if (!CQ.w[0])
1762       CQ.w[1]++;
1770       fast_get_BID64_check_OF (sign_x ^ sign_y, diff_expon, CQ.w[0], rnd_mode,
1787       get_BID64_UF (sign_x ^ sign_y, diff_expon, CQ.w[0], CA4.w[1] | CA4.w[0], rmode, pfpsf);