arm/insts/fplib.cc

1211     int a_sgn, a_exp, b_sgn, b_exp, x_sgn, x_exp;
1247     x_exp = a_exp;
1262     x_mnt = fp16_normalise(x_mnt, &x_exp);
1264     return fp16_round(x_sgn, x_exp + FP16_EXP_BITS - 3, x_mnt << 1,
1271     int a_sgn, a_exp, b_sgn, b_exp, x_sgn, x_exp;
1307     x_exp = a_exp;
1322     x_mnt = fp32_normalise(x_mnt, &x_exp);
1324     return fp32_round(x_sgn, x_exp + FP32_EXP_BITS - 3, x_mnt << 1,
1331     int a_sgn, a_exp, b_sgn, b_exp, x_sgn, x_exp;
1367     x_exp = a_exp;
1382     x_mnt = fp64_normalise(x_mnt, &x_exp);
1384     return fp64_round(x_sgn, x_exp + FP64_EXP_BITS - 3, x_mnt << 1,
1391     int a_sgn, a_exp, b_sgn, b_exp, x_sgn, x_exp;
1415     x_exp = a_exp + b_exp - FP16_EXP_BIAS + 2 * FP16_EXP_BITS + 1;
1417     x_mnt = fp32_normalise(x_mnt, &x_exp);
1422     return fp16_round(x_sgn, x_exp, x_mnt, mode, flags);
1428     int a_sgn, a_exp, b_sgn, b_exp, x_sgn, x_exp;
1452     x_exp = a_exp + b_exp - FP32_EXP_BIAS + 2 * FP32_EXP_BITS + 1;
1454     x_mnt = fp64_normalise(x_mnt, &x_exp);
1459     return fp32_round(x_sgn, x_exp, x_mnt, mode, flags);
1465     int a_sgn, a_exp, b_sgn, b_exp, x_sgn, x_exp;
1489     x_exp = a_exp + b_exp - FP64_EXP_BIAS + 2 * FP64_EXP_BITS + 1;
1491     fp128_normalise(&x0_mnt, &x1_mnt, &x_exp);
1496     return fp64_round(x_sgn, x_exp, x0_mnt, mode, flags);
1503     int a_sgn, a_exp, b_sgn, b_exp, c_sgn, c_exp, x_sgn, x_exp, y_sgn, y_exp;
1542     x_exp = a_exp + 2 * FP16_EXP_BITS - 3;
1550         y_exp = x_exp;
1554     if (x_exp >= y_exp) {
1555         y_mnt = (lsr32(y_mnt, x_exp - y_exp) |
1556                  !!(y_mnt & (lsl32(1, x_exp - y_exp) - 1)));
1557         y_exp = x_exp;
1559         x_mnt = (lsr32(x_mnt, y_exp - x_exp) |
1560                  !!(x_mnt & (lsl32(1, y_exp - x_exp) - 1)));
1561         x_exp = y_exp;
1578     x_mnt = fp32_normalise(x_mnt, &x_exp);
1581     return fp16_round(x_sgn, x_exp + scale, x_mnt, mode, flags);
1588     int a_sgn, a_exp, b_sgn, b_exp, c_sgn, c_exp, x_sgn, x_exp, y_sgn, y_exp;
1627     x_exp = a_exp + 2 * FP32_EXP_BITS - 3;
1635         y_exp = x_exp;
1639     if (x_exp >= y_exp) {
1640         y_mnt = (lsr64(y_mnt, x_exp - y_exp) |
1641                  !!(y_mnt & (lsl64(1, x_exp - y_exp) - 1)));
1642         y_exp = x_exp;
1644         x_mnt = (lsr64(x_mnt, y_exp - x_exp) |
1645                  !!(x_mnt & (lsl64(1, y_exp - x_exp) - 1)));
1646         x_exp = y_exp;
1663     x_mnt = fp64_normalise(x_mnt, &x_exp);
1666     return fp32_round(x_sgn, x_exp + scale, x_mnt, mode, flags);
1673     int a_sgn, a_exp, b_sgn, b_exp, c_sgn, c_exp, x_sgn, x_exp, y_sgn, y_exp;
1712     x_exp = a_exp + FP64_EXP_BITS;
1721         y_exp = x_exp;
1725     if (x_exp >= y_exp) {
1728                x_exp - y_exp < 128 ? 128 - (x_exp - y_exp) : 0);
1729         lsr128(&y0_mnt, &y1_mnt, y0_mnt, y1_mnt, x_exp - y_exp);
1731         y_exp = x_exp;
1735                y_exp - x_exp < 128 ? 128 - (y_exp - x_exp) : 0);
1736         lsr128(&x0_mnt, &x1_mnt, x0_mnt, x1_mnt, y_exp - x_exp);
1738         x_exp = y_exp;
1755     fp128_normalise(&x0_mnt, &x1_mnt, &x_exp);
1758     return fp64_round(x_sgn, x_exp + scale, x0_mnt, mode, flags);
1764     int a_sgn, a_exp, b_sgn, b_exp, x_sgn, x_exp;
1791     x_exp = a_exp - b_exp + (FP16_EXP_BIAS + FP16_BITS + 2 * FP16_EXP_BITS - 3);
1797     x_mnt = fp32_normalise(x_mnt, &x_exp);
1800     return fp16_round(x_sgn, x_exp, x_mnt, mode, flags);
1806     int a_sgn, a_exp, b_sgn, b_exp, x_sgn, x_exp;
1833     x_exp = a_exp - b_exp + (FP32_EXP_BIAS + FP32_BITS + 2 * FP32_EXP_BITS - 3);
1839     x_mnt = fp64_normalise(x_mnt, &x_exp);
1842     return fp32_round(x_sgn, x_exp, x_mnt, mode, flags);
1848     int a_sgn, a_exp, b_sgn, b_exp, x_sgn, x_exp, c;
1883     x_exp = a_exp - b_exp + FP64_EXP_BIAS + 8;
1895     x_mnt = fp64_normalise(x_mnt, &x_exp);
1897     return fp64_round(x_sgn, x_exp, x_mnt << 1 | !!c, mode, flags);
2028     int a_sgn, a_exp, x_sgn, x_exp;
2064     x_exp = (a_exp + 27) >> 1;
2072     x_mnt = fp16_normalise(x_mnt, &x_exp);
2074     return fp16_round(x_sgn, x_exp, x_mnt << 1 | (t1 != t0), mode, flags);
2080     int a_sgn, a_exp, x_sgn, x_exp;
2119     x_exp = (a_exp + 147) >> 1;
2127     x_mnt = fp32_normalise(x_mnt, &x_exp);
2129     return fp32_round(x_sgn, x_exp, x_mnt << 1 | (t1 != t0), mode, flags);
2135     int a_sgn, a_exp, x_sgn, x_exp, c;
2188     x_exp = (a_exp + 1053) >> 1;
2197     x_mnt = fp64_normalise(x_mnt, &x_exp);
2199     return fp64_round(x_sgn, x_exp, x_mnt << 1 | !!c, mode, flags);
4832     int x_exp = FP16_EXP_BIAS + FP64_BITS - 1 - fbits;
4841     x_mnt = fp64_normalise(x_mnt, &x_exp);
4845     return fp16_round(x_sgn, x_exp, x_mnt, mode, flags);
4852     int x_exp = FP32_EXP_BIAS + FP64_BITS - 1 - fbits;
4861     x_mnt = fp64_normalise(x_mnt, &x_exp);
4865     return fp32_round(x_sgn, x_exp, x_mnt, mode, flags);
4872     int x_exp = FP64_EXP_BIAS + FP64_BITS - 1 - fbits;
4880     x_mnt = fp64_normalise(x_mnt, &x_exp);
4882     return fp64_round(x_sgn, x_exp, x_mnt << 1, mode, flags);