← Back to branch summary

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

« back to all changes in this revision

Viewing changes to math/s_csqrtl.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_csqrtl.c
1
1
/* Complex square root of long double value.
2
 
   Copyright (C) 1997, 1998, 2005, 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
   Based on an algorithm by Stephen L. Moshier <moshier@world.std.com>.
5
5
   Contributed by Ulrich Drepper <drepper@cygnus.com>, 1997.
15
15
   Lesser General Public License for more details.
16
16
 
17
17
   You should have received a copy of the GNU Lesser General Public
18
 
   License along with the GNU C Library; if not, write to the Free
19
 
   Software Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA
20
 
   02111-1307 USA.  */
 
18
   License along with the GNU C Library; if not, see
 
19
   <http://www.gnu.org/licenses/>.  */
21
20
 
22
21
#include <complex.h>
23
22
#include <math.h>
24
23
#include <math_private.h>
25
 
 
 
24
#include <float.h>
26
25
 
27
26
__complex__ long double
28
27
__csqrtl (__complex__ long double x)
84
83
      else
85
84
        {
86
85
          long double d, r, s;
 
86
          int scale = 0;
 
87
 
 
88
          if (fabsl (__real__ x) > LDBL_MAX / 2.0L
 
89
              || fabsl (__imag__ x) > LDBL_MAX / 2.0L)
 
90
            {
 
91
              scale = 1;
 
92
              __real__ x = __scalbnl (__real__ x, -2 * scale);
 
93
              __imag__ x = __scalbnl (__imag__ x, -2 * scale);
 
94
            }
 
95
          else if (fabsl (__real__ x) < LDBL_MIN
 
96
                   && fabsl (__imag__ x) < LDBL_MIN)
 
97
            {
 
98
              scale = -(LDBL_MANT_DIG / 2);
 
99
              __real__ x = __scalbnl (__real__ x, -2 * scale);
 
100
              __imag__ x = __scalbnl (__imag__ x, -2 * scale);
 
101
            }
87
102
 
88
103
          d = __ieee754_hypotl (__real__ x, __imag__ x);
89
104
          /* Use the identity   2  Re res  Im res = Im x
99
114
              r = fabsl ((0.5L * __imag__ x) / s);
100
115
            }
101
116
 
 
117
          if (scale)
 
118
            {
 
119
              r = __scalbnl (r, scale);
 
120
              s = __scalbnl (s, scale);
 
121
            }
 
122
 
102
123
          __real__ res = r;
103
124
          __imag__ res = __copysignl (s, __imag__ x);
104
125
        }