int fsqrts(void *frD, void *frB) { FP_DECL_D(B); FP_DECL_D(R); FP_DECL_EX; #ifdef DEBUG printk("%s: %p %p %p %p\n", __func__, frD, frB); #endif FP_UNPACK_DP(B, frB); #ifdef DEBUG printk("B: %ld %lu %lu %ld (%ld)\n", B_s, B_f1, B_f0, B_e, B_c); #endif if (B_s && B_c != FP_CLS_ZERO) FP_SET_EXCEPTION(EFLAG_VXSQRT); if (B_c == FP_CLS_NAN) FP_SET_EXCEPTION(EFLAG_VXSNAN); FP_SQRT_D(R, B); #ifdef DEBUG printk("R: %ld %lu %lu %ld (%ld)\n", R_s, R_f1, R_f0, R_e, R_c); #endif __FP_PACK_DS(frD, R); return FP_CUR_EXCEPTIONS; }
int fdivs(void *frD, void *frA, void *frB) { FP_DECL_D(A); FP_DECL_D(B); FP_DECL_D(R); FP_DECL_EX; #ifdef DEBUG printk("%s: %p %p %p\n", __func__, frD, frA, frB); #endif FP_UNPACK_DP(A, frA); FP_UNPACK_DP(B, frB); #ifdef DEBUG printk("A: %ld %lu %lu %ld (%ld)\n", A_s, A_f1, A_f0, A_e, A_c); printk("B: %ld %lu %lu %ld (%ld)\n", B_s, B_f1, B_f0, B_e, B_c); #endif if (A_c == FP_CLS_ZERO && B_c == FP_CLS_ZERO) { FP_SET_EXCEPTION(EFLAG_VXZDZ); #ifdef DEBUG printk("%s: FPSCR_VXZDZ raised\n", __func__); #endif } if (A_c == FP_CLS_INF && B_c == FP_CLS_INF) { FP_SET_EXCEPTION(EFLAG_VXIDI); #ifdef DEBUG printk("%s: FPSCR_VXIDI raised\n", __func__); #endif } if (B_c == FP_CLS_ZERO && A_c != FP_CLS_ZERO) { FP_SET_EXCEPTION(EFLAG_DIVZERO); if (__FPU_TRAP_P(EFLAG_DIVZERO)) return FP_CUR_EXCEPTIONS; } FP_DIV_D(R, A, B); #ifdef DEBUG printk("D: %ld %lu %lu %ld (%ld)\n", R_s, R_f1, R_f0, R_e, R_c); #endif __FP_PACK_DS(frD, R); return FP_CUR_EXCEPTIONS; }
int fnmsub(void *frD, void *frA, void *frB, void *frC) { FP_DECL_D(R); FP_DECL_D(A); FP_DECL_D(B); FP_DECL_D(C); FP_DECL_D(T); FP_DECL_EX; #ifdef DEBUG printk("%s: %p %p %p %p\n", __func__, frD, frA, frB, frC); #endif FP_UNPACK_DP(A, frA); FP_UNPACK_DP(B, frB); FP_UNPACK_DP(C, frC); #ifdef DEBUG printk("A: %ld %lu %lu %ld (%ld)\n", A_s, A_f1, A_f0, A_e, A_c); printk("B: %ld %lu %lu %ld (%ld)\n", B_s, B_f1, B_f0, B_e, B_c); printk("C: %ld %lu %lu %ld (%ld)\n", C_s, C_f1, C_f0, C_e, C_c); #endif if ((A_c == FP_CLS_INF && C_c == FP_CLS_ZERO) || (A_c == FP_CLS_ZERO && C_c == FP_CLS_INF)) FP_SET_EXCEPTION(EFLAG_VXIMZ); FP_MUL_D(T, A, C); if (B_c != FP_CLS_NAN) B_s ^= 1; if (T_s != B_s && T_c == FP_CLS_INF && B_c == FP_CLS_INF) FP_SET_EXCEPTION(EFLAG_VXISI); FP_ADD_D(R, T, B); if (R_c != FP_CLS_NAN) R_s ^= 1; #ifdef DEBUG printk("D: %ld %lu %lu %ld (%ld)\n", R_s, R_f1, R_f0, R_e, R_c); #endif __FP_PACK_D(frD, R); return FP_CUR_EXCEPTIONS; }
static long internal_equality (long al, long ah, long bl, long bh, long neq) { FP_DECL_EX; FP_DECL_Q(A); FP_DECL_Q(B); long r; FP_UNPACK_RAW_Q(A, a); FP_UNPACK_RAW_Q(B, b); if ((A_e == _FP_EXPMAX_Q && !_FP_FRAC_ZEROP_2(A)) || (B_e == _FP_EXPMAX_Q && !_FP_FRAC_ZEROP_2(B))) { /* EQ and NE signal invalid operation only if either operand is SNaN. */ if (FP_ISSIGNAN_Q(A) || FP_ISSIGNAN_Q(B)) { FP_SET_EXCEPTION(FP_EX_INVALID); FP_HANDLE_EXCEPTIONS; } return -1; } r = (A_e == B_e && _FP_FRAC_EQ_2 (A, B) && (A_s == B_s || (!A_e && _FP_FRAC_ZEROP_2(A)))); r ^= neq; return r; }
CMPtype __c6xabi_ltd(DFtype a, DFtype b) { FP_DECL_EX; FP_DECL_D(A); FP_DECL_D(B); CMPtype r; FP_UNPACK_RAW_D(A, a); FP_UNPACK_RAW_D(B, b); FP_CMP_D(r, A, B, 2); if (r == 2 && (FP_ISSIGNAN_D(A) || FP_ISSIGNAN_D(B))) FP_SET_EXCEPTION(FP_EX_INVALID); FP_HANDLE_EXCEPTIONS; return r < 0; }
int _Q_fne(const long double a, const long double b) { FP_DECL_EX; FP_DECL_Q(A); FP_DECL_Q(B); int r; FP_UNPACK_RAW_Q(A, a); FP_UNPACK_RAW_Q(B, b); FP_CMP_EQ_Q(r, A, B); if (r && (FP_ISSIGNAN_Q(A) || FP_ISSIGNAN_Q(B))) FP_SET_EXCEPTION(FP_EX_INVALID); FP_HANDLE_EXCEPTIONS; return r; }
int __eqsf2(SFtype a, SFtype b) { FP_DECL_EX; FP_DECL_S(A); FP_DECL_S(B); int r; FP_UNPACK_RAW_S(A, a); FP_UNPACK_RAW_S(B, b); FP_CMP_EQ_S(r, A, B); if (r && (FP_ISSIGNAN_S(A) || FP_ISSIGNAN_S(B))) FP_SET_EXCEPTION(FP_EX_INVALID); FP_HANDLE_EXCEPTIONS; return r; }
CMPtype __letf2(TFtype a, TFtype b) { FP_DECL_EX; FP_DECL_Q(A); FP_DECL_Q(B); CMPtype r; FP_UNPACK_RAW_Q(A, a); FP_UNPACK_RAW_Q(B, b); FP_CMP_Q(r, A, B, 2); if (r == 2 && (FP_ISSIGNAN_Q(A) || FP_ISSIGNAN_Q(B))) FP_SET_EXCEPTION(FP_EX_INVALID); FP_HANDLE_EXCEPTIONS; return r; }
int __gedf2(double a, double b) { FP_DECL_EX; FP_DECL_D(A); FP_DECL_D(B); int r; FP_UNPACK_RAW_D(A, a); FP_UNPACK_RAW_D(B, b); FP_CMP_D(r, A, B, -2); if (r == -2 && (FP_ISSIGNAN_D(A) || FP_ISSIGNAN_D(B))) FP_SET_EXCEPTION(FP_EX_INVALID); FP_HANDLE_EXCEPTIONS; return r; }
int _Q_flt(const long double a, const long double b) { FP_DECL_EX; FP_DECL_Q(A); FP_DECL_Q(B); int r; FP_UNPACK_RAW_Q(A, a); FP_UNPACK_RAW_Q(B, b); FP_CMP_Q(r, B, A, 3); if (r == 3) FP_SET_EXCEPTION(FP_EX_INVALID); FP_HANDLE_EXCEPTIONS; return (r == 1); }
CMPtype __eqdf2(DFtype a, DFtype b) { FP_DECL_EX; FP_DECL_D(A); FP_DECL_D(B); CMPtype r; FP_UNPACK_RAW_D(A, a); FP_UNPACK_RAW_D(B, b); FP_CMP_EQ_D(r, A, B); if (r && (FP_ISSIGNAN_D(A) || FP_ISSIGNAN_D(B))) FP_SET_EXCEPTION(FP_EX_INVALID); FP_HANDLE_EXCEPTIONS; return r; }
CMPtype __c6xabi_eqf(SFtype a, SFtype b) { FP_DECL_EX; FP_DECL_S(A); FP_DECL_S(B); CMPtype r; FP_UNPACK_RAW_S(A, a); FP_UNPACK_RAW_S(B, b); FP_CMP_EQ_S(r, A, B); if (r && (FP_ISSIGNAN_S(A) || FP_ISSIGNAN_S(B))) FP_SET_EXCEPTION(FP_EX_INVALID); FP_HANDLE_EXCEPTIONS; return !r; }
int _Q_cmp(const long double a, const long double b) { FP_DECL_EX; FP_DECL_Q(A); FP_DECL_Q(B); int r; FP_UNPACK_RAW_Q(A, a); FP_UNPACK_RAW_Q(B, b); FP_CMP_Q(r, B, A, 3); if (r == -1) r = 2; if (r == 3 && (FP_ISSIGNAN_Q(A) || FP_ISSIGNAN_Q(B))) FP_SET_EXCEPTION(FP_EX_INVALID); FP_HANDLE_EXCEPTIONS; return r; }
CMPtype __gedf2 (DFtype a, DFtype b) { FP_DECL_EX; FP_DECL_D (A); FP_DECL_D (B); CMPtype r; FP_INIT_EXCEPTIONS; FP_UNPACK_RAW_D (A, a); FP_UNPACK_RAW_D (B, b); FP_CMP_D (r, A, B, -2); if (r == -2) FP_SET_EXCEPTION (FP_EX_INVALID); FP_HANDLE_EXCEPTIONS; return r; }
void fdivs(void *ft, void *fa, void *fb) { FP_DECL_S(A); FP_DECL_S(B); FP_DECL_S(R); FP_DECL_EX; FP_UNPACK_SP(A, fa); FP_UNPACK_SP(B, fb); if (B_c == FP_CLS_ZERO && A_c != FP_CLS_ZERO) FP_SET_EXCEPTION(FP_EX_DIVZERO); FP_DIV_S(R, A, B); FP_PACK_SP(ft, R); __FPU_FPCSR |= FP_CUR_EXCEPTIONS; }
int fmuls(void *frD, void *frA, void *frB) { FP_DECL_D(A); FP_DECL_D(B); FP_DECL_D(R); FP_DECL_EX; #ifdef DEBUG printk("%s: %p %p %p\n", __func__, frD, frA, frB); #endif FP_UNPACK_DP(A, frA); FP_UNPACK_DP(B, frB); #ifdef DEBUG printk("A: %ld %lu %lu %ld (%ld) [%08lx.%08lx %lx]\n", A_s, A_f1, A_f0, A_e, A_c, A_f1, A_f0, A_e + 1023); printk("B: %ld %lu %lu %ld (%ld) [%08lx.%08lx %lx]\n", B_s, B_f1, B_f0, B_e, B_c, B_f1, B_f0, B_e + 1023); #endif if ((A_c == FP_CLS_INF && B_c == FP_CLS_ZERO) || (A_c == FP_CLS_ZERO && B_c == FP_CLS_INF)) FP_SET_EXCEPTION(EFLAG_VXIMZ); FP_MUL_D(R, A, B); #ifdef DEBUG printk("D: %ld %lu %lu %ld (%ld) [%08lx.%08lx %lx]\n", R_s, R_f1, R_f0, R_e, R_c, R_f1, R_f0, R_e + 1023); #endif __FP_PACK_DS(frD, R); return FP_CUR_EXCEPTIONS; }
int fsubs(void *frD, void *frA, void *frB) { FP_DECL_D(A); FP_DECL_D(B); FP_DECL_D(R); FP_DECL_EX; #ifdef DEBUG printk("%s: %p %p %p\n", __func__, frD, frA, frB); #endif FP_UNPACK_DP(A, frA); FP_UNPACK_DP(B, frB); #ifdef DEBUG printk("A: %ld %lu %lu %ld (%ld)\n", A_s, A_f1, A_f0, A_e, A_c); printk("B: %ld %lu %lu %ld (%ld)\n", B_s, B_f1, B_f0, B_e, B_c); #endif if (B_c != FP_CLS_NAN) B_s ^= 1; if (A_s != B_s && A_c == FP_CLS_INF && B_c == FP_CLS_INF) FP_SET_EXCEPTION(EFLAG_VXISI); FP_ADD_D(R, A, B); #ifdef DEBUG printk("D: %ld %lu %lu %ld (%ld)\n", R_s, R_f1, R_f0, R_e, R_c); #endif __FP_PACK_DS(frD, R); return FP_CUR_EXCEPTIONS; }