d/freebsd-dev
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Revert the last change, the conversion from long double to double can raise

Browse Source 
unwanted underflow exceptions.

Pointed out by:	das
This commit is contained in:
stefanf 2005-04-28 19:45:55 +00:00
parent e055b330dc
commit b663127cfb
3 changed files with 12 additions and 12 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
lib/msun/src/s_ceill.c
View File

@ -50,7 +50,7 @@ static char rcsid[] = "$FreeBSD$";
} while (0)
#endif
static const double huge = 1.0e300;
static const long double huge = 1.0e300;
long double
ceill(long double x)
@ -60,7 +60,7 @@ ceill(long double x)
if (e < MANH_SIZE - 1) {
if (e < 0) { /* raise inexact if x != 0 */
if (huge + (double)x > 0.0)
if (huge + x > 0.0)
if (u.bits.exp > 0 ||
(u.bits.manh | u.bits.manl) != 0)
u.e = u.bits.sign ? 0.0 : 1.0;
@ -76,7 +76,7 @@ ceill(long double x)
#endif
INC_MANH(u, 1llu << (MANH_SIZE - e - 1));
}
if (huge + (double)x > 0.0) { /* raise inexact flag */
if (huge + x > 0.0) { /* raise inexact flag */
u.bits.manh &= ~m;
u.bits.manl = 0;
}
@ -95,7 +95,7 @@ ceill(long double x)
INC_MANH(u, 1);
}
}
if (huge + (double)x > 0.0) /* raise inexact flag */
if (huge + x > 0.0) /* raise inexact flag */
u.bits.manl &= ~m;
}
return (u.e);

8
lib/msun/src/s_floorl.c
View File

@ -50,7 +50,7 @@ static char rcsid[] = "$FreeBSD$";
} while (0)
#endif
static const double huge = 1.0e300;
static const long double huge = 1.0e300;
long double
floorl(long double x)
@ -60,7 +60,7 @@ floorl(long double x)
if (e < MANH_SIZE - 1) {
if (e < 0) { /* raise inexact if x != 0 */
if (huge + (double)x > 0.0)
if (huge + x > 0.0)
if (u.bits.exp > 0 ||
(u.bits.manh | u.bits.manl) != 0)
u.e = u.bits.sign ? -1.0 : 0.0;
@ -76,7 +76,7 @@ floorl(long double x)
#endif
INC_MANH(u, 1llu << (MANH_SIZE - e - 1));
}
if (huge + (double)x > 0.0) { /* raise inexact flag */
if (huge + x > 0.0) { /* raise inexact flag */
u.bits.manh &= ~m;
u.bits.manl = 0;
}
@ -95,7 +95,7 @@ floorl(long double x)
INC_MANH(u, 1);
}
}
if (huge + (double)x > 0.0) /* raise inexact flag */
if (huge + x > 0.0) /* raise inexact flag */
u.bits.manl &= ~m;
}
return (u.e);

8
lib/msun/src/s_truncl.c
View File

@ -36,7 +36,7 @@ static char rcsid[] = "$FreeBSD$";
#define MANH_SIZE LDBL_MANH_SIZE
#endif
static const double huge = 1.0e300;
static const long double huge = 1.0e300;
long double
truncl(long double x)
@ -46,13 +46,13 @@ truncl(long double x)
if (e < MANH_SIZE - 1) {
if (e < 0) { /* raise inexact if x != 0 */
if (huge + (double)x > 0.0)
if (huge + x > 0.0)
u.e = 0.0;
} else {
uint64_t m = ((1llu << MANH_SIZE) - 1) >> (e + 1);
if (((u.bits.manh & m) | u.bits.manl) == 0)
return (x); /* x is integral */
if (huge + (double)x > 0.0) { /* raise inexact flag */
if (huge + x > 0.0) { /* raise inexact flag */
u.bits.manh &= ~m;
u.bits.manl = 0;
}
@ -61,7 +61,7 @@ truncl(long double x)
uint64_t m = (uint64_t)-1 >> (64 - LDBL_MANT_DIG + e + 1);
if ((u.bits.manl & m) == 0)
return (x); /* x is integral */
if (huge + (double)x > 0.0) /* raise inexact flag */
if (huge + x > 0.0) /* raise inexact flag */
u.bits.manl &= ~m;
}
return (u.e);