drew
/
acadia-newlib
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

* mingwex/math/round.c: Rewrite. 

* mingwex/math/roundf.c: Rewrite.
	* mingwex/math/roundl.c: Rewrite.
This commit is contained in:
Danny Smith 2004-03-29 08:22:20 +00:00
parent 20e0dcfb3e
commit 521a0aa9a8
4 changed files with 21 additions and 79 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
winsup/mingw/ChangeLog
View File

@ -1,3 +1,9 @@
2004-03-29 Danny Smith <dannysmith@users.sourceforge.net>
* mingwex/math/round.c: Rewrite.
* mingwex/math/roundf.c: Rewrite.
* mingwex/math/roundl.c: Rewrite.
2004-03-25 Danny Smith <dannysmith@users.sourceforge.net>
* include/_mingw.h (__MSVCRT_VERSION__): Define default as 0x0600.

31
winsup/mingw/mingwex/math/round.c
View File

@ -1,29 +1,8 @@
#include <fenv.h>
#include <math.h>
double
round (double x) {
double retval;
unsigned short saved_cw, _cw;
__asm__ (
"fnstcw %0;"
: "=m" (saved_cw)
); /* save control word */
_cw = ~(FE_TONEAREST | FE_DOWNWARD | FE_UPWARD | FE_TOWARDZERO)
| (x > 0.0 ? FE_UPWARD : FE_DOWNWARD); /* round away from zero */
__asm__ (
"fldcw %0;"
:
: "m" (_cw)
); /* load the rounding control */
__asm__ (
"frndint;"
: "=t" (retval)
: "0" (x)
); /* do the rounding */
__asm__ (
"fldcw %0;"
:
: "m" (saved_cw)
); /* restore control word */
return retval;
round (double x)
{
/* Add +/- 0.5 then then round towards zero. */
return trunc ( x + (x >= 0.0 ? 0.5 : -0.5));
}

31
winsup/mingw/mingwex/math/roundf.c
View File

@ -1,29 +1,8 @@
#include <fenv.h>
#include <math.h>
float
roundf (float x) {
double retval;
unsigned short saved_cw, _cw;
__asm__ (
"fnstcw %0;"
: "=m" (saved_cw)
); /* save control word */
_cw = ~(FE_TONEAREST | FE_DOWNWARD | FE_UPWARD | FE_TOWARDZERO)
| (x > 0.0 ? FE_UPWARD : FE_DOWNWARD); /* round away from zero */
__asm__ (
"fldcw %0;"
:
: "m" (_cw)
); /* load the rounding control */
__asm__ (
"frndint;"
: "=t" (retval)
: "0" (x)
); /* do the rounding */
__asm__ (
"fldcw %0;"
:
: "m" (saved_cw)
); /* restore control word */
return retval;
roundf (float x)
{
/* Add +/- 0.5 then then round towards zero. */
return truncf ( x + (x >= 0.0F ? 0.5F : -0.5F));
}

32
winsup/mingw/mingwex/math/roundl.c
View File

@ -1,30 +1,8 @@
#include <fenv.h>
#include <math.h>
long double
roundl (long double x) {
long double retval;
unsigned short saved_cw, _cw;
__asm__ (
"fnstcw %0;"
: "=m" (saved_cw)
); /* save control word */
_cw = ~(FE_TONEAREST | FE_DOWNWARD | FE_UPWARD | FE_TOWARDZERO)
| (x > 0.0 ? FE_UPWARD : FE_DOWNWARD); /* round away from zero */
__asm__ (
"fldcw %0;"
:
: "m" (_cw)
); /* load the rounding control */
__asm__ (
"frndint;"
: "=t" (retval)
: "0" (x)
); /* do the rounding */
__asm__ (
"fldcw %0;"
:
: "m" (saved_cw)
); /* restore control word */
return retval;
roundl (long double x)
{
/* Add +/- 0.5 then then round towards zero. */
return truncl ( x + (x >= 0.0L ? 0.5L : -0.5L));
}