x86-64: Add vector log1p/log1pf implementation to libmvec

[thirdparty/glibc.git] / math / w_powf.c

diff --git a/math/w_powf.c b/math/w_powf.c
index 2e786e2f4e644d37acc71bd68e7575538baa573f..aea9058ebe3dd36ef99e74b252afb9e5db9120a0 100644 (file)
--- a/math/w_powf.c
+++ b/math/w_powf.c
@@ -1,74 +1,8 @@
-/* Copyright (C) 2011-2016 Free Software Foundation, Inc.
-   This file is part of the GNU C Library.
-   Contributed by Ulrich Drepper <drepper@gmail.com>, 2011.
-
-   The GNU C Library is free software; you can redistribute it and/or
-   modify it under the terms of the GNU Lesser General Public
-   License as published by the Free Software Foundation; either
-   version 2.1 of the License, or (at your option) any later version.
-
-   The GNU C Library is distributed in the hope that it will be useful,
-   but WITHOUT ANY WARRANTY; without even the implied warranty of
-   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
-   Lesser General Public License for more details.
-
-   You should have received a copy of the GNU Lesser General Public
-   License along with the GNU C Library; if not, see
-   <http://www.gnu.org/licenses/>.  */
-
-#include <math.h>
-#include <math_private.h>
-
-
-/* wrapper powf */
-float
-__powf (float x, float y)
-{
-  float z = __ieee754_powf (x, y);
-  if (__glibc_unlikely (!isfinite (z)))
-    {
-      if (_LIB_VERSION != _IEEE_)
-       {
-         if (isnan (x))
-           {
-             if (y == 0.0f)
-               /* pow(NaN,0.0) */
-               return __kernel_standard_f (x, y, 142);
-           }
-         else if (isfinite (x) && isfinite (y))
-           {
-             if (isnan (z))
-               /* pow neg**non-int */
-               return __kernel_standard_f (x, y, 124);
-             else if (x == 0.0f && y < 0.0f)
-               {
-                 if (signbit (x) && signbit (z))
-                   /* pow(-0.0,negative) */
-                   return __kernel_standard_f (x, y, 123);
-                 else
-                   /* pow(+0.0,negative) */
-                   return __kernel_standard_f (x, y, 143);
-               }
-             else
-               /* pow overflow */
-               return __kernel_standard_f (x, y, 121);
-           }
-       }
-    }
-  else if (__builtin_expect (z == 0.0f, 0) && isfinite (x) && isfinite (y)
-          && _LIB_VERSION != _IEEE_)
-    {
-      if (x == 0.0f)
-       {
-         if (y == 0.0f)
-           /* pow(0.0,0.0) */
-           return __kernel_standard_f (x, y, 120);
-       }
-      else
-       /* pow underflow */
-       return __kernel_standard_f (x, y, 122);
-    }
-
-  return z;
-}
-weak_alias (__powf, powf)
+#include <math-type-macros-float.h>
+#undef __USE_WRAPPER_TEMPLATE
+#define __USE_WRAPPER_TEMPLATE 1
+#undef declare_mgen_alias
+#define declare_mgen_alias(a, b)
+#include <w_pow_template.c>
+versioned_symbol (libm, __powf, powf, GLIBC_2_27);
+libm_alias_float_other (__pow, pow)