123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112 |
- #if ! defined USE_LONG_DOUBLE
- # include <config.h>
- #endif
- #include <math.h>
- #include <float.h>
- #undef MIN
- #ifdef USE_LONG_DOUBLE
- # define FUNC ceill
- # define DOUBLE long double
- # define MANT_DIG LDBL_MANT_DIG
- # define MIN LDBL_MIN
- # define L_(literal) literal##L
- #elif ! defined USE_FLOAT
- # define FUNC ceil
- # define DOUBLE double
- # define MANT_DIG DBL_MANT_DIG
- # define MIN DBL_MIN
- # define L_(literal) literal
- #else
- # define FUNC ceilf
- # define DOUBLE float
- # define MANT_DIG FLT_MANT_DIG
- # define MIN FLT_MIN
- # define L_(literal) literal##f
- #endif
- #if defined __hpux || defined __sgi || defined __ICC
- # define MINUS_ZERO (-MIN * MIN)
- #else
- # define MINUS_ZERO L_(-0.0)
- #endif
- static const DOUBLE TWO_MANT_DIG =
-
- (DOUBLE) (1U << ((MANT_DIG - 1) / 5))
- * (DOUBLE) (1U << ((MANT_DIG - 1 + 1) / 5))
- * (DOUBLE) (1U << ((MANT_DIG - 1 + 2) / 5))
- * (DOUBLE) (1U << ((MANT_DIG - 1 + 3) / 5))
- * (DOUBLE) (1U << ((MANT_DIG - 1 + 4) / 5));
- DOUBLE
- FUNC (DOUBLE x)
- {
-
- volatile DOUBLE y = x;
- volatile DOUBLE z = y;
- if (z > L_(0.0))
- {
-
- if (z < TWO_MANT_DIG)
- {
-
- z += TWO_MANT_DIG;
- z -= TWO_MANT_DIG;
-
- if (z < y)
- z += L_(1.0);
- }
- }
- else if (z < L_(0.0))
- {
-
- if (z > L_(-1.0))
- z = MINUS_ZERO;
-
- else if (z > - TWO_MANT_DIG)
- {
-
- z -= TWO_MANT_DIG;
- z += TWO_MANT_DIG;
-
- if (z < y)
- z += L_(1.0);
- }
- }
- return z;
- }
|