Add implementations for fmin, fminf, and fminl. Testing infrastructure update is splitted to https://reviews.llvm.org/D83931.

4 files changed, 261 insertions, 0 deletions

diff --git a/libc/test/src/math/CMakeLists.txt b/libc/test/src/math/CMakeLists.txt
index 6c6cb9036e53..a6898b370621 100644
--- a/libc/test/src/math/CMakeLists.txt
+++ b/libc/test/src/math/CMakeLists.txt
@@ -437,3 +437,39 @@ add_math_unittest(
     libc.src.math.modfl
     libc.utils.FPUtil.fputil
 )
+
+add_math_unittest(
+  fminf_test
+  SUITE
+    libc_math_unittests
+  SRCS
+    fminf_test.cpp
+  DEPENDS
+    libc.include.math
+    libc.src.math.fminf
+    libc.utils.FPUtil.fputil
+)
+
+add_math_unittest(
+  fmin_test
+  SUITE
+    libc_math_unittests
+  SRCS
+    fmin_test.cpp
+  DEPENDS
+    libc.include.math
+    libc.src.math.fmin
+    libc.utils.FPUtil.fputil
+)
+
+add_math_unittest(
+  fminl_test
+  SUITE
+    libc_math_unittests
+  SRCS
+    fminl_test.cpp
+  DEPENDS
+    libc.include.math
+    libc.src.math.fminl
+    libc.utils.FPUtil.fputil
+)
\ No newline at end of file
diff --git a/libc/test/src/math/fmin_test.cpp b/libc/test/src/math/fmin_test.cpp
new file mode 100644
index 000000000000..c1467e4cbb90
--- /dev/null
+++ b/libc/test/src/math/fmin_test.cpp
@@ -0,0 +1,75 @@
+//===-- Unittests for fmin -----------------------------------------------===//
+//
+// Part of the LLVM Project, under the Apache License v2.0 with LLVM Exceptions.
+// See https://llvm.org/LICENSE.txt for license information.
+// SPDX-License-Identifier: Apache-2.0 WITH LLVM-exception
+//
+//===---------------------------------------------------------------------===//
+
+#include <vector>
+
+#include "include/math.h"
+#include "src/math/fmin.h"
+#include "utils/FPUtil/FPBits.h"
+#include "utils/UnitTest/Test.h"
+
+using FPBits = __llvm_libc::fputil::FPBits<double>;
+
+double nan = static_cast<double>(FPBits::buildNaN(1));
+double inf = static_cast<double>(FPBits::inf());
+double negInf = static_cast<double>(FPBits::negInf());
+
+TEST(FminTest, NaNArg) {
+  EXPECT_EQ(inf, __llvm_libc::fmin(nan, inf));
+  EXPECT_EQ(negInf, __llvm_libc::fmin(negInf, nan));
+  EXPECT_EQ(0.0, __llvm_libc::fmin(nan, 0.0));
+  EXPECT_EQ(-0.0, __llvm_libc::fmin(-0.0, nan));
+  EXPECT_EQ(-1.2345, __llvm_libc::fmin(nan, -1.2345));
+  EXPECT_EQ(1.2345, __llvm_libc::fmin(1.2345, nan));
+  EXPECT_NE(isnan(__llvm_libc::fmin(nan, nan)), 0);
+}
+
+TEST(FminTest, InfArg) {
+  EXPECT_EQ(negInf, __llvm_libc::fmin(negInf, inf));
+  EXPECT_EQ(0.0, __llvm_libc::fmin(inf, 0.0));
+  EXPECT_EQ(-0.0, __llvm_libc::fmin(-0.0, inf));
+  EXPECT_EQ(1.2345, __llvm_libc::fmin(inf, 1.2345));
+  EXPECT_EQ(-1.2345, __llvm_libc::fmin(-1.2345, inf));
+}
+
+TEST(FminTest, NegInfArg) {
+  EXPECT_EQ(negInf, __llvm_libc::fmin(inf, negInf));
+  EXPECT_EQ(negInf, __llvm_libc::fmin(negInf, 0.0));
+  EXPECT_EQ(negInf, __llvm_libc::fmin(-0.0, negInf));
+  EXPECT_EQ(negInf, __llvm_libc::fmin(negInf, -1.2345));
+  EXPECT_EQ(negInf, __llvm_libc::fmin(1.2345, negInf));
+}
+
+TEST(FminTest, BothZero) {
+  EXPECT_EQ(0.0, __llvm_libc::fmin(0.0, 0.0));
+  EXPECT_EQ(-0.0, __llvm_libc::fmin(-0.0, 0.0));
+  EXPECT_EQ(-0.0, __llvm_libc::fmin(0.0, -0.0));
+  EXPECT_EQ(-0.0, __llvm_libc::fmin(-0.0, -0.0));
+}
+
+TEST(FminTest, InFloatRange) {
+  using UIntType = FPBits::UIntType;
+  constexpr UIntType count = 10000000;
+  constexpr UIntType step = UIntType(-1) / count;
+  for (UIntType i = 0, v = 0, w = UIntType(-1); i <= count;
+       ++i, v += step, w -= step) {
+    double x = FPBits(v), y = FPBits(w);
+    if (isnan(x) || isinf(x))
+      continue;
+    if (isnan(y) || isinf(y))
+      continue;
+    if ((x == 0) && (y == 0))
+      continue;
+
+    if (x < y) {
+      ASSERT_EQ(x, __llvm_libc::fmin(x, y));
+    } else {
+      ASSERT_EQ(y, __llvm_libc::fmin(x, y));
+    }
+  }
+}
diff --git a/libc/test/src/math/fminf_test.cpp b/libc/test/src/math/fminf_test.cpp
new file mode 100644
index 000000000000..8d03d92f7237
--- /dev/null
+++ b/libc/test/src/math/fminf_test.cpp
@@ -0,0 +1,75 @@
+//===-- Unittests for fminf ----------------------------------------------===//
+//
+// Part of the LLVM Project, under the Apache License v2.0 with LLVM Exceptions.
+// See https://llvm.org/LICENSE.txt for license information.
+// SPDX-License-Identifier: Apache-2.0 WITH LLVM-exception
+//
+//===---------------------------------------------------------------------===//
+
+#include <vector>
+
+#include "include/math.h"
+#include "src/math/fminf.h"
+#include "utils/FPUtil/FPBits.h"
+#include "utils/UnitTest/Test.h"
+
+using FPBits = __llvm_libc::fputil::FPBits<float>;
+
+float nan = static_cast<float>(FPBits::buildNaN(1));
+float inf = static_cast<float>(FPBits::inf());
+float negInf = static_cast<float>(FPBits::negInf());
+
+TEST(FminfTest, NaNArg) {
+  EXPECT_EQ(inf, __llvm_libc::fminf(nan, inf));
+  EXPECT_EQ(negInf, __llvm_libc::fminf(negInf, nan));
+  EXPECT_EQ(0.0f, __llvm_libc::fminf(nan, 0.0f));
+  EXPECT_EQ(-0.0f, __llvm_libc::fminf(-0.0f, nan));
+  EXPECT_EQ(-1.2345f, __llvm_libc::fminf(nan, -1.2345f));
+  EXPECT_EQ(1.2345f, __llvm_libc::fminf(1.2345f, nan));
+  EXPECT_NE(isnan(__llvm_libc::fminf(nan, nan)), 0);
+}
+
+TEST(FminfTest, InfArg) {
+  EXPECT_EQ(negInf, __llvm_libc::fminf(negInf, inf));
+  EXPECT_EQ(0.0f, __llvm_libc::fminf(inf, 0.0f));
+  EXPECT_EQ(-0.0f, __llvm_libc::fminf(-0.0f, inf));
+  EXPECT_EQ(1.2345f, __llvm_libc::fminf(inf, 1.2345f));
+  EXPECT_EQ(-1.2345f, __llvm_libc::fminf(-1.2345f, inf));
+}
+
+TEST(FminfTest, NegInfArg) {
+  EXPECT_EQ(negInf, __llvm_libc::fminf(inf, negInf));
+  EXPECT_EQ(negInf, __llvm_libc::fminf(negInf, 0.0f));
+  EXPECT_EQ(negInf, __llvm_libc::fminf(-0.0f, negInf));
+  EXPECT_EQ(negInf, __llvm_libc::fminf(negInf, -1.2345f));
+  EXPECT_EQ(negInf, __llvm_libc::fminf(1.2345f, negInf));
+}
+
+TEST(FminfTest, BothZero) {
+  EXPECT_EQ(0.0f, __llvm_libc::fminf(0.0f, 0.0f));
+  EXPECT_EQ(-0.0f, __llvm_libc::fminf(-0.0f, 0.0f));
+  EXPECT_EQ(-0.0f, __llvm_libc::fminf(0.0f, -0.0f));
+  EXPECT_EQ(-0.0f, __llvm_libc::fminf(-0.0f, -0.0f));
+}
+
+TEST(FminfTest, InFloatRange) {
+  using UIntType = FPBits::UIntType;
+  constexpr UIntType count = 10000000;
+  constexpr UIntType step = UIntType(-1) / count;
+  for (UIntType i = 0, v = 0, w = UIntType(-1); i <= count;
+       ++i, v += step, w -= step) {
+    float x = FPBits(v), y = FPBits(w);
+    if (isnan(x) || isinf(x))
+      continue;
+    if (isnan(y) || isinf(y))
+      continue;
+    if ((x == 0) && (y == 0))
+      continue;
+
+    if (x < y) {
+      ASSERT_EQ(x, __llvm_libc::fminf(x, y));
+    } else {
+      ASSERT_EQ(y, __llvm_libc::fminf(x, y));
+    }
+  }
+}
diff --git a/libc/test/src/math/fminl_test.cpp b/libc/test/src/math/fminl_test.cpp
new file mode 100644
index 000000000000..519c1f33fba6
--- /dev/null
+++ b/libc/test/src/math/fminl_test.cpp
@@ -0,0 +1,75 @@
+//===-- Unittests for fmin -----------------------------------------------===//
+//
+// Part of the LLVM Project, under the Apache License v2.0 with LLVM Exceptions.
+// See https://llvm.org/LICENSE.txt for license information.
+// SPDX-License-Identifier: Apache-2.0 WITH LLVM-exception
+//
+//===---------------------------------------------------------------------===//
+
+#include <vector>
+
+#include "include/math.h"
+#include "src/math/fminl.h"
+#include "utils/FPUtil/FPBits.h"
+#include "utils/UnitTest/Test.h"
+
+using FPBits = __llvm_libc::fputil::FPBits<long double>;
+
+long double nan = static_cast<long double>(FPBits::buildNaN(1));
+long double inf = static_cast<long double>(FPBits::inf());
+long double negInf = static_cast<long double>(FPBits::negInf());
+
+TEST(FminlTest, NaNArg) {
+  EXPECT_EQ(inf, __llvm_libc::fminl(nan, inf));
+  EXPECT_EQ(negInf, __llvm_libc::fminl(negInf, nan));
+  EXPECT_EQ(0.0L, __llvm_libc::fminl(nan, 0.0L));
+  EXPECT_EQ(-0.0L, __llvm_libc::fminl(-0.0L, nan));
+  EXPECT_EQ(-1.2345L, __llvm_libc::fminl(nan, -1.2345L));
+  EXPECT_EQ(1.2345L, __llvm_libc::fminl(1.2345L, nan));
+  EXPECT_NE(isnan(__llvm_libc::fminl(nan, nan)), 0);
+}
+
+TEST(FminlTest, InfArg) {
+  EXPECT_EQ(negInf, __llvm_libc::fminl(negInf, inf));
+  EXPECT_EQ(0.0L, __llvm_libc::fminl(inf, 0.0L));
+  EXPECT_EQ(-0.0L, __llvm_libc::fminl(-0.0L, inf));
+  EXPECT_EQ(1.2345L, __llvm_libc::fminl(inf, 1.2345L));
+  EXPECT_EQ(-1.2345L, __llvm_libc::fminl(-1.2345L, inf));
+}
+
+TEST(FminlTest, NegInfArg) {
+  EXPECT_EQ(negInf, __llvm_libc::fminl(inf, negInf));
+  EXPECT_EQ(negInf, __llvm_libc::fminl(negInf, 0.0L));
+  EXPECT_EQ(negInf, __llvm_libc::fminl(-0.0L, negInf));
+  EXPECT_EQ(negInf, __llvm_libc::fminl(negInf, -1.2345L));
+  EXPECT_EQ(negInf, __llvm_libc::fminl(1.2345L, negInf));
+}
+
+TEST(FminlTest, BothZero) {
+  EXPECT_EQ(0.0L, __llvm_libc::fminl(0.0L, 0.0L));
+  EXPECT_EQ(-0.0L, __llvm_libc::fminl(-0.0L, 0.0L));
+  EXPECT_EQ(-0.0L, __llvm_libc::fminl(0.0L, -0.0L));
+  EXPECT_EQ(-0.0L, __llvm_libc::fminl(-0.0L, -0.0L));
+}
+
+TEST(FminlTest, InLongDoubleRange) {
+  using UIntType = FPBits::UIntType;
+  constexpr UIntType count = 10000000;
+  constexpr UIntType step = UIntType(-1) / count;
+  for (UIntType i = 0, v = 0, w = UIntType(-1); i <= count;
+       ++i, v += step, w -= step) {
+    long double x = FPBits(v), y = FPBits(w);
+    if (isnan(x) || isinf(x))
+      continue;
+    if (isnan(y) || isinf(y))
+      continue;
+    if ((x == 0) && (y == 0))
+      continue;
+
+    if (x < y) {
+      ASSERT_EQ(x, __llvm_libc::fminl(x, y));
+    } else {
+      ASSERT_EQ(y, __llvm_libc::fminl(x, y));
+    }
+  }
+}