aboutsummaryrefslogtreecommitdiff
diff options
context:
space:
mode:
authorAdhemerval Zanella <adhemerval.zanella@linaro.org>2023-12-26 16:43:54 -0300
committerAdhemerval Zanella <adhemerval.zanella@linaro.org>2023-12-27 10:22:36 -0300
commit28e3562609428955f34960abe83d530a116d32cf (patch)
tree10ac5734f50fe941e7ee7ddf4778e905e895484e
parent96af11c2077c95860317169303c3ce25bcc0f43f (diff)
downloadglibc-28e3562609428955f34960abe83d530a116d32cf.zip
glibc-28e3562609428955f34960abe83d530a116d32cf.tar.gz
glibc-28e3562609428955f34960abe83d530a116d32cf.tar.bz2
mips: Implement floor with hardware floating-point rounding instruction
-rw-r--r--sysdeps/mips/fpu/round_to_integer.h4
-rw-r--r--sysdeps/mips/fpu/s_floor.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 7ac6783..a2146e4 100644
--- a/sysdeps/mips/fpu/round_to_integer.h
+++ b/sysdeps/mips/fpu/round_to_integer.h
@@ -25,6 +25,7 @@
enum round_mode
{
CEIL,
+ FLOOR,
};
static inline double
@@ -42,6 +43,9 @@ round_to_integer_double (enum round_mode mode, double x)
case CEIL:
asm ("ceil.l.d %0, %0" : "+f" (r));
break;
+ case FLOOR:
+ asm ("floor.l.d %0, %0" : "+f" (r));
+ break;
}
libc_fesetenv (&fe);
diff --git a/sysdeps/mips/fpu/s_floor.c b/sysdeps/mips/fpu/s_floor.c
new file mode 100644
index 0000000..52b60c7
--- /dev/null
+++ b/sysdeps/mips/fpu/s_floor.c
@@ -0,0 +1,36 @@
+/* Largest integral value not greater than argument. 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
+__floor (double x)
+{
+ return round_to_integer_double (FLOOR, x);
+}
+libm_alias_double (__floor, floor)
+#else
+# include <sysdeps/ieee754/dbl-64/s_floor.c>
+#endif