1 //===-- Unittests for frexpl ----------------------------------------------===//
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/frexpl.h"
11 #include "utils/FPUtil/BasicOperations.h"
12 #include "utils/FPUtil/FPBits.h"
13 #include "utils/UnitTest/Test.h"
14
15 using FPBits = __llvm_libc::fputil::FPBits<long double>;
16
TEST(FrexplTest,SpecialNumbers)17 TEST(FrexplTest, SpecialNumbers) {
18 int exponent;
19
20 EXPECT_TRUE(FPBits::inf() == __llvm_libc::frexpl(FPBits::inf(), &exponent));
21 EXPECT_TRUE(FPBits::negInf() ==
22 __llvm_libc::frexpl(FPBits::negInf(), &exponent));
23
24 EXPECT_TRUE(FPBits::zero() == __llvm_libc::frexpl(FPBits::zero(), &exponent));
25 EXPECT_EQ(exponent, 0);
26
27 EXPECT_TRUE(FPBits::negZero() ==
28 __llvm_libc::frexpl(FPBits::negZero(), &exponent));
29 EXPECT_EQ(exponent, 0);
30
31 EXPECT_TRUE(
32 FPBits(__llvm_libc::frexpl(FPBits::buildNaN(1), &exponent)).isNaN());
33 }
34
TEST(FrexplTest,PowersOfTwo)35 TEST(FrexplTest, PowersOfTwo) {
36 int exponent;
37
38 EXPECT_TRUE(0.5l == __llvm_libc::frexpl(1.0l, &exponent));
39 EXPECT_EQ(exponent, 1);
40 EXPECT_TRUE(-0.5l == __llvm_libc::frexpl(-1.0l, &exponent));
41 EXPECT_EQ(exponent, 1);
42
43 EXPECT_TRUE(0.5l == __llvm_libc::frexpl(2.0l, &exponent));
44 EXPECT_EQ(exponent, 2);
45 EXPECT_TRUE(-0.5l == __llvm_libc::frexpl(-2.0l, &exponent));
46 EXPECT_EQ(exponent, 2);
47
48 EXPECT_TRUE(0.5l == __llvm_libc::frexpl(4.0l, &exponent));
49 EXPECT_EQ(exponent, 3);
50 EXPECT_TRUE(-0.5l == __llvm_libc::frexpl(-4.0l, &exponent));
51 EXPECT_EQ(exponent, 3);
52
53 EXPECT_TRUE(0.5l == __llvm_libc::frexpl(8.0l, &exponent));
54 EXPECT_EQ(exponent, 4);
55 EXPECT_TRUE(-0.5l == __llvm_libc::frexpl(-8.0l, &exponent));
56 EXPECT_EQ(exponent, 4);
57
58 EXPECT_TRUE(0.5l == __llvm_libc::frexpl(16.0l, &exponent));
59 EXPECT_EQ(exponent, 5);
60 EXPECT_TRUE(-0.5l == __llvm_libc::frexpl(-16.0l, &exponent));
61 EXPECT_EQ(exponent, 5);
62
63 EXPECT_TRUE(0.5l == __llvm_libc::frexpl(32.0l, &exponent));
64 EXPECT_EQ(exponent, 6);
65 EXPECT_TRUE(-0.5l == __llvm_libc::frexpl(-32.0l, &exponent));
66 EXPECT_EQ(exponent, 6);
67 }
68
TEST(FrexplTest,SomeIntegers)69 TEST(FrexplTest, SomeIntegers) {
70 int exponent;
71
72 EXPECT_TRUE(0.75l == __llvm_libc::frexpl(24.0l, &exponent));
73 EXPECT_EQ(exponent, 5);
74 EXPECT_TRUE(-0.75l == __llvm_libc::frexpl(-24.0l, &exponent));
75 EXPECT_EQ(exponent, 5);
76
77 EXPECT_TRUE(0.625l == __llvm_libc::frexpl(40.0l, &exponent));
78 EXPECT_EQ(exponent, 6);
79 EXPECT_TRUE(-0.625l == __llvm_libc::frexpl(-40.0l, &exponent));
80 EXPECT_EQ(exponent, 6);
81
82 EXPECT_TRUE(0.78125l == __llvm_libc::frexpl(800.0l, &exponent));
83 EXPECT_EQ(exponent, 10);
84 EXPECT_TRUE(-0.78125l == __llvm_libc::frexpl(-800.0l, &exponent));
85 EXPECT_EQ(exponent, 10);
86 }
87
TEST(FrexplTest,LongDoubleRange)88 TEST(FrexplTest, LongDoubleRange) {
89 using UIntType = FPBits::UIntType;
90 constexpr UIntType count = 10000000;
91 constexpr UIntType step = UIntType(-1) / count;
92 for (UIntType i = 0, v = 0; i <= count; ++i, v += step) {
93 long double x = FPBits(v);
94 if (isnan(x) || isinf(x) || x == 0.0l)
95 continue;
96
97 int exponent;
98 long double frac = __llvm_libc::frexpl(x, &exponent);
99
100 ASSERT_TRUE(__llvm_libc::fputil::abs(frac) < 1.0l);
101 ASSERT_TRUE(__llvm_libc::fputil::abs(frac) >= 0.5l);
102 }
103 }
104