xref: /openbmc/linux/arch/mips/math-emu/ieee754int.h (revision 293d5b43)
1 /*
2  * IEEE754 floating point
3  * common internal header file
4  */
5 /*
6  * MIPS floating point support
7  * Copyright (C) 1994-2000 Algorithmics Ltd.
8  *
9  *  This program is free software; you can distribute it and/or modify it
10  *  under the terms of the GNU General Public License (Version 2) as
11  *  published by the Free Software Foundation.
12  *
13  *  This program is distributed in the hope it will be useful, but WITHOUT
14  *  ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
15  *  FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
16  *  for more details.
17  *
18  *  You should have received a copy of the GNU General Public License along
19  *  with this program; if not, write to the Free Software Foundation, Inc.,
20  *  51 Franklin St, Fifth Floor, Boston, MA  02110-1301  USA.
21  */
22 #ifndef __IEEE754INT_H
23 #define __IEEE754INT_H
24 
25 #include "ieee754.h"
26 
27 #define CLPAIR(x, y)	((x)*6+(y))
28 
29 static inline void ieee754_clearcx(void)
30 {
31 	ieee754_csr.cx = 0;
32 }
33 
34 static inline void ieee754_setcx(const unsigned int flags)
35 {
36 	ieee754_csr.cx |= flags;
37 	ieee754_csr.sx |= flags;
38 }
39 
40 static inline int ieee754_setandtestcx(const unsigned int x)
41 {
42 	ieee754_setcx(x);
43 
44 	return ieee754_csr.mx & x;
45 }
46 
47 static inline int ieee754_class_nan(int xc)
48 {
49 	return xc >= IEEE754_CLASS_SNAN;
50 }
51 
52 #define COMPXSP \
53 	unsigned xm; int xe; int xs __maybe_unused; int xc
54 
55 #define COMPYSP \
56 	unsigned ym; int ye; int ys; int yc
57 
58 #define COMPZSP \
59 	unsigned zm; int ze; int zs; int zc
60 
61 #define EXPLODESP(v, vc, vs, ve, vm)					\
62 {									\
63 	vs = SPSIGN(v);							\
64 	ve = SPBEXP(v);							\
65 	vm = SPMANT(v);							\
66 	if (ve == SP_EMAX+1+SP_EBIAS) {					\
67 		if (vm == 0)						\
68 			vc = IEEE754_CLASS_INF;				\
69 		else if (ieee754_csr.nan2008 ^ !(vm & SP_MBIT(SP_FBITS - 1))) \
70 			vc = IEEE754_CLASS_QNAN;			\
71 		else							\
72 			vc = IEEE754_CLASS_SNAN;			\
73 	} else if (ve == SP_EMIN-1+SP_EBIAS) {				\
74 		if (vm) {						\
75 			ve = SP_EMIN;					\
76 			vc = IEEE754_CLASS_DNORM;			\
77 		} else							\
78 			vc = IEEE754_CLASS_ZERO;			\
79 	} else {							\
80 		ve -= SP_EBIAS;						\
81 		vm |= SP_HIDDEN_BIT;					\
82 		vc = IEEE754_CLASS_NORM;				\
83 	}								\
84 }
85 #define EXPLODEXSP EXPLODESP(x, xc, xs, xe, xm)
86 #define EXPLODEYSP EXPLODESP(y, yc, ys, ye, ym)
87 #define EXPLODEZSP EXPLODESP(z, zc, zs, ze, zm)
88 
89 
90 #define COMPXDP \
91 	u64 xm; int xe; int xs __maybe_unused; int xc
92 
93 #define COMPYDP \
94 	u64 ym; int ye; int ys; int yc
95 
96 #define COMPZDP \
97 	u64 zm; int ze; int zs; int zc
98 
99 #define EXPLODEDP(v, vc, vs, ve, vm)					\
100 {									\
101 	vm = DPMANT(v);							\
102 	vs = DPSIGN(v);							\
103 	ve = DPBEXP(v);							\
104 	if (ve == DP_EMAX+1+DP_EBIAS) {					\
105 		if (vm == 0)						\
106 			vc = IEEE754_CLASS_INF;				\
107 		else if (ieee754_csr.nan2008 ^ !(vm & DP_MBIT(DP_FBITS - 1))) \
108 			vc = IEEE754_CLASS_QNAN;			\
109 		else							\
110 			vc = IEEE754_CLASS_SNAN;			\
111 	} else if (ve == DP_EMIN-1+DP_EBIAS) {				\
112 		if (vm) {						\
113 			ve = DP_EMIN;					\
114 			vc = IEEE754_CLASS_DNORM;			\
115 		} else							\
116 			vc = IEEE754_CLASS_ZERO;			\
117 	} else {							\
118 		ve -= DP_EBIAS;						\
119 		vm |= DP_HIDDEN_BIT;					\
120 		vc = IEEE754_CLASS_NORM;				\
121 	}								\
122 }
123 #define EXPLODEXDP EXPLODEDP(x, xc, xs, xe, xm)
124 #define EXPLODEYDP EXPLODEDP(y, yc, ys, ye, ym)
125 #define EXPLODEZDP EXPLODEDP(z, zc, zs, ze, zm)
126 
127 #define FLUSHDP(v, vc, vs, ve, vm)					\
128 	if (vc==IEEE754_CLASS_DNORM) {					\
129 		if (ieee754_csr.nod) {					\
130 			ieee754_setcx(IEEE754_INEXACT);			\
131 			vc = IEEE754_CLASS_ZERO;			\
132 			ve = DP_EMIN-1+DP_EBIAS;			\
133 			vm = 0;						\
134 			v = ieee754dp_zero(vs);				\
135 		}							\
136 	}
137 
138 #define FLUSHSP(v, vc, vs, ve, vm)					\
139 	if (vc==IEEE754_CLASS_DNORM) {					\
140 		if (ieee754_csr.nod) {					\
141 			ieee754_setcx(IEEE754_INEXACT);			\
142 			vc = IEEE754_CLASS_ZERO;			\
143 			ve = SP_EMIN-1+SP_EBIAS;			\
144 			vm = 0;						\
145 			v = ieee754sp_zero(vs);				\
146 		}							\
147 	}
148 
149 #define FLUSHXDP FLUSHDP(x, xc, xs, xe, xm)
150 #define FLUSHYDP FLUSHDP(y, yc, ys, ye, ym)
151 #define FLUSHZDP FLUSHDP(z, zc, zs, ze, zm)
152 #define FLUSHXSP FLUSHSP(x, xc, xs, xe, xm)
153 #define FLUSHYSP FLUSHSP(y, yc, ys, ye, ym)
154 #define FLUSHZSP FLUSHSP(z, zc, zs, ze, zm)
155 
156 #endif /* __IEEE754INT_H  */
157