[MATH-EMU]: Fix underflow exception reporting.
[linux-2.6.git] / include / math-emu / soft-fp.h
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
100 #ifndef FP_TRAPPING_EXCEPTIONS
101 #define FP_TRAPPING_EXCPETIONS 0
102 #endif
103
104 #define FP_SET_EXCEPTION(ex)                            \
105   _fex |= (ex)
106   
107 #define FP_UNSET_EXCEPTION(ex)                          \
108   _fex &= ~(ex)
109
110 #define FP_CUR_EXCEPTIONS                               \
111   (_fex)
112
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__ */