projects
/
musl
/ blobdiff
commit
grep
author
committer
pickaxe
?
search:
re
summary
|
shortlog
|
log
|
commit
|
commitdiff
|
tree
raw
|
inline
| side by side
optimize malloc0
[musl]
/
src
/
math
/
truncl.c
diff --git
a/src/math/truncl.c
b/src/math/truncl.c
index
3eedb08
..
f07b193
100644
(file)
--- a/
src/math/truncl.c
+++ b/
src/math/truncl.c
@@
-6,11
+6,9
@@
long double truncl(long double x)
return trunc(x);
}
#elif (LDBL_MANT_DIG == 64 || LDBL_MANT_DIG == 113) && LDBL_MAX_EXP == 16384
return trunc(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 truncl(long double x)
{
union ldshape u = {x};
long double truncl(long double x)
{
union ldshape u = {x};
@@
-27,7
+25,7
@@
long double truncl(long double x)
/* y = int(|x|) - |x|, where int(|x|) is an integer neighbor of |x| */
if (s)
x = -x;
/* y = int(|x|) - |x|, where int(|x|) is an integer neighbor of |x| */
if (s)
x = -x;
- y = x +
TOINT - TOINT
- x;
+ y = x +
toint - toint
- x;
if (y > 0)
y -= 1;
x += y;
if (y > 0)
y -= 1;
x += y;