arm/insts/fplib.cc

85 #define FP64_MANT_BITS (FP64_BITS - FP64_EXP_BITS - 1)
89 #define FP64_EXP(x) ((x) >> FP64_MANT_BITS & ((1ULL << FP64_EXP_BITS) - 1))
93 #define FP64_MANT(x) ((x) & ((1ULL << FP64_MANT_BITS) - 1))
314     return sgn << (FP64_BITS - 1) | exp << FP64_MANT_BITS | FP64_MANT(mnt);
386     return fp64_pack(0, FP64_EXP_INF, 1ULL << (FP64_MANT_BITS - 1));
438         *mnt |= 1ULL << FP64_MANT_BITS;
481     return fp64_is_NaN(exp, mnt) && !(mnt >> (FP64_MANT_BITS - 1) & 1);
499     return exp == FP64_EXP_INF && (mnt >> (FP64_MANT_BITS - 1) & 1);
543     if (!(a >> (FP64_MANT_BITS - 1) & 1)) {
545         a |= 1ULL << (FP64_MANT_BITS - 1);
871     uint64_t int_mant; // mantissa for result, less than (2 << FP64_MANT_BITS)
883     mnt = (4ULL << FP64_MANT_BITS | mnt >> (FP64_EXP_BITS - 1) |
905         if (int_mant == 1ULL << FP64_MANT_BITS) {
909         if (int_mant == 2ULL << FP64_MANT_BITS) {
2522                      op >> (FP64_MANT_BITS - FP16_MANT_BITS));
2538                      op >> (FP64_MANT_BITS - FP32_MANT_BITS));
2545                      1ULL << (FP64_MANT_BITS - 1) |
2546                      (uint64_t)op << (FP64_MANT_BITS - FP16_MANT_BITS));
2553                      1ULL << (FP64_MANT_BITS - 1) |
2554                      (uint64_t)op << (FP64_MANT_BITS - FP32_MANT_BITS));
2572     return fp64_pack(sgn, FP64_EXP_BIAS, 1ULL << (FP64_MANT_BITS - 1));
2590     return fp64_pack(sgn, FP64_EXP_BIAS + 1, 1ULL << (FP64_MANT_BITS - 1));
2683         if (!(mnt >> (FP64_MANT_BITS - 1) & 1) || alt_hp) {
2699                         mnt >> (FP64_MANT_BITS - FP16_BITS) |
2700                         !!(mnt & ((1ULL << (FP64_MANT_BITS - FP16_BITS)) - 1)),
2766         if (!(mnt >> (FP64_MANT_BITS - 1) & 1)) {
2776                         mnt >> (FP64_MANT_BITS - FP32_BITS) |
2777                         !!(mnt & ((1ULL << (FP64_MANT_BITS - FP32_BITS)) - 1)),
2816                            (uint64_t)mnt << (FP64_MANT_BITS - FP16_BITS + 1));
2854                            (uint64_t)mnt << (FP64_MANT_BITS - FP32_BITS + 1));
3108     return ((((op >> 6) & ((1 << FP64_EXP_BITS) - 1)) << FP64_MANT_BITS) |
3127     return fp64_pack(sgn, mnt >> FP64_MANT_BITS ? exp : 0, mnt);
3158     if (!((uint64_t)~(*op1 << 1) >> FP64_MANT_BITS) &&
3159         (uint64_t)~(*op2 << 1) >> FP64_MANT_BITS)
3161     if (!((uint64_t)~(*op2 << 1) >> FP64_MANT_BITS) &&
3162         (uint64_t)~(*op1 << 1) >> FP64_MANT_BITS)
3612                            mnt << (FP64_MANT_BITS - 8));
3852     } else if (!((uint64_t)(op << 1) >> (FP64_MANT_BITS - 1))) {
3881         fraction <<= FP64_MANT_BITS - 8;
4201     int expint = FP64_EXP_BIAS + FP64_MANT_BITS;
4539         (uint64_t)FP64_EXP_BIAS << FP64_MANT_BITS; // 1.0
4643                          (uint64_t)mnt << (FP64_MANT_BITS - FP16_MANT_BITS),
4675                          (uint64_t)mnt << (FP64_MANT_BITS - FP16_MANT_BITS),
4704                          (uint64_t)mnt << (FP64_MANT_BITS - FP32_MANT_BITS),
4763                          (uint64_t)mnt << (FP64_MANT_BITS - FP16_MANT_BITS),
4793                          (uint64_t)mnt << (FP64_MANT_BITS - FP32_MANT_BITS),