|
0
|
1 #define _FP_W_TYPE_SIZE 64
|
|
|
2 #define _FP_W_TYPE unsigned long
|
|
|
3 #define _FP_WS_TYPE signed long
|
|
|
4 #define _FP_I_TYPE long
|
|
|
5
|
|
|
6 typedef int TItype __attribute__ ((mode (TI)));
|
|
|
7 typedef unsigned int UTItype __attribute__ ((mode (TI)));
|
|
|
8
|
|
|
9 #define TI_BITS (__CHAR_BIT__ * (int)sizeof(TItype))
|
|
|
10
|
|
|
11 /* The type of the result of a floating point comparison. This must
|
|
|
12 match `__libgcc_cmp_return__' in GCC for the target. */
|
|
|
13 typedef int __gcc_CMPtype __attribute__ ((mode (__libgcc_cmp_return__)));
|
|
|
14 #define CMPtype __gcc_CMPtype
|
|
|
15
|
|
|
16 #define _FP_MUL_MEAT_Q(R,X,Y) \
|
|
|
17 _FP_MUL_MEAT_2_wide(_FP_WFRACBITS_Q,R,X,Y,umul_ppmm)
|
|
|
18
|
|
|
19 #define _FP_DIV_MEAT_Q(R,X,Y) _FP_DIV_MEAT_2_udiv(Q,R,X,Y)
|
|
|
20
|
|
|
21 #define _FP_NANFRAC_S _FP_QNANBIT_S
|
|
|
22 #define _FP_NANFRAC_D _FP_QNANBIT_D
|
|
|
23 #define _FP_NANFRAC_E _FP_QNANBIT_E, 0
|
|
|
24 #define _FP_NANFRAC_Q _FP_QNANBIT_Q, 0
|
|
|
25 #define _FP_NANSIGN_S 1
|
|
|
26 #define _FP_NANSIGN_D 1
|
|
|
27 #define _FP_NANSIGN_E 1
|
|
|
28 #define _FP_NANSIGN_Q 1
|
|
|
29
|
|
|
30 #define _FP_KEEPNANFRACP 1
|
|
|
31
|
|
|
32 /* Here is something Intel misdesigned: the specs don't define
|
|
|
33 the case where we have two NaNs with same mantissas, but
|
|
|
34 different sign. Different operations pick up different NaNs. */
|
|
|
35 #define _FP_CHOOSENAN(fs, wc, R, X, Y, OP) \
|
|
|
36 do { \
|
|
|
37 if (_FP_FRAC_GT_##wc(X, Y) \
|
|
|
38 || (_FP_FRAC_EQ_##wc(X,Y) && (OP == '+' || OP == '*'))) \
|
|
|
39 { \
|
|
|
40 R##_s = X##_s; \
|
|
|
41 _FP_FRAC_COPY_##wc(R,X); \
|
|
|
42 } \
|
|
|
43 else \
|
|
|
44 { \
|
|
|
45 R##_s = Y##_s; \
|
|
|
46 _FP_FRAC_COPY_##wc(R,Y); \
|
|
|
47 } \
|
|
|
48 R##_c = FP_CLS_NAN; \
|
|
|
49 } while (0)
|
|
|
50
|
|
|
51 #define FP_EX_INVALID 0x01
|
|
|
52 #define FP_EX_DENORM 0x02
|
|
|
53 #define FP_EX_DIVZERO 0x04
|
|
|
54 #define FP_EX_OVERFLOW 0x08
|
|
|
55 #define FP_EX_UNDERFLOW 0x10
|
|
|
56 #define FP_EX_INEXACT 0x20
|
|
|
57
|
|
|
58 struct fenv
|
|
|
59 {
|
|
|
60 unsigned short int __control_word;
|
|
|
61 unsigned short int __unused1;
|
|
|
62 unsigned short int __status_word;
|
|
|
63 unsigned short int __unused2;
|
|
|
64 unsigned short int __tags;
|
|
|
65 unsigned short int __unused3;
|
|
|
66 unsigned int __eip;
|
|
|
67 unsigned short int __cs_selector;
|
|
|
68 unsigned int __opcode:11;
|
|
|
69 unsigned int __unused4:5;
|
|
|
70 unsigned int __data_offset;
|
|
|
71 unsigned short int __data_selector;
|
|
|
72 unsigned short int __unused5;
|
|
|
73 };
|
|
|
74
|
|
|
75 #define FP_HANDLE_EXCEPTIONS \
|
|
|
76 do { \
|
|
|
77 if (_fex & FP_EX_INVALID) \
|
|
|
78 { \
|
|
|
79 float f = 0.0; \
|
|
|
80 __asm__ __volatile__ ("divss %0, %0 " : : "x" (f)); \
|
|
|
81 } \
|
|
|
82 if (_fex & FP_EX_DIVZERO) \
|
|
|
83 { \
|
|
|
84 float f = 1.0, g = 0.0; \
|
|
|
85 __asm__ __volatile__ ("divss %1, %0" : : "x" (f), "x" (g)); \
|
|
|
86 } \
|
|
|
87 if (_fex & FP_EX_OVERFLOW) \
|
|
|
88 { \
|
|
|
89 struct fenv temp; \
|
|
|
90 __asm__ __volatile__ ("fnstenv %0" : "=m" (temp)); \
|
|
|
91 temp.__status_word |= FP_EX_OVERFLOW; \
|
|
|
92 __asm__ __volatile__ ("fldenv %0" : : "m" (temp)); \
|
|
|
93 __asm__ __volatile__ ("fwait"); \
|
|
|
94 } \
|
|
|
95 if (_fex & FP_EX_UNDERFLOW) \
|
|
|
96 { \
|
|
|
97 struct fenv temp; \
|
|
|
98 __asm__ __volatile__ ("fnstenv %0" : "=m" (temp)); \
|
|
|
99 temp.__status_word |= FP_EX_UNDERFLOW; \
|
|
|
100 __asm__ __volatile__ ("fldenv %0" : : "m" (temp)); \
|
|
|
101 __asm__ __volatile__ ("fwait"); \
|
|
|
102 } \
|
|
|
103 if (_fex & FP_EX_INEXACT) \
|
|
|
104 { \
|
|
|
105 struct fenv temp; \
|
|
|
106 __asm__ __volatile__ ("fnstenv %0" : "=m" (temp)); \
|
|
|
107 temp.__status_word |= FP_EX_INEXACT; \
|
|
|
108 __asm__ __volatile__ ("fldenv %0" : : "m" (temp)); \
|
|
|
109 __asm__ __volatile__ ("fwait"); \
|
|
|
110 } \
|
|
|
111 } while (0)
|
|
|
112
|
|
|
113 #define FP_RND_NEAREST 0
|
|
|
114 #define FP_RND_ZERO 0xc00
|
|
|
115 #define FP_RND_PINF 0x800
|
|
|
116 #define FP_RND_MINF 0x400
|
|
|
117
|
|
|
118 #define _FP_DECL_EX \
|
|
|
119 unsigned short _fcw __attribute__ ((unused)) = FP_RND_NEAREST
|
|
|
120
|
|
|
121 #define FP_INIT_ROUNDMODE \
|
|
|
122 do { \
|
|
|
123 __asm__ ("fnstcw %0" : "=m" (_fcw)); \
|
|
|
124 } while (0)
|
|
|
125
|
|
|
126 #define FP_ROUNDMODE (_fcw & 0xc00)
|
|
|
127
|
|
|
128 #define __LITTLE_ENDIAN 1234
|
|
|
129 #define __BIG_ENDIAN 4321
|
|
|
130
|
|
|
131 #define __BYTE_ORDER __LITTLE_ENDIAN
|
|
|
132
|
|
|
133 /* Define ALIASNAME as a strong alias for NAME. */
|
|
|
134 #if defined __MACH__
|
|
|
135 /* Mach-O doesn't support aliasing. If these functions ever return
|
|
|
136 anything but CMPtype we need to revisit this... */
|
|
|
137 #define strong_alias(name, aliasname) \
|
|
|
138 CMPtype aliasname (TFtype a, TFtype b) { return name(a, b); }
|
|
|
139 #else
|
|
|
140 # define strong_alias(name, aliasname) _strong_alias(name, aliasname)
|
|
|
141 # define _strong_alias(name, aliasname) \
|
|
|
142 extern __typeof (name) aliasname __attribute__ ((alias (#name)));
|
|
|
143 #endif
|
