ieee754int.h 3.5 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149
  1. /* SPDX-License-Identifier: GPL-2.0-only */
  2. /*
  3. * IEEE754 floating point
  4. * common internal header file
  5. */
  6. /*
  7. * MIPS floating point support
  8. * Copyright (C) 1994-2000 Algorithmics Ltd.
  9. */
  10. #ifndef __IEEE754INT_H
  11. #define __IEEE754INT_H
  12. #include "ieee754.h"
  13. #define CLPAIR(x, y) ((x)*6+(y))
  14. enum maddf_flags {
  15. MADDF_NEGATE_PRODUCT = 1 << 0,
  16. MADDF_NEGATE_ADDITION = 1 << 1,
  17. };
  18. static inline void ieee754_clearcx(void)
  19. {
  20. ieee754_csr.cx = 0;
  21. }
  22. static inline void ieee754_setcx(const unsigned int flags)
  23. {
  24. ieee754_csr.cx |= flags;
  25. ieee754_csr.sx |= flags;
  26. }
  27. static inline int ieee754_setandtestcx(const unsigned int x)
  28. {
  29. ieee754_setcx(x);
  30. return ieee754_csr.mx & x;
  31. }
  32. static inline int ieee754_class_nan(int xc)
  33. {
  34. return xc >= IEEE754_CLASS_SNAN;
  35. }
  36. #define COMPXSP \
  37. unsigned int xm; int xe; int xs __maybe_unused; int xc
  38. #define COMPYSP \
  39. unsigned int ym; int ye; int ys; int yc
  40. #define COMPZSP \
  41. unsigned int zm; int ze; int zs; int zc
  42. #define EXPLODESP(v, vc, vs, ve, vm) \
  43. { \
  44. vs = SPSIGN(v); \
  45. ve = SPBEXP(v); \
  46. vm = SPMANT(v); \
  47. if (ve == SP_EMAX+1+SP_EBIAS) { \
  48. if (vm == 0) \
  49. vc = IEEE754_CLASS_INF; \
  50. else if (ieee754_csr.nan2008 ^ !(vm & SP_MBIT(SP_FBITS - 1))) \
  51. vc = IEEE754_CLASS_QNAN; \
  52. else \
  53. vc = IEEE754_CLASS_SNAN; \
  54. } else if (ve == SP_EMIN-1+SP_EBIAS) { \
  55. if (vm) { \
  56. ve = SP_EMIN; \
  57. vc = IEEE754_CLASS_DNORM; \
  58. } else \
  59. vc = IEEE754_CLASS_ZERO; \
  60. } else { \
  61. ve -= SP_EBIAS; \
  62. vm |= SP_HIDDEN_BIT; \
  63. vc = IEEE754_CLASS_NORM; \
  64. } \
  65. }
  66. #define EXPLODEXSP EXPLODESP(x, xc, xs, xe, xm)
  67. #define EXPLODEYSP EXPLODESP(y, yc, ys, ye, ym)
  68. #define EXPLODEZSP EXPLODESP(z, zc, zs, ze, zm)
  69. #define COMPXDP \
  70. u64 xm; int xe; int xs __maybe_unused; int xc
  71. #define COMPYDP \
  72. u64 ym; int ye; int ys; int yc
  73. #define COMPZDP \
  74. u64 zm; int ze; int zs; int zc
  75. #define EXPLODEDP(v, vc, vs, ve, vm) \
  76. { \
  77. vm = DPMANT(v); \
  78. vs = DPSIGN(v); \
  79. ve = DPBEXP(v); \
  80. if (ve == DP_EMAX+1+DP_EBIAS) { \
  81. if (vm == 0) \
  82. vc = IEEE754_CLASS_INF; \
  83. else if (ieee754_csr.nan2008 ^ !(vm & DP_MBIT(DP_FBITS - 1))) \
  84. vc = IEEE754_CLASS_QNAN; \
  85. else \
  86. vc = IEEE754_CLASS_SNAN; \
  87. } else if (ve == DP_EMIN-1+DP_EBIAS) { \
  88. if (vm) { \
  89. ve = DP_EMIN; \
  90. vc = IEEE754_CLASS_DNORM; \
  91. } else \
  92. vc = IEEE754_CLASS_ZERO; \
  93. } else { \
  94. ve -= DP_EBIAS; \
  95. vm |= DP_HIDDEN_BIT; \
  96. vc = IEEE754_CLASS_NORM; \
  97. } \
  98. }
  99. #define EXPLODEXDP EXPLODEDP(x, xc, xs, xe, xm)
  100. #define EXPLODEYDP EXPLODEDP(y, yc, ys, ye, ym)
  101. #define EXPLODEZDP EXPLODEDP(z, zc, zs, ze, zm)
  102. #define FLUSHDP(v, vc, vs, ve, vm) \
  103. if (vc==IEEE754_CLASS_DNORM) { \
  104. if (ieee754_csr.nod) { \
  105. ieee754_setcx(IEEE754_INEXACT); \
  106. vc = IEEE754_CLASS_ZERO; \
  107. ve = DP_EMIN-1+DP_EBIAS; \
  108. vm = 0; \
  109. v = ieee754dp_zero(vs); \
  110. } \
  111. }
  112. #define FLUSHSP(v, vc, vs, ve, vm) \
  113. if (vc==IEEE754_CLASS_DNORM) { \
  114. if (ieee754_csr.nod) { \
  115. ieee754_setcx(IEEE754_INEXACT); \
  116. vc = IEEE754_CLASS_ZERO; \
  117. ve = SP_EMIN-1+SP_EBIAS; \
  118. vm = 0; \
  119. v = ieee754sp_zero(vs); \
  120. } \
  121. }
  122. #define FLUSHXDP FLUSHDP(x, xc, xs, xe, xm)
  123. #define FLUSHYDP FLUSHDP(y, yc, ys, ye, ym)
  124. #define FLUSHZDP FLUSHDP(z, zc, zs, ze, zm)
  125. #define FLUSHXSP FLUSHSP(x, xc, xs, xe, xm)
  126. #define FLUSHYSP FLUSHSP(y, yc, ys, ye, ym)
  127. #define FLUSHZSP FLUSHSP(z, zc, zs, ze, zm)
  128. #endif /* __IEEE754INT_H */