xref: /freebsd/lib/libc/powerpc64/softfloat/softfloat.h (revision 2a63c3be158216222d89a073dcbd6a72ee4aab5a)
1*840b91ccSNathan Whitehorn /*	$NetBSD: softfloat.h,v 1.6 2002/05/12 13:12:46 bjh21 Exp $	*/
2*840b91ccSNathan Whitehorn 
3*840b91ccSNathan Whitehorn /* This is a derivative work. */
4*840b91ccSNathan Whitehorn 
5*840b91ccSNathan Whitehorn /*
6*840b91ccSNathan Whitehorn ===============================================================================
7*840b91ccSNathan Whitehorn 
8*840b91ccSNathan Whitehorn This C header file is part of the SoftFloat IEC/IEEE Floating-point
9*840b91ccSNathan Whitehorn Arithmetic Package, Release 2a.
10*840b91ccSNathan Whitehorn 
11*840b91ccSNathan Whitehorn Written by John R. Hauser.  This work was made possible in part by the
12*840b91ccSNathan Whitehorn International Computer Science Institute, located at Suite 600, 1947 Center
13*840b91ccSNathan Whitehorn Street, Berkeley, California 94704.  Funding was partially provided by the
14*840b91ccSNathan Whitehorn National Science Foundation under grant MIP-9311980.  The original version
15*840b91ccSNathan Whitehorn of this code was written as part of a project to build a fixed-point vector
16*840b91ccSNathan Whitehorn processor in collaboration with the University of California at Berkeley,
17*840b91ccSNathan Whitehorn overseen by Profs. Nelson Morgan and John Wawrzynek.  More information
18*840b91ccSNathan Whitehorn is available through the Web page `http://HTTP.CS.Berkeley.EDU/~jhauser/
19*840b91ccSNathan Whitehorn arithmetic/SoftFloat.html'.
20*840b91ccSNathan Whitehorn 
21*840b91ccSNathan Whitehorn THIS SOFTWARE IS DISTRIBUTED AS IS, FOR FREE.  Although reasonable effort
22*840b91ccSNathan Whitehorn has been made to avoid it, THIS SOFTWARE MAY CONTAIN FAULTS THAT WILL AT
23*840b91ccSNathan Whitehorn TIMES RESULT IN INCORRECT BEHAVIOR.  USE OF THIS SOFTWARE IS RESTRICTED TO
24*840b91ccSNathan Whitehorn PERSONS AND ORGANIZATIONS WHO CAN AND WILL TAKE FULL RESPONSIBILITY FOR ANY
25*840b91ccSNathan Whitehorn AND ALL LOSSES, COSTS, OR OTHER PROBLEMS ARISING FROM ITS USE.
26*840b91ccSNathan Whitehorn 
27*840b91ccSNathan Whitehorn Derivative works are acceptable, even for commercial purposes, so long as
28*840b91ccSNathan Whitehorn (1) they include prominent notice that the work is derivative, and (2) they
29*840b91ccSNathan Whitehorn include prominent notice akin to these four paragraphs for those parts of
30*840b91ccSNathan Whitehorn this code that are retained.
31*840b91ccSNathan Whitehorn 
32*840b91ccSNathan Whitehorn ===============================================================================
33*840b91ccSNathan Whitehorn */
34*840b91ccSNathan Whitehorn 
35*840b91ccSNathan Whitehorn /*
36*840b91ccSNathan Whitehorn -------------------------------------------------------------------------------
37*840b91ccSNathan Whitehorn The macro `FLOATX80' must be defined to enable the extended double-precision
38*840b91ccSNathan Whitehorn floating-point format `floatx80'.  If this macro is not defined, the
39*840b91ccSNathan Whitehorn `floatx80' type will not be defined, and none of the functions that either
40*840b91ccSNathan Whitehorn input or output the `floatx80' type will be defined.  The same applies to
41*840b91ccSNathan Whitehorn the `FLOAT128' macro and the quadruple-precision format `float128'.
42*840b91ccSNathan Whitehorn -------------------------------------------------------------------------------
43*840b91ccSNathan Whitehorn */
44*840b91ccSNathan Whitehorn /* #define FLOATX80 */
45*840b91ccSNathan Whitehorn /* #define FLOAT128 */
46*840b91ccSNathan Whitehorn 
47*840b91ccSNathan Whitehorn #include <machine/ieeefp.h>
48*840b91ccSNathan Whitehorn 
49*840b91ccSNathan Whitehorn /*
50*840b91ccSNathan Whitehorn -------------------------------------------------------------------------------
51*840b91ccSNathan Whitehorn Software IEC/IEEE floating-point types.
52*840b91ccSNathan Whitehorn -------------------------------------------------------------------------------
53*840b91ccSNathan Whitehorn */
54*840b91ccSNathan Whitehorn typedef unsigned int float32;
55*840b91ccSNathan Whitehorn typedef unsigned long long float64;
56*840b91ccSNathan Whitehorn #ifdef FLOATX80
57*840b91ccSNathan Whitehorn typedef struct {
58*840b91ccSNathan Whitehorn     unsigned short high;
59*840b91ccSNathan Whitehorn     unsigned long long low;
60*840b91ccSNathan Whitehorn } floatx80;
61*840b91ccSNathan Whitehorn #endif
62*840b91ccSNathan Whitehorn #ifdef FLOAT128
63*840b91ccSNathan Whitehorn typedef struct {
64*840b91ccSNathan Whitehorn     unsigned long long high, low;
65*840b91ccSNathan Whitehorn } float128;
66*840b91ccSNathan Whitehorn #endif
67*840b91ccSNathan Whitehorn 
68*840b91ccSNathan Whitehorn /*
69*840b91ccSNathan Whitehorn -------------------------------------------------------------------------------
70*840b91ccSNathan Whitehorn Software IEC/IEEE floating-point underflow tininess-detection mode.
71*840b91ccSNathan Whitehorn -------------------------------------------------------------------------------
72*840b91ccSNathan Whitehorn */
73*840b91ccSNathan Whitehorn #ifndef SOFTFLOAT_FOR_GCC
74*840b91ccSNathan Whitehorn extern int8 float_detect_tininess;
75*840b91ccSNathan Whitehorn #endif
76*840b91ccSNathan Whitehorn enum {
77*840b91ccSNathan Whitehorn     float_tininess_after_rounding  = 0,
78*840b91ccSNathan Whitehorn     float_tininess_before_rounding = 1
79*840b91ccSNathan Whitehorn };
80*840b91ccSNathan Whitehorn 
81*840b91ccSNathan Whitehorn /*
82*840b91ccSNathan Whitehorn -------------------------------------------------------------------------------
83*840b91ccSNathan Whitehorn Software IEC/IEEE floating-point rounding mode.
84*840b91ccSNathan Whitehorn -------------------------------------------------------------------------------
85*840b91ccSNathan Whitehorn */
86*840b91ccSNathan Whitehorn extern fp_rnd_t float_rounding_mode;
87*840b91ccSNathan Whitehorn enum {
88*840b91ccSNathan Whitehorn     float_round_nearest_even = FP_RN,
89*840b91ccSNathan Whitehorn     float_round_to_zero      = FP_RZ,
90*840b91ccSNathan Whitehorn     float_round_down         = FP_RM,
91*840b91ccSNathan Whitehorn     float_round_up           = FP_RP
92*840b91ccSNathan Whitehorn };
93*840b91ccSNathan Whitehorn 
94*840b91ccSNathan Whitehorn /*
95*840b91ccSNathan Whitehorn -------------------------------------------------------------------------------
96*840b91ccSNathan Whitehorn Software IEC/IEEE floating-point exception flags.
97*840b91ccSNathan Whitehorn -------------------------------------------------------------------------------
98*840b91ccSNathan Whitehorn */
99*840b91ccSNathan Whitehorn typedef fp_except_t fp_except;
100*840b91ccSNathan Whitehorn 
101*840b91ccSNathan Whitehorn extern fp_except float_exception_flags;
102*840b91ccSNathan Whitehorn extern fp_except float_exception_mask;
103*840b91ccSNathan Whitehorn enum {
104*840b91ccSNathan Whitehorn     float_flag_inexact   = FP_X_IMP,
105*840b91ccSNathan Whitehorn     float_flag_underflow = FP_X_UFL,
106*840b91ccSNathan Whitehorn     float_flag_overflow  = FP_X_OFL,
107*840b91ccSNathan Whitehorn     float_flag_divbyzero = FP_X_DZ,
108*840b91ccSNathan Whitehorn     float_flag_invalid   = FP_X_INV
109*840b91ccSNathan Whitehorn };
110*840b91ccSNathan Whitehorn 
111*840b91ccSNathan Whitehorn /*
112*840b91ccSNathan Whitehorn -------------------------------------------------------------------------------
113*840b91ccSNathan Whitehorn Routine to raise any or all of the software IEC/IEEE floating-point
114*840b91ccSNathan Whitehorn exception flags.
115*840b91ccSNathan Whitehorn -------------------------------------------------------------------------------
116*840b91ccSNathan Whitehorn */
117*840b91ccSNathan Whitehorn void float_raise( fp_except );
118*840b91ccSNathan Whitehorn 
119*840b91ccSNathan Whitehorn /*
120*840b91ccSNathan Whitehorn -------------------------------------------------------------------------------
121*840b91ccSNathan Whitehorn Software IEC/IEEE integer-to-floating-point conversion routines.
122*840b91ccSNathan Whitehorn -------------------------------------------------------------------------------
123*840b91ccSNathan Whitehorn */
124*840b91ccSNathan Whitehorn float32 int32_to_float32( int );
125*840b91ccSNathan Whitehorn float64 int32_to_float64( int );
126*840b91ccSNathan Whitehorn #ifdef FLOATX80
127*840b91ccSNathan Whitehorn floatx80 int32_to_floatx80( int );
128*840b91ccSNathan Whitehorn #endif
129*840b91ccSNathan Whitehorn #ifdef FLOAT128
130*840b91ccSNathan Whitehorn float128 int32_to_float128( int );
131*840b91ccSNathan Whitehorn #endif
132*840b91ccSNathan Whitehorn float32 int64_to_float32( long long );
133*840b91ccSNathan Whitehorn float64 int64_to_float64( long long );
134*840b91ccSNathan Whitehorn #ifdef FLOATX80
135*840b91ccSNathan Whitehorn floatx80 int64_to_floatx80( long long );
136*840b91ccSNathan Whitehorn #endif
137*840b91ccSNathan Whitehorn #ifdef FLOAT128
138*840b91ccSNathan Whitehorn float128 int64_to_float128( long long );
139*840b91ccSNathan Whitehorn #endif
140*840b91ccSNathan Whitehorn 
141*840b91ccSNathan Whitehorn /*
142*840b91ccSNathan Whitehorn -------------------------------------------------------------------------------
143*840b91ccSNathan Whitehorn Software IEC/IEEE single-precision conversion routines.
144*840b91ccSNathan Whitehorn -------------------------------------------------------------------------------
145*840b91ccSNathan Whitehorn */
146*840b91ccSNathan Whitehorn int float32_to_int32( float32 );
147*840b91ccSNathan Whitehorn int float32_to_int32_round_to_zero( float32 );
148*840b91ccSNathan Whitehorn unsigned int float32_to_uint32_round_to_zero( float32 );
149*840b91ccSNathan Whitehorn long long float32_to_int64( float32 );
150*840b91ccSNathan Whitehorn long long float32_to_int64_round_to_zero( float32 );
151*840b91ccSNathan Whitehorn float64 float32_to_float64( float32 );
152*840b91ccSNathan Whitehorn #ifdef FLOATX80
153*840b91ccSNathan Whitehorn floatx80 float32_to_floatx80( float32 );
154*840b91ccSNathan Whitehorn #endif
155*840b91ccSNathan Whitehorn #ifdef FLOAT128
156*840b91ccSNathan Whitehorn float128 float32_to_float128( float32 );
157*840b91ccSNathan Whitehorn #endif
158*840b91ccSNathan Whitehorn 
159*840b91ccSNathan Whitehorn /*
160*840b91ccSNathan Whitehorn -------------------------------------------------------------------------------
161*840b91ccSNathan Whitehorn Software IEC/IEEE single-precision operations.
162*840b91ccSNathan Whitehorn -------------------------------------------------------------------------------
163*840b91ccSNathan Whitehorn */
164*840b91ccSNathan Whitehorn float32 float32_round_to_int( float32 );
165*840b91ccSNathan Whitehorn float32 float32_add( float32, float32 );
166*840b91ccSNathan Whitehorn float32 float32_sub( float32, float32 );
167*840b91ccSNathan Whitehorn float32 float32_mul( float32, float32 );
168*840b91ccSNathan Whitehorn float32 float32_div( float32, float32 );
169*840b91ccSNathan Whitehorn float32 float32_rem( float32, float32 );
170*840b91ccSNathan Whitehorn float32 float32_sqrt( float32 );
171*840b91ccSNathan Whitehorn int float32_eq( float32, float32 );
172*840b91ccSNathan Whitehorn int float32_le( float32, float32 );
173*840b91ccSNathan Whitehorn int float32_lt( float32, float32 );
174*840b91ccSNathan Whitehorn int float32_eq_signaling( float32, float32 );
175*840b91ccSNathan Whitehorn int float32_le_quiet( float32, float32 );
176*840b91ccSNathan Whitehorn int float32_lt_quiet( float32, float32 );
177*840b91ccSNathan Whitehorn #ifndef SOFTFLOAT_FOR_GCC
178*840b91ccSNathan Whitehorn int float32_is_signaling_nan( float32 );
179*840b91ccSNathan Whitehorn #endif
180*840b91ccSNathan Whitehorn 
181*840b91ccSNathan Whitehorn /*
182*840b91ccSNathan Whitehorn -------------------------------------------------------------------------------
183*840b91ccSNathan Whitehorn Software IEC/IEEE double-precision conversion routines.
184*840b91ccSNathan Whitehorn -------------------------------------------------------------------------------
185*840b91ccSNathan Whitehorn */
186*840b91ccSNathan Whitehorn int float64_to_int32( float64 );
187*840b91ccSNathan Whitehorn int float64_to_int32_round_to_zero( float64 );
188*840b91ccSNathan Whitehorn unsigned int float64_to_uint32_round_to_zero( float64 );
189*840b91ccSNathan Whitehorn long long float64_to_int64( float64 );
190*840b91ccSNathan Whitehorn long long float64_to_int64_round_to_zero( float64 );
191*840b91ccSNathan Whitehorn float32 float64_to_float32( float64 );
192*840b91ccSNathan Whitehorn #ifdef FLOATX80
193*840b91ccSNathan Whitehorn floatx80 float64_to_floatx80( float64 );
194*840b91ccSNathan Whitehorn #endif
195*840b91ccSNathan Whitehorn #ifdef FLOAT128
196*840b91ccSNathan Whitehorn float128 float64_to_float128( float64 );
197*840b91ccSNathan Whitehorn #endif
198*840b91ccSNathan Whitehorn 
199*840b91ccSNathan Whitehorn /*
200*840b91ccSNathan Whitehorn -------------------------------------------------------------------------------
201*840b91ccSNathan Whitehorn Software IEC/IEEE double-precision operations.
202*840b91ccSNathan Whitehorn -------------------------------------------------------------------------------
203*840b91ccSNathan Whitehorn */
204*840b91ccSNathan Whitehorn float64 float64_round_to_int( float64 );
205*840b91ccSNathan Whitehorn float64 float64_add( float64, float64 );
206*840b91ccSNathan Whitehorn float64 float64_sub( float64, float64 );
207*840b91ccSNathan Whitehorn float64 float64_mul( float64, float64 );
208*840b91ccSNathan Whitehorn float64 float64_div( float64, float64 );
209*840b91ccSNathan Whitehorn float64 float64_rem( float64, float64 );
210*840b91ccSNathan Whitehorn float64 float64_sqrt( float64 );
211*840b91ccSNathan Whitehorn int float64_eq( float64, float64 );
212*840b91ccSNathan Whitehorn int float64_le( float64, float64 );
213*840b91ccSNathan Whitehorn int float64_lt( float64, float64 );
214*840b91ccSNathan Whitehorn int float64_eq_signaling( float64, float64 );
215*840b91ccSNathan Whitehorn int float64_le_quiet( float64, float64 );
216*840b91ccSNathan Whitehorn int float64_lt_quiet( float64, float64 );
217*840b91ccSNathan Whitehorn #ifndef SOFTFLOAT_FOR_GCC
218*840b91ccSNathan Whitehorn int float64_is_signaling_nan( float64 );
219*840b91ccSNathan Whitehorn #endif
220*840b91ccSNathan Whitehorn 
221*840b91ccSNathan Whitehorn #ifdef FLOATX80
222*840b91ccSNathan Whitehorn 
223*840b91ccSNathan Whitehorn /*
224*840b91ccSNathan Whitehorn -------------------------------------------------------------------------------
225*840b91ccSNathan Whitehorn Software IEC/IEEE extended double-precision conversion routines.
226*840b91ccSNathan Whitehorn -------------------------------------------------------------------------------
227*840b91ccSNathan Whitehorn */
228*840b91ccSNathan Whitehorn int floatx80_to_int32( floatx80 );
229*840b91ccSNathan Whitehorn int floatx80_to_int32_round_to_zero( floatx80 );
230*840b91ccSNathan Whitehorn long long floatx80_to_int64( floatx80 );
231*840b91ccSNathan Whitehorn long long floatx80_to_int64_round_to_zero( floatx80 );
232*840b91ccSNathan Whitehorn float32 floatx80_to_float32( floatx80 );
233*840b91ccSNathan Whitehorn float64 floatx80_to_float64( floatx80 );
234*840b91ccSNathan Whitehorn #ifdef FLOAT128
235*840b91ccSNathan Whitehorn float128 floatx80_to_float128( floatx80 );
236*840b91ccSNathan Whitehorn #endif
237*840b91ccSNathan Whitehorn 
238*840b91ccSNathan Whitehorn /*
239*840b91ccSNathan Whitehorn -------------------------------------------------------------------------------
240*840b91ccSNathan Whitehorn Software IEC/IEEE extended double-precision rounding precision.  Valid
241*840b91ccSNathan Whitehorn values are 32, 64, and 80.
242*840b91ccSNathan Whitehorn -------------------------------------------------------------------------------
243*840b91ccSNathan Whitehorn */
244*840b91ccSNathan Whitehorn extern int floatx80_rounding_precision;
245*840b91ccSNathan Whitehorn 
246*840b91ccSNathan Whitehorn /*
247*840b91ccSNathan Whitehorn -------------------------------------------------------------------------------
248*840b91ccSNathan Whitehorn Software IEC/IEEE extended double-precision operations.
249*840b91ccSNathan Whitehorn -------------------------------------------------------------------------------
250*840b91ccSNathan Whitehorn */
251*840b91ccSNathan Whitehorn floatx80 floatx80_round_to_int( floatx80 );
252*840b91ccSNathan Whitehorn floatx80 floatx80_add( floatx80, floatx80 );
253*840b91ccSNathan Whitehorn floatx80 floatx80_sub( floatx80, floatx80 );
254*840b91ccSNathan Whitehorn floatx80 floatx80_mul( floatx80, floatx80 );
255*840b91ccSNathan Whitehorn floatx80 floatx80_div( floatx80, floatx80 );
256*840b91ccSNathan Whitehorn floatx80 floatx80_rem( floatx80, floatx80 );
257*840b91ccSNathan Whitehorn floatx80 floatx80_sqrt( floatx80 );
258*840b91ccSNathan Whitehorn int floatx80_eq( floatx80, floatx80 );
259*840b91ccSNathan Whitehorn int floatx80_le( floatx80, floatx80 );
260*840b91ccSNathan Whitehorn int floatx80_lt( floatx80, floatx80 );
261*840b91ccSNathan Whitehorn int floatx80_eq_signaling( floatx80, floatx80 );
262*840b91ccSNathan Whitehorn int floatx80_le_quiet( floatx80, floatx80 );
263*840b91ccSNathan Whitehorn int floatx80_lt_quiet( floatx80, floatx80 );
264*840b91ccSNathan Whitehorn int floatx80_is_signaling_nan( floatx80 );
265*840b91ccSNathan Whitehorn 
266*840b91ccSNathan Whitehorn #endif
267*840b91ccSNathan Whitehorn 
268*840b91ccSNathan Whitehorn #ifdef FLOAT128
269*840b91ccSNathan Whitehorn 
270*840b91ccSNathan Whitehorn /*
271*840b91ccSNathan Whitehorn -------------------------------------------------------------------------------
272*840b91ccSNathan Whitehorn Software IEC/IEEE quadruple-precision conversion routines.
273*840b91ccSNathan Whitehorn -------------------------------------------------------------------------------
274*840b91ccSNathan Whitehorn */
275*840b91ccSNathan Whitehorn int float128_to_int32( float128 );
276*840b91ccSNathan Whitehorn int float128_to_int32_round_to_zero( float128 );
277*840b91ccSNathan Whitehorn long long float128_to_int64( float128 );
278*840b91ccSNathan Whitehorn long long float128_to_int64_round_to_zero( float128 );
279*840b91ccSNathan Whitehorn float32 float128_to_float32( float128 );
280*840b91ccSNathan Whitehorn float64 float128_to_float64( float128 );
281*840b91ccSNathan Whitehorn #ifdef FLOATX80
282*840b91ccSNathan Whitehorn floatx80 float128_to_floatx80( float128 );
283*840b91ccSNathan Whitehorn #endif
284*840b91ccSNathan Whitehorn 
285*840b91ccSNathan Whitehorn /*
286*840b91ccSNathan Whitehorn -------------------------------------------------------------------------------
287*840b91ccSNathan Whitehorn Software IEC/IEEE quadruple-precision operations.
288*840b91ccSNathan Whitehorn -------------------------------------------------------------------------------
289*840b91ccSNathan Whitehorn */
290*840b91ccSNathan Whitehorn float128 float128_round_to_int( float128 );
291*840b91ccSNathan Whitehorn float128 float128_add( float128, float128 );
292*840b91ccSNathan Whitehorn float128 float128_sub( float128, float128 );
293*840b91ccSNathan Whitehorn float128 float128_mul( float128, float128 );
294*840b91ccSNathan Whitehorn float128 float128_div( float128, float128 );
295*840b91ccSNathan Whitehorn float128 float128_rem( float128, float128 );
296*840b91ccSNathan Whitehorn float128 float128_sqrt( float128 );
297*840b91ccSNathan Whitehorn int float128_eq( float128, float128 );
298*840b91ccSNathan Whitehorn int float128_le( float128, float128 );
299*840b91ccSNathan Whitehorn int float128_lt( float128, float128 );
300*840b91ccSNathan Whitehorn int float128_eq_signaling( float128, float128 );
301*840b91ccSNathan Whitehorn int float128_le_quiet( float128, float128 );
302*840b91ccSNathan Whitehorn int float128_lt_quiet( float128, float128 );
303*840b91ccSNathan Whitehorn int float128_is_signaling_nan( float128 );
304*840b91ccSNathan Whitehorn 
305*840b91ccSNathan Whitehorn #endif
306*840b91ccSNathan Whitehorn 
307