1 //===-- Unittests for ceill -----------------------------------------------===//
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 "include/math.h"
10 #include "src/math/ceill.h"
11 #include "utils/FPUtil/FPBits.h"
12 #include "utils/MPFRWrapper/MPFRUtils.h"
13 #include "utils/UnitTest/Test.h"
14 
15 using FPBits = __llvm_libc::fputil::FPBits<long double>;
16 
17 namespace mpfr = __llvm_libc::testing::mpfr;
18 
19 // Zero tolerance; As in, exact match with MPFR result.
20 static constexpr mpfr::Tolerance tolerance{mpfr::Tolerance::doublePrecision, 0,
21                                            0};
22 TEST(CeillTest, SpecialNumbers) {
23   ASSERT_TRUE(FPBits::zero() == __llvm_libc::ceill(FPBits::zero()));
24   ASSERT_TRUE(FPBits::negZero() == __llvm_libc::ceill(FPBits::negZero()));
25 
26   ASSERT_TRUE(FPBits::inf() == __llvm_libc::ceill(FPBits::inf()));
27   ASSERT_TRUE(FPBits::negInf() == __llvm_libc::ceill(FPBits::negInf()));
28 
29   long double nan = FPBits::buildNaN(1);
30   ASSERT_NE(isnan(nan), 0);
31   ASSERT_NE(isnan(__llvm_libc::ceill(nan)), 0);
32 }
33 
34 TEST(CeillTest, RoundedNumbers) {
35   ASSERT_TRUE(FPBits(1.0l) == __llvm_libc::ceill(1.0l));
36   ASSERT_TRUE(FPBits(-1.0l) == __llvm_libc::ceill(-1.0l));
37   ASSERT_TRUE(FPBits(10.0l) == __llvm_libc::ceill(10.0l));
38   ASSERT_TRUE(FPBits(-10.0l) == __llvm_libc::ceill(-10.0l));
39   ASSERT_TRUE(FPBits(1234.0l) == __llvm_libc::ceill(1234.0l));
40   ASSERT_TRUE(FPBits(-1234.0l) == __llvm_libc::ceill(-1234.0l));
41 }
42 
43 TEST(CeillTest, Fractions) {
44   ASSERT_TRUE(FPBits(2.0l) == __llvm_libc::ceill(1.3l));
45   ASSERT_TRUE(FPBits(-1.0l) == __llvm_libc::ceill(-1.3l));
46   ASSERT_TRUE(FPBits(2.0l) == __llvm_libc::ceill(1.5l));
47   ASSERT_TRUE(FPBits(-1.0l) == __llvm_libc::ceill(-1.5l));
48   ASSERT_TRUE(FPBits(2.0l) == __llvm_libc::ceill(1.75l));
49   ASSERT_TRUE(FPBits(-1.0l) == __llvm_libc::ceill(-1.75l));
50   ASSERT_TRUE(FPBits(11.0l) == __llvm_libc::ceill(10.32l));
51   ASSERT_TRUE(FPBits(-10.0l) == __llvm_libc::ceill(-10.32l));
52   ASSERT_TRUE(FPBits(11.0l) == __llvm_libc::ceill(10.65l));
53   ASSERT_TRUE(FPBits(-10.0l) == __llvm_libc::ceill(-10.65l));
54   ASSERT_TRUE(FPBits(1235.0l) == __llvm_libc::ceill(1234.18l));
55   ASSERT_TRUE(FPBits(-1234.0l) == __llvm_libc::ceill(-1234.18l));
56   ASSERT_TRUE(FPBits(1235.0l) == __llvm_libc::ceill(1234.96l));
57   ASSERT_TRUE(FPBits(-1234.0l) == __llvm_libc::ceill(-1234.96l));
58 }
59 
60 TEST(CeillTest, InLongDoubleRange) {
61   using UIntType = FPBits::UIntType;
62   constexpr UIntType count = 10000000;
63   constexpr UIntType step = UIntType(-1) / count;
64   for (UIntType i = 0, v = 0; i <= count; ++i, v += step) {
65     long double x = FPBits(v);
66     if (isnan(x) || isinf(x))
67       continue;
68 
69     ASSERT_MPFR_MATCH(mpfr::Operation::Ceil, x, __llvm_libc::ceill(x),
70                       tolerance);
71   }
72 }
73