ieee754.h 9.8 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311
  1. /* SPDX-License-Identifier: GPL-2.0-only */
  2. /*
  3. * MIPS floating point support
  4. * Copyright (C) 1994-2000 Algorithmics Ltd.
  5. *
  6. * Nov 7, 2000
  7. * Modification to allow integration with Linux kernel
  8. *
  9. * Kevin D. Kissell, kevink@mips.com and Carsten Langgard, carstenl@mips.com
  10. * Copyright (C) 2000 MIPS Technologies, Inc. All rights reserved.
  11. */
  12. #ifndef __ARCH_MIPS_MATH_EMU_IEEE754_H
  13. #define __ARCH_MIPS_MATH_EMU_IEEE754_H
  14. #include <linux/compiler.h>
  15. #include <asm/byteorder.h>
  16. #include <linux/kernel.h>
  17. #include <linux/types.h>
  18. #include <linux/sched.h>
  19. #include <asm/bitfield.h>
  20. union ieee754dp {
  21. struct {
  22. __BITFIELD_FIELD(unsigned int sign:1,
  23. __BITFIELD_FIELD(unsigned int bexp:11,
  24. __BITFIELD_FIELD(u64 mant:52,
  25. ;)))
  26. };
  27. u64 bits;
  28. };
  29. union ieee754sp {
  30. struct {
  31. __BITFIELD_FIELD(unsigned sign:1,
  32. __BITFIELD_FIELD(unsigned bexp:8,
  33. __BITFIELD_FIELD(unsigned mant:23,
  34. ;)))
  35. };
  36. u32 bits;
  37. };
  38. /*
  39. * single precision (often aka float)
  40. */
  41. int ieee754sp_class(union ieee754sp x);
  42. union ieee754sp ieee754sp_abs(union ieee754sp x);
  43. union ieee754sp ieee754sp_neg(union ieee754sp x);
  44. union ieee754sp ieee754sp_add(union ieee754sp x, union ieee754sp y);
  45. union ieee754sp ieee754sp_sub(union ieee754sp x, union ieee754sp y);
  46. union ieee754sp ieee754sp_mul(union ieee754sp x, union ieee754sp y);
  47. union ieee754sp ieee754sp_div(union ieee754sp x, union ieee754sp y);
  48. union ieee754sp ieee754sp_fint(int x);
  49. union ieee754sp ieee754sp_flong(s64 x);
  50. union ieee754sp ieee754sp_fdp(union ieee754dp x);
  51. union ieee754sp ieee754sp_rint(union ieee754sp x);
  52. int ieee754sp_tint(union ieee754sp x);
  53. s64 ieee754sp_tlong(union ieee754sp x);
  54. int ieee754sp_cmp(union ieee754sp x, union ieee754sp y, int cop, int sig);
  55. union ieee754sp ieee754sp_sqrt(union ieee754sp x);
  56. union ieee754sp ieee754sp_maddf(union ieee754sp z, union ieee754sp x,
  57. union ieee754sp y);
  58. union ieee754sp ieee754sp_msubf(union ieee754sp z, union ieee754sp x,
  59. union ieee754sp y);
  60. union ieee754sp ieee754sp_madd(union ieee754sp z, union ieee754sp x,
  61. union ieee754sp y);
  62. union ieee754sp ieee754sp_msub(union ieee754sp z, union ieee754sp x,
  63. union ieee754sp y);
  64. union ieee754sp ieee754sp_nmadd(union ieee754sp z, union ieee754sp x,
  65. union ieee754sp y);
  66. union ieee754sp ieee754sp_nmsub(union ieee754sp z, union ieee754sp x,
  67. union ieee754sp y);
  68. int ieee754sp_2008class(union ieee754sp x);
  69. union ieee754sp ieee754sp_fmin(union ieee754sp x, union ieee754sp y);
  70. union ieee754sp ieee754sp_fmina(union ieee754sp x, union ieee754sp y);
  71. union ieee754sp ieee754sp_fmax(union ieee754sp x, union ieee754sp y);
  72. union ieee754sp ieee754sp_fmaxa(union ieee754sp x, union ieee754sp y);
  73. /*
  74. * double precision (often aka double)
  75. */
  76. int ieee754dp_class(union ieee754dp x);
  77. union ieee754dp ieee754dp_add(union ieee754dp x, union ieee754dp y);
  78. union ieee754dp ieee754dp_sub(union ieee754dp x, union ieee754dp y);
  79. union ieee754dp ieee754dp_mul(union ieee754dp x, union ieee754dp y);
  80. union ieee754dp ieee754dp_div(union ieee754dp x, union ieee754dp y);
  81. union ieee754dp ieee754dp_abs(union ieee754dp x);
  82. union ieee754dp ieee754dp_neg(union ieee754dp x);
  83. union ieee754dp ieee754dp_fint(int x);
  84. union ieee754dp ieee754dp_flong(s64 x);
  85. union ieee754dp ieee754dp_fsp(union ieee754sp x);
  86. union ieee754dp ieee754dp_rint(union ieee754dp x);
  87. int ieee754dp_tint(union ieee754dp x);
  88. s64 ieee754dp_tlong(union ieee754dp x);
  89. int ieee754dp_cmp(union ieee754dp x, union ieee754dp y, int cop, int sig);
  90. union ieee754dp ieee754dp_sqrt(union ieee754dp x);
  91. union ieee754dp ieee754dp_maddf(union ieee754dp z, union ieee754dp x,
  92. union ieee754dp y);
  93. union ieee754dp ieee754dp_msubf(union ieee754dp z, union ieee754dp x,
  94. union ieee754dp y);
  95. union ieee754dp ieee754dp_madd(union ieee754dp z, union ieee754dp x,
  96. union ieee754dp y);
  97. union ieee754dp ieee754dp_msub(union ieee754dp z, union ieee754dp x,
  98. union ieee754dp y);
  99. union ieee754dp ieee754dp_nmadd(union ieee754dp z, union ieee754dp x,
  100. union ieee754dp y);
  101. union ieee754dp ieee754dp_nmsub(union ieee754dp z, union ieee754dp x,
  102. union ieee754dp y);
  103. int ieee754dp_2008class(union ieee754dp x);
  104. union ieee754dp ieee754dp_fmin(union ieee754dp x, union ieee754dp y);
  105. union ieee754dp ieee754dp_fmina(union ieee754dp x, union ieee754dp y);
  106. union ieee754dp ieee754dp_fmax(union ieee754dp x, union ieee754dp y);
  107. union ieee754dp ieee754dp_fmaxa(union ieee754dp x, union ieee754dp y);
  108. /* 5 types of floating point number
  109. */
  110. enum {
  111. IEEE754_CLASS_NORM = 0x00,
  112. IEEE754_CLASS_ZERO = 0x01,
  113. IEEE754_CLASS_DNORM = 0x02,
  114. IEEE754_CLASS_INF = 0x03,
  115. IEEE754_CLASS_SNAN = 0x04,
  116. IEEE754_CLASS_QNAN = 0x05,
  117. };
  118. /* exception numbers */
  119. #define IEEE754_INEXACT 0x01
  120. #define IEEE754_UNDERFLOW 0x02
  121. #define IEEE754_OVERFLOW 0x04
  122. #define IEEE754_ZERO_DIVIDE 0x08
  123. #define IEEE754_INVALID_OPERATION 0x10
  124. /* cmp operators
  125. */
  126. #define IEEE754_CLT 0x01
  127. #define IEEE754_CEQ 0x02
  128. #define IEEE754_CGT 0x04
  129. #define IEEE754_CUN 0x08
  130. /*
  131. * The control status register
  132. */
  133. struct _ieee754_csr {
  134. __BITFIELD_FIELD(unsigned fcc:7, /* condition[7:1] */
  135. __BITFIELD_FIELD(unsigned nod:1, /* set 1 for no denormals */
  136. __BITFIELD_FIELD(unsigned c:1, /* condition[0] */
  137. __BITFIELD_FIELD(unsigned pad0:3,
  138. __BITFIELD_FIELD(unsigned abs2008:1, /* IEEE 754-2008 ABS/NEG.fmt */
  139. __BITFIELD_FIELD(unsigned nan2008:1, /* IEEE 754-2008 NaN mode */
  140. __BITFIELD_FIELD(unsigned cx:6, /* exceptions this operation */
  141. __BITFIELD_FIELD(unsigned mx:5, /* exception enable mask */
  142. __BITFIELD_FIELD(unsigned sx:5, /* exceptions total */
  143. __BITFIELD_FIELD(unsigned rm:2, /* current rounding mode */
  144. ;))))))))))
  145. };
  146. #define ieee754_csr (*(struct _ieee754_csr *)(&current->thread.fpu.fcr31))
  147. static inline unsigned int ieee754_getrm(void)
  148. {
  149. return (ieee754_csr.rm);
  150. }
  151. static inline unsigned int ieee754_setrm(unsigned int rm)
  152. {
  153. return (ieee754_csr.rm = rm);
  154. }
  155. /*
  156. * get current exceptions
  157. */
  158. static inline unsigned int ieee754_getcx(void)
  159. {
  160. return (ieee754_csr.cx);
  161. }
  162. /* test for current exception condition
  163. */
  164. static inline int ieee754_cxtest(unsigned int n)
  165. {
  166. return (ieee754_csr.cx & n);
  167. }
  168. /*
  169. * get sticky exceptions
  170. */
  171. static inline unsigned int ieee754_getsx(void)
  172. {
  173. return (ieee754_csr.sx);
  174. }
  175. /* clear sticky conditions
  176. */
  177. static inline unsigned int ieee754_clrsx(void)
  178. {
  179. return (ieee754_csr.sx = 0);
  180. }
  181. /* test for sticky exception condition
  182. */
  183. static inline int ieee754_sxtest(unsigned int n)
  184. {
  185. return (ieee754_csr.sx & n);
  186. }
  187. /* debugging */
  188. union ieee754sp ieee754sp_dump(char *s, union ieee754sp x);
  189. union ieee754dp ieee754dp_dump(char *s, union ieee754dp x);
  190. #define IEEE754_SPCVAL_PZERO 0 /* +0.0 */
  191. #define IEEE754_SPCVAL_NZERO 1 /* -0.0 */
  192. #define IEEE754_SPCVAL_PONE 2 /* +1.0 */
  193. #define IEEE754_SPCVAL_NONE 3 /* -1.0 */
  194. #define IEEE754_SPCVAL_PTEN 4 /* +10.0 */
  195. #define IEEE754_SPCVAL_NTEN 5 /* -10.0 */
  196. #define IEEE754_SPCVAL_PINFINITY 6 /* +inf */
  197. #define IEEE754_SPCVAL_NINFINITY 7 /* -inf */
  198. #define IEEE754_SPCVAL_INDEF_LEG 8 /* legacy quiet NaN */
  199. #define IEEE754_SPCVAL_INDEF_2008 9 /* IEEE 754-2008 quiet NaN */
  200. #define IEEE754_SPCVAL_PMAX 10 /* +max norm */
  201. #define IEEE754_SPCVAL_NMAX 11 /* -max norm */
  202. #define IEEE754_SPCVAL_PMIN 12 /* +min norm */
  203. #define IEEE754_SPCVAL_NMIN 13 /* -min norm */
  204. #define IEEE754_SPCVAL_PMIND 14 /* +min denorm */
  205. #define IEEE754_SPCVAL_NMIND 15 /* -min denorm */
  206. #define IEEE754_SPCVAL_P1E31 16 /* + 1.0e31 */
  207. #define IEEE754_SPCVAL_P1E63 17 /* + 1.0e63 */
  208. extern const union ieee754dp __ieee754dp_spcvals[];
  209. extern const union ieee754sp __ieee754sp_spcvals[];
  210. #define ieee754dp_spcvals ((const union ieee754dp *)__ieee754dp_spcvals)
  211. #define ieee754sp_spcvals ((const union ieee754sp *)__ieee754sp_spcvals)
  212. /*
  213. * Return infinity with given sign
  214. */
  215. #define ieee754dp_inf(sn) (ieee754dp_spcvals[IEEE754_SPCVAL_PINFINITY+(sn)])
  216. #define ieee754dp_zero(sn) (ieee754dp_spcvals[IEEE754_SPCVAL_PZERO+(sn)])
  217. #define ieee754dp_one(sn) (ieee754dp_spcvals[IEEE754_SPCVAL_PONE+(sn)])
  218. #define ieee754dp_ten(sn) (ieee754dp_spcvals[IEEE754_SPCVAL_PTEN+(sn)])
  219. #define ieee754dp_indef() (ieee754dp_spcvals[IEEE754_SPCVAL_INDEF_LEG + \
  220. ieee754_csr.nan2008])
  221. #define ieee754dp_max(sn) (ieee754dp_spcvals[IEEE754_SPCVAL_PMAX+(sn)])
  222. #define ieee754dp_min(sn) (ieee754dp_spcvals[IEEE754_SPCVAL_PMIN+(sn)])
  223. #define ieee754dp_mind(sn) (ieee754dp_spcvals[IEEE754_SPCVAL_PMIND+(sn)])
  224. #define ieee754dp_1e31() (ieee754dp_spcvals[IEEE754_SPCVAL_P1E31])
  225. #define ieee754dp_1e63() (ieee754dp_spcvals[IEEE754_SPCVAL_P1E63])
  226. #define ieee754sp_inf(sn) (ieee754sp_spcvals[IEEE754_SPCVAL_PINFINITY+(sn)])
  227. #define ieee754sp_zero(sn) (ieee754sp_spcvals[IEEE754_SPCVAL_PZERO+(sn)])
  228. #define ieee754sp_one(sn) (ieee754sp_spcvals[IEEE754_SPCVAL_PONE+(sn)])
  229. #define ieee754sp_ten(sn) (ieee754sp_spcvals[IEEE754_SPCVAL_PTEN+(sn)])
  230. #define ieee754sp_indef() (ieee754sp_spcvals[IEEE754_SPCVAL_INDEF_LEG + \
  231. ieee754_csr.nan2008])
  232. #define ieee754sp_max(sn) (ieee754sp_spcvals[IEEE754_SPCVAL_PMAX+(sn)])
  233. #define ieee754sp_min(sn) (ieee754sp_spcvals[IEEE754_SPCVAL_PMIN+(sn)])
  234. #define ieee754sp_mind(sn) (ieee754sp_spcvals[IEEE754_SPCVAL_PMIND+(sn)])
  235. #define ieee754sp_1e31() (ieee754sp_spcvals[IEEE754_SPCVAL_P1E31])
  236. #define ieee754sp_1e63() (ieee754sp_spcvals[IEEE754_SPCVAL_P1E63])
  237. /*
  238. * Indefinite integer value
  239. */
  240. static inline int ieee754si_indef(void)
  241. {
  242. return ieee754_csr.nan2008 ? 0 : INT_MAX;
  243. }
  244. static inline s64 ieee754di_indef(void)
  245. {
  246. return ieee754_csr.nan2008 ? 0 : S64_MAX;
  247. }
  248. /*
  249. * Overflow integer value
  250. */
  251. static inline int ieee754si_overflow(int xs)
  252. {
  253. return ieee754_csr.nan2008 && xs ? INT_MIN : INT_MAX;
  254. }
  255. static inline s64 ieee754di_overflow(int xs)
  256. {
  257. return ieee754_csr.nan2008 && xs ? S64_MIN : S64_MAX;
  258. }
  259. /* result types for xctx.rt */
  260. #define IEEE754_RT_SP 0
  261. #define IEEE754_RT_DP 1
  262. #define IEEE754_RT_XP 2
  263. #define IEEE754_RT_SI 3
  264. #define IEEE754_RT_DI 4
  265. /* compat */
  266. #define ieee754dp_fix(x) ieee754dp_tint(x)
  267. #define ieee754sp_fix(x) ieee754sp_tint(x)
  268. #endif /* __ARCH_MIPS_MATH_EMU_IEEE754_H */