reg_constant.c 3.7 KB

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