This is the mail archive of the glibc-cvs@sourceware.org mailing list for the glibc project.


Index Nav: [Date Index] [Subject Index] [Author Index] [Thread Index]
Message Nav: [Date Prev] [Date Next] [Thread Prev] [Thread Next]
Other format: [Raw text]

GNU C Library master sources branch master updated. glibc-2.28.9000-43-g43cfdf8


This is an automated email from the git hooks/post-receive script. It was
generated because a ref change was pushed to the repository containing
the project "GNU C Library master sources".

The branch, master has been updated
       via  43cfdf8f4811c63bca8a2185d82a6d92977a125b (commit)
      from  690652882b499defb3d950dfeff8fe421d13cab5 (commit)

Those revisions listed above that are new to this repository have
not appeared on any other notification email; so we list those
revisions in full, below.

- Log -----------------------------------------------------------------
http://sourceware.org/git/gitweb.cgi?p=glibc.git;a=commitdiff;h=43cfdf8f4811c63bca8a2185d82a6d92977a125b

commit 43cfdf8f4811c63bca8a2185d82a6d92977a125b
Author: Szabolcs Nagy <szabolcs.nagy@arm.com>
Date:   Wed Jul 4 12:29:29 2018 +0100

    Clean up converttoint handling and document the semantics
    
    This patch currently only affects aarch64.
    
    The roundtoint and converttoint internal functions are only called with small
    values, so 32 bit result is enough for converttoint and it is a signed int
    conversion so the return type is changed to int32_t.
    
    The original idea was to help the compiler keeping the result in uint64_t,
    then it's clear that no sign extension is needed and there is no accidental
    undefined or implementation defined signed int arithmetics.
    
    But it turns out gcc does a good job with inlining so changing the type has
    no overhead and the semantics of the conversion is less surprising this way.
    Since we want to allow the asuint64 (x + 0x1.8p52) style conversion, the top
    bits were never usable and the existing code ensures that only the bottom
    32 bits of the conversion result are used.
    
    On aarch64 the neon intrinsics (which round ties to even) are changed to
    round and lround (which round ties away from zero) this does not affect the
    results in a significant way, but more portable (relies on round and lround
    being inlined which works with -fno-math-errno).
    
    The TOINT_SHIFT and TOINT_RINT macros were removed, only keep separate code
    paths for TOINT_INTRINSICS and !TOINT_INTRINSICS.
    
    	* sysdeps/aarch64/fpu/math_private.h (roundtoint): Use round.
    	(converttoint): Use lround.
    	* sysdeps/ieee754/flt-32/math_config.h (roundtoint): Declare and
    	document the semantics when TOINT_INTRINSICS is set.
    	(converttoint): Likewise.
    	(TOINT_RINT): Remove.
    	(TOINT_SHIFT): Remove.
    	* sysdeps/ieee754/flt-32/e_expf.c (__expf): Remove the TOINT_RINT code
    	path.

diff --git a/ChangeLog b/ChangeLog
index 48fd7eb..54c4b5a 100644
--- a/ChangeLog
+++ b/ChangeLog
@@ -1,3 +1,16 @@
+2018-08-10  Wilco Dijkstra  <wdijkstr@arm.com>
+	    Szabolcs Nagy  <szabolcs.nagy@arm.com>
+
+	* sysdeps/aarch64/fpu/math_private.h (roundtoint): Use round.
+	(converttoint): Use lround.
+	* sysdeps/ieee754/flt-32/math_config.h (roundtoint): Declare and
+	document the semantics when TOINT_INTRINSICS is set.
+	(converttoint): Likewise.
+	(TOINT_RINT): Remove.
+	(TOINT_SHIFT): Remove.
+	* sysdeps/ieee754/flt-32/e_expf.c (__expf): Remove the TOINT_RINT code
+	path.
+
 2018-08-10  Florian Weimer  <fweimer@redhat.com>
 
 	[BZ #23497]
diff --git a/sysdeps/aarch64/fpu/math_private.h b/sysdeps/aarch64/fpu/math_private.h
index fcd02c0..d2e0abc 100644
--- a/sysdeps/aarch64/fpu/math_private.h
+++ b/sysdeps/aarch64/fpu/math_private.h
@@ -21,6 +21,8 @@
 
 #include <fenv.h>
 #include <fpu_control.h>
+#include <stdint.h>
+#include <math.h>
 
 static __always_inline void
 libc_feholdexcept_aarch64 (fenv_t *envp)
@@ -298,25 +300,20 @@ libc_feresetround_noex_aarch64_ctx (struct rm_ctx *ctx)
 #define libc_feresetround_noexf_ctx	libc_feresetround_noex_aarch64_ctx
 #define libc_feresetround_noexl_ctx	libc_feresetround_noex_aarch64_ctx
 
-/* Hack: only include the large arm_neon.h when needed.  */
-#ifdef _MATH_CONFIG_H
-# include <arm_neon.h>
-
-/* ACLE intrinsics for frintn and fcvtns instructions.  */
-# define TOINT_INTRINSICS 1
+/* Use inline round and lround instructions.  */
+#define TOINT_INTRINSICS 1
 
 static inline double_t
 roundtoint (double_t x)
 {
-  return vget_lane_f64 (vrndn_f64 (vld1_f64 (&x)), 0);
+  return round (x);
 }
 
-static inline uint64_t
+static inline int32_t
 converttoint (double_t x)
 {
-  return vcvtnd_s64_f64 (x);
+  return lround (x);
 }
-#endif
 
 #include_next <math_private.h>
 
diff --git a/sysdeps/ieee754/flt-32/e_expf.c b/sysdeps/ieee754/flt-32/e_expf.c
index f2238bf..384a586 100644
--- a/sysdeps/ieee754/flt-32/e_expf.c
+++ b/sysdeps/ieee754/flt-32/e_expf.c
@@ -85,10 +85,7 @@ __expf (float x)
 #if TOINT_INTRINSICS
   kd = roundtoint (z);
   ki = converttoint (z);
-#elif TOINT_RINT
-  kd = rint (z);
-  ki = (long) kd;
-#elif TOINT_SHIFT
+#else
 # define SHIFT __exp2f_data.shift
   kd = math_narrow_eval ((double) (z + SHIFT)); /* Needs to be double.  */
   ki = asuint64 (kd);
diff --git a/sysdeps/ieee754/flt-32/math_config.h b/sysdeps/ieee754/flt-32/math_config.h
index 9c4ef30..8ca7532 100644
--- a/sysdeps/ieee754/flt-32/math_config.h
+++ b/sysdeps/ieee754/flt-32/math_config.h
@@ -38,13 +38,23 @@
 #endif
 
 #ifndef TOINT_INTRINSICS
+/* When set, the roundtoint and converttoint functions are provided with
+   the semantics documented below.  */
 # define TOINT_INTRINSICS 0
 #endif
-#ifndef TOINT_RINT
-# define TOINT_RINT 0
-#endif
-#ifndef TOINT_SHIFT
-# define TOINT_SHIFT 1
+
+#if TOINT_INTRINSICS
+/* Round x to nearest int in all rounding modes, ties have to be rounded
+   consistently with converttoint so the results match.  If the result
+   would be outside of [-2^31, 2^31-1] then the semantics is unspecified.  */
+static inline double_t
+roundtoint (double_t x);
+
+/* Convert x to nearest int in all rounding modes, ties have to be rounded
+   consistently with roundtoint.  If the result is not representible in an
+   int32_t then the semantics is unspecified.  */
+static inline int32_t
+converttoint (double_t x);
 #endif
 
 static inline uint32_t

-----------------------------------------------------------------------

Summary of changes:
 ChangeLog                            |   13 +++++++++++++
 sysdeps/aarch64/fpu/math_private.h   |   17 +++++++----------
 sysdeps/ieee754/flt-32/e_expf.c      |    5 +----
 sysdeps/ieee754/flt-32/math_config.h |   20 +++++++++++++++-----
 4 files changed, 36 insertions(+), 19 deletions(-)


hooks/post-receive
-- 
GNU C Library master sources


Index Nav: [Date Index] [Subject Index] [Author Index] [Thread Index]
Message Nav: [Date Prev] [Date Next] [Thread Prev] [Thread Next]