sysdeps/libm-ieee754/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  */
  20
  21 #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
  31 {
  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(__finite(y)&&y<0.0)
  48               if (signbit (x) && signbit (z))
  49                 return __kernel_standard(x,y,223); /* pow(-0.0,negative) */
  50               else
  51                 return __kernel_standard(x,y,243); /* pow(+0.0,negative) */
  52             return z;
  53         }
  54         if(!__finitel(z)) {
  55             if(__finitel(x)&&__finitel(y)) {
  56                 if(__isnanl(z))
  57                     return __kernel_standard(x,y,224); /* pow neg**non-int */
  58                 else
  59                     return __kernel_standard(x,y,221); /* pow overflow */
  60             }
  61         }
  62         if(z==0.0&&__finitel(x)&&__finitel(y))
  63             return __kernel_standard(x,y,222); /* pow underflow */
  64         return z;
  65 #endif
  66 }
  67 weak_alias (__powl, powl)