^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 1) // SPDX-License-Identifier: GPL-2.0
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 2) /*---------------------------------------------------------------------------+
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 3) | reg_constant.c |
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 4) | |
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 5) | All of the constant FPU_REGs |
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 6) | |
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 7) | Copyright (C) 1992,1993,1994,1997 |
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 8) | W. Metzenthen, 22 Parker St, Ormond, Vic 3163, |
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 9) | Australia. E-mail billm@suburbia.net |
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 10) | |
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 11) | |
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 12) +---------------------------------------------------------------------------*/
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 13)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 14) #include "fpu_system.h"
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 15) #include "fpu_emu.h"
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 16) #include "status_w.h"
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 17) #include "reg_constant.h"
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 18) #include "control_w.h"
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 19)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 20) #define MAKE_REG(s, e, l, h) { l, h, \
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 21) (u16)((EXTENDED_Ebias+(e)) | ((SIGN_##s != 0)*0x8000)) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 22)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 23) FPU_REG const CONST_1 = MAKE_REG(POS, 0, 0x00000000, 0x80000000);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 24) #if 0
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 25) FPU_REG const CONST_2 = MAKE_REG(POS, 1, 0x00000000, 0x80000000);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 26) FPU_REG const CONST_HALF = MAKE_REG(POS, -1, 0x00000000, 0x80000000);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 27) #endif /* 0 */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 28) static FPU_REG const CONST_L2T = MAKE_REG(POS, 1, 0xcd1b8afe, 0xd49a784b);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 29) static FPU_REG const CONST_L2E = MAKE_REG(POS, 0, 0x5c17f0bc, 0xb8aa3b29);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 30) FPU_REG const CONST_PI = MAKE_REG(POS, 1, 0x2168c235, 0xc90fdaa2);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 31) FPU_REG const CONST_PI2 = MAKE_REG(POS, 0, 0x2168c235, 0xc90fdaa2);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 32) FPU_REG const CONST_PI4 = MAKE_REG(POS, -1, 0x2168c235, 0xc90fdaa2);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 33) static FPU_REG const CONST_LG2 = MAKE_REG(POS, -2, 0xfbcff799, 0x9a209a84);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 34) static FPU_REG const CONST_LN2 = MAKE_REG(POS, -1, 0xd1cf79ac, 0xb17217f7);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 35)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 36) /* Extra bits to take pi/2 to more than 128 bits precision. */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 37) FPU_REG const CONST_PI2extra = MAKE_REG(NEG, -66,
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 38) 0xfc8f8cbb, 0xece675d1);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 39)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 40) /* Only the sign (and tag) is used in internal zeroes */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 41) FPU_REG const CONST_Z = MAKE_REG(POS, EXP_UNDER, 0x0, 0x0);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 42)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 43) /* Only the sign and significand (and tag) are used in internal NaNs */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 44) /* The 80486 never generates one of these
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 45) FPU_REG const CONST_SNAN = MAKE_REG(POS, EXP_OVER, 0x00000001, 0x80000000);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 46) */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 47) /* This is the real indefinite QNaN */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 48) FPU_REG const CONST_QNaN = MAKE_REG(NEG, EXP_OVER, 0x00000000, 0xC0000000);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 49)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 50) /* Only the sign (and tag) is used in internal infinities */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 51) FPU_REG const CONST_INF = MAKE_REG(POS, EXP_OVER, 0x00000000, 0x80000000);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 52)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 53) static void fld_const(FPU_REG const * c, int adj, u_char tag)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 54) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 55) FPU_REG *st_new_ptr;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 56)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 57) if (STACK_OVERFLOW) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 58) FPU_stack_overflow();
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 59) return;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 60) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 61) push();
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 62) reg_copy(c, st_new_ptr);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 63) st_new_ptr->sigl += adj; /* For all our fldxxx constants, we don't need to
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 64) borrow or carry. */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 65) FPU_settag0(tag);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 66) clear_C1();
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 67) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 68)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 69) /* A fast way to find out whether x is one of RC_DOWN or RC_CHOP
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 70) (and not one of RC_RND or RC_UP).
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 71) */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 72) #define DOWN_OR_CHOP(x) (x & RC_DOWN)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 73)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 74) static void fld1(int rc)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 75) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 76) fld_const(&CONST_1, 0, TAG_Valid);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 77) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 78)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 79) static void fldl2t(int rc)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 80) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 81) fld_const(&CONST_L2T, (rc == RC_UP) ? 1 : 0, TAG_Valid);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 82) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 83)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 84) static void fldl2e(int rc)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 85) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 86) fld_const(&CONST_L2E, DOWN_OR_CHOP(rc) ? -1 : 0, TAG_Valid);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 87) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 88)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 89) static void fldpi(int rc)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 90) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 91) fld_const(&CONST_PI, DOWN_OR_CHOP(rc) ? -1 : 0, TAG_Valid);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 92) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 93)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 94) static void fldlg2(int rc)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 95) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 96) fld_const(&CONST_LG2, DOWN_OR_CHOP(rc) ? -1 : 0, TAG_Valid);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 97) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 98)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 99) static void fldln2(int rc)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 100) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 101) fld_const(&CONST_LN2, DOWN_OR_CHOP(rc) ? -1 : 0, TAG_Valid);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 102) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 103)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 104) static void fldz(int rc)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 105) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 106) fld_const(&CONST_Z, 0, TAG_Zero);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 107) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 108)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 109) typedef void (*FUNC_RC) (int);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 110)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 111) static FUNC_RC constants_table[] = {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 112) fld1, fldl2t, fldl2e, fldpi, fldlg2, fldln2, fldz, (FUNC_RC) FPU_illegal
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 113) };
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 114)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 115) void fconst(void)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 116) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 117) (constants_table[FPU_rm]) (control_word & CW_RC);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 118) }