Fix the inaccuracy of j0f/j1f/y0f/y1f [BZ #14469, #14470, #14471, #14472]

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

diff --git a/math/w_tgamma_compat.c b/math/w_tgamma_compat.c
index f843475f107859ca56e2df9ebbb91530001e536a..910d2fe4904e9b04952aba51e64cbed800cbb936 100644 (file)
--- a/math/w_tgamma_compat.c
+++ b/math/w_tgamma_compat.c
@@ -19,6 +19,7 @@
 #include <math.h>
 #include <math_private.h>
 #include <math-svid-compat.h>
+#include <libm-alias-double.h>
 
 #if LIBM_SVID_COMPAT
 double
@@ -32,7 +33,7 @@ __tgamma(double x)
           && _LIB_VERSION != _IEEE_) {
          if (x == 0.0)
            return __kernel_standard(x,x,50); /* tgamma pole */
-         else if(__floor(x)==x&&x<0.0)
+         else if(floor(x)==x&&x<0.0)
            return __kernel_standard(x,x,41); /* tgamma domain */
          else if (y == 0)
            __set_errno (ERANGE); /* tgamma underflow */
@@ -41,9 +42,5 @@ __tgamma(double x)
        }
        return local_signgam < 0 ? -y : y;
 }
-weak_alias (__tgamma, tgamma)
-# ifdef NO_LONG_DOUBLE
-strong_alias (__tgamma, __tgammal)
-weak_alias (__tgamma, tgammal)
-# endif
+libm_alias_double (__tgamma, tgamma)
 #endif