1 /* Header file for fp-bit.c.  */
2 /* Copyright (C) 2000-2019 Free Software Foundation, Inc.
3 
4 This file is part of GCC.
5 
6 GCC is free software; you can redistribute it and/or modify it under
7 the terms of the GNU General Public License as published by the Free
8 Software Foundation; either version 3, or (at your option) any later
9 version.
10 
11 GCC is distributed in the hope that it will be useful, but WITHOUT ANY
12 WARRANTY; without even the implied warranty of MERCHANTABILITY or
13 FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
14 for more details.
15 
16 Under Section 7 of GPL version 3, you are granted additional
17 permissions described in the GCC Runtime Library Exception, version
18 3.1, as published by the Free Software Foundation.
19 
20 You should have received a copy of the GNU General Public License and
21 a copy of the GCC Runtime Library Exception along with this program;
22 see the files COPYING3 and COPYING.RUNTIME respectively.  If not, see
23 <http://www.gnu.org/licenses/>.  */
24 
25 #ifndef GCC_FP_BIT_H
26 #define GCC_FP_BIT_H
27 
28 /* Defining FINE_GRAINED_LIBRARIES allows one to select which routines
29    from this file are compiled via additional -D options.
30 
31    This avoids the need to pull in the entire fp emulation library
32    when only a small number of functions are needed.
33 
34    If FINE_GRAINED_LIBRARIES is not defined, then compile every
35    suitable routine.  */
36 #ifndef FINE_GRAINED_LIBRARIES
37 #define L_pack_df
38 #define L_unpack_df
39 #define L_pack_sf
40 #define L_unpack_sf
41 #define L_addsub_sf
42 #define L_addsub_df
43 #define L_mul_sf
44 #define L_mul_df
45 #define L_div_sf
46 #define L_div_df
47 #define L_fpcmp_parts_sf
48 #define L_fpcmp_parts_df
49 #define L_compare_sf
50 #define L_compare_df
51 #define L_eq_sf
52 #define L_eq_df
53 #define L_ne_sf
54 #define L_ne_df
55 #define L_gt_sf
56 #define L_gt_df
57 #define L_ge_sf
58 #define L_ge_df
59 #define L_lt_sf
60 #define L_lt_df
61 #define L_le_sf
62 #define L_le_df
63 #define L_unord_sf
64 #define L_unord_df
65 #define L_usi_to_sf
66 #define L_usi_to_df
67 #define L_si_to_sf
68 #define L_si_to_df
69 #define L_sf_to_si
70 #define L_df_to_si
71 #define L_f_to_usi
72 #define L_df_to_usi
73 #define L_negate_sf
74 #define L_negate_df
75 #define L_make_sf
76 #define L_make_df
77 #define L_sf_to_df
78 #define L_df_to_sf
79 #ifdef FLOAT
80 #define L_thenan_sf
81 #else
82 #define L_thenan_df
83 #endif
84 #endif /* ! FINE_GRAINED_LIBRARIES */
85 
86 #if __LDBL_MANT_DIG__ == 113 || __LDBL_MANT_DIG__ == 106
87 # if defined(TFLOAT) || defined(L_sf_to_tf) || defined(L_df_to_tf)
88 #  define TMODES
89 # endif
90 #endif
91 
92 typedef float SFtype __attribute__ ((mode (SF)));
93 typedef float DFtype __attribute__ ((mode (DF)));
94 #ifdef TMODES
95 typedef float TFtype __attribute__ ((mode (TF)));
96 #endif
97 
98 typedef int HItype __attribute__ ((mode (HI)));
99 typedef int SItype __attribute__ ((mode (SI)));
100 typedef int DItype __attribute__ ((mode (DI)));
101 #ifdef TMODES
102 typedef int TItype __attribute__ ((mode (TI)));
103 #endif
104 
105 /* The type of the result of a floating point comparison.  This must
106    match `__libgcc_cmp_return__' in GCC for the target.  */
107 #ifndef CMPtype
108 typedef int CMPtype __attribute__ ((mode (__libgcc_cmp_return__)));
109 #endif
110 
111 typedef unsigned int UHItype __attribute__ ((mode (HI)));
112 typedef unsigned int USItype __attribute__ ((mode (SI)));
113 typedef unsigned int UDItype __attribute__ ((mode (DI)));
114 #ifdef TMODES
115 typedef unsigned int UTItype __attribute__ ((mode (TI)));
116 #endif
117 
118 #define MAX_USI_INT  (~(USItype)0)
119 #define MAX_SI_INT   ((SItype) (MAX_USI_INT >> 1))
120 #define BITS_PER_SI  (4 * __CHAR_BIT__)
121 #ifdef TMODES
122 #define MAX_UDI_INT  (~(UDItype)0)
123 #define MAX_DI_INT   ((DItype) (MAX_UDI_INT >> 1))
124 #define BITS_PER_DI  (8 * __CHAR_BIT__)
125 #endif
126 
127 #ifdef FLOAT_ONLY
128 #define NO_DI_MODE
129 #endif
130 
131 #if __BYTE_ORDER__ == __ORDER_LITTLE_ENDIAN__
132 #define FLOAT_BIT_ORDER_MISMATCH
133 #endif
134 
135 #if __BYTE_ORDER__ != __FLOAT_WORD_ORDER__
136 #define FLOAT_WORD_ORDER_MISMATCH
137 #endif
138 
139 #ifdef TFLOAT
140 # ifndef TMODES
141 #  error "TFLOAT requires long double to have 113 bits of mantissa"
142 # endif
143 
144 #	define PREFIXFPDP tp
145 #	define PREFIXSFDF tf
146 #	define NGARDS 10L /* Is this right? */
147 #	define GARDROUND 0x1ff
148 #	define GARDMASK  0x3ff
149 #	define GARDMSB   0x200
150 #	define FRAC_NBITS 128
151 
152 # if __LDBL_MANT_DIG__ == 113 /* IEEE quad */
153 #	define EXPBITS 15
154 #	define EXPBIAS 16383
155 #	define EXPMAX (0x7fff)
156 #	define QUIET_NAN ((TItype)0x8 << 108)
157 #	define FRACHIGH  ((TItype)0x8 << 124)
158 #	define FRACHIGH2 ((TItype)0xc << 124)
159 #	define FRACBITS 112
160 # endif
161 
162 # if __LDBL_MANT_DIG__ == 106 /* IBM extended (double+double) */
163 #	define EXPBITS 11
164 #	define EXPBIAS 1023
165 #	define EXPMAX (0x7ff)
166 #	define QUIET_NAN ((TItype)0x8 << (48 + 64))
167 #	define FRACHIGH  ((TItype)0x8 << 124)
168 #	define FRACHIGH2 ((TItype)0xc << 124)
169 #	define FRACBITS 105
170 #	define HALFFRACBITS 52
171 #	define HALFSHIFT 64
172 # endif
173 
174 #	define pack_d __pack_t
175 #	define unpack_d __unpack_t
176 #	define __fpcmp_parts __fpcmp_parts_t
177 	typedef UTItype fractype;
178 	typedef UDItype halffractype;
179 	typedef USItype qrtrfractype;
180 #define qrtrfractype qrtrfractype
181 	typedef TFtype FLO_type;
182 	typedef TItype intfrac;
183 #elif defined FLOAT
184 #	define NGARDS    7L
185 #	define GARDROUND 0x3f
186 #	define GARDMASK  0x7f
187 #	define GARDMSB   0x40
188 #	define EXPBITS 8
189 #	define EXPBIAS 127
190 #	define FRACBITS 23
191 #	define EXPMAX (0xff)
192 #	define QUIET_NAN 0x400000L
193 #	define FRAC_NBITS 32
194 #	define FRACHIGH  0x80000000L
195 #	define FRACHIGH2 0xc0000000L
196 #	define pack_d __pack_f
197 #	define unpack_d __unpack_f
198 #	define __fpcmp_parts __fpcmp_parts_f
199 	typedef USItype fractype;
200 	typedef UHItype halffractype;
201 	typedef SFtype FLO_type;
202 	typedef SItype intfrac;
203 
204 #else
205 #	define PREFIXFPDP dp
206 #	define PREFIXSFDF df
207 #	define NGARDS 8L
208 #	define GARDROUND 0x7f
209 #	define GARDMASK  0xff
210 #	define GARDMSB   0x80
211 #	define EXPBITS 11
212 #	define EXPBIAS 1023
213 #	define FRACBITS 52
214 #	define EXPMAX (0x7ff)
215 #	define QUIET_NAN 0x8000000000000LL
216 #	define FRAC_NBITS 64
217 #	define FRACHIGH  0x8000000000000000LL
218 #	define FRACHIGH2 0xc000000000000000LL
219 #	define pack_d __pack_d
220 #	define unpack_d __unpack_d
221 #	define __fpcmp_parts __fpcmp_parts_d
222 	typedef UDItype fractype;
223 	typedef USItype halffractype;
224 	typedef DFtype FLO_type;
225 	typedef DItype intfrac;
226 #endif /* FLOAT */
227 
228 #ifdef TFLOAT
229 #	define add 		__addtf3
230 #	define sub 		__subtf3
231 #	define multiply 	__multf3
232 #	define divide 		__divtf3
233 #	define compare 		__cmptf2
234 #	define _eq_f2 		__eqtf2
235 #	define _ne_f2 		__netf2
236 #	define _gt_f2 		__gttf2
237 #	define _ge_f2 		__getf2
238 #	define _lt_f2 		__lttf2
239 #	define _le_f2 		__letf2
240 #	define _unord_f2	__unordtf2
241 #	define usi_to_float 	__floatunsitf
242 #	define si_to_float 	__floatsitf
243 #	define float_to_si 	__fixtfsi
244 #	define float_to_usi 	__fixunstfsi
245 #	define negate 		__negtf2
246 #	define tf_to_sf		__trunctfsf2
247 #	define tf_to_df		__trunctfdf2
248 #elif defined FLOAT
249 #	define add 		__addsf3
250 #	define sub 		__subsf3
251 #	define multiply 	__mulsf3
252 #	define divide 		__divsf3
253 #	define compare 		__cmpsf2
254 #	define _eq_f2 		__eqsf2
255 #	define _ne_f2 		__nesf2
256 #	define _gt_f2 		__gtsf2
257 #	define _ge_f2 		__gesf2
258 #	define _lt_f2 		__ltsf2
259 #	define _le_f2 		__lesf2
260 #	define _unord_f2	__unordsf2
261 #	define usi_to_float 	__floatunsisf
262 #	define si_to_float 	__floatsisf
263 #	define float_to_si 	__fixsfsi
264 #	define float_to_usi 	__fixunssfsi
265 #	define negate 		__negsf2
266 #	define sf_to_df		__extendsfdf2
267 #	define sf_to_tf		__extendsftf2
268 #else
269 #	define add 		__adddf3
270 #	define sub 		__subdf3
271 #	define multiply 	__muldf3
272 #	define divide 		__divdf3
273 #	define compare 		__cmpdf2
274 #	define _eq_f2 		__eqdf2
275 #	define _ne_f2 		__nedf2
276 #	define _gt_f2 		__gtdf2
277 #	define _ge_f2 		__gedf2
278 #	define _lt_f2 		__ltdf2
279 #	define _le_f2 		__ledf2
280 #	define _unord_f2	__unorddf2
281 #	define usi_to_float 	__floatunsidf
282 #	define si_to_float 	__floatsidf
283 #	define float_to_si 	__fixdfsi
284 #	define float_to_usi 	__fixunsdfsi
285 #	define negate 		__negdf2
286 #	define df_to_sf		__truncdfsf2
287 #	define df_to_tf		__extenddftf2
288 #endif /* FLOAT */
289 
290 #ifndef INLINE
291 #define INLINE __inline__
292 #endif
293 
294 /* Preserve the sticky-bit when shifting fractions to the right.  */
295 #define LSHIFT(a, s) { a = (a >> s) | !!(a & (((fractype) 1 << s) - 1)); }
296 
297 /* numeric parameters */
298 /* F_D_BITOFF is the number of bits offset between the MSB of the mantissa
299    of a float and of a double. Assumes there are only two float types.
300    (double::FRAC_BITS+double::NGARDS-(float::FRAC_BITS+float::NGARDS))
301  */
302 #define F_D_BITOFF (52+8-(23+7))
303 
304 #ifdef TMODES
305 # define F_T_BITOFF (__LDBL_MANT_DIG__-1+10-(23+7))
306 # define D_T_BITOFF (__LDBL_MANT_DIG__-1+10-(52+8))
307 #endif
308 
309 
310 #define NORMAL_EXPMIN (-(EXPBIAS)+1)
311 #define IMPLICIT_1 ((fractype)1<<(FRACBITS+NGARDS))
312 #define IMPLICIT_2 ((fractype)1<<(FRACBITS+1+NGARDS))
313 
314 /* common types */
315 
316 typedef enum
317 {
318   CLASS_SNAN,
319   CLASS_QNAN,
320   CLASS_ZERO,
321   CLASS_NUMBER,
322   CLASS_INFINITY
323 } fp_class_type;
324 
325 typedef struct
326 {
327 #ifdef SMALL_MACHINE
328   char class;
329   unsigned char sign;
330   short normal_exp;
331 #else
332   fp_class_type class;
333   unsigned int sign;
334   int normal_exp;
335 #endif
336 
337   union
338     {
339       fractype ll;
340       halffractype l[2];
341     } fraction;
342 } fp_number_type;
343 
344 typedef union
345 {
346   FLO_type value;
347   fractype value_raw;
348 
349 #ifndef FLOAT
350 # ifdef qrtrfractype
351   qrtrfractype qwords[4];
352 # else
353   halffractype words[2];
354 # endif
355 #endif
356 
357 #ifdef FLOAT_BIT_ORDER_MISMATCH
358   struct
359     {
360       fractype fraction:FRACBITS __attribute__ ((packed));
361       unsigned int exp:EXPBITS __attribute__ ((packed));
362       unsigned int sign:1 __attribute__ ((packed));
363     }
364   bits;
365 #endif
366 
367 #ifdef _DEBUG_BITFLOAT
368   struct
369     {
370       unsigned int sign:1 __attribute__ ((packed));
371       unsigned int exp:EXPBITS __attribute__ ((packed));
372       fractype fraction:FRACBITS __attribute__ ((packed));
373     }
374   bits_big_endian;
375 
376   struct
377     {
378       fractype fraction:FRACBITS __attribute__ ((packed));
379       unsigned int exp:EXPBITS __attribute__ ((packed));
380       unsigned int sign:1 __attribute__ ((packed));
381     }
382   bits_little_endian;
383 #endif
384 }
385 FLO_union_type;
386 
387 /* Prototypes.  */
388 
389 #if defined(L_pack_df) || defined(L_pack_sf) || defined(L_pack_tf)
390 extern FLO_type pack_d (const fp_number_type *);
391 #endif
392 
393 extern void unpack_d (FLO_union_type *, fp_number_type *);
394 
395 #if defined(L_addsub_sf) || defined(L_addsub_df) || defined(L_addsub_tf)
396 extern FLO_type add (FLO_type, FLO_type);
397 extern FLO_type sub (FLO_type, FLO_type);
398 #endif
399 
400 #if defined(L_mul_sf) || defined(L_mul_df) || defined(L_mul_tf)
401 extern FLO_type multiply (FLO_type, FLO_type);
402 #endif
403 
404 #if defined(L_div_sf) || defined(L_div_df) || defined(L_div_tf)
405 extern FLO_type divide (FLO_type, FLO_type);
406 #endif
407 
408 extern int __fpcmp_parts (fp_number_type *, fp_number_type *);
409 
410 #if defined(L_compare_sf) || defined(L_compare_df) || defined(L_compare_tf)
411 extern CMPtype compare (FLO_type, FLO_type);
412 #endif
413 
414 #if defined(L_eq_sf) || defined(L_eq_df) || defined(L_eq_tf)
415 extern CMPtype _eq_f2 (FLO_type, FLO_type);
416 #endif
417 
418 #if defined(L_ne_sf) || defined(L_ne_df) || defined(L_ne_tf)
419 extern CMPtype _ne_f2 (FLO_type, FLO_type);
420 #endif
421 
422 #if defined(L_gt_sf) || defined(L_gt_df) || defined(L_gt_tf)
423 extern CMPtype _gt_f2 (FLO_type, FLO_type);
424 #endif
425 
426 #if defined(L_ge_sf) || defined(L_ge_df) || defined(L_ge_tf)
427 extern CMPtype _ge_f2 (FLO_type, FLO_type);
428 #endif
429 
430 #if defined(L_lt_sf) || defined(L_lt_df) || defined(L_lt_tf)
431 extern CMPtype _lt_f2 (FLO_type, FLO_type);
432 #endif
433 
434 #if defined(L_le_sf) || defined(L_le_df) || defined(L_le_tf)
435 extern CMPtype _le_f2 (FLO_type, FLO_type);
436 #endif
437 
438 #if defined(L_unord_sf) || defined(L_unord_df) || defined(L_unord_tf)
439 extern CMPtype _unord_f2 (FLO_type, FLO_type);
440 #endif
441 
442 #if defined(L_si_to_sf) || defined(L_si_to_df) || defined(L_si_to_tf)
443 extern FLO_type si_to_float (SItype);
444 #endif
445 
446 #if defined(L_sf_to_si) || defined(L_df_to_si) || defined(L_tf_to_si)
447 extern SItype float_to_si (FLO_type);
448 #endif
449 
450 #if defined(L_tf_to_usi)
451 extern USItype float_to_usi (FLO_type);
452 #endif
453 
454 #if defined(L_usi_to_sf) || defined(L_usi_to_df) || defined(L_usi_to_tf)
455 extern FLO_type usi_to_float (USItype);
456 #endif
457 
458 #if defined(L_negate_sf) || defined(L_negate_df) || defined(L_negate_tf)
459 extern FLO_type negate (FLO_type);
460 #endif
461 
462 #ifdef FLOAT
463 #if defined(L_make_sf)
464 extern SFtype __make_fp (fp_class_type, unsigned int, int, USItype);
465 #endif
466 #ifndef FLOAT_ONLY
467 extern DFtype __make_dp (fp_class_type, unsigned int, int, UDItype);
468 #if defined(L_sf_to_df)
469 extern DFtype sf_to_df (SFtype);
470 #endif
471 #if defined(L_sf_to_tf) && defined(TMODES)
472 extern TFtype sf_to_tf (SFtype);
473 #endif
474 #endif /* ! FLOAT_ONLY */
475 #endif /* FLOAT */
476 
477 #ifndef FLOAT
478 extern SFtype __make_fp (fp_class_type, unsigned int, int, USItype);
479 #if defined(L_make_df)
480 extern DFtype __make_dp (fp_class_type, unsigned int, int, UDItype);
481 #endif
482 #if defined(L_df_to_sf)
483 extern SFtype df_to_sf (DFtype);
484 #endif
485 #if defined(L_df_to_tf) && defined(TMODES)
486 extern TFtype df_to_tf (DFtype);
487 #endif
488 #endif /* ! FLOAT */
489 
490 #ifdef TMODES
491 extern DFtype __make_dp (fp_class_type, unsigned int, int, UDItype);
492 extern TFtype __make_tp (fp_class_type, unsigned int, int, UTItype);
493 #ifdef TFLOAT
494 #if defined(L_tf_to_sf)
495 extern SFtype tf_to_sf (TFtype);
496 #endif
497 #if defined(L_tf_to_df)
498 extern DFtype tf_to_df (TFtype);
499 #endif
500 #if defined(L_di_to_tf)
501 extern TFtype di_to_df (DItype);
502 #endif
503 #endif /* TFLOAT */
504 #endif /* TMODES */
505 
506 #endif /* ! GCC_FP_BIT_H */
507