← Back to branch summary

~ubuntu-branches/ubuntu/utopic/eglibc/utopic

« back to all changes in this revision

Viewing changes to sysdeps/ieee754/dbl-64/s_logb.c

  • Committer: Package Import Robot
  • Author(s): Adam Conrad
  • Date: 2012-10-26 05:14:58 UTC
  • mfrom: (1.5.1) (4.4.22 experimental)
  • Revision ID: package-import@ubuntu.com-20121026051458-oryotr4i03ob5pab
Tags: 2.16-0ubuntu1
* Merge with unreleased 2.16 in Debian experimental, remaining changes:
  - Drop the Breaks line from libc6, which refers to a Debian transition
  - Remove the libc6 recommends on libc6-i686, which we don't build
  - Enable libc6{,-dev}-armel on armhf and libc6{-dev}-armhf on armel
  - Ship update-locale and validlocale in /usr/sbin in libc-bin
  - Don't build locales or locales-all in Ubuntu, we rely on langpacks
  - Heavily mangle the way we do service restarting on major upgrades
  - Use different MIN_KERNEL_SUPPORTED versions than Debian, due to
    buildd needs.  This should be universally bumped to 3.2.0 once all
    our buildds (including the PPA guests) are running precise kernels
  - Build i386 variants as -march=i686, build amd64 with -O3, and build
    ppc64 variants (both 64-bit and 32-bit) with -O3 -fno-tree-vectorize
  - Re-enable unsubmitted-ldconfig-cache-abi.diff and rebuild the cache
    on upgrades from previous versions that used a different constant
  - debian/patches/any/local-CVE-2012-3406.diff: switch to malloc when
    array grows too large to handle via alloca extension (CVE-2012-3406)
  - Build generic i386/i686 flavour with -mno-tls-direct-seg-refs
* Changes added/dropped with this merge while reducing our delta:
  - Stop building glibc docs from the eglibc source, and instead make
    the glibc-docs stub have a hard dependency on glibc-doc-reference
  - Remove outdated conflicts against ancient versions of ia32-libs
  - Drop the tzdata dependency from libc6, it's in required and minimal
  - Use gcc-4.7/g++-4.7 by default on all our supported architectures
  - Save our historical changelog as changelog.ubuntu in the source
  - Drop nscd's libaudit build-dep for now, as libaudit is in universe
  - Drop the unnecessary Breaks from libc6 to locales and locales-all
  - Ship xen's ld.so.conf.d snippet as /etc/ld.so.conf.d/libc6-xen.conf
* Disable hard failures on the test suite for the first upload to raring

Show diffs side-by-side

added added

removed removed

Lines of Context:
sysdeps/ieee754/dbl-64/s_logb.c
10
10
 * ====================================================
11
11
 */
12
12
 
13
 
#if defined(LIBM_SCCS) && !defined(lint)
14
 
static char rcsid[] = "$NetBSD: s_logb.c,v 1.8 1995/05/10 20:47:50 jtc Exp $";
15
 
#endif
16
 
 
17
13
/*
18
14
 * double logb(x)
19
15
 * IEEE 754 logb. Included to pass IEEE test suite. Not recommend.
20
16
 * Use ilogb instead.
21
17
 */
22
18
 
23
 
#include "math.h"
24
 
#include "math_private.h"
 
19
#include <math.h>
 
20
#include <math_private.h>
25
21
 
26
 
#ifdef __STDC__
27
 
        double __logb(double x)
28
 
#else
29
 
        double __logb(x)
30
 
        double x;
31
 
#endif
 
22
double
 
23
__logb (double x)
32
24
{
33
 
        int32_t lx,ix;
34
 
        EXTRACT_WORDS(ix,lx,x);
35
 
        ix &= 0x7fffffff;                       /* high |x| */
36
 
        if((ix|lx)==0) return -1.0/fabs(x);
37
 
        if(ix>=0x7ff00000) return x*x;
38
 
        if((ix>>=20)==0)                        /* IEEE 754 logb */
39
 
                return -1022.0;
40
 
        else
41
 
                return (double) (ix-1023);
 
25
  int32_t lx, ix, rix;
 
26
  double ret;
 
27
 
 
28
  EXTRACT_WORDS (ix, lx, x);
 
29
  ix &= 0x7fffffff;             /* high |x| */
 
30
  if ((ix | lx) == 0)
 
31
    return -1.0 / fabs (x);
 
32
  if (ix >= 0x7ff00000)
 
33
    return x * x;
 
34
  if (__builtin_expect ((rix = ix >> 20) == 0, 0))
 
35
    {
 
36
      /* POSIX specifies that denormal number is treated as
 
37
         though it were normalized.  */
 
38
      int ma;
 
39
      if (ix == 0)
 
40
        ma = __builtin_clz (lx) + 32;
 
41
      else
 
42
        ma = __builtin_clz (ix);
 
43
      rix -= ma - 12;
 
44
    }
 
45
  ret = (double) (rix - 1023);
 
46
  /* The test is to avoid logb_downward (0.0) == -0.0.  */
 
47
  return ret == -0.0 ? 0.0 : ret;
42
48
}
43
49
weak_alias (__logb, logb)
44
50
#ifdef NO_LONG_DOUBLE
45
 
strong_alias (__logb, __logbl)
46
 
weak_alias (__logb, logbl)
 
51
strong_alias (__logb, __logbl) weak_alias (__logb, logbl)
47
52
#endif