arm/insts/fplib.cc

556     uint16_t b_mnt = FP16_MANT(b);
561     if (fp16_is_signalling_NaN(b_exp, b_mnt))
567     if (fp16_is_NaN(b_exp, b_mnt))
579     uint32_t b_mnt = FP32_MANT(b);
584     if (fp32_is_signalling_NaN(b_exp, b_mnt))
590     if (fp32_is_NaN(b_exp, b_mnt))
602     uint64_t b_mnt = FP64_MANT(b);
607     if (fp64_is_signalling_NaN(b_exp, b_mnt))
613     if (fp64_is_NaN(b_exp, b_mnt))
625     uint16_t b_mnt = FP16_MANT(b);
632     if (fp16_is_signalling_NaN(b_exp, b_mnt))
640     if (fp16_is_NaN(b_exp, b_mnt))
654     uint32_t b_mnt = FP32_MANT(b);
661     if (fp32_is_signalling_NaN(b_exp, b_mnt))
669     if (fp32_is_NaN(b_exp, b_mnt))
683     uint64_t b_mnt = FP64_MANT(b);
690     if (fp64_is_signalling_NaN(b_exp, b_mnt))
698     if (fp64_is_NaN(b_exp, b_mnt))
948     uint16_t a_mnt, b_mnt;
951     fp16_unpack(&b_sgn, &b_exp, &b_mnt, b, mode, flags);
954         fp16_is_NaN(b_exp, b_mnt)) {
956             fp16_is_signalling_NaN(b_exp, b_mnt))
960     return a == b || (!a_mnt && !b_mnt);
967     uint16_t a_mnt, b_mnt;
970     fp16_unpack(&b_sgn, &b_exp, &b_mnt, b, mode, flags);
973         fp16_is_NaN(b_exp, b_mnt)) {
977     if (!a_mnt && !b_mnt)
983     if (a_mnt != b_mnt)
984         return a_sgn ^ (a_mnt > b_mnt);
992     uint16_t a_mnt, b_mnt;
995     fp16_unpack(&b_sgn, &b_exp, &b_mnt, b, mode, flags);
998         fp16_is_NaN(b_exp, b_mnt)) {
1002     if (!a_mnt && !b_mnt)
1008     if (a_mnt != b_mnt)
1009         return a_sgn ^ (a_mnt > b_mnt);
1017     uint16_t a_mnt, b_mnt;
1020     fp16_unpack(&b_sgn, &b_exp, &b_mnt, b, mode, flags);
1023         fp16_is_NaN(b_exp, b_mnt)) {
1025             fp16_is_signalling_NaN(b_exp, b_mnt))
1036     uint32_t a_mnt, b_mnt;
1039     fp32_unpack(&b_sgn, &b_exp, &b_mnt, b, mode, flags);
1042         fp32_is_NaN(b_exp, b_mnt)) {
1044             fp32_is_signalling_NaN(b_exp, b_mnt))
1048     return a == b || (!a_mnt && !b_mnt);
1055     uint32_t a_mnt, b_mnt;
1058     fp32_unpack(&b_sgn, &b_exp, &b_mnt, b, mode, flags);
1061         fp32_is_NaN(b_exp, b_mnt)) {
1065     if (!a_mnt && !b_mnt)
1071     if (a_mnt != b_mnt)
1072         return a_sgn ^ (a_mnt > b_mnt);
1080     uint32_t a_mnt, b_mnt;
1083     fp32_unpack(&b_sgn, &b_exp, &b_mnt, b, mode, flags);
1086         fp32_is_NaN(b_exp, b_mnt)) {
1090     if (!a_mnt && !b_mnt)
1096     if (a_mnt != b_mnt)
1097         return a_sgn ^ (a_mnt > b_mnt);
1105     uint32_t a_mnt, b_mnt;
1108     fp32_unpack(&b_sgn, &b_exp, &b_mnt, b, mode, flags);
1111         fp32_is_NaN(b_exp, b_mnt)) {
1113             fp32_is_signalling_NaN(b_exp, b_mnt))
1124     uint64_t a_mnt, b_mnt;
1127     fp64_unpack(&b_sgn, &b_exp, &b_mnt, b, mode, flags);
1130         fp64_is_NaN(b_exp, b_mnt)) {
1132             fp64_is_signalling_NaN(b_exp, b_mnt))
1136     return a == b || (!a_mnt && !b_mnt);
1143     uint64_t a_mnt, b_mnt;
1146     fp64_unpack(&b_sgn, &b_exp, &b_mnt, b, mode, flags);
1149         fp64_is_NaN(b_exp, b_mnt)) {
1153     if (!a_mnt && !b_mnt)
1159     if (a_mnt != b_mnt)
1160         return a_sgn ^ (a_mnt > b_mnt);
1168     uint64_t a_mnt, b_mnt;
1171     fp64_unpack(&b_sgn, &b_exp, &b_mnt, b, mode, flags);
1174         fp64_is_NaN(b_exp, b_mnt)) {
1178     if (!a_mnt && !b_mnt)
1184     if (a_mnt != b_mnt)
1185         return a_sgn ^ (a_mnt > b_mnt);
1193     uint64_t a_mnt, b_mnt;
1196     fp64_unpack(&b_sgn, &b_exp, &b_mnt, b, mode, flags);
1199         fp64_is_NaN(b_exp, b_mnt)) {
1201             fp64_is_signalling_NaN(b_exp, b_mnt))
1212     uint16_t a_mnt, b_mnt, x, x_mnt;
1215     fp16_unpack(&b_sgn, &b_exp, &b_mnt, b, mode, flags);
1231     } else if (!a_mnt && !b_mnt && a_sgn == b_sgn) {
1236     b_mnt <<= 3;
1238         b_mnt = (lsr16(b_mnt, a_exp - b_exp) |
1239                  !!(b_mnt & (lsl16(1, a_exp - b_exp) - 1)));
1249         x_mnt = a_mnt + b_mnt;
1250     } else if (a_mnt >= b_mnt) {
1251         x_mnt = a_mnt - b_mnt;
1254         x_mnt = b_mnt - a_mnt;
1272     uint32_t a_mnt, b_mnt, x, x_mnt;
1275     fp32_unpack(&b_sgn, &b_exp, &b_mnt, b, mode, flags);
1291     } else if (!a_mnt && !b_mnt && a_sgn == b_sgn) {
1296     b_mnt <<= 3;
1298         b_mnt = (lsr32(b_mnt, a_exp - b_exp) |
1299                  !!(b_mnt & (lsl32(1, a_exp - b_exp) - 1)));
1309         x_mnt = a_mnt + b_mnt;
1310     } else if (a_mnt >= b_mnt) {
1311         x_mnt = a_mnt - b_mnt;
1314         x_mnt = b_mnt - a_mnt;
1332     uint64_t a_mnt, b_mnt, x, x_mnt;
1335     fp64_unpack(&b_sgn, &b_exp, &b_mnt, b, mode, flags);
1351     } else if (!a_mnt && !b_mnt && a_sgn == b_sgn) {
1356     b_mnt <<= 3;
1358         b_mnt = (lsr64(b_mnt, a_exp - b_exp) |
1359                  !!(b_mnt & (lsl64(1, a_exp - b_exp) - 1)));
1369         x_mnt = a_mnt + b_mnt;
1370     } else if (a_mnt >= b_mnt) {
1371         x_mnt = a_mnt - b_mnt;
1374         x_mnt = b_mnt - a_mnt;
1392     uint16_t a_mnt, b_mnt, x;
1396     fp16_unpack(&b_sgn, &b_exp, &b_mnt, b, mode, flags);
1403     if ((a_exp == FP16_EXP_INF && !b_mnt) ||
1409     } else if (!a_mnt || !b_mnt) {
1416     x_mnt = (uint32_t)a_mnt * b_mnt;
1429     uint32_t a_mnt, b_mnt, x;
1433     fp32_unpack(&b_sgn, &b_exp, &b_mnt, b, mode, flags);
1440     if ((a_exp == FP32_EXP_INF && !b_mnt) ||
1446     } else if (!a_mnt || !b_mnt) {
1453     x_mnt = (uint64_t)a_mnt * b_mnt;
1466     uint64_t a_mnt, b_mnt, x;
1470     fp64_unpack(&b_sgn, &b_exp, &b_mnt, b, mode, flags);
1477     if ((a_exp == FP64_EXP_INF && !b_mnt) ||
1483     } else if (!a_mnt || !b_mnt) {
1490     mul62x62(&x0_mnt, &x1_mnt, a_mnt, b_mnt);
1504     uint16_t a_mnt, b_mnt, c_mnt, x;
1508     fp16_unpack(&b_sgn, &b_exp, &b_mnt, b, mode, flags);
1515         ((!b_mnt && fp16_is_infinity(c_exp, c_mnt)) ||
1516          (!c_mnt && fp16_is_infinity(b_exp, b_mnt)))) {
1527         (c_exp == FP16_EXP_INF && !b_mnt) ||
1538     if (!a_mnt && (!b_mnt || !c_mnt) && a_sgn == (b_sgn ^ c_sgn))
1548     y_mnt = (uint32_t)b_mnt * c_mnt << 3;
1589     uint32_t a_mnt, b_mnt, c_mnt, x;
1593     fp32_unpack(&b_sgn, &b_exp, &b_mnt, b, mode, flags);
1600         ((!b_mnt && fp32_is_infinity(c_exp, c_mnt)) ||
1601          (!c_mnt && fp32_is_infinity(b_exp, b_mnt)))) {
1612         (c_exp == FP32_EXP_INF && !b_mnt) ||
1623     if (!a_mnt && (!b_mnt || !c_mnt) && a_sgn == (b_sgn ^ c_sgn))
1633     y_mnt = (uint64_t)b_mnt * c_mnt << 3;
1674     uint64_t a_mnt, b_mnt, c_mnt, x;
1678     fp64_unpack(&b_sgn, &b_exp, &b_mnt, b, mode, flags);
1685         ((!b_mnt && fp64_is_infinity(c_exp, c_mnt)) ||
1686          (!c_mnt && fp64_is_infinity(b_exp, b_mnt)))) {
1697         (c_exp == FP64_EXP_INF && !b_mnt) ||
1708     if (!a_mnt && (!b_mnt || !c_mnt) && a_sgn == (b_sgn ^ c_sgn))
1719     mul62x62(&y0_mnt, &y1_mnt, b_mnt, c_mnt << 3);
1765     uint16_t a_mnt, b_mnt, x;
1769     fp16_unpack(&b_sgn, &b_exp, &b_mnt, b, mode, flags);
1776         (!a_mnt && !b_mnt)) {
1780     if (a_exp == FP16_EXP_INF || !b_mnt) {
1792     x_mnt = ((uint32_t)a_mnt << (FP16_MANT_BITS - FP16_EXP_BITS + 3)) / b_mnt;
1793     x_mnt |= (x_mnt * b_mnt !=
1807     uint32_t a_mnt, b_mnt, x;
1811     fp32_unpack(&b_sgn, &b_exp, &b_mnt, b, mode, flags);
1818         (!a_mnt && !b_mnt)) {
1822     if (a_exp == FP32_EXP_INF || !b_mnt) {
1834     x_mnt = ((uint64_t)a_mnt << (FP32_MANT_BITS - FP32_EXP_BITS + 3)) / b_mnt;
1835     x_mnt |= (x_mnt * b_mnt !=
1849     uint64_t a_mnt, b_mnt, x, x_mnt, x0_mnt, x1_mnt;
1852     fp64_unpack(&b_sgn, &b_exp, &b_mnt, b, mode, flags);
1859         (!a_mnt && !b_mnt)) {
1863     if (a_exp == FP64_EXP_INF || !b_mnt) {
1873     b_mnt = fp64_normalise(b_mnt, &b_exp);
1874     x_mnt = ~(uint64_t)0 / (b_mnt >> 31);
1875     mul64x32(&x0_mnt, &x1_mnt, b_mnt, x_mnt);
1889     mul62x62(&x0_mnt, &x1_mnt, b_mnt >> 2, x_mnt + 1);