Commit | Line | Data |
---|---|---|
1da177e4 LT |
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_OVERFLOW | |
55 | #define FP_EX_OVERFLOW 0 | |
56 | #endif | |
57 | #ifndef FP_EX_UNDERFLOW | |
58 | #define FP_EX_UNDERFLOW | |
59 | #endif | |
60 | #ifndef FP_EX_DIVZERO | |
61 | #define FP_EX_DIVZERO 0 | |
62 | #endif | |
63 | #ifndef FP_EX_INEXACT | |
64 | #define FP_EX_INEXACT 0 | |
65 | #endif | |
66 | #ifndef FP_EX_DENORM | |
67 | #define FP_EX_DENORM 0 | |
68 | #endif | |
69 | ||
70 | #ifdef _FP_DECL_EX | |
71 | #define FP_DECL_EX \ | |
72 | int _fex = 0; \ | |
73 | _FP_DECL_EX | |
74 | #else | |
75 | #define FP_DECL_EX int _fex = 0 | |
76 | #endif | |
77 | ||
78 | #ifndef FP_INIT_ROUNDMODE | |
79 | #define FP_INIT_ROUNDMODE do {} while (0) | |
80 | #endif | |
81 | ||
82 | #ifndef FP_HANDLE_EXCEPTIONS | |
83 | #define FP_HANDLE_EXCEPTIONS do {} while (0) | |
84 | #endif | |
85 | ||
86 | /* By default we never flush denormal input operands to signed zero. */ | |
87 | #ifndef FP_DENORM_ZERO | |
88 | #define FP_DENORM_ZERO 0 | |
89 | #endif | |
90 | ||
91 | #ifndef FP_INHIBIT_RESULTS | |
92 | /* By default we write the results always. | |
93 | * sfp-machine may override this and e.g. | |
94 | * check if some exceptions are unmasked | |
95 | * and inhibit it in such a case. | |
96 | */ | |
97 | #define FP_INHIBIT_RESULTS 0 | |
98 | #endif | |
99 | ||
40584961 | 100 | #ifndef FP_TRAPPING_EXCEPTIONS |
2a677896 | 101 | #define FP_TRAPPING_EXCEPTIONS 0 |
40584961 DM |
102 | #endif |
103 | ||
1da177e4 LT |
104 | #define FP_SET_EXCEPTION(ex) \ |
105 | _fex |= (ex) | |
106 | ||
107 | #define FP_UNSET_EXCEPTION(ex) \ | |
108 | _fex &= ~(ex) | |
109 | ||
40584961 DM |
110 | #define FP_CUR_EXCEPTIONS \ |
111 | (_fex) | |
112 | ||
1da177e4 LT |
113 | #define FP_CLEAR_EXCEPTIONS \ |
114 | _fex = 0 | |
115 | ||
116 | #define _FP_ROUND_NEAREST(wc, X) \ | |
117 | do { \ | |
118 | if ((_FP_FRAC_LOW_##wc(X) & 15) != _FP_WORK_ROUND) \ | |
119 | _FP_FRAC_ADDI_##wc(X, _FP_WORK_ROUND); \ | |
120 | } while (0) | |
121 | ||
122 | #define _FP_ROUND_ZERO(wc, X) 0 | |
123 | ||
124 | #define _FP_ROUND_PINF(wc, X) \ | |
125 | do { \ | |
126 | if (!X##_s && (_FP_FRAC_LOW_##wc(X) & 7)) \ | |
127 | _FP_FRAC_ADDI_##wc(X, _FP_WORK_LSB); \ | |
128 | } while (0) | |
129 | ||
130 | #define _FP_ROUND_MINF(wc, X) \ | |
131 | do { \ | |
132 | if (X##_s && (_FP_FRAC_LOW_##wc(X) & 7)) \ | |
133 | _FP_FRAC_ADDI_##wc(X, _FP_WORK_LSB); \ | |
134 | } while (0) | |
135 | ||
136 | #define _FP_ROUND(wc, X) \ | |
137 | do { \ | |
138 | if (_FP_FRAC_LOW_##wc(X) & 7) \ | |
139 | FP_SET_EXCEPTION(FP_EX_INEXACT); \ | |
140 | switch (FP_ROUNDMODE) \ | |
141 | { \ | |
142 | case FP_RND_NEAREST: \ | |
143 | _FP_ROUND_NEAREST(wc,X); \ | |
144 | break; \ | |
145 | case FP_RND_ZERO: \ | |
146 | _FP_ROUND_ZERO(wc,X); \ | |
147 | break; \ | |
148 | case FP_RND_PINF: \ | |
149 | _FP_ROUND_PINF(wc,X); \ | |
150 | break; \ | |
151 | case FP_RND_MINF: \ | |
152 | _FP_ROUND_MINF(wc,X); \ | |
153 | break; \ | |
154 | } \ | |
155 | } while (0) | |
156 | ||
157 | #define FP_CLS_NORMAL 0 | |
158 | #define FP_CLS_ZERO 1 | |
159 | #define FP_CLS_INF 2 | |
160 | #define FP_CLS_NAN 3 | |
161 | ||
162 | #define _FP_CLS_COMBINE(x,y) (((x) << 2) | (y)) | |
163 | ||
164 | #include <math-emu/op-1.h> | |
165 | #include <math-emu/op-2.h> | |
166 | #include <math-emu/op-4.h> | |
167 | #include <math-emu/op-8.h> | |
168 | #include <math-emu/op-common.h> | |
169 | ||
170 | /* Sigh. Silly things longlong.h needs. */ | |
171 | #define UWtype _FP_W_TYPE | |
172 | #define W_TYPE_SIZE _FP_W_TYPE_SIZE | |
173 | ||
174 | typedef int SItype __attribute__((mode(SI))); | |
175 | typedef int DItype __attribute__((mode(DI))); | |
176 | typedef unsigned int USItype __attribute__((mode(SI))); | |
177 | typedef unsigned int UDItype __attribute__((mode(DI))); | |
178 | #if _FP_W_TYPE_SIZE == 32 | |
179 | typedef unsigned int UHWtype __attribute__((mode(HI))); | |
180 | #elif _FP_W_TYPE_SIZE == 64 | |
181 | typedef USItype UHWtype; | |
182 | #endif | |
183 | ||
184 | #ifndef umul_ppmm | |
185 | #include <stdlib/longlong.h> | |
186 | #endif | |
187 | ||
188 | #endif /* __MATH_EMU_SOFT_FP_H__ */ |