]> git.ipfire.org Git - thirdparty/glibc.git/commitdiff
Use ceil functions not __ceil functions in glibc libm.
authorJoseph Myers <joseph@codesourcery.com>
Mon, 17 Sep 2018 20:42:06 +0000 (20:42 +0000)
committerJoseph Myers <joseph@codesourcery.com>
Mon, 17 Sep 2018 20:42:06 +0000 (20:42 +0000)
Continuing the move to use, within libm, public names for libm
functions that can be inlined as built-in functions on many
architectures, this patch moves calls to __ceil functions to call the
corresponding ceil names instead, with asm redirection to __ceil when
the calls are not inlined.

Tested for x86_64, and with build-many-glibcs.py.

* include/math.h [!_ISOMAC && !(__FINITE_MATH_ONLY__ &&
__FINITE_MATH_ONLY__ > 0) && !NO_MATH_REDIRECT] (ceil): Redirect
using MATH_REDIRECT.
* sysdeps/aarch64/fpu/s_ceil.c: Define NO_MATH_REDIRECT before
header inclusion.
* sysdeps/aarch64/fpu/s_ceilf.c: Likewise.
* sysdeps/ieee754/dbl-64/s_ceil.c: Likewise.
* sysdeps/ieee754/dbl-64/wordsize-64/s_ceil.c: Likewise.
* sysdeps/ieee754/float128/s_ceilf128.c: Likewise.
* sysdeps/ieee754/flt-32/s_ceilf.c: Likewise.
* sysdeps/ieee754/ldbl-128/s_ceill.c: Likewise.
* sysdeps/ieee754/ldbl-128ibm/s_ceill.c: Likewise.
* sysdeps/m68k/m680x0/fpu/s_ceil_template.c: Likewise.
* sysdeps/powerpc/powerpc32/power4/fpu/multiarch/s_ceil.c: Likewise.
* sysdeps/powerpc/powerpc32/power4/fpu/multiarch/s_ceilf.c: Likewise.
* sysdeps/powerpc/powerpc64/fpu/multiarch/s_ceil.c: Likewise.
* sysdeps/powerpc/powerpc64/fpu/multiarch/s_ceilf.c: Likewise.
* sysdeps/riscv/rv64/rvd/s_ceil.c: Likewise.
* sysdeps/riscv/rvf/s_ceilf.c: Likewise.
* sysdeps/sparc/sparc64/fpu/multiarch/s_ceil.c: Likewise.
* sysdeps/sparc/sparc64/fpu/multiarch/s_ceilf.c: Likewise.
* sysdeps/x86_64/fpu/multiarch/s_ceil.c: Likewise.
* sysdeps/x86_64/fpu/multiarch/s_ceilf.c: Likewise.
* sysdeps/powerpc/fpu/math_private.h [_ARCH_PWR5X] (__ceil):
Remove macro.
* sysdeps/ieee754/dbl-64/e_gamma_r.c (gamma_positive): Use ceil
functions instead of __ceil variants.
* sysdeps/ieee754/flt-32/e_gammaf_r.c (gammaf_positive): Likewise.
* sysdeps/ieee754/ldbl-128/e_gammal_r.c (gammal_positive):
Likewise.
* sysdeps/ieee754/ldbl-128ibm/e_gammal_r.c (gammal_positive):
Likewise.
* sysdeps/ieee754/ldbl-128ibm/s_truncl.c (__truncl): Likewise.
* sysdeps/ieee754/ldbl-96/e_gammal_r.c (gammal_positive):
Likewise.
* sysdeps/powerpc/power5+/fpu/s_modf.c (__modf): Likewise.
* sysdeps/powerpc/power5+/fpu/s_modff.c (__modff): Likewise.

30 files changed:
ChangeLog
include/math.h
sysdeps/aarch64/fpu/s_ceil.c
sysdeps/aarch64/fpu/s_ceilf.c
sysdeps/ieee754/dbl-64/e_gamma_r.c
sysdeps/ieee754/dbl-64/s_ceil.c
sysdeps/ieee754/dbl-64/wordsize-64/s_ceil.c
sysdeps/ieee754/float128/s_ceilf128.c
sysdeps/ieee754/flt-32/e_gammaf_r.c
sysdeps/ieee754/flt-32/s_ceilf.c
sysdeps/ieee754/ldbl-128/e_gammal_r.c
sysdeps/ieee754/ldbl-128/s_ceill.c
sysdeps/ieee754/ldbl-128ibm/e_gammal_r.c
sysdeps/ieee754/ldbl-128ibm/s_ceill.c
sysdeps/ieee754/ldbl-128ibm/s_truncl.c
sysdeps/ieee754/ldbl-96/e_gammal_r.c
sysdeps/m68k/m680x0/fpu/s_ceil_template.c
sysdeps/powerpc/fpu/math_private.h
sysdeps/powerpc/power5+/fpu/s_modf.c
sysdeps/powerpc/power5+/fpu/s_modff.c
sysdeps/powerpc/powerpc32/power4/fpu/multiarch/s_ceil.c
sysdeps/powerpc/powerpc32/power4/fpu/multiarch/s_ceilf.c
sysdeps/powerpc/powerpc64/fpu/multiarch/s_ceil.c
sysdeps/powerpc/powerpc64/fpu/multiarch/s_ceilf.c
sysdeps/riscv/rv64/rvd/s_ceil.c
sysdeps/riscv/rvf/s_ceilf.c
sysdeps/sparc/sparc64/fpu/multiarch/s_ceil.c
sysdeps/sparc/sparc64/fpu/multiarch/s_ceilf.c
sysdeps/x86_64/fpu/multiarch/s_ceil.c
sysdeps/x86_64/fpu/multiarch/s_ceilf.c

index 1532d35b9400a253da505d6f2f1059895c3e6080..b2a3a1f686f9d4d330f4b4f479c72a0e4778846e 100644 (file)
--- a/ChangeLog
+++ b/ChangeLog
@@ -1,5 +1,43 @@
 2018-09-17  Joseph Myers  <joseph@codesourcery.com>
 
+       * include/math.h [!_ISOMAC && !(__FINITE_MATH_ONLY__ &&
+       __FINITE_MATH_ONLY__ > 0) && !NO_MATH_REDIRECT] (ceil): Redirect
+       using MATH_REDIRECT.
+       * sysdeps/aarch64/fpu/s_ceil.c: Define NO_MATH_REDIRECT before
+       header inclusion.
+       * sysdeps/aarch64/fpu/s_ceilf.c: Likewise.
+       * sysdeps/ieee754/dbl-64/s_ceil.c: Likewise.
+       * sysdeps/ieee754/dbl-64/wordsize-64/s_ceil.c: Likewise.
+       * sysdeps/ieee754/float128/s_ceilf128.c: Likewise.
+       * sysdeps/ieee754/flt-32/s_ceilf.c: Likewise.
+       * sysdeps/ieee754/ldbl-128/s_ceill.c: Likewise.
+       * sysdeps/ieee754/ldbl-128ibm/s_ceill.c: Likewise.
+       * sysdeps/m68k/m680x0/fpu/s_ceil_template.c: Likewise.
+       * sysdeps/powerpc/powerpc32/power4/fpu/multiarch/s_ceil.c: Likewise.
+       * sysdeps/powerpc/powerpc32/power4/fpu/multiarch/s_ceilf.c: Likewise.
+       * sysdeps/powerpc/powerpc64/fpu/multiarch/s_ceil.c: Likewise.
+       * sysdeps/powerpc/powerpc64/fpu/multiarch/s_ceilf.c: Likewise.
+       * sysdeps/riscv/rv64/rvd/s_ceil.c: Likewise.
+       * sysdeps/riscv/rvf/s_ceilf.c: Likewise.
+       * sysdeps/sparc/sparc64/fpu/multiarch/s_ceil.c: Likewise.
+       * sysdeps/sparc/sparc64/fpu/multiarch/s_ceilf.c: Likewise.
+       * sysdeps/x86_64/fpu/multiarch/s_ceil.c: Likewise.
+       * sysdeps/x86_64/fpu/multiarch/s_ceilf.c: Likewise.
+       * sysdeps/powerpc/fpu/math_private.h [_ARCH_PWR5X] (__ceil):
+       Remove macro.
+       * sysdeps/ieee754/dbl-64/e_gamma_r.c (gamma_positive): Use ceil
+       functions instead of __ceil variants.
+       * sysdeps/ieee754/flt-32/e_gammaf_r.c (gammaf_positive): Likewise.
+       * sysdeps/ieee754/ldbl-128/e_gammal_r.c (gammal_positive):
+       Likewise.
+       * sysdeps/ieee754/ldbl-128ibm/e_gammal_r.c (gammal_positive):
+       Likewise.
+       * sysdeps/ieee754/ldbl-128ibm/s_truncl.c (__truncl): Likewise.
+       * sysdeps/ieee754/ldbl-96/e_gammal_r.c (gammal_positive):
+       Likewise.
+       * sysdeps/powerpc/power5+/fpu/s_modf.c (__modf): Likewise.
+       * sysdeps/powerpc/power5+/fpu/s_modff.c (__modff): Likewise.
+
        [BZ #21286]
        * sysdeps/unix/sysv/linux/bits/siginfo-consts.h (SI_DETHREAD): New
        constant.
index f93ecaf6e63f59f30ebbc7e1b9fd183c9f383c10..e96f9712575d289a3ff62696b34860c7e3e7ae6e 100644 (file)
@@ -101,6 +101,7 @@ fabsf128 (_Float128 x)
 #   endif
 #   define MATH_REDIRECT_UNARY_ARGS(TYPE) TYPE
 MATH_REDIRECT (sqrt, "__ieee754_", MATH_REDIRECT_UNARY_ARGS)
+MATH_REDIRECT (ceil, "__", MATH_REDIRECT_UNARY_ARGS)
 MATH_REDIRECT (floor, "__", MATH_REDIRECT_UNARY_ARGS)
 MATH_REDIRECT (rint, "__", MATH_REDIRECT_UNARY_ARGS)
 #  endif
index 4e3b5fdc4dde4bd560fe3a7e7b66672f4fbc4a19..b374da5b2f4564a0820571ca2d0157d83c6c2921 100644 (file)
@@ -16,6 +16,7 @@
    License along with the GNU C Library; if not, see
    <http://www.gnu.org/licenses/>.  */
 
+#define NO_MATH_REDIRECT
 #include <math.h>
 #include <libm-alias-double.h>
 
index d01dadc2460b966c16fd3ec39d359946dd42fc91..417395440e1f975ab43628b947c1aae205680878 100644 (file)
@@ -16,6 +16,7 @@
    License along with the GNU C Library; if not, see
    <http://www.gnu.org/licenses/>.  */
 
+#define NO_MATH_REDIRECT
 #include <math.h>
 #include <libm-alias-float.h>
 
index 9f7af6f0c9318bbad80511689c176799cbffb9af..823ffdbced5215af54562d92eb88108206e182bf 100644 (file)
@@ -61,7 +61,7 @@ gamma_positive (double x, int *exp2_adj)
     {
       /* Adjust into the range for using exp (lgamma).  */
       *exp2_adj = 0;
-      double n = __ceil (x - 1.5);
+      double n = ceil (x - 1.5);
       double x_adj = x - n;
       double eps;
       double prod = __gamma_product (x_adj, 0, n, &eps);
@@ -78,7 +78,7 @@ gamma_positive (double x, int *exp2_adj)
        {
          /* Adjust into the range for applying Stirling's
             approximation.  */
-         double n = __ceil (12.0 - x);
+         double n = ceil (12.0 - x);
          x_adj = math_narrow_eval (x + n);
          x_eps = (x - (x_adj - n));
          prod = __gamma_product (x_adj - n, x_eps, n, &eps);
index 5a7434c737812fbd4bff3a23f5ac4dffb5f62745..4c96286c69fc6cc6f545bca932a7053615849f15 100644 (file)
@@ -17,6 +17,7 @@
  *     Bit twiddling.
  */
 
+#define NO_MATH_REDIRECT
 #include <math.h>
 #include <math_private.h>
 #include <libm-alias-double.h>
index b99829d2b0d71c69df75936d0f9afe4c24d7330b..4bb93d06331ccabd18b9527614c44fcf04fdd842 100644 (file)
@@ -17,6 +17,7 @@
  *     Bit twiddling.
  */
 
+#define NO_MATH_REDIRECT
 #include <math.h>
 #include <math_private.h>
 #include <libm-alias-double.h>
index 0af15f5f03ec75b28f12ffb243d781c04e9f79bd..34949ac98ffcc59cc7c79dc578aafc6ed4a5e300 100644 (file)
@@ -1,2 +1,3 @@
+#define NO_MATH_REDIRECT
 #include <float128_private.h>
 #include "../ldbl-128/s_ceill.c"
index 28da9c46d27e1f2ff4253b4fa5462c2c574d4726..68f8f44c6533f8338d03e47a2090e9b294a1811d 100644 (file)
@@ -71,7 +71,7 @@ gammaf_positive (float x, int *exp2_adj)
        {
          /* Adjust into the range for applying Stirling's
             approximation.  */
-         float n = __ceilf (4.0f - x);
+         float n = ceilf (4.0f - x);
          x_adj = math_narrow_eval (x + n);
          x_eps = (x - (x_adj - n));
          prod = __gamma_productf (x_adj - n, x_eps, n, &eps);
index f289ec23415044c90494443d0d31fc4bf3379d31..25cba0807c8e0c8694d493d4d3076c735d9e5a18 100644 (file)
@@ -13,6 +13,7 @@
  * ====================================================
  */
 
+#define NO_MATH_REDIRECT
 #include <math.h>
 #include <math_private.h>
 #include <libm-alias-float.h>
index 9fd0dc046dc6a20d8966fc13c66422806d1f68cb..f620463ef99f375378eea52e59557cd746fe6a4e 100644 (file)
@@ -69,7 +69,7 @@ gammal_positive (_Float128 x, int *exp2_adj)
     {
       /* Adjust into the range for using exp (lgamma).  */
       *exp2_adj = 0;
-      _Float128 n = __ceill (x - L(1.5));
+      _Float128 n = ceill (x - L(1.5));
       _Float128 x_adj = x - n;
       _Float128 eps;
       _Float128 prod = __gamma_productl (x_adj, 0, n, &eps);
@@ -86,7 +86,7 @@ gammal_positive (_Float128 x, int *exp2_adj)
        {
          /* Adjust into the range for applying Stirling's
             approximation.  */
-         _Float128 n = __ceill (24 - x);
+         _Float128 n = ceill (24 - x);
          x_adj = x + n;
          x_eps = (x - (x_adj - n));
          prod = __gamma_productl (x_adj - n, x_eps, n, &eps);
index e6aba5f2afd23d3732ce35ad99d3d8cd3decbb52..2ec55de25a27a74ab6edf6d2e4c6125ebcdab44c 100644 (file)
@@ -24,6 +24,7 @@ static char rcsid[] = "$NetBSD: $";
  *     Bit twiddling.
  */
 
+#define NO_MATH_REDIRECT
 #include <math.h>
 #include <math_private.h>
 #include <libm-alias-ldouble.h>
index 8906da48461d667a1e9f13b4377123cc4b099121..c1e0efa5e12df352d933fb570dffe59fffb3fd60 100644 (file)
@@ -68,7 +68,7 @@ gammal_positive (long double x, int *exp2_adj)
     {
       /* Adjust into the range for using exp (lgamma).  */
       *exp2_adj = 0;
-      long double n = __ceill (x - 1.5L);
+      long double n = ceill (x - 1.5L);
       long double x_adj = x - n;
       long double eps;
       long double prod = __gamma_productl (x_adj, 0, n, &eps);
@@ -85,7 +85,7 @@ gammal_positive (long double x, int *exp2_adj)
        {
          /* Adjust into the range for applying Stirling's
             approximation.  */
-         long double n = __ceill (23.0L - x);
+         long double n = ceill (23.0L - x);
          x_adj = x + n;
          x_eps = (x - (x_adj - n));
          prod = __gamma_productl (x_adj - n, x_eps, n, &eps);
index 7dcff021c45486b6d4447e8c9ea028ae5b1be789..800835411a2273ac2ab221b1a071b23a59527e5f 100644 (file)
@@ -17,6 +17,7 @@
    License along with the GNU C Library; if not, see
    <http://www.gnu.org/licenses/>.  */
 
+#define NO_MATH_REDIRECT
 #include <math.h>
 #include <math_private.h>
 #include <math_ldbl_opt.h>
index a366c42599e82a2cfb53a087cd19c86e042948bb..96465678838b3d00687a8f000cf9bb769d222194 100644 (file)
@@ -47,7 +47,7 @@ __truncl (long double x)
       else
        {
          /* The high part is a nonzero integer.  */
-         lo = xh > 0 ? floor (xl) : __ceil (xl);
+         lo = xh > 0 ? floor (xl) : ceil (xl);
          xh = hi;
          xl = lo;
          ldbl_canonicalize_int (&xh, &xl);
index e88a137abbc8b14105df5f1c2fe9c10ecb5611e9..addb60d748fd7ed0d2850b3ddaa2d74b3520c183 100644 (file)
@@ -62,7 +62,7 @@ gammal_positive (long double x, int *exp2_adj)
     {
       /* Adjust into the range for using exp (lgamma).  */
       *exp2_adj = 0;
-      long double n = __ceill (x - 1.5L);
+      long double n = ceill (x - 1.5L);
       long double x_adj = x - n;
       long double eps;
       long double prod = __gamma_productl (x_adj, 0, n, &eps);
@@ -79,7 +79,7 @@ gammal_positive (long double x, int *exp2_adj)
        {
          /* Adjust into the range for applying Stirling's
             approximation.  */
-         long double n = __ceill (13.0L - x);
+         long double n = ceill (13.0L - x);
          x_adj = x + n;
          x_eps = (x - (x_adj - n));
          prod = __gamma_productl (x_adj - n, x_eps, n, &eps);
index a4d032d90c20d287f093f4ac821f72b4312f3648..be32b19eb6cf868ca2741cec51119090aa97dd5f 100644 (file)
@@ -16,6 +16,7 @@
    License along with the GNU C Library.  If not, see
    <http://www.gnu.org/licenses/>.  */
 
+#define NO_MATH_REDIRECT
 #include <math.h>
 
 FLOAT
index dc6dbfa474c97b97012000a9a9693ff4e57e6a0e..b71278ade46cd43bc72dd26170feac80e0279a4b 100644 (file)
@@ -77,26 +77,6 @@ __ieee754_sqrtf128 (_Float128 __x)
      __z; })
 # endif
 
-# ifndef __ceil
-#  define __ceil(x)                    \
-    ({ double __z;                     \
-     __asm __volatile (                        \
-       "       frip %0,%1\n"           \
-               : "=f" (__z)            \
-               : "f" (x));             \
-     __z; })
-# endif
-# ifndef __ceilf
-#  define __ceilf(x)                   \
-    ({ float __z;                      \
-     __asm __volatile (                        \
-       "       frip %0,%1\n"           \
-       "       frsp %0,%0\n"           \
-               : "=f" (__z)            \
-               : "f" (x));             \
-     __z; })
-# endif
-
 #endif /* defined _ARCH_PWR5X */
 
 #endif /* _PPC_MATH_PRIVATE_H_ */
index 3d3b3e7866dc8fec61956bd7198d9ea884f41191..20c828c778a0d485f782e358b87aae713588a701 100644 (file)
@@ -41,7 +41,7 @@ __modf (double x, double *iptr)
     }
   else
     {
-      *iptr = __ceil (x);
+      *iptr = ceil (x);
       return __copysign (x - *iptr, x);
     }
 }
index de44218982835a252c1d0e369c3099985e699b86..807bfe1e03f3409529ab7541ff28535c9a49e3e7 100644 (file)
@@ -40,7 +40,7 @@ __modff (float x, float *iptr)
     }
   else
     {
-      *iptr = __ceilf (x);
+      *iptr = ceilf (x);
       return __copysignf (x - *iptr, x);
     }
 }
index 5f3c0713fbb9a85f9dcc442a551a917326c5fd61..ae0f628759b9dbdc449084e0b71e8cb5827016a8 100644 (file)
@@ -16,6 +16,7 @@
    License along with the GNU C Library; if not, see
    <http://www.gnu.org/licenses/>.  */
 
+#define NO_MATH_REDIRECT
 #include <math.h>
 #include <math_ldbl_opt.h>
 #include <shlib-compat.h>
index 1f10ed65853ea161c82d03a172d055cc6c4c4283..4dcf2b0dfe83f797df178a5fb16e215a4920e193 100644 (file)
@@ -16,6 +16,7 @@
    License along with the GNU C Library; if not, see
    <http://www.gnu.org/licenses/>.  */
 
+#define NO_MATH_REDIRECT
 #include <math.h>
 #include <math_ldbl_opt.h>
 #include <shlib-compat.h>
index 5cde4eb46f26ec82fe82cabe8619cd2a4df4bd58..2ea9bb878e025c134b36be5b356e5d10cc1f6e62 100644 (file)
@@ -16,6 +16,7 @@
    License along with the GNU C Library; if not, see
    <http://www.gnu.org/licenses/>.  */
 
+#define NO_MATH_REDIRECT
 #include <math.h>
 #include <math_ldbl_opt.h>
 #include <shlib-compat.h>
index 18697e52b3ac16166b1d2debb7b31cb1f1fc6148..b7a7a0afda092bc25d6a5eb18a68f3c93bb32eb9 100644 (file)
@@ -16,6 +16,7 @@
    License along with the GNU C Library; if not, see
    <http://www.gnu.org/licenses/>.  */
 
+#define NO_MATH_REDIRECT
 #include <math.h>
 #include <math_ldbl_opt.h>
 #include <shlib-compat.h>
index 30caa968690888fcbd8833bfad1f37ea1965bfbc..697f067722da51a4672c529b7d0460270fcb253d 100644 (file)
@@ -16,6 +16,7 @@
    License along with the GNU C Library.  If not, see
    <http://www.gnu.org/licenses/>.  */
 
+#define NO_MATH_REDIRECT
 #include <math.h>
 #include <math_private.h>
 #include <fenv_private.h>
index a56f3626105f7bc69a5c47409237591b4e6b1032..11ca3aba6e8779a62beb878744bf29d2f54833d5 100644 (file)
@@ -16,6 +16,7 @@
    License along with the GNU C Library.  If not, see
    <http://www.gnu.org/licenses/>.  */
 
+#define NO_MATH_REDIRECT
 #include <math.h>
 #include <math_private.h>
 #include <fenv_private.h>
index 3d34622d43472a36a89831f6f7aa7f7c5068d33d..728bb9265785d0ba4bcf8008d217ba5848bbfaf1 100644 (file)
@@ -16,6 +16,7 @@
    License along with the GNU C Library; if not, see
    <http://www.gnu.org/licenses/>.  */
 
+#define NO_MATH_REDIRECT
 #include <sparc-ifunc.h>
 #include <math.h>
 #include <libm-alias-double.h>
index 18d2eb565d0f313873f266bc1849eaa724055e62..b1149173fe00f94dbfdbc716c32351a272daca86 100644 (file)
@@ -16,6 +16,7 @@
    License along with the GNU C Library; if not, see
    <http://www.gnu.org/licenses/>.  */
 
+#define NO_MATH_REDIRECT
 #include <sparc-ifunc.h>
 #include <math.h>
 #include <libm-alias-float.h>
index 070fcdddeaa186bd6ce635b8288b14b66f5304c3..e846dbeafa58c9e484199a5e11a90c0c78625973 100644 (file)
@@ -16,6 +16,7 @@
    License along with the GNU C Library; if not, see
    <http://www.gnu.org/licenses/>.  */
 
+#define NO_MATH_REDIRECT
 #include <libm-alias-double.h>
 
 #define ceil __redirect_ceil
index db0c6c4bc30711ad9c1609c64cbb873ad825552a..62f3ee7081cde6c514904a8536b0b00c11cb992f 100644 (file)
@@ -16,6 +16,7 @@
    License along with the GNU C Library; if not, see
    <http://www.gnu.org/licenses/>.  */
 
+#define NO_MATH_REDIRECT
 #include <libm-alias-float.h>
 
 #define ceilf __redirect_ceilf