← Back to branch summary

~ubuntu-branches/ubuntu/precise/eglibc/precise-201308281639

« back to all changes in this revision

Viewing changes to math/w_powl.c

  • Committer: Package Import Robot
  • Author(s): Matthias Klose
  • Date: 2012-02-08 01:58:09 UTC
  • mfrom: (1.5.3) (288.1.12 precise)
  • Revision ID: package-import@ubuntu.com-20120208015809-ulscst7uteq3e22z
Tags: 2.15~pre6-0ubuntu10
Merge from Debian (r5151, 2.13-26).

Show diffs side-by-side

added added

removed removed

Lines of Context:
math/w_powl.c
1
 
/* w_powl.c -- long double version of w_pow.c.
2
 
 * Conversion to long double by Ulrich Drepper,
3
 
 * Cygnus Support, drepper@cygnus.com.
4
 
 */
5
 
 
6
 
/*
7
 
 * ====================================================
8
 
 * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
9
 
 *
10
 
 * Developed at SunPro, a Sun Microsystems, Inc. business.
11
 
 * Permission to use, copy, modify, and distribute this
12
 
 * software is freely granted, provided that this notice
13
 
 * is preserved.
14
 
 * ====================================================
15
 
 */
16
 
 
17
 
/*
18
 
 * wrapper powl(x,y) return x**y
19
 
 */
 
1
/* Copyright (C) 2011 Free Software Foundation, Inc.
 
2
   This file is part of the GNU C Library.
 
3
   Contributed by Ulrich Drepper <drepper@gmail.com>, 2011.
 
4
 
 
5
   The GNU C Library is free software; you can redistribute it and/or
 
6
   modify it under the terms of the GNU Lesser General Public
 
7
   License as published by the Free Software Foundation; either
 
8
   version 2.1 of the License, or (at your option) any later version.
 
9
 
 
10
   The GNU C Library is distributed in the hope that it will be useful,
 
11
   but WITHOUT ANY WARRANTY; without even the implied warranty of
 
12
   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
 
13
   Lesser General Public License for more details.
 
14
 
 
15
   You should have received a copy of the GNU Lesser General Public
 
16
   License along with the GNU C Library; if not, write to the Free
 
17
   Software Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA
 
18
   02111-1307 USA.  */
20
19
 
21
20
#include <math.h>
22
 
#include "math_private.h"
23
 
 
24
 
 
25
 
#ifdef __STDC__
26
 
        long double __powl(long double x, long double y)/* wrapper powl */
27
 
#else
28
 
        long double __powl(x,y)                 /* wrapper powl */
29
 
        long double x,y;
30
 
#endif
 
21
#include <math_private.h>
 
22
 
 
23
 
 
24
/* wrapper powl */
 
25
long double
 
26
__powl (long double x, long double y)
31
27
{
32
 
#ifdef _IEEE_LIBM
33
 
        return  __ieee754_powl(x,y);
34
 
#else
35
 
        long double z;
36
 
        z=__ieee754_powl(x,y);
37
 
        if(_LIB_VERSION == _IEEE_|| __isnanl(y)) return z;
38
 
        if(__isnanl(x)) {
39
 
            if(y==0.0)
40
 
                return __kernel_standard(x,y,242); /* pow(NaN,0.0) */
41
 
            else
42
 
                return z;
43
 
        }
44
 
        if(x==0.0) {
45
 
            if(y==0.0)
46
 
                return __kernel_standard(x,y,220); /* pow(0.0,0.0) */
47
 
            if(__finitel(y)&&y<0.0) {
48
 
              if (signbit (x) && signbit (z))
49
 
                return __kernel_standard(x,y,223); /* pow(-0.0,negative) */
 
28
  long double z = __ieee754_powl (x, y);
 
29
  if (__builtin_expect (!__finitel (z), 0))
 
30
    {
 
31
      if (_LIB_VERSION != _IEEE_)
 
32
        {
 
33
          if (__isnanl (x))
 
34
            {
 
35
              if (y == 0.0L)
 
36
                /* pow(NaN,0.0) */
 
37
                return __kernel_standard (x, y, 242);
 
38
            }
 
39
          else if (__finitel (x) && __finitel (y))
 
40
            {
 
41
              if (__isnanl (z))
 
42
                /* pow neg**non-int */
 
43
                return __kernel_standard (x, y, 224);
 
44
              else if (x == 0.0L && y < 0.0L)
 
45
                {
 
46
                  if (signbit (x) && signbit (z))
 
47
                    /* pow(-0.0,negative) */
 
48
                    return __kernel_standard (x, y, 223);
 
49
                  else
 
50
                    /* pow(+0.0,negative) */
 
51
                    return __kernel_standard (x, y, 243);
 
52
                }
50
53
              else
51
 
                return __kernel_standard(x,y,243); /* pow(+0.0,negative) */
52
 
            }
53
 
            return z;
54
 
        }
55
 
        if(!__finitel(z)) {
56
 
            if(__finitel(x)&&__finitel(y)) {
57
 
                if(__isnanl(z))
58
 
                    return __kernel_standard(x,y,224); /* pow neg**non-int */
59
 
                else
60
 
                    return __kernel_standard(x,y,221); /* pow overflow */
61
 
            }
62
 
        }
63
 
        if(z==0.0&&__finitel(x)&&__finitel(y))
64
 
            return __kernel_standard(x,y,222); /* pow underflow */
65
 
        return z;
66
 
#endif
 
54
                /* pow overflow */
 
55
                return __kernel_standard (x, y, 221);
 
56
            }
 
57
        }
 
58
    }
 
59
  else if (__builtin_expect (z == 0.0L, 0) && __finitel (x) && __finitel (y)
 
60
           && _LIB_VERSION != _IEEE_)
 
61
    {
 
62
      if (x == 0.0L)
 
63
        {
 
64
          if (y == 0.0L)
 
65
            /* pow(0.0,0.0) */
 
66
            return __kernel_standard (x, y, 220);
 
67
        }
 
68
      else
 
69
        /* pow underflow */
 
70
        return __kernel_standard (x, y, 222);
 
71
    }
 
72
 
 
73
  return z;
67
74
}
68
75
weak_alias (__powl, powl)