diff options
| author | Ulrich Drepper <drepper@gmail.com> | 2012-01-10 22:26:22 -0500 |
|---|---|---|
| committer | Ulrich Drepper <drepper@gmail.com> | 2012-01-10 22:26:22 -0500 |
| commit | daa891c0e8900d0f03c16fde9c09f8d6f9a26a7d (patch) | |
| tree | 3b1bf170ebd79ce294d3a09e2c71d133f4bcdd73 /sysdeps/ieee754/ldbl-96 | |
| parent | e58ef0f20437e121a06e0eb3d6c6c10c3630746e (diff) | |
| download | glibc-daa891c0e8900d0f03c16fde9c09f8d6f9a26a7d.tar glibc-daa891c0e8900d0f03c16fde9c09f8d6f9a26a7d.tar.gz glibc-daa891c0e8900d0f03c16fde9c09f8d6f9a26a7d.tar.bz2 glibc-daa891c0e8900d0f03c16fde9c09f8d6f9a26a7d.zip | |
Optimize s_modf
Diffstat (limited to 'sysdeps/ieee754/ldbl-96')
| -rw-r--r-- | sysdeps/ieee754/ldbl-96/s_modfl.c | 20 |
1 files changed, 4 insertions, 16 deletions
diff --git a/sysdeps/ieee754/ldbl-96/s_modfl.c b/sysdeps/ieee754/ldbl-96/s_modfl.c index 99c0752efa..c7659119f3 100644 --- a/sysdeps/ieee754/ldbl-96/s_modfl.c +++ b/sysdeps/ieee754/ldbl-96/s_modfl.c @@ -14,10 +14,6 @@ * ==================================================== */ -#if defined(LIBM_SCCS) && !defined(lint) -static char rcsid[] = "$NetBSD: $"; -#endif - /* * modfl(long double x, long double *iptr) * return fraction part of x, and return x's integral part in *iptr. @@ -31,18 +27,10 @@ static char rcsid[] = "$NetBSD: $"; #include "math.h" #include "math_private.h" -#ifdef __STDC__ static const long double one = 1.0; -#else -static long double one = 1.0; -#endif -#ifdef __STDC__ - long double __modfl(long double x, long double *iptr) -#else - long double __modfl(x, iptr) - long double x,*iptr; -#endif +long double +__modfl(long double x, long double *iptr) { int32_t i0,i1,j0; u_int32_t i,se; @@ -50,7 +38,7 @@ static long double one = 1.0; j0 = (se&0x7fff)-0x3fff; /* exponent of x */ if(j0<32) { /* integer part in high x */ if(j0<0) { /* |x|<1 */ - SET_LDOUBLE_WORDS(*iptr,se&0x8000,0,0); /* *iptr = +-0 */ + SET_LDOUBLE_WORDS(*iptr,se&0x8000,0,0); /* *iptr = +-0 */ return x; } else { i = (0x7fffffff)>>j0; @@ -63,7 +51,7 @@ static long double one = 1.0; return x - *iptr; } } - } else if (j0>63) { /* no fraction part */ + } else if (__builtin_expect(j0>63, 0)) { /* no fraction part */ *iptr = x*one; /* We must handle NaNs separately. */ if (j0 == 0x4000 && ((i0 & 0x7fffffff) | i1)) |