arm/insts/fplib.cc

553     int a_exp = FP16_EXP(a);
559     if (fp16_is_signalling_NaN(a_exp, a_mnt))
565     if (fp16_is_NaN(a_exp, a_mnt))
576     int a_exp = FP32_EXP(a);
582     if (fp32_is_signalling_NaN(a_exp, a_mnt))
588     if (fp32_is_NaN(a_exp, a_mnt))
599     int a_exp = FP64_EXP(a);
605     if (fp64_is_signalling_NaN(a_exp, a_mnt))
611     if (fp64_is_NaN(a_exp, a_mnt))
622     int a_exp = FP16_EXP(a);
630     if (fp16_is_signalling_NaN(a_exp, a_mnt))
638     if (fp16_is_NaN(a_exp, a_mnt))
651     int a_exp = FP32_EXP(a);
659     if (fp32_is_signalling_NaN(a_exp, a_mnt))
667     if (fp32_is_NaN(a_exp, a_mnt))
680     int a_exp = FP64_EXP(a);
688     if (fp64_is_signalling_NaN(a_exp, a_mnt))
696     if (fp64_is_NaN(a_exp, a_mnt))
947     int a_sgn, a_exp, b_sgn, b_exp;
950     fp16_unpack(&a_sgn, &a_exp, &a_mnt, a, mode, flags);
953     if (fp16_is_NaN(a_exp, a_mnt) ||
955         if (fp16_is_signalling_NaN(a_exp, a_mnt) ||
966     int a_sgn, a_exp, b_sgn, b_exp;
969     fp16_unpack(&a_sgn, &a_exp, &a_mnt, a, mode, flags);
972     if (fp16_is_NaN(a_exp, a_mnt) ||
981     if (a_exp != b_exp)
982         return a_sgn ^ (a_exp > b_exp);
991     int a_sgn, a_exp, b_sgn, b_exp;
994     fp16_unpack(&a_sgn, &a_exp, &a_mnt, a, mode, flags);
997     if (fp16_is_NaN(a_exp, a_mnt) ||
1006     if (a_exp != b_exp)
1007         return a_sgn ^ (a_exp > b_exp);
1016     int a_sgn, a_exp, b_sgn, b_exp;
1019     fp16_unpack(&a_sgn, &a_exp, &a_mnt, a, mode, flags);
1022     if (fp16_is_NaN(a_exp, a_mnt) ||
1024         if (fp16_is_signalling_NaN(a_exp, a_mnt) ||
1035     int a_sgn, a_exp, b_sgn, b_exp;
1038     fp32_unpack(&a_sgn, &a_exp, &a_mnt, a, mode, flags);
1041     if (fp32_is_NaN(a_exp, a_mnt) ||
1043         if (fp32_is_signalling_NaN(a_exp, a_mnt) ||
1054     int a_sgn, a_exp, b_sgn, b_exp;
1057     fp32_unpack(&a_sgn, &a_exp, &a_mnt, a, mode, flags);
1060     if (fp32_is_NaN(a_exp, a_mnt) ||
1069     if (a_exp != b_exp)
1070         return a_sgn ^ (a_exp > b_exp);
1079     int a_sgn, a_exp, b_sgn, b_exp;
1082     fp32_unpack(&a_sgn, &a_exp, &a_mnt, a, mode, flags);
1085     if (fp32_is_NaN(a_exp, a_mnt) ||
1094     if (a_exp != b_exp)
1095         return a_sgn ^ (a_exp > b_exp);
1104     int a_sgn, a_exp, b_sgn, b_exp;
1107     fp32_unpack(&a_sgn, &a_exp, &a_mnt, a, mode, flags);
1110     if (fp32_is_NaN(a_exp, a_mnt) ||
1112         if (fp32_is_signalling_NaN(a_exp, a_mnt) ||
1123     int a_sgn, a_exp, b_sgn, b_exp;
1126     fp64_unpack(&a_sgn, &a_exp, &a_mnt, a, mode, flags);
1129     if (fp64_is_NaN(a_exp, a_mnt) ||
1131         if (fp64_is_signalling_NaN(a_exp, a_mnt) ||
1142     int a_sgn, a_exp, b_sgn, b_exp;
1145     fp64_unpack(&a_sgn, &a_exp, &a_mnt, a, mode, flags);
1148     if (fp64_is_NaN(a_exp, a_mnt) ||
1157     if (a_exp != b_exp)
1158         return a_sgn ^ (a_exp > b_exp);
1167     int a_sgn, a_exp, b_sgn, b_exp;
1170     fp64_unpack(&a_sgn, &a_exp, &a_mnt, a, mode, flags);
1173     if (fp64_is_NaN(a_exp, a_mnt) ||
1182     if (a_exp != b_exp)
1183         return a_sgn ^ (a_exp > b_exp);
1192     int a_sgn, a_exp, b_sgn, b_exp;
1195     fp64_unpack(&a_sgn, &a_exp, &a_mnt, a, mode, flags);
1198     if (fp64_is_NaN(a_exp, a_mnt) ||
1200         if (fp64_is_signalling_NaN(a_exp, a_mnt) ||
1211     int a_sgn, a_exp, b_sgn, b_exp, x_sgn, x_exp;
1214     fp16_unpack(&a_sgn, &a_exp, &a_mnt, a, mode, flags);
1224     if (a_exp == FP16_EXP_INF && b_exp == FP16_EXP_INF && a_sgn != b_sgn) {
1227     } else if (a_exp == FP16_EXP_INF) {
1237     if (a_exp >= b_exp) {
1238         b_mnt = (lsr16(b_mnt, a_exp - b_exp) |
1239                  !!(b_mnt & (lsl16(1, a_exp - b_exp) - 1)));
1240         b_exp = a_exp;
1242         a_mnt = (lsr16(a_mnt, b_exp - a_exp) |
1243                  !!(a_mnt & (lsl16(1, b_exp - a_exp) - 1)));
1244         a_exp = b_exp;
1247     x_exp = a_exp;
1271     int a_sgn, a_exp, b_sgn, b_exp, x_sgn, x_exp;
1274     fp32_unpack(&a_sgn, &a_exp, &a_mnt, a, mode, flags);
1284     if (a_exp == FP32_EXP_INF && b_exp == FP32_EXP_INF && a_sgn != b_sgn) {
1287     } else if (a_exp == FP32_EXP_INF) {
1297     if (a_exp >= b_exp) {
1298         b_mnt = (lsr32(b_mnt, a_exp - b_exp) |
1299                  !!(b_mnt & (lsl32(1, a_exp - b_exp) - 1)));
1300         b_exp = a_exp;
1302         a_mnt = (lsr32(a_mnt, b_exp - a_exp) |
1303                  !!(a_mnt & (lsl32(1, b_exp - a_exp) - 1)));
1304         a_exp = b_exp;
1307     x_exp = a_exp;
1331     int a_sgn, a_exp, b_sgn, b_exp, x_sgn, x_exp;
1334     fp64_unpack(&a_sgn, &a_exp, &a_mnt, a, mode, flags);
1344     if (a_exp == FP64_EXP_INF && b_exp == FP64_EXP_INF && a_sgn != b_sgn) {
1347     } else if (a_exp == FP64_EXP_INF) {
1357     if (a_exp >= b_exp) {
1358         b_mnt = (lsr64(b_mnt, a_exp - b_exp) |
1359                  !!(b_mnt & (lsl64(1, a_exp - b_exp) - 1)));
1360         b_exp = a_exp;
1362         a_mnt = (lsr64(a_mnt, b_exp - a_exp) |
1363                  !!(a_mnt & (lsl64(1, b_exp - a_exp) - 1)));
1364         a_exp = b_exp;
1367     x_exp = a_exp;
1391     int a_sgn, a_exp, b_sgn, b_exp, x_sgn, x_exp;
1395     fp16_unpack(&a_sgn, &a_exp, &a_mnt, a, mode, flags);
1403     if ((a_exp == FP16_EXP_INF && !b_mnt) ||
1407     } else if (a_exp == FP16_EXP_INF || b_exp == FP16_EXP_INF) {
1415     x_exp = a_exp + b_exp - FP16_EXP_BIAS + 2 * FP16_EXP_BITS + 1;
1428     int a_sgn, a_exp, b_sgn, b_exp, x_sgn, x_exp;
1432     fp32_unpack(&a_sgn, &a_exp, &a_mnt, a, mode, flags);
1440     if ((a_exp == FP32_EXP_INF && !b_mnt) ||
1444     } else if (a_exp == FP32_EXP_INF || b_exp == FP32_EXP_INF) {
1452     x_exp = a_exp + b_exp - FP32_EXP_BIAS + 2 * FP32_EXP_BITS + 1;
1465     int a_sgn, a_exp, b_sgn, b_exp, x_sgn, x_exp;
1469     fp64_unpack(&a_sgn, &a_exp, &a_mnt, a, mode, flags);
1477     if ((a_exp == FP64_EXP_INF && !b_mnt) ||
1481     } else if (a_exp == FP64_EXP_INF || b_exp == FP64_EXP_INF) {
1489     x_exp = a_exp + b_exp - FP64_EXP_BIAS + 2 * FP64_EXP_BITS + 1;
1503     int a_sgn, a_exp, b_sgn, b_exp, c_sgn, c_exp, x_sgn, x_exp, y_sgn, y_exp;
1507     fp16_unpack(&a_sgn, &a_exp, &a_mnt, a, mode, flags);
1514     if (fp16_is_quiet_NaN(a_exp, a_mnt) &&
1528         (a_exp == FP16_EXP_INF &&
1534     if (a_exp == FP16_EXP_INF)
1542     x_exp = a_exp + 2 * FP16_EXP_BITS - 3;
1588     int a_sgn, a_exp, b_sgn, b_exp, c_sgn, c_exp, x_sgn, x_exp, y_sgn, y_exp;
1592     fp32_unpack(&a_sgn, &a_exp, &a_mnt, a, mode, flags);
1599     if (fp32_is_quiet_NaN(a_exp, a_mnt) &&
1613         (a_exp == FP32_EXP_INF &&
1619     if (a_exp == FP32_EXP_INF)
1627     x_exp = a_exp + 2 * FP32_EXP_BITS - 3;
1673     int a_sgn, a_exp, b_sgn, b_exp, c_sgn, c_exp, x_sgn, x_exp, y_sgn, y_exp;
1677     fp64_unpack(&a_sgn, &a_exp, &a_mnt, a, mode, flags);
1684     if (fp64_is_quiet_NaN(a_exp, a_mnt) &&
1698         (a_exp == FP64_EXP_INF &&
1704     if (a_exp == FP64_EXP_INF)
1712     x_exp = a_exp + FP64_EXP_BITS;
1764     int a_sgn, a_exp, b_sgn, b_exp, x_sgn, x_exp;
1768     fp16_unpack(&a_sgn, &a_exp, &a_mnt, a, mode, flags);
1775     if ((a_exp == FP16_EXP_INF && b_exp == FP16_EXP_INF) ||
1780     if (a_exp == FP16_EXP_INF || !b_mnt) {
1781         if (a_exp != FP16_EXP_INF)
1789     a_mnt = fp16_normalise(a_mnt, &a_exp);
1791     x_exp = a_exp - b_exp + (FP16_EXP_BIAS + FP16_BITS + 2 * FP16_EXP_BITS - 3);
1806     int a_sgn, a_exp, b_sgn, b_exp, x_sgn, x_exp;
1810     fp32_unpack(&a_sgn, &a_exp, &a_mnt, a, mode, flags);
1817     if ((a_exp == FP32_EXP_INF && b_exp == FP32_EXP_INF) ||
1822     if (a_exp == FP32_EXP_INF || !b_mnt) {
1823         if (a_exp != FP32_EXP_INF)
1831     a_mnt = fp32_normalise(a_mnt, &a_exp);
1833     x_exp = a_exp - b_exp + (FP32_EXP_BIAS + FP32_BITS + 2 * FP32_EXP_BITS - 3);
1848     int a_sgn, a_exp, b_sgn, b_exp, x_sgn, x_exp, c;
1851     fp64_unpack(&a_sgn, &a_exp, &a_mnt, a, mode, flags);
1858     if ((a_exp == FP64_EXP_INF && b_exp == FP64_EXP_INF) ||
1863     if (a_exp == FP64_EXP_INF || !b_mnt) {
1864         if (a_exp != FP64_EXP_INF)
1872     a_mnt = fp64_normalise(a_mnt, &a_exp);
1883     x_exp = a_exp - b_exp + FP64_EXP_BIAS + 8;
1926     int a_sgn, a_exp;
1929     fp16_unpack(&a_sgn, &a_exp, &a_mnt, a, mode, flags);
1932     if (fp16_is_NaN(a_exp, a_mnt)) {
1942     if (a_exp == FP16_EXP_INF) {
1948     a_exp += b;
1951     a_mnt = fp16_normalise(a_mnt, &a_exp);
1953     return fp16_round(a_sgn, a_exp + FP16_EXP_BITS - 3, a_mnt << 1,
1960     int a_sgn, a_exp;
1963     fp32_unpack(&a_sgn, &a_exp, &a_mnt, a, mode, flags);
1966     if (fp32_is_NaN(a_exp, a_mnt)) {
1976     if (a_exp == FP32_EXP_INF) {
1982     a_exp += b;
1985     a_mnt = fp32_normalise(a_mnt, &a_exp);
1987     return fp32_round(a_sgn, a_exp + FP32_EXP_BITS - 3, a_mnt << 1,
1994     int a_sgn, a_exp;
1997     fp64_unpack(&a_sgn, &a_exp, &a_mnt, a, mode, flags);
2000     if (fp64_is_NaN(a_exp, a_mnt)) {
2010     if (a_exp == FP64_EXP_INF) {
2016     a_exp += b;
2019     a_mnt = fp64_normalise(a_mnt, &a_exp);
2021     return fp64_round(a_sgn, a_exp + FP64_EXP_BITS - 3, a_mnt << 1,
2028     int a_sgn, a_exp, x_sgn, x_exp;
2032     fp16_unpack(&a_sgn, &a_exp, &a_mnt, a, mode, flags);
2035     if (fp16_is_NaN(a_exp, a_mnt))
2041     if (a_exp == FP16_EXP_INF && !a_sgn)
2048     a_mnt = fp16_normalise(a_mnt, &a_exp);
2049     if (a_exp & 1) {
2050         ++a_exp;
2064     x_exp = (a_exp + 27) >> 1;
2080     int a_sgn, a_exp, x_sgn, x_exp;
2084     fp32_unpack(&a_sgn, &a_exp, &a_mnt, a, mode, flags);
2087     if (fp32_is_NaN(a_exp, a_mnt))
2093     if (a_exp == FP32_EXP_INF && !a_sgn)
2100     a_mnt = fp32_normalise(a_mnt, &a_exp);
2101     if (!(a_exp & 1)) {
2102         ++a_exp;
2119     x_exp = (a_exp + 147) >> 1;
2135     int a_sgn, a_exp, x_sgn, x_exp, c;
2139     fp64_unpack(&a_sgn, &a_exp, &a_mnt, a, mode, flags);
2142     if (fp64_is_NaN(a_exp, a_mnt))
2148     if (a_exp == FP64_EXP_INF && !a_sgn)
2155     a_mnt = fp64_normalise(a_mnt, &a_exp);
2156     if (a_exp & 1) {
2157         ++a_exp;
2188     x_exp = (a_exp + 1053) >> 1;