about summary refs log tree commit diff
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.gz
glibc-78768e9002314aa427986fe6fb7aad14df877c5e.tar.xz
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