1 /* Native implementation of soft float functions. Only a single status
2 context is supported */
6 void set_float_rounding_mode(int val STATUS_PARAM
)
8 STATUS(float_rounding_mode
) = val
;
9 #if defined(_BSD) && !defined(__APPLE__)
11 #elif defined(__arm__)
19 void set_floatx80_rounding_precision(int val STATUS_PARAM
)
21 STATUS(floatx80_rounding_precision
) = val
;
26 #define lrint(d) ((int32_t)rint(d))
27 #define llrint(d) ((int64_t)rint(d))
30 #if defined(__powerpc__)
32 /* correct (but slow) PowerPC rint() (glibc version is incorrect) */
33 double qemu_rint(double x
)
35 double y
= 4503599627370496.0;
46 #define rint qemu_rint
49 /*----------------------------------------------------------------------------
50 | Software IEC/IEEE integer-to-floating-point conversion routines.
51 *----------------------------------------------------------------------------*/
52 float32
int32_to_float32(int v STATUS_PARAM
)
57 float64
int32_to_float64(int v STATUS_PARAM
)
63 floatx80
int32_to_floatx80(int v STATUS_PARAM
)
68 float32
int64_to_float32( int64_t v STATUS_PARAM
)
72 float64
int64_to_float64( int64_t v STATUS_PARAM
)
77 floatx80
int64_to_floatx80( int64_t v STATUS_PARAM
)
83 /*----------------------------------------------------------------------------
84 | Software IEC/IEEE single-precision conversion routines.
85 *----------------------------------------------------------------------------*/
86 int float32_to_int32( float32 a STATUS_PARAM
)
90 int float32_to_int32_round_to_zero( float32 a STATUS_PARAM
)
94 int64_t float32_to_int64( float32 a STATUS_PARAM
)
99 int64_t float32_to_int64_round_to_zero( float32 a STATUS_PARAM
)
104 float64
float32_to_float64( float32 a STATUS_PARAM
)
109 floatx80
float32_to_floatx80( float32 a STATUS_PARAM
)
115 /*----------------------------------------------------------------------------
116 | Software IEC/IEEE single-precision operations.
117 *----------------------------------------------------------------------------*/
118 float32
float32_round_to_int( float32 a STATUS_PARAM
)
123 float32
float32_rem( float32 a
, float32 b STATUS_PARAM
)
125 return remainderf(a
, b
);
128 float32
float32_sqrt( float32 a STATUS_PARAM
)
132 char float32_compare( float32 a
, float32 b STATUS_PARAM
)
144 char float32_compare_quiet( float32 a
, float32 b STATUS_PARAM
)
150 } else if (isgreater(a
, b
)) {
156 char float32_is_signaling_nan( float32 a1
)
162 return ( ( ( a
>>22 ) & 0x1FF ) == 0x1FE ) && ( a
& 0x003FFFFF );
165 /*----------------------------------------------------------------------------
166 | Software IEC/IEEE double-precision conversion routines.
167 *----------------------------------------------------------------------------*/
168 int float64_to_int32( float64 a STATUS_PARAM
)
172 int float64_to_int32_round_to_zero( float64 a STATUS_PARAM
)
176 int64_t float64_to_int64( float64 a STATUS_PARAM
)
180 int64_t float64_to_int64_round_to_zero( float64 a STATUS_PARAM
)
184 float32
float64_to_float32( float64 a STATUS_PARAM
)
189 floatx80
float64_to_floatx80( float64 a STATUS_PARAM
)
195 float128
float64_to_float128( float64 a STATUS_PARAM
)
201 /*----------------------------------------------------------------------------
202 | Software IEC/IEEE double-precision operations.
203 *----------------------------------------------------------------------------*/
204 float64
float64_round_to_int( float64 a STATUS_PARAM
)
207 switch(STATUS(float_rounding_mode
)) {
209 case float_round_nearest_even
:
210 asm("rndd %0, %1" : "=f" (a
) : "f"(a
));
212 case float_round_down
:
213 asm("rnddm %0, %1" : "=f" (a
) : "f"(a
));
216 asm("rnddp %0, %1" : "=f" (a
) : "f"(a
));
218 case float_round_to_zero
:
219 asm("rnddz %0, %1" : "=f" (a
) : "f"(a
));
227 float64
float64_rem( float64 a
, float64 b STATUS_PARAM
)
229 return remainder(a
, b
);
232 float64
float64_sqrt( float64 a STATUS_PARAM
)
236 char float64_compare( float64 a
, float64 b STATUS_PARAM
)
248 char float64_compare_quiet( float64 a
, float64 b STATUS_PARAM
)
254 } else if (isgreater(a
, b
)) {
260 char float64_is_signaling_nan( float64 a1
)
267 ( ( ( a
>>51 ) & 0xFFF ) == 0xFFE )
268 && ( a
& LIT64( 0x0007FFFFFFFFFFFF ) );
274 /*----------------------------------------------------------------------------
275 | Software IEC/IEEE extended double-precision conversion routines.
276 *----------------------------------------------------------------------------*/
277 int floatx80_to_int32( floatx80 a STATUS_PARAM
)
281 int floatx80_to_int32_round_to_zero( floatx80 a STATUS_PARAM
)
285 int64_t floatx80_to_int64( floatx80 a STATUS_PARAM
)
289 int64_t floatx80_to_int64_round_to_zero( floatx80 a STATUS_PARAM
)
293 float32
floatx80_to_float32( floatx80 a STATUS_PARAM
)
297 float64
floatx80_to_float64( floatx80 a STATUS_PARAM
)
302 /*----------------------------------------------------------------------------
303 | Software IEC/IEEE extended double-precision operations.
304 *----------------------------------------------------------------------------*/
305 floatx80
floatx80_round_to_int( floatx80 a STATUS_PARAM
)
309 floatx80
floatx80_rem( floatx80 a
, floatx80 b STATUS_PARAM
)
311 return remainderl(a
, b
);
313 floatx80
floatx80_sqrt( floatx80 a STATUS_PARAM
)
317 char floatx80_compare( floatx80 a
, floatx80 b STATUS_PARAM
)
329 char floatx80_compare_quiet( floatx80 a
, floatx80 b STATUS_PARAM
)
335 } else if (isgreater(a
, b
)) {
341 char floatx80_is_signaling_nan( floatx80 a1
)
345 return ( ( u
.i
.high
& 0x7FFF ) == 0x7FFF ) && (bits64
) ( u
.i
.low
<<1 );