Macintosh-Tools/macemu
1
0
Fork 0
mirror of https://github.com/kanjitalk755/macemu.git synced 2024-11-26 10:49:21 +00:00
Code Issues Projects Releases Wiki Activity

add generic roundf() from glibc (for IRIX/mips with older libm)

Browse Source
This commit is contained in:
gbeauche 2006-01-17 23:09:18 +00:00
parent 5ee0a6741b
commit 2e7f194ff2
2 changed files with 55 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

52
SheepShaver/src/kpx_cpu/src/mathlib/mathlib.cpp
View File

@ -51,6 +51,14 @@ do { \
(i) = gf_u.word; \
} while (0)
// Set a float from a 32 bit int
#define SET_FLOAT_WORD(d,i) \
do { \
ieee_float_shape_type sf_u; \
sf_u.word = (i); \
(d) = sf_u.value; \
} while (0)
/**
* Representation of an IEEE 754 double
@ -184,3 +192,47 @@ int mathlib_signbitl(long double x)
{
unimplemented("signbitl");
}
/**
* 7.12.9.6 The round functions
**/
float mathlib_roundf(float x)
{
int32 i0, j0;
static const float huge = 1.0e30;
GET_FLOAT_WORD (i0, x);
j0 = ((i0 >> 23) & 0xff) - 0x7f;
if (j0 < 23) {
if (j0 < 0) {
if (huge + x > 0.0F) {
i0 &= 0x80000000;
if (j0 == -1)
i0 |= 0x3f800000;
}
}
else {
u_int32_t i = 0x007fffff >> j0;
if ((i0 & i) == 0)
/* X is integral. */
return x;
if (huge + x > 0.0F) {
/* Raise inexact if x != 0. */
i0 += 0x00400000 >> j0;
i0 &= ~i;
}
}
}
else {
if (j0 == 0x80)
/* Inf or NaN. */
return x + x;
else
return x;
}
SET_FLOAT_WORD (x, i0);
return x;
}

3
SheepShaver/src/kpx_cpu/src/mathlib/mathlib.hpp
View File

@ -134,6 +134,9 @@ extern "C" float roundf(float x);
#ifdef HAVE_ROUND
extern "C" double round(double x);
#define roundf(x) (float)round(x)
#else
extern float mathlib_roundf(float);
#define roundf(x) mathlib_roundf(x)
#endif
#endif