aboutsummaryrefslogtreecommitdiff
diff options
context:
space:
mode:
authorAdhemerval Zanella <adhemerval.zanella@linaro.org>2023-12-26 17:03:07 -0300
committerAdhemerval Zanella <adhemerval.zanella@linaro.org>2023-12-27 10:22:36 -0300
commit78768e9002314aa427986fe6fb7aad14df877c5e (patch)
treea12b13a617d840f89e8f06e40aff382d70553cb4
parentb046c66f0874fa75d1bc859851f7d2e28481905f (diff)
downloadglibc-78768e9002314aa427986fe6fb7aad14df877c5e.tar
glibc-78768e9002314aa427986fe6fb7aad14df877c5e.tar.gz
glibc-78768e9002314aa427986fe6fb7aad14df877c5e.tar.bz2
glibc-78768e9002314aa427986fe6fb7aad14df877c5e.zip
mips: Implement roundeven with hardware floating-point rounding instruction
-rw-r--r--sysdeps/mips/fpu/round_to_integer.h4
-rw-r--r--sysdeps/mips/fpu/s_roundeven.c36
2 files changed, 40 insertions, 0 deletions
diff --git a/sysdeps/mips/fpu/round_to_integer.h b/sysdeps/mips/fpu/round_to_integer.h
index 78975ce3eb..b2ba56b458 100644
--- a/sysdeps/mips/fpu/round_to_integer.h
+++ b/sysdeps/mips/fpu/round_to_integer.h
@@ -27,6 +27,7 @@ enum round_mode
CEIL,
FLOOR,
TRUNC,
+ ROUNDEVEN,
};
static inline double
@@ -50,6 +51,9 @@ round_to_integer_double (enum round_mode mode, double x)
case TRUNC:
asm ("trunc.l.d %0, %0" : "+f" (r));
break;
+ case ROUNDEVEN:
+ asm ("round.l.d %0, %0" : "+f" (r));
+ break;
}
libc_fesetenv (&fe);
diff --git a/sysdeps/mips/fpu/s_roundeven.c b/sysdeps/mips/fpu/s_roundeven.c
new file mode 100644
index 0000000000..55cb315bf1
--- /dev/null
+++ b/sysdeps/mips/fpu/s_roundeven.c
@@ -0,0 +1,36 @@
+/* Round towards the nearest even number. MIPS version.
+ Copyright (C) 2023 Free Software Foundation, Inc.
+ This file is part of the GNU C Library.
+
+ 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
+ <https://www.gnu.org/licenses/>. */
+
+#if ((__mips_fpr == 64) \
+ && (__mips_hard_float == 1) \
+ && ((__mips == 32 && __mips_isa_rev > 1) || __mips == 64))
+
+#define NO_MATH_REDIRECT
+#include <math.h>
+#include <round_to_integer.h>
+#include <libm-alias-double.h>
+
+double
+__roundeven (double x)
+{
+ return round_to_integer_double (ROUNDEVEN, x);
+}
+libm_alias_double (__roundeven, roundeven)
+#else
+# include <sysdeps/ieee754/dbl-64/s_roundeven.c>
+#endif