projects
/
musl
/ blobdiff
commit
grep
author
committer
pickaxe
?
search:
re
summary
|
shortlog
|
log
|
commit
|
commitdiff
|
tree
raw
|
inline
| side by side
math: use the rounding idiom consistently
[musl]
/
src
/
math
/
rintl.c
diff --git
a/src/math/rintl.c
b/src/math/rintl.c
index
2672507
..
374327d
100644
(file)
--- a/
src/math/rintl.c
+++ b/
src/math/rintl.c
@@
-6,11
+6,9
@@
long double rintl(long double x)
return rint(x);
}
#elif (LDBL_MANT_DIG == 64 || LDBL_MANT_DIG == 113) && LDBL_MAX_EXP == 16384
return rint(x);
}
#elif (LDBL_MANT_DIG == 64 || LDBL_MANT_DIG == 113) && LDBL_MAX_EXP == 16384
-#if LDBL_MANT_DIG == 64
-#define TOINT 0x1p63
-#elif LDBL_MANT_DIG == 113
-#define TOINT 0x1p112
-#endif
+
+static const long double toint = 1/LDBL_EPSILON;
+
long double rintl(long double x)
{
union ldshape u = {x};
long double rintl(long double x)
{
union ldshape u = {x};
@@
-21,9
+19,9
@@
long double rintl(long double x)
if (e >= 0x3fff+LDBL_MANT_DIG-1)
return x;
if (s)
if (e >= 0x3fff+LDBL_MANT_DIG-1)
return x;
if (s)
- y = x -
TOINT + TOINT
;
+ y = x -
toint + toint
;
else
else
- y = x +
TOINT - TOINT
;
+ y = x +
toint - toint
;
if (y == 0)
return 0*x;
return y;
if (y == 0)
return 0*x;
return y;