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