|
 |
7c0489 |
From cf5397eb5b33bab37c16bcb2d1bbddbce1a27de2 Mon Sep 17 00:00:00 2001
|
|
|
7c0489 |
From: Stefan Liebler <stli@linux.ibm.com>
|
|
|
7c0489 |
Date: Wed, 11 Dec 2019 15:09:33 +0100
|
|
|
7c0489 |
Subject: [PATCH 25/28] S390: Implement roundtoint and converttoint and define
|
|
|
7c0489 |
TOINT_INTRINSICS.
|
|
|
7c0489 |
|
|
|
7c0489 |
This patch implements roundtoint and convertoint for s390
|
|
|
7c0489 |
by using the load-fp-integer and convert-to-fixed instructions.
|
|
|
7c0489 |
Both functions are using "round to nearest with ties away from zero"
|
|
|
7c0489 |
rounding mode and do not raise inexact exceptions.
|
|
|
7c0489 |
|
|
|
7c0489 |
(cherry picked from commit 2763d3145a326aa9afa613fe9e1b444cf912a883)
|
|
|
7c0489 |
---
|
|
|
7c0489 |
sysdeps/s390/fpu/math_private.h | 53 +++++++++++++++++++++++++++++++++
|
|
|
7c0489 |
1 file changed, 53 insertions(+)
|
|
|
7c0489 |
create mode 100644 sysdeps/s390/fpu/math_private.h
|
|
|
7c0489 |
|
|
|
7c0489 |
diff --git a/sysdeps/s390/fpu/math_private.h b/sysdeps/s390/fpu/math_private.h
|
|
|
7c0489 |
new file mode 100644
|
|
|
7c0489 |
index 0000000000..a1ae91a87c
|
|
|
7c0489 |
--- /dev/null
|
|
|
7c0489 |
+++ b/sysdeps/s390/fpu/math_private.h
|
|
|
7c0489 |
@@ -0,0 +1,53 @@
|
|
|
7c0489 |
+/* Configure optimized libm functions. S390 version.
|
|
|
7c0489 |
+ Copyright (C) 2019 Free Software Foundation, Inc.
|
|
|
7c0489 |
+ This file is part of the GNU C Library.
|
|
|
7c0489 |
+
|
|
|
7c0489 |
+ The GNU C Library is free software; you can redistribute it and/or
|
|
|
7c0489 |
+ modify it under the terms of the GNU Lesser General Public
|
|
|
7c0489 |
+ License as published by the Free Software Foundation; either
|
|
|
7c0489 |
+ version 2.1 of the License, or (at your option) any later version.
|
|
|
7c0489 |
+
|
|
|
7c0489 |
+ The GNU C Library is distributed in the hope that it will be useful,
|
|
|
7c0489 |
+ but WITHOUT ANY WARRANTY; without even the implied warranty of
|
|
|
7c0489 |
+ MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
|
|
|
7c0489 |
+ Lesser General Public License for more details.
|
|
|
7c0489 |
+
|
|
|
7c0489 |
+ You should have received a copy of the GNU Lesser General Public
|
|
|
7c0489 |
+ License along with the GNU C Library; if not, see
|
|
|
7c0489 |
+ <https://www.gnu.org/licenses/>. */
|
|
|
7c0489 |
+
|
|
|
7c0489 |
+#ifndef S390_MATH_PRIVATE_H
|
|
|
7c0489 |
+#define S390_MATH_PRIVATE_H 1
|
|
|
7c0489 |
+
|
|
|
7c0489 |
+#include <stdint.h>
|
|
|
7c0489 |
+#include <math.h>
|
|
|
7c0489 |
+
|
|
|
7c0489 |
+#ifdef HAVE_S390_MIN_Z196_ZARCH_ASM_SUPPORT
|
|
|
7c0489 |
+# define TOINT_INTRINSICS 1
|
|
|
7c0489 |
+
|
|
|
7c0489 |
+static inline double_t
|
|
|
7c0489 |
+roundtoint (double_t x)
|
|
|
7c0489 |
+{
|
|
|
7c0489 |
+ double_t y;
|
|
|
7c0489 |
+ /* The z196 zarch "load fp integer" (fidbra) instruction is rounding
|
|
|
7c0489 |
+ x to the nearest integer with ties away from zero (M3-field: 1)
|
|
|
7c0489 |
+ where inexact exceptions are suppressed (M4-field: 4). */
|
|
|
7c0489 |
+ __asm__ ("fidbra %0,1,%1,4" : "=f" (y) : "f" (x));
|
|
|
7c0489 |
+ return y;
|
|
|
7c0489 |
+}
|
|
|
7c0489 |
+
|
|
|
7c0489 |
+static inline int32_t
|
|
|
7c0489 |
+converttoint (double_t x)
|
|
|
7c0489 |
+{
|
|
|
7c0489 |
+ int32_t y;
|
|
|
7c0489 |
+ /* The z196 zarch "convert to fixed" (cfdbra) instruction is rounding
|
|
|
7c0489 |
+ x to the nearest integer with ties away from zero (M3-field: 1)
|
|
|
7c0489 |
+ where inexact exceptions are suppressed (M4-field: 4). */
|
|
|
7c0489 |
+ __asm__ ("cfdbra %0,1,%1,4" : "=d" (y) : "f" (x) : "cc");
|
|
|
7c0489 |
+ return y;
|
|
|
7c0489 |
+}
|
|
|
7c0489 |
+#endif
|
|
|
7c0489 |
+
|
|
|
7c0489 |
+#include_next <math_private.h>
|
|
|
7c0489 |
+
|
|
|
7c0489 |
+#endif
|
|
|
7c0489 |
--
|
|
|
7c0489 |
2.18.2
|
|
|
7c0489 |
|