diff options
Diffstat (limited to 'sysdeps/ieee754/dbl-64/wordsize-64/s_ceil.c')
-rw-r--r-- | sysdeps/ieee754/dbl-64/wordsize-64/s_ceil.c | 16 |
1 files changed, 4 insertions, 12 deletions
diff --git a/sysdeps/ieee754/dbl-64/wordsize-64/s_ceil.c b/sysdeps/ieee754/dbl-64/wordsize-64/s_ceil.c index c687525d4b..b99829d2b0 100644 --- a/sysdeps/ieee754/dbl-64/wordsize-64/s_ceil.c +++ b/sysdeps/ieee754/dbl-64/wordsize-64/s_ceil.c @@ -15,14 +15,11 @@ * Return x rounded toward -inf to integral value * Method: * Bit twiddling. - * Exception: - * Inexact flag raised if x not equal to ceil(x). */ #include <math.h> #include <math_private.h> - -static const double huge = 1.0e300; +#include <libm-alias-double.h> double __ceil(double x) @@ -32,14 +29,13 @@ __ceil(double x) EXTRACT_WORDS64(i0,x); j0 = ((i0>>52)&0x7ff)-0x3ff; if(j0<=51) { - if(j0<0) { /* raise inexact if x != 0 */ - math_force_eval(huge+x);/* return 0*sign(x) if |x|<1 */ + if(j0<0) { + /* return 0*sign(x) if |x|<1 */ if(i0<0) {i0=INT64_C(0x8000000000000000);} else if(i0!=0) { i0=INT64_C(0x3ff0000000000000);} } else { i = INT64_C(0x000fffffffffffff)>>j0; if((i0&i)==0) return x; /* x is integral */ - math_force_eval(huge+x); /* raise inexact flag */ if(i0>0) i0 += UINT64_C(0x0010000000000000)>>j0; i0 &= (~i); } @@ -51,9 +47,5 @@ __ceil(double x) return x; } #ifndef __ceil -weak_alias (__ceil, ceil) -# ifdef NO_LONG_DOUBLE -strong_alias (__ceil, __ceill) -weak_alias (__ceil, ceill) -# endif +libm_alias_double (__ceil, ceil) #endif |