← Back to branch summary

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

« back to all changes in this revision

Viewing changes to math/s_cexpf.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_cexpf.c
1
1
/* Return value of complex exponential function for float complex value.
2
 
   Copyright (C) 1997, 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 <fenv.h>
23
22
#include <math.h>
24
23
#include <math_private.h>
25
 
 
 
24
#include <float.h>
26
25
 
27
26
__complex__ float
28
27
__cexpf (__complex__ float x)
37
36
      if (__builtin_expect (icls >= FP_ZERO, 1))
38
37
        {
39
38
          /* Imaginary part is finite.  */
40
 
          float exp_val = __ieee754_expf (__real__ x);
 
39
          const int t = (int) ((FLT_MAX_EXP - 1) * M_LN2);
41
40
          float sinix, cosix;
42
41
 
43
42
          __sincosf (__imag__ x, &sinix, &cosix);
44
43
 
45
 
          if (isfinite (exp_val))
46
 
            {
 
44
          if (__real__ x > t)
 
45
            {
 
46
              float exp_t = __ieee754_expf (t);
 
47
              __real__ x -= t;
 
48
              sinix *= exp_t;
 
49
              cosix *= exp_t;
 
50
              if (__real__ x > t)
 
51
                {
 
52
                  __real__ x -= t;
 
53
                  sinix *= exp_t;
 
54
                  cosix *= exp_t;
 
55
                }
 
56
            }
 
57
          if (__real__ x > t)
 
58
            {
 
59
              /* Overflow (original real part of x > 3t).  */
 
60
              __real__ retval = FLT_MAX * cosix;
 
61
              __imag__ retval = FLT_MAX * sinix;
 
62
            }
 
63
          else
 
64
            {
 
65
              float exp_val = __ieee754_expf (__real__ x);
47
66
              __real__ retval = exp_val * cosix;
48
67
              __imag__ retval = exp_val * sinix;
49
68
            }
50
 
          else
51
 
            {
52
 
              __real__ retval = __copysignf (exp_val, cosix);
53
 
              __imag__ retval = __copysignf (exp_val, sinix);
54
 
            }
55
69
        }
56
70
      else
57
71
        {