← Back to branch summary

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

« back to all changes in this revision

Viewing changes to math/s_clogl.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:
math/s_clogl.c
1
1
/* Compute complex natural logarithm.
2
 
   Copyright (C) 1997, 1998, 2011 Free Software Foundation, Inc.
 
2
   Copyright (C) 1997-2012 Free Software Foundation, Inc.
3
3
   This file is part of the GNU C Library.
4
4
   Contributed by Ulrich Drepper <drepper@cygnus.com>, 1997.
5
5
 
14
14
   Lesser General Public License for more details.
15
15
 
16
16
   You should have received a copy of the GNU Lesser General Public
17
 
   License along with the GNU C Library; if not, write to the Free
18
 
   Software Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA
19
 
   02111-1307 USA.  */
 
17
   License along with the GNU C Library; if not, see
 
18
   <http://www.gnu.org/licenses/>.  */
20
19
 
21
20
#include <complex.h>
22
21
#include <math.h>
23
22
#include <math_private.h>
24
 
 
 
23
#include <float.h>
25
24
 
26
25
__complex__ long double
27
26
__clogl (__complex__ long double x)
41
40
  else if (__builtin_expect (rcls != FP_NAN && icls != FP_NAN, 1))
42
41
    {
43
42
      /* Neither real nor imaginary part is NaN.  */
44
 
      __real__ result = __ieee754_logl (__ieee754_hypotl (__real__ x,
45
 
                                                          __imag__ x));
 
43
      long double d;
 
44
      int scale = 0;
 
45
 
 
46
      if (fabsl (__real__ x) > LDBL_MAX / 2.0L
 
47
          || fabsl (__imag__ x) > LDBL_MAX / 2.0L)
 
48
        {
 
49
          scale = -1;
 
50
          __real__ x = __scalbnl (__real__ x, scale);
 
51
          __imag__ x = __scalbnl (__imag__ x, scale);
 
52
        }
 
53
      else if (fabsl (__real__ x) < LDBL_MIN
 
54
               && fabsl (__imag__ x) < LDBL_MIN)
 
55
        {
 
56
          scale = LDBL_MANT_DIG;
 
57
          __real__ x = __scalbnl (__real__ x, scale);
 
58
          __imag__ x = __scalbnl (__imag__ x, scale);
 
59
        }
 
60
 
 
61
      d = __ieee754_hypotl (__real__ x, __imag__ x);
 
62
 
 
63
      __real__ result = __ieee754_logl (d) - scale * M_LN2l;
46
64
      __imag__ result = __ieee754_atan2l (__imag__ x, __real__ x);
47
65
    }
48
66
  else