]>
Commit | Line | Data |
---|---|---|
d876f532 UD |
1 | /* Software floating-point emulation. |
2 | Definitions for IEEE Quad Precision. | |
fe0b1e85 | 3 | Copyright (C) 1997,1998,1999,2006 Free Software Foundation, Inc. |
d876f532 UD |
4 | This file is part of the GNU C Library. |
5 | Contributed by Richard Henderson (rth@cygnus.com), | |
6 | Jakub Jelinek (jj@ultra.linux.cz), | |
7 | David S. Miller (davem@redhat.com) and | |
8 | Peter Maydell (pmaydell@chiark.greenend.org.uk). | |
9 | ||
10 | The GNU C Library is free software; you can redistribute it and/or | |
41bdb6e2 AJ |
11 | modify it under the terms of the GNU Lesser General Public |
12 | License as published by the Free Software Foundation; either | |
13 | version 2.1 of the License, or (at your option) any later version. | |
d876f532 UD |
14 | |
15 | The GNU C Library is distributed in the hope that it will be useful, | |
16 | but WITHOUT ANY WARRANTY; without even the implied warranty of | |
17 | MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU | |
41bdb6e2 | 18 | Lesser General Public License for more details. |
d876f532 | 19 | |
41bdb6e2 AJ |
20 | You should have received a copy of the GNU Lesser General Public |
21 | License along with the GNU C Library; if not, write to the Free | |
a334319f UD |
22 | Software Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA |
23 | 02111-1307 USA. */ | |
d876f532 UD |
24 | |
25 | #if _FP_W_TYPE_SIZE < 32 | |
26 | #error "Here's a nickel, kid. Go buy yourself a real computer." | |
27 | #endif | |
28 | ||
29 | #if _FP_W_TYPE_SIZE < 64 | |
30 | #define _FP_FRACTBITS_Q (4*_FP_W_TYPE_SIZE) | |
31 | #else | |
32 | #define _FP_FRACTBITS_Q (2*_FP_W_TYPE_SIZE) | |
33 | #endif | |
34 | ||
35 | #define _FP_FRACBITS_Q 113 | |
36 | #define _FP_FRACXBITS_Q (_FP_FRACTBITS_Q - _FP_FRACBITS_Q) | |
37 | #define _FP_WFRACBITS_Q (_FP_WORKBITS + _FP_FRACBITS_Q) | |
38 | #define _FP_WFRACXBITS_Q (_FP_FRACTBITS_Q - _FP_WFRACBITS_Q) | |
39 | #define _FP_EXPBITS_Q 15 | |
40 | #define _FP_EXPBIAS_Q 16383 | |
41 | #define _FP_EXPMAX_Q 32767 | |
42 | ||
43 | #define _FP_QNANBIT_Q \ | |
44 | ((_FP_W_TYPE)1 << (_FP_FRACBITS_Q-2) % _FP_W_TYPE_SIZE) | |
fe0b1e85 RM |
45 | #define _FP_QNANBIT_SH_Q \ |
46 | ((_FP_W_TYPE)1 << (_FP_FRACBITS_Q-2+_FP_WORKBITS) % _FP_W_TYPE_SIZE) | |
d876f532 UD |
47 | #define _FP_IMPLBIT_Q \ |
48 | ((_FP_W_TYPE)1 << (_FP_FRACBITS_Q-1) % _FP_W_TYPE_SIZE) | |
fe0b1e85 RM |
49 | #define _FP_IMPLBIT_SH_Q \ |
50 | ((_FP_W_TYPE)1 << (_FP_FRACBITS_Q-1+_FP_WORKBITS) % _FP_W_TYPE_SIZE) | |
d876f532 UD |
51 | #define _FP_OVERFLOW_Q \ |
52 | ((_FP_W_TYPE)1 << (_FP_WFRACBITS_Q % _FP_W_TYPE_SIZE)) | |
53 | ||
fe0b1e85 RM |
54 | typedef float TFtype __attribute__((mode(TF))); |
55 | ||
d876f532 UD |
56 | #if _FP_W_TYPE_SIZE < 64 |
57 | ||
58 | union _FP_UNION_Q | |
59 | { | |
fe0b1e85 | 60 | TFtype flt; |
d876f532 UD |
61 | struct |
62 | { | |
63 | #if __BYTE_ORDER == __BIG_ENDIAN | |
64 | unsigned sign : 1; | |
65 | unsigned exp : _FP_EXPBITS_Q; | |
66 | unsigned long frac3 : _FP_FRACBITS_Q - (_FP_IMPLBIT_Q != 0)-(_FP_W_TYPE_SIZE * 3); | |
67 | unsigned long frac2 : _FP_W_TYPE_SIZE; | |
68 | unsigned long frac1 : _FP_W_TYPE_SIZE; | |
69 | unsigned long frac0 : _FP_W_TYPE_SIZE; | |
70 | #else | |
71 | unsigned long frac0 : _FP_W_TYPE_SIZE; | |
72 | unsigned long frac1 : _FP_W_TYPE_SIZE; | |
73 | unsigned long frac2 : _FP_W_TYPE_SIZE; | |
74 | unsigned long frac3 : _FP_FRACBITS_Q - (_FP_IMPLBIT_Q != 0)-(_FP_W_TYPE_SIZE * 3); | |
75 | unsigned exp : _FP_EXPBITS_Q; | |
76 | unsigned sign : 1; | |
77 | #endif /* not bigendian */ | |
78 | } bits __attribute__((packed)); | |
79 | }; | |
80 | ||
81 | ||
82 | #define FP_DECL_Q(X) _FP_DECL(4,X) | |
83 | #define FP_UNPACK_RAW_Q(X,val) _FP_UNPACK_RAW_4(Q,X,val) | |
84 | #define FP_UNPACK_RAW_QP(X,val) _FP_UNPACK_RAW_4_P(Q,X,val) | |
85 | #define FP_PACK_RAW_Q(val,X) _FP_PACK_RAW_4(Q,val,X) | |
86 | #define FP_PACK_RAW_QP(val,X) \ | |
87 | do { \ | |
88 | if (!FP_INHIBIT_RESULTS) \ | |
89 | _FP_PACK_RAW_4_P(Q,val,X); \ | |
90 | } while (0) | |
91 | ||
92 | #define FP_UNPACK_Q(X,val) \ | |
93 | do { \ | |
94 | _FP_UNPACK_RAW_4(Q,X,val); \ | |
95 | _FP_UNPACK_CANONICAL(Q,4,X); \ | |
96 | } while (0) | |
97 | ||
98 | #define FP_UNPACK_QP(X,val) \ | |
99 | do { \ | |
100 | _FP_UNPACK_RAW_4_P(Q,X,val); \ | |
101 | _FP_UNPACK_CANONICAL(Q,4,X); \ | |
102 | } while (0) | |
103 | ||
fe0b1e85 RM |
104 | #define FP_UNPACK_SEMIRAW_Q(X,val) \ |
105 | do { \ | |
106 | _FP_UNPACK_RAW_4(Q,X,val); \ | |
107 | _FP_UNPACK_SEMIRAW(Q,4,X); \ | |
108 | } while (0) | |
109 | ||
110 | #define FP_UNPACK_SEMIRAW_QP(X,val) \ | |
111 | do { \ | |
112 | _FP_UNPACK_RAW_4_P(Q,X,val); \ | |
113 | _FP_UNPACK_SEMIRAW(Q,4,X); \ | |
114 | } while (0) | |
115 | ||
d876f532 UD |
116 | #define FP_PACK_Q(val,X) \ |
117 | do { \ | |
118 | _FP_PACK_CANONICAL(Q,4,X); \ | |
119 | _FP_PACK_RAW_4(Q,val,X); \ | |
120 | } while (0) | |
121 | ||
122 | #define FP_PACK_QP(val,X) \ | |
123 | do { \ | |
124 | _FP_PACK_CANONICAL(Q,4,X); \ | |
125 | if (!FP_INHIBIT_RESULTS) \ | |
126 | _FP_PACK_RAW_4_P(Q,val,X); \ | |
127 | } while (0) | |
128 | ||
fe0b1e85 RM |
129 | #define FP_PACK_SEMIRAW_Q(val,X) \ |
130 | do { \ | |
131 | _FP_PACK_SEMIRAW(Q,4,X); \ | |
132 | _FP_PACK_RAW_4(Q,val,X); \ | |
133 | } while (0) | |
134 | ||
135 | #define FP_PACK_SEMIRAW_QP(val,X) \ | |
136 | do { \ | |
137 | _FP_PACK_SEMIRAW(Q,4,X); \ | |
138 | if (!FP_INHIBIT_RESULTS) \ | |
139 | _FP_PACK_RAW_4_P(Q,val,X); \ | |
140 | } while (0) | |
141 | ||
d876f532 UD |
142 | #define FP_ISSIGNAN_Q(X) _FP_ISSIGNAN(Q,4,X) |
143 | #define FP_NEG_Q(R,X) _FP_NEG(Q,4,R,X) | |
144 | #define FP_ADD_Q(R,X,Y) _FP_ADD(Q,4,R,X,Y) | |
145 | #define FP_SUB_Q(R,X,Y) _FP_SUB(Q,4,R,X,Y) | |
146 | #define FP_MUL_Q(R,X,Y) _FP_MUL(Q,4,R,X,Y) | |
147 | #define FP_DIV_Q(R,X,Y) _FP_DIV(Q,4,R,X,Y) | |
148 | #define FP_SQRT_Q(R,X) _FP_SQRT(Q,4,R,X) | |
149 | #define _FP_SQRT_MEAT_Q(R,S,T,X,Q) _FP_SQRT_MEAT_4(R,S,T,X,Q) | |
150 | ||
151 | #define FP_CMP_Q(r,X,Y,un) _FP_CMP(Q,4,r,X,Y,un) | |
152 | #define FP_CMP_EQ_Q(r,X,Y) _FP_CMP_EQ(Q,4,r,X,Y) | |
e7b8c7bc | 153 | #define FP_CMP_UNORD_Q(r,X,Y) _FP_CMP_UNORD(Q,4,r,X,Y) |
d876f532 UD |
154 | |
155 | #define FP_TO_INT_Q(r,X,rsz,rsg) _FP_TO_INT(Q,4,r,X,rsz,rsg) | |
156 | #define FP_FROM_INT_Q(X,r,rs,rt) _FP_FROM_INT(Q,4,X,r,rs,rt) | |
157 | ||
158 | #define _FP_FRAC_HIGH_Q(X) _FP_FRAC_HIGH_4(X) | |
159 | #define _FP_FRAC_HIGH_RAW_Q(X) _FP_FRAC_HIGH_4(X) | |
160 | ||
161 | #else /* not _FP_W_TYPE_SIZE < 64 */ | |
162 | union _FP_UNION_Q | |
163 | { | |
fe0b1e85 | 164 | TFtype flt /* __attribute__((mode(TF))) */ ; |
a598e3fc RH |
165 | struct { |
166 | _FP_W_TYPE a, b; | |
167 | } longs; | |
d876f532 UD |
168 | struct { |
169 | #if __BYTE_ORDER == __BIG_ENDIAN | |
170 | unsigned sign : 1; | |
171 | unsigned exp : _FP_EXPBITS_Q; | |
172 | unsigned long frac1 : _FP_FRACBITS_Q-(_FP_IMPLBIT_Q != 0)-_FP_W_TYPE_SIZE; | |
173 | unsigned long frac0 : _FP_W_TYPE_SIZE; | |
174 | #else | |
175 | unsigned long frac0 : _FP_W_TYPE_SIZE; | |
176 | unsigned long frac1 : _FP_FRACBITS_Q-(_FP_IMPLBIT_Q != 0)-_FP_W_TYPE_SIZE; | |
177 | unsigned exp : _FP_EXPBITS_Q; | |
178 | unsigned sign : 1; | |
179 | #endif | |
180 | } bits; | |
181 | }; | |
182 | ||
183 | #define FP_DECL_Q(X) _FP_DECL(2,X) | |
184 | #define FP_UNPACK_RAW_Q(X,val) _FP_UNPACK_RAW_2(Q,X,val) | |
185 | #define FP_UNPACK_RAW_QP(X,val) _FP_UNPACK_RAW_2_P(Q,X,val) | |
186 | #define FP_PACK_RAW_Q(val,X) _FP_PACK_RAW_2(Q,val,X) | |
187 | #define FP_PACK_RAW_QP(val,X) \ | |
188 | do { \ | |
189 | if (!FP_INHIBIT_RESULTS) \ | |
190 | _FP_PACK_RAW_2_P(Q,val,X); \ | |
191 | } while (0) | |
192 | ||
193 | #define FP_UNPACK_Q(X,val) \ | |
194 | do { \ | |
195 | _FP_UNPACK_RAW_2(Q,X,val); \ | |
196 | _FP_UNPACK_CANONICAL(Q,2,X); \ | |
197 | } while (0) | |
198 | ||
199 | #define FP_UNPACK_QP(X,val) \ | |
200 | do { \ | |
201 | _FP_UNPACK_RAW_2_P(Q,X,val); \ | |
202 | _FP_UNPACK_CANONICAL(Q,2,X); \ | |
203 | } while (0) | |
204 | ||
fe0b1e85 RM |
205 | #define FP_UNPACK_SEMIRAW_Q(X,val) \ |
206 | do { \ | |
207 | _FP_UNPACK_RAW_2(Q,X,val); \ | |
208 | _FP_UNPACK_SEMIRAW(Q,2,X); \ | |
209 | } while (0) | |
210 | ||
211 | #define FP_UNPACK_SEMIRAW_QP(X,val) \ | |
212 | do { \ | |
213 | _FP_UNPACK_RAW_2_P(Q,X,val); \ | |
214 | _FP_UNPACK_SEMIRAW(Q,2,X); \ | |
215 | } while (0) | |
216 | ||
d876f532 UD |
217 | #define FP_PACK_Q(val,X) \ |
218 | do { \ | |
219 | _FP_PACK_CANONICAL(Q,2,X); \ | |
220 | _FP_PACK_RAW_2(Q,val,X); \ | |
221 | } while (0) | |
222 | ||
223 | #define FP_PACK_QP(val,X) \ | |
224 | do { \ | |
225 | _FP_PACK_CANONICAL(Q,2,X); \ | |
226 | if (!FP_INHIBIT_RESULTS) \ | |
227 | _FP_PACK_RAW_2_P(Q,val,X); \ | |
228 | } while (0) | |
229 | ||
fe0b1e85 RM |
230 | #define FP_PACK_SEMIRAW_Q(val,X) \ |
231 | do { \ | |
232 | _FP_PACK_SEMIRAW(Q,2,X); \ | |
233 | _FP_PACK_RAW_2(Q,val,X); \ | |
234 | } while (0) | |
235 | ||
236 | #define FP_PACK_SEMIRAW_QP(val,X) \ | |
237 | do { \ | |
238 | _FP_PACK_SEMIRAW(Q,2,X); \ | |
239 | if (!FP_INHIBIT_RESULTS) \ | |
240 | _FP_PACK_RAW_2_P(Q,val,X); \ | |
241 | } while (0) | |
242 | ||
d876f532 UD |
243 | #define FP_ISSIGNAN_Q(X) _FP_ISSIGNAN(Q,2,X) |
244 | #define FP_NEG_Q(R,X) _FP_NEG(Q,2,R,X) | |
245 | #define FP_ADD_Q(R,X,Y) _FP_ADD(Q,2,R,X,Y) | |
246 | #define FP_SUB_Q(R,X,Y) _FP_SUB(Q,2,R,X,Y) | |
247 | #define FP_MUL_Q(R,X,Y) _FP_MUL(Q,2,R,X,Y) | |
248 | #define FP_DIV_Q(R,X,Y) _FP_DIV(Q,2,R,X,Y) | |
249 | #define FP_SQRT_Q(R,X) _FP_SQRT(Q,2,R,X) | |
250 | #define _FP_SQRT_MEAT_Q(R,S,T,X,Q) _FP_SQRT_MEAT_2(R,S,T,X,Q) | |
251 | ||
252 | #define FP_CMP_Q(r,X,Y,un) _FP_CMP(Q,2,r,X,Y,un) | |
253 | #define FP_CMP_EQ_Q(r,X,Y) _FP_CMP_EQ(Q,2,r,X,Y) | |
1e832e37 | 254 | #define FP_CMP_UNORD_Q(r,X,Y) _FP_CMP_UNORD(Q,2,r,X,Y) |
d876f532 UD |
255 | |
256 | #define FP_TO_INT_Q(r,X,rsz,rsg) _FP_TO_INT(Q,2,r,X,rsz,rsg) | |
257 | #define FP_FROM_INT_Q(X,r,rs,rt) _FP_FROM_INT(Q,2,X,r,rs,rt) | |
258 | ||
259 | #define _FP_FRAC_HIGH_Q(X) _FP_FRAC_HIGH_2(X) | |
260 | #define _FP_FRAC_HIGH_RAW_Q(X) _FP_FRAC_HIGH_2(X) | |
261 | ||
262 | #endif /* not _FP_W_TYPE_SIZE < 64 */ |