← Back to branch summary

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

« back to all changes in this revision

Viewing changes to math/s_csinf.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_csinf.c
1
1
/* Complex sine function for float.
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
 
 
25
23
#include <math_private.h>
26
 
 
 
24
#include <float.h>
27
25
 
28
26
__complex__ float
29
27
__csinf (__complex__ float x)
41
39
      if (__builtin_expect (rcls >= FP_ZERO, 1))
42
40
        {
43
41
          /* Real part is finite.  */
44
 
          float sinh_val = __ieee754_sinhf (__imag__ x);
45
 
          float cosh_val = __ieee754_coshf (__imag__ x);
 
42
          const int t = (int) ((FLT_MAX_EXP - 1) * M_LN2);
46
43
          float sinix, cosix;
47
44
 
48
45
          __sincosf (__real__ x, &sinix, &cosix);
49
46
 
50
 
          __real__ retval = cosh_val * sinix;
51
 
          __imag__ retval = sinh_val * cosix;
 
47
          if (fabsf (__imag__ x) > t)
 
48
            {
 
49
              float exp_t = __ieee754_expf (t);
 
50
              float ix = fabsf (__imag__ x);
 
51
              if (signbit (__imag__ x))
 
52
                cosix = -cosix;
 
53
              ix -= t;
 
54
              sinix *= exp_t / 2.0f;
 
55
              cosix *= exp_t / 2.0f;
 
56
              if (ix > t)
 
57
                {
 
58
                  ix -= t;
 
59
                  sinix *= exp_t;
 
60
                  cosix *= exp_t;
 
61
                }
 
62
              if (ix > t)
 
63
                {
 
64
                  /* Overflow (original imaginary part of x > 3t).  */
 
65
                  __real__ retval = FLT_MAX * sinix;
 
66
                  __imag__ retval = FLT_MAX * cosix;
 
67
                }
 
68
              else
 
69
                {
 
70
                  float exp_val = __ieee754_expf (ix);
 
71
                  __real__ retval = exp_val * sinix;
 
72
                  __imag__ retval = exp_val * cosix;
 
73
                }
 
74
            }
 
75
          else
 
76
            {
 
77
              __real__ retval = __ieee754_coshf (__imag__ x) * sinix;
 
78
              __imag__ retval = __ieee754_sinhf (__imag__ x) * cosix;
 
79
            }
52
80
 
53
81
          if (negate)
54
82
            __real__ retval = -__real__ retval;