X-Git-Url: http://nsz.repo.hu/git/?a=blobdiff_plain;f=src%2Fmath%2Fmodfl.c;h=a47b1924f75d1125b712e3ccf5b0f39af27fd8f2;hb=89c5016cbacc8d2ea6c631ce12fae99cee48b62b;hp=f736bba4d95cc017bae7f65078b9f3a8c1ecd9a4;hpb=b9cd1d4d5073eba4c3afbe7666fc2d9c77b42c7c;p=musl

diff --git a/src/math/modfl.c b/src/math/modfl.c
index f736bba4..a47b1924 100644
--- a/src/math/modfl.c
+++ b/src/math/modfl.c
@@ -3,14 +3,17 @@
 #if LDBL_MANT_DIG == 53 && LDBL_MAX_EXP == 1024
 long double modfl(long double x, long double *iptr)
 {
-	return modf(x, (double *)iptr);
+	double d;
+	long double r;
+
+	r = modf(x, &d);
+	*iptr = d;
+	return r;
 }
 #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 modfl(long double x, long double *iptr)
 {
 	union ldshape u = {x};
@@ -35,7 +38,7 @@ long double modfl(long double x, long double *iptr)
 
 	/* raises spurious inexact */
 	absx = s ? -x : x;
-	y = absx + TOINT - TOINT - absx;
+	y = absx + toint - toint - absx;
 	if (y == 0) {
 		*iptr = x;
 		return s ? -0.0 : 0.0;