config/libbid/bid128_div.c

39      UINT128 CX, CY, T128, CQ, CR, CA, TP128, Qh, res;
52 valid_y = unpack_BID128_value (&sign_y, &exponent_y, &CY, y);
93   if ((!CY.w[0]) && !(CY.w[1] & 0x0001ffffffffffffull)) {
123     res.w[1] = CY.w[1] & QUIET_MASK64;
124     res.w[0] = CY.w[0];
148 if (__unsigned_compare_gt_128 (CY, CX)) {
149   // CX < CY
154   // fx ~ CX,   fy ~ CY
156   fy.d = (float) CY.w[1] * f64.d + (float) CY.w[0];
169   if (__unsigned_compare_gt_128 (CY, CA))
180   // get CQ = CX/CY
181   __div_128_by_128 (&CQ, &CR, CX, CY);
213 __div_256_by_128 (&CQ, &CA4, CY);
231   // check whether CX, CY are short
232   if (!CX.w[1] && !CY.w[1] && (CX.w[0] <= 1024) && (CY.w[0] <= 1024)) {
233     i = (int) CY.w[0] - 1;
389   // 2*CA4 - CY
392   __sub_borrow_out (CA4r.w[0], carry64, CA4r.w[0], CY.w[0]);
393   CA4r.w[1] = CA4r.w[1] - CY.w[1] - carry64;
404   // 2*CA4 - CY
407   __sub_borrow_out (CA4r.w[0], carry64, CA4r.w[0], CY.w[0]);
408   CA4r.w[1] = CA4r.w[1] - CY.w[1] - carry64;
423     // 2*CA4 - CY
426     __sub_borrow_out (CA4r.w[0], carry64, CA4r.w[0], CY.w[0]);
427     CA4r.w[1] = CA4r.w[1] - CY.w[1] - carry64;
436     // 2*CA4 - CY
439     __sub_borrow_out (CA4r.w[0], carry64, CA4r.w[0], CY.w[0]);
440     CA4r.w[1] = CA4r.w[1] - CY.w[1] - carry64;
490      UINT128 CX, CY, T128, CQ, CR, CA, TP128, Qh, res;
503 valid_y = unpack_BID64 (&sign_y, &exponent_y, &CY.w[0], y);
547     if(!CY.w[0]) {
573 CY.w[1] = 0;
583   res.w[0] = (CY.w[0] & 0x0003ffffffffffffull);
585   res.w[1] |= ((CY.w[0]) & 0xfc00000000000000ull);
609 if (__unsigned_compare_gt_128 (CY, CX)) {
610   // CX < CY
615   // fx ~ CX,   fy ~ CY
617   fy.d = (float) CY.w[1] * f64.d + (float) CY.w[0];
630   if (__unsigned_compare_gt_128 (CY, CA))
641   // get CQ = CX/CY
642   __div_128_by_128 (&CQ, &CR, CX, CY);
674 __div_256_by_128 (&CQ, &CA4, CY);
693     // check whether CX, CY are short
694     if (!CX.w[1] && !CY.w[1] && (CX.w[0] <= 1024) && (CY.w[0] <= 1024)) {
695       i = (int) CY.w[0] - 1;
853   // 2*CA4 - CY
856   __sub_borrow_out (CA4r.w[0], carry64, CA4r.w[0], CY.w[0]);
857   CA4r.w[1] = CA4r.w[1] - CY.w[1] - carry64;
868   // 2*CA4 - CY
871   __sub_borrow_out (CA4r.w[0], carry64, CA4r.w[0], CY.w[0]);
872   CA4r.w[1] = CA4r.w[1] - CY.w[1] - carry64;
887     // 2*CA4 - CY
890     __sub_borrow_out (CA4r.w[0], carry64, CA4r.w[0], CY.w[0]);
891     CA4r.w[1] = CA4r.w[1] - CY.w[1] - carry64;
900     // 2*CA4 - CY
903     __sub_borrow_out (CA4r.w[0], carry64, CA4r.w[0], CY.w[0]);
904     CA4r.w[1] = CA4r.w[1] - CY.w[1] - carry64;
949      UINT128 CX, CY, T128, CQ, CR, CA, TP128, Qh, res;
962 valid_y = unpack_BID128_value (&sign_y, &exponent_y, &CY, y);
1006   if ((!CY.w[0]) && !(CY.w[1] & 0x0001ffffffffffffull)) {
1038     res.w[1] = CY.w[1] & QUIET_MASK64;
1039     res.w[0] = CY.w[0];
1063 if (__unsigned_compare_gt_128 (CY, CX)) {
1064   // CX < CY
1069   // fx ~ CX,   fy ~ CY
1071   fy.d = (float) CY.w[1] * f64.d + (float) CY.w[0];
1084   if (__unsigned_compare_gt_128 (CY, CA))
1095   // get CQ = CX/CY
1096   __div_128_by_128 (&CQ, &CR, CX, CY);
1128 __div_256_by_128 (&CQ, &CA4, CY);
1147     // check whether CX, CY are short
1148     if (!CX.w[1] && !CY.w[1] && (CX.w[0] <= 1024) && (CY.w[0] <= 1024)) {
1149       i = (int) CY.w[0] - 1;
1310   // 2*CA4 - CY
1313   __sub_borrow_out (CA4r.w[0], carry64, CA4r.w[0], CY.w[0]);
1314   CA4r.w[1] = CA4r.w[1] - CY.w[1] - carry64;
1325   // 2*CA4 - CY
1328   __sub_borrow_out (CA4r.w[0], carry64, CA4r.w[0], CY.w[0]);
1329   CA4r.w[1] = CA4r.w[1] - CY.w[1] - carry64;
1344     // 2*CA4 - CY
1347     __sub_borrow_out (CA4r.w[0], carry64, CA4r.w[0], CY.w[0]);
1348     CA4r.w[1] = CA4r.w[1] - CY.w[1] - carry64;
1357     // 2*CA4 - CY
1360     __sub_borrow_out (CA4r.w[0], carry64, CA4r.w[0], CY.w[0]);
1361     CA4r.w[1] = CA4r.w[1] - CY.w[1] - carry64;
1406      UINT128 CX, CY, T128, CQ, CR, CA, TP128, Qh, res;
1419 valid_y = unpack_BID64 (&sign_y, &exponent_y, &CY.w[0], y);
1459 	if (!CY.w[0]) {
1480 CY.w[1] = 0;
1490   res.w[0] = (CY.w[0] & 0x0003ffffffffffffull);
1492   res.w[1] |= ((CY.w[0]) & 0xfc00000000000000ull);
1515 if (__unsigned_compare_gt_128 (CY, CX)) {
1516   // CX < CY
1521   // fx ~ CX,   fy ~ CY
1523   fy.d = (float) CY.w[1] * f64.d + (float) CY.w[0];
1536   if (__unsigned_compare_gt_128 (CY, CA))
1547   // get CQ = CX/CY
1548   __div_128_by_128 (&CQ, &CR, CX, CY);
1580 __div_256_by_128 (&CQ, &CA4, CY);
1599     // check whether CX, CY are short
1600     if (!CX.w[1] && !CY.w[1] && (CX.w[0] <= 1024) && (CY.w[0] <= 1024)) {
1601       i = (int) CY.w[0] - 1;
1759   // 2*CA4 - CY
1762   __sub_borrow_out (CA4r.w[0], carry64, CA4r.w[0], CY.w[0]);
1763   CA4r.w[1] = CA4r.w[1] - CY.w[1] - carry64;
1774   // 2*CA4 - CY
1777   __sub_borrow_out (CA4r.w[0], carry64, CA4r.w[0], CY.w[0]);
1778   CA4r.w[1] = CA4r.w[1] - CY.w[1] - carry64;
1793     // 2*CA4 - CY
1796     __sub_borrow_out (CA4r.w[0], carry64, CA4r.w[0], CY.w[0]);
1797     CA4r.w[1] = CA4r.w[1] - CY.w[1] - carry64;
1806     // 2*CA4 - CY
1809     __sub_borrow_out (CA4r.w[0], carry64, CA4r.w[0], CY.w[0]);
1810     CA4r.w[1] = CA4r.w[1] - CY.w[1] - carry64;