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

* mingwex/math/lround.c: Rewrite. 

* mingwex/math/lroundf.c: Rewrite.
	* mingwex/math/lroundl.c: Rewrite.
	* mingwex/math/llround.c: Rewrite.
	* mingwex/math/llroundf.c: Rewrite.
	* mingwex/math/llroundl.c: Rewrite.
This commit is contained in:
Danny Smith 2004-04-22 04:38:24 +00:00
parent b7ede86cfe
commit 1a4b7623f4
7 changed files with 101 additions and 116 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

9
winsup/mingw/ChangeLog
View File

@ -1,3 +1,12 @@
2004-04-22 Danny Smith <dannysmith@users.sourceforge.net>
* mingwex/math/lround.c: Rewrite.
* mingwex/math/lroundf.c: Rewrite.
* mingwex/math/lroundl.c: Rewrite.
* mingwex/math/llround.c: Rewrite.
* mingwex/math/llroundf.c: Rewrite.
* mingwex/math/llroundl.c: Rewrite.
2004-04-20 Earnie Boyd <earnie@users.sf.net> 2004-04-20 Earnie Boyd <earnie@users.sf.net>
* CONTRIBUTORS: New file. * CONTRIBUTORS: New file.

35
winsup/mingw/mingwex/math/llround.c
View File

@ -1,24 +1,19 @@
#include <fenv.h>
#include <math.h> #include <math.h>
#include <limits.h>
#include <errno.h>
long long long long
llround (double x) { llround (double x)
long long retval; {
unsigned short saved_cw, _cw; /* Add +/- 0.5, then round towards zero. */
__asm__ ( double tmp = trunc (x + (x >= 0.0 ? 0.5 : -0.5));
"fnstcw %0;" : "=m" (saved_cw) if (!isfinite (tmp)
); /* save control word */ || tmp > (double)LONG_LONG_MAX
_cw = ~(FE_TONEAREST | FE_DOWNWARD | FE_UPWARD | FE_TOWARDZERO) || tmp < (double)LONG_LONG_MIN)
| (x > 0.0 ? FE_UPWARD : FE_DOWNWARD); /* round away from zero */ {
__asm__ ( errno = ERANGE;
"fldcw %0;" : : "m" (_cw) /* Undefined behaviour, so we could return anything. */
); /* load the rounding control */ /* return tmp > 0.0 ? LONG_LONG_MAX : LONG_LONG_MIN; */
__asm__ __volatile__ ( }
"fistpll %0" : "=m" (retval) : "t" (x) : "st" return (long long)tmp;
);
__asm__ (
"fldcw %0;" : : "m" (saved_cw)
); /* restore control word */
return retval;
} }

34
winsup/mingw/mingwex/math/llroundf.c
View File

@ -1,23 +1,19 @@
#include <fenv.h>
#include <math.h> #include <math.h>
#include <limits.h>
#include <errno.h>
long long long long
llroundf (float x) { llroundf (float x)
long long retval; {
unsigned short saved_cw, _cw; /* Add +/- 0.5, then round towards zero. */
__asm__ ( float tmp = truncf (x + (x >= 0.0F ? 0.5F : -0.5F));
"fnstcw %0;" : "=m" (saved_cw) if (!isfinite (tmp)
); /* save control word */ || tmp > (float)LONG_LONG_MAX
_cw = ~(FE_TONEAREST | FE_DOWNWARD | FE_UPWARD | FE_TOWARDZERO) || tmp < (float)LONG_LONG_MIN)
| (x > 0.0 ? FE_UPWARD : FE_DOWNWARD); /* round away from zero */ {
__asm__ ( errno = ERANGE;
"fldcw %0;" : : "m" (_cw) /* Undefined behaviour, so we could return anything. */
); /* load the rounding control */ /* return tmp > 0.0F ? LONG_LONG_MAX : LONG_LONG_MIN; */
__asm__ __volatile__ ( }
"fistpll %0" : "=m" (retval) : "t" (x) : "st" return (long long)tmp;
);
__asm__ (
"fldcw %0;" : : "m" (saved_cw)
); /* restore control word */
return retval;
} }

33
winsup/mingw/mingwex/math/llroundl.c
View File

@ -1,22 +1,19 @@
#include <fenv.h>
#include <math.h> #include <math.h>
#include <limits.h>
#include <errno.h>
long long long long
llroundl (long double x) { llroundl (long double x)
long long retval; {
unsigned short saved_cw, _cw; /* Add +/- 0.5, then round towards zero. */
__asm__ ( long double tmp = truncl (x + (x >= 0.0L ? 0.5L : -0.5L));
"fnstcw %0;" : "=m" (saved_cw) if (!isfinite (tmp)
); /* save control word */ || tmp > (long double)LONG_LONG_MAX
_cw = ~(FE_TONEAREST | FE_DOWNWARD | FE_UPWARD | FE_TOWARDZERO) || tmp < (long double)LONG_LONG_MIN)
| (x > 0.0 ? FE_UPWARD : FE_DOWNWARD); /* round away from zero */ {
__asm__ ( errno = ERANGE;
"fldcw %0;" : : "m" (_cw) /* Undefined behaviour, so we could return anything. */
); /* load the rounding control */ /* return tmp > 0.0L ? LONG_LONG_MAX : LONG_LONG_MIN; */
__asm__ __volatile__ ( }
"fistpll %0" : "=m" (retval) : "t" (x) : "st"); return (long long)tmp;
__asm__ (
"fldcw %0;" : : "m" (saved_cw)
); /* restore control word */
return retval;
} }

35
winsup/mingw/mingwex/math/lround.c
View File

@ -1,24 +1,19 @@
#include <fenv.h>
#include <math.h> #include <math.h>
#include <limits.h>
#include <errno.h>
long long
lround (double x) { lround (double x)
long retval; {
unsigned short saved_cw, _cw; /* Add +/- 0.5 then then round towards zero. */
__asm__ ( double tmp = trunc (x + (x >= 0.0 ? 0.5 : -0.5));
"fnstcw %0;" : "=m" (saved_cw) if (!isfinite (tmp)
); /* save control word */ || tmp > (double)LONG_MAX
_cw = ~(FE_TONEAREST | FE_DOWNWARD | FE_UPWARD | FE_TOWARDZERO) || tmp < (double)LONG_MIN)
| (x > 0.0 ? FE_UPWARD : FE_DOWNWARD); /* round away from zero */ {
__asm__ ( errno = ERANGE;
"fldcw %0;" : : "m" (_cw) /* Undefined behaviour, so we could return anything. */
); /* load the rounding control */ /* return tmp > 0.0 ? LONG_MAX : LONG_MIN; */
__asm__ __volatile__ ( }
"fistpl %0" : "=m" (retval) : "t" (x) : "st" return (long)tmp;
);
__asm__ (
"fldcw %0;" : : "m" (saved_cw)
); /* restore control word */
return retval;
} }

34
winsup/mingw/mingwex/math/lroundf.c
View File

@ -1,23 +1,19 @@
#include <fenv.h>
#include <math.h> #include <math.h>
#include <limits.h>
#include <errno.h>
long long
lroundf (float x) { lroundf (float x)
long retval; {
unsigned short saved_cw, _cw; /* Add +/- 0.5, then round towards zero. */
__asm__ ( float tmp = truncf (x + (x >= 0.0F ? 0.5F : -0.5F));
"fnstcw %0;" : "=m" (saved_cw) if (!isfinite (tmp)
); /* save control word */ || tmp > (float)LONG_MAX
_cw = ~(FE_TONEAREST | FE_DOWNWARD | FE_UPWARD | FE_TOWARDZERO) || tmp < (float)LONG_MIN)
| (x > 0.0 ? FE_UPWARD : FE_DOWNWARD); /* round away from zero */ {
__asm__ ( errno = ERANGE;
"fldcw %0;" : : "m" (_cw) /* Undefined behaviour, so we could return anything. */
); /* load the rounding control */ /* return tmp > 0.0F ? LONG_MAX : LONG_MIN; */
__asm__ __volatile__ ( }
"fistpl %0" : "=m" (retval) : "t" (x) : "st" return (long)tmp;
);
__asm__ (
"fldcw %0;" : : "m" (saved_cw)
); /* restore control word */
return retval;
} }

33
winsup/mingw/mingwex/math/lroundl.c
View File

@ -1,22 +1,19 @@
#include <fenv.h>
#include <math.h> #include <math.h>
#include <limits.h>
#include <errno.h>
long long
lroundl (long double x) { lroundl (long double x)
long retval; {
unsigned short saved_cw, _cw; /* Add +/- 0.5, then round towards zero. */
__asm__ ( long double tmp = truncl (x + (x >= 0.0L ? 0.5L : -0.5L));
"fnstcw %0;" : "=m" (saved_cw) if (!isfinite (tmp)
); /* save control word */ || tmp > (long double)LONG_MAX
_cw = ~(FE_TONEAREST | FE_DOWNWARD | FE_UPWARD | FE_TOWARDZERO) || tmp < (long double)LONG_MIN)
| (x > 0.0 ? FE_UPWARD : FE_DOWNWARD); /* round away from zero */ {
__asm__ ( errno = ERANGE;
"fldcw %0;" : : "m" (_cw) /* Undefined behaviour, so we could return anything. */
); /* load the rounding control */ /* return tmp > 0.0L ? LONG_MAX : LONG_MIN; */
__asm__ __volatile__ ( }
"fistpl %0" : "=m" (retval) : "t" (x) : "st"); return (long)tmp;
__asm__ (
"fldcw %0;" : : "m" (saved_cw)
); /* restore control word */
return retval;
} }