1 | /* Software floating-point emulation. |
2 | Copyright (C) 1997,1998,1999 Free Software Foundation, Inc. |
3 | This file is part of the GNU C Library. |
4 | Contributed by Richard Henderson (rth@cygnus.com), |
5 | Jakub Jelinek (jj@ultra.linux.cz), |
6 | David S. Miller (davem@redhat.com) and |
7 | Peter Maydell (pmaydell@chiark.greenend.org.uk). |
8 | |
9 | The GNU C Library is free software; you can redistribute it and/or |
10 | modify it under the terms of the GNU Library General Public License as |
11 | published by the Free Software Foundation; either version 2 of the |
12 | License, or (at your option) any later version. |
13 | |
14 | The GNU C Library is distributed in the hope that it will be useful, |
15 | but WITHOUT ANY WARRANTY; without even the implied warranty of |
16 | MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU |
17 | Library General Public License for more details. |
18 | |
19 | You should have received a copy of the GNU Library General Public |
20 | License along with the GNU C Library; see the file COPYING.LIB. If |
21 | not, write to the Free Software Foundation, Inc., |
22 | 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA. */ |
23 | |
24 | #ifndef __MATH_EMU_SOFT_FP_H__ |
25 | #define __MATH_EMU_SOFT_FP_H__ |
26 | |
27 | #include <asm/sfp-machine.h> |
28 | |
29 | /* Allow sfp-machine to have its own byte order definitions. */ |
30 | #ifndef __BYTE_ORDER |
31 | #include <endian.h> |
32 | #endif |
33 | |
34 | #define _FP_WORKBITS 3 |
35 | #define _FP_WORK_LSB ((_FP_W_TYPE)1 << 3) |
36 | #define _FP_WORK_ROUND ((_FP_W_TYPE)1 << 2) |
37 | #define _FP_WORK_GUARD ((_FP_W_TYPE)1 << 1) |
38 | #define _FP_WORK_STICKY ((_FP_W_TYPE)1 << 0) |
39 | |
40 | #ifndef FP_RND_NEAREST |
41 | # define FP_RND_NEAREST 0 |
42 | # define FP_RND_ZERO 1 |
43 | # define FP_RND_PINF 2 |
44 | # define FP_RND_MINF 3 |
45 | #ifndef FP_ROUNDMODE |
46 | # define FP_ROUNDMODE FP_RND_NEAREST |
47 | #endif |
48 | #endif |
49 | |
50 | /* By default don't care about exceptions. */ |
51 | #ifndef FP_EX_INVALID |
52 | #define FP_EX_INVALID 0 |
53 | #endif |
54 | #ifndef FP_EX_INVALID_SNAN |
55 | #define FP_EX_INVALID_SNAN 0 |
56 | #endif |
57 | /* inf - inf */ |
58 | #ifndef FP_EX_INVALID_ISI |
59 | #define FP_EX_INVALID_ISI 0 |
60 | #endif |
61 | /* inf / inf */ |
62 | #ifndef FP_EX_INVALID_IDI |
63 | #define FP_EX_INVALID_IDI 0 |
64 | #endif |
65 | /* 0 / 0 */ |
66 | #ifndef FP_EX_INVALID_ZDZ |
67 | #define FP_EX_INVALID_ZDZ 0 |
68 | #endif |
69 | /* inf * 0 */ |
70 | #ifndef FP_EX_INVALID_IMZ |
71 | #define FP_EX_INVALID_IMZ 0 |
72 | #endif |
73 | #ifndef FP_EX_OVERFLOW |
74 | #define FP_EX_OVERFLOW 0 |
75 | #endif |
76 | #ifndef FP_EX_UNDERFLOW |
77 | #define FP_EX_UNDERFLOW |
78 | #endif |
79 | #ifndef FP_EX_DIVZERO |
80 | #define FP_EX_DIVZERO 0 |
81 | #endif |
82 | #ifndef FP_EX_INEXACT |
83 | #define FP_EX_INEXACT 0 |
84 | #endif |
85 | #ifndef FP_EX_DENORM |
86 | #define FP_EX_DENORM 0 |
87 | #endif |
88 | |
89 | #ifdef _FP_DECL_EX |
90 | #define FP_DECL_EX \ |
91 | int _fex = 0; \ |
92 | _FP_DECL_EX |
93 | #else |
94 | #define FP_DECL_EX int _fex = 0 |
95 | #endif |
96 | |
97 | #ifndef FP_INIT_ROUNDMODE |
98 | #define FP_INIT_ROUNDMODE do {} while (0) |
99 | #endif |
100 | |
101 | #ifndef FP_HANDLE_EXCEPTIONS |
102 | #define FP_HANDLE_EXCEPTIONS do {} while (0) |
103 | #endif |
104 | |
105 | /* By default we never flush denormal input operands to signed zero. */ |
106 | #ifndef FP_DENORM_ZERO |
107 | #define FP_DENORM_ZERO 0 |
108 | #endif |
109 | |
110 | #ifndef FP_INHIBIT_RESULTS |
111 | /* By default we write the results always. |
112 | * sfp-machine may override this and e.g. |
113 | * check if some exceptions are unmasked |
114 | * and inhibit it in such a case. |
115 | */ |
116 | #define FP_INHIBIT_RESULTS 0 |
117 | #endif |
118 | |
119 | #ifndef FP_TRAPPING_EXCEPTIONS |
120 | #define FP_TRAPPING_EXCEPTIONS 0 |
121 | #endif |
122 | |
123 | #define FP_SET_EXCEPTION(ex) \ |
124 | _fex |= (ex) |
125 | |
126 | #define FP_UNSET_EXCEPTION(ex) \ |
127 | _fex &= ~(ex) |
128 | |
129 | #define FP_CUR_EXCEPTIONS \ |
130 | (_fex) |
131 | |
132 | #define FP_CLEAR_EXCEPTIONS \ |
133 | _fex = 0 |
134 | |
135 | #define _FP_ROUND_NEAREST(wc, X) \ |
136 | do { \ |
137 | if ((_FP_FRAC_LOW_##wc(X) & 15) != _FP_WORK_ROUND) \ |
138 | _FP_FRAC_ADDI_##wc(X, _FP_WORK_ROUND); \ |
139 | } while (0) |
140 | |
141 | #define _FP_ROUND_ZERO(wc, X) (void)0 |
142 | |
143 | #define _FP_ROUND_PINF(wc, X) \ |
144 | do { \ |
145 | if (!X##_s && (_FP_FRAC_LOW_##wc(X) & 7)) \ |
146 | _FP_FRAC_ADDI_##wc(X, _FP_WORK_LSB); \ |
147 | } while (0) |
148 | |
149 | #define _FP_ROUND_MINF(wc, X) \ |
150 | do { \ |
151 | if (X##_s && (_FP_FRAC_LOW_##wc(X) & 7)) \ |
152 | _FP_FRAC_ADDI_##wc(X, _FP_WORK_LSB); \ |
153 | } while (0) |
154 | |
155 | #define _FP_ROUND(wc, X) \ |
156 | do { \ |
157 | if (_FP_FRAC_LOW_##wc(X) & 7) \ |
158 | FP_SET_EXCEPTION(FP_EX_INEXACT); \ |
159 | switch (FP_ROUNDMODE) \ |
160 | { \ |
161 | case FP_RND_NEAREST: \ |
162 | _FP_ROUND_NEAREST(wc,X); \ |
163 | break; \ |
164 | case FP_RND_ZERO: \ |
165 | _FP_ROUND_ZERO(wc,X); \ |
166 | break; \ |
167 | case FP_RND_PINF: \ |
168 | _FP_ROUND_PINF(wc,X); \ |
169 | break; \ |
170 | case FP_RND_MINF: \ |
171 | _FP_ROUND_MINF(wc,X); \ |
172 | break; \ |
173 | } \ |
174 | } while (0) |
175 | |
176 | #define FP_CLS_NORMAL 0 |
177 | #define FP_CLS_ZERO 1 |
178 | #define FP_CLS_INF 2 |
179 | #define FP_CLS_NAN 3 |
180 | |
181 | #define _FP_CLS_COMBINE(x,y) (((x) << 2) | (y)) |
182 | |
183 | #include <math-emu/op-1.h> |
184 | #include <math-emu/op-2.h> |
185 | #include <math-emu/op-4.h> |
186 | #include <math-emu/op-8.h> |
187 | #include <math-emu/op-common.h> |
188 | |
189 | /* Sigh. Silly things longlong.h needs. */ |
190 | #define UWtype _FP_W_TYPE |
191 | #define W_TYPE_SIZE _FP_W_TYPE_SIZE |
192 | |
193 | typedef int SItype __attribute__((mode(SI))); |
194 | typedef int DItype __attribute__((mode(DI))); |
195 | typedef unsigned int USItype __attribute__((mode(SI))); |
196 | typedef unsigned int UDItype __attribute__((mode(DI))); |
197 | #if _FP_W_TYPE_SIZE == 32 |
198 | typedef unsigned int UHWtype __attribute__((mode(HI))); |
199 | #elif _FP_W_TYPE_SIZE == 64 |
200 | typedef USItype UHWtype; |
201 | #endif |
202 | |
203 | #ifndef umul_ppmm |
204 | #include <stdlib/longlong.h> |
205 | #endif |
206 | |
207 | #endif /* __MATH_EMU_SOFT_FP_H__ */ |
208 | |