1 //===-- Utility class to test floor[f|l] ------------------------*- C++ -*-===//
2 //
3 // Part of the LLVM Project, under the Apache License v2.0 with LLVM Exceptions.
4 // See https://llvm.org/LICENSE.txt for license information.
5 // SPDX-License-Identifier: Apache-2.0 WITH LLVM-exception
6 //
7 //===----------------------------------------------------------------------===//
8 
9 #include "src/__support/FPUtil/BasicOperations.h"
10 #include "src/__support/FPUtil/NearestIntegerOperations.h"
11 #include "utils/MPFRWrapper/MPFRUtils.h"
12 #include "utils/UnitTest/FPMatcher.h"
13 #include "utils/UnitTest/Test.h"
14 
15 #include <math.h>
16 
17 namespace mpfr = __llvm_libc::testing::mpfr;
18 
19 template <typename T> class ModfTest : public __llvm_libc::testing::Test {
20 
21   DECLARE_SPECIAL_CONSTANTS(T)
22 
23 public:
24   typedef T (*ModfFunc)(T, T *);
25 
testSpecialNumbers(ModfFunc func)26   void testSpecialNumbers(ModfFunc func) {
27     T integral;
28 
29     EXPECT_FP_EQ(zero, func(zero, &integral));
30     EXPECT_FP_EQ(integral, zero);
31     EXPECT_FP_EQ(neg_zero, func(neg_zero, &integral));
32     EXPECT_FP_EQ(integral, neg_zero);
33 
34     EXPECT_FP_EQ(zero, func(inf, &integral));
35     EXPECT_FP_EQ(inf, integral);
36     EXPECT_FP_EQ(neg_zero, func(neg_inf, &integral));
37     EXPECT_FP_EQ(neg_inf, integral);
38 
39     EXPECT_FP_EQ(aNaN, func(aNaN, &integral));
40   }
41 
testIntegers(ModfFunc func)42   void testIntegers(ModfFunc func) {
43     T integral;
44 
45     EXPECT_FP_EQ(T(0.0), func(T(1.0), &integral));
46     EXPECT_FP_EQ(T(1.0), integral);
47 
48     EXPECT_FP_EQ(T(-0.0), func(T(-1.0), &integral));
49     EXPECT_FP_EQ(T(-1.0), integral);
50 
51     EXPECT_FP_EQ(T(0.0), func(T(10.0), &integral));
52     EXPECT_FP_EQ(T(10.0), integral);
53 
54     EXPECT_FP_EQ(T(-0.0), func(T(-10.0), &integral));
55     EXPECT_FP_EQ(T(-10.0), integral);
56 
57     EXPECT_FP_EQ(T(0.0), func(T(12345.0), &integral));
58     EXPECT_FP_EQ(T(12345.0), integral);
59 
60     EXPECT_FP_EQ(T(-0.0), func(T(-12345.0), &integral));
61     EXPECT_FP_EQ(T(-12345.0), integral);
62   }
63 
testFractions(ModfFunc func)64   void testFractions(ModfFunc func) {
65     T integral;
66 
67     EXPECT_FP_EQ(T(0.5), func(T(1.5), &integral));
68     EXPECT_FP_EQ(integral, T(1.0));
69 
70     EXPECT_FP_EQ(T(-0.5), func(T(-1.5), &integral));
71     EXPECT_FP_EQ(integral, T(-1.0));
72 
73     EXPECT_FP_EQ(T(0.75), func(T(10.75), &integral));
74     EXPECT_FP_EQ(integral, T(10.0));
75 
76     EXPECT_FP_EQ(T(-0.75), func(T(-10.75), &integral));
77     EXPECT_FP_EQ(integral, T(-10.0));
78 
79     EXPECT_FP_EQ(T(0.125), func(T(100.125), &integral));
80     EXPECT_FP_EQ(integral, T(100.0));
81 
82     EXPECT_FP_EQ(T(-0.125), func(T(-100.125), &integral));
83     EXPECT_FP_EQ(integral, T(-100.0));
84   }
85 
testRange(ModfFunc func)86   void testRange(ModfFunc func) {
87     constexpr UIntType COUNT = 10000000;
88     constexpr UIntType STEP = UIntType(-1) / COUNT;
89     for (UIntType i = 0, v = 0; i <= COUNT; ++i, v += STEP) {
90       T x = T(FPBits(v));
91       if (isnan(x) || isinf(x) || x == T(0.0))
92         continue;
93 
94       T integral;
95       T frac = func(x, &integral);
96       ASSERT_TRUE(__llvm_libc::fputil::abs(frac) < 1.0l);
97       ASSERT_TRUE(__llvm_libc::fputil::trunc(x) == integral);
98       ASSERT_TRUE(integral + frac == x);
99     }
100   }
101 };
102 
103 #define LIST_MODF_TESTS(T, func)                                               \
104   using LlvmLibcModfTest = ModfTest<T>;                                        \
105   TEST_F(LlvmLibcModfTest, SpecialNumbers) { testSpecialNumbers(&func); }      \
106   TEST_F(LlvmLibcModfTest, RoundedNubmers) { testIntegers(&func); }            \
107   TEST_F(LlvmLibcModfTest, Fractions) { testFractions(&func); }                \
108   TEST_F(LlvmLibcModfTest, Range) { testRange(&func); }
109