[glibc/azanella/mips-hw-fp-round] mips: Implement floor with hardware floating-point rounding instruction

[Adhemerval Zanella <azanella@sourceware.org>]

https://sourceware.org/git/gitweb.cgi?p=glibc.git;h=f5b0332d9200c3ca70e8a29d7a7b6d0519aa00da

commit f5b0332d9200c3ca70e8a29d7a7b6d0519aa00da
Author: Adhemerval Zanella <adhemerval.zanella@linaro.org>
Date:   Tue Dec 26 16:43:54 2023 -0300

    mips: Implement floor with hardware floating-point rounding instruction

Diff:
---
 sysdeps/mips/fpu/round_to_integer.h |  4 ++++
 sysdeps/mips/fpu/s_floor.c          | 36 ++++++++++++++++++++++++++++++++++++
 2 files changed, 40 insertions(+)

diff --git a/sysdeps/mips/fpu/round_to_integer.h b/sysdeps/mips/fpu/round_to_integer.h
index 7ac6783e54..a2146e4a71 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 0000000000..52b60c7ebc
--- /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