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.
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.
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.
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
21 #include <math_private.h>
26 __powf (float x, float y)
28 float z = __ieee754_powf (x, y);
29 if (__builtin_expect (!__finitef (z), 0))
31 if (_LIB_VERSION != _IEEE_)
37 return __kernel_standard_f (x, y, 142);
39 else if (__finitef (x) && __finitef (y))
42 /* pow neg**non-int */
43 return __kernel_standard_f (x, y, 124);
44 else if (x == 0.0f && y < 0.0f)
46 if (signbit (x) && signbit (z))
47 /* pow(-0.0,negative) */
48 return __kernel_standard_f (x, y, 123);
50 /* pow(+0.0,negative) */
51 return __kernel_standard_f (x, y, 143);
55 return __kernel_standard_f (x, y, 121);
59 else if (__builtin_expect (z == 0.0f, 0) && __finitef (x) && __finitef (y)
60 && _LIB_VERSION != _IEEE_)
66 return __kernel_standard_f (x, y, 120);
70 return __kernel_standard_f (x, y, 122);
75 weak_alias (__powf, powf)