1 //===-- Unittests for tanf ------------------------------------------------===//
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 "hdr/math_macros.h"
10 #include "src/__support/FPUtil/FPBits.h"
11 #include "src/errno/libc_errno.h"
12 #include "src/math/tanf.h"
13 #include "test/UnitTest/FPMatcher.h"
14 #include "test/UnitTest/Test.h"
15 #include "test/src/math/sdcomp26094.h"
16 #include "utils/MPFRWrapper/MPFRUtils.h"
17
18 #include <stdint.h>
19
20 using LlvmLibcTanfTest = LIBC_NAMESPACE::testing::FPTest<float>;
21
22 using LIBC_NAMESPACE::testing::SDCOMP26094_VALUES;
23
24 namespace mpfr = LIBC_NAMESPACE::testing::mpfr;
25
TEST_F(LlvmLibcTanfTest,SpecialNumbers)26 TEST_F(LlvmLibcTanfTest, SpecialNumbers) {
27 LIBC_NAMESPACE::libc_errno = 0;
28
29 EXPECT_FP_EQ(aNaN, LIBC_NAMESPACE::tanf(aNaN));
30 EXPECT_MATH_ERRNO(0);
31
32 EXPECT_FP_EQ(0.0f, LIBC_NAMESPACE::tanf(0.0f));
33 EXPECT_MATH_ERRNO(0);
34
35 EXPECT_FP_EQ(-0.0f, LIBC_NAMESPACE::tanf(-0.0f));
36 EXPECT_MATH_ERRNO(0);
37
38 EXPECT_FP_EQ(aNaN, LIBC_NAMESPACE::tanf(inf));
39 EXPECT_MATH_ERRNO(EDOM);
40
41 EXPECT_FP_EQ(aNaN, LIBC_NAMESPACE::tanf(neg_inf));
42 EXPECT_MATH_ERRNO(EDOM);
43 }
44
TEST_F(LlvmLibcTanfTest,InFloatRange)45 TEST_F(LlvmLibcTanfTest, InFloatRange) {
46 constexpr uint32_t COUNT = 100'000;
47 constexpr uint32_t STEP = UINT32_MAX / COUNT;
48 for (uint32_t i = 0, v = 0; i <= COUNT; ++i, v += STEP) {
49 float x = FPBits(v).get_val();
50 if (FPBits(v).is_nan() || FPBits(v).is_inf())
51 continue;
52 ASSERT_MPFR_MATCH_ALL_ROUNDING(mpfr::Operation::Tan, x,
53 LIBC_NAMESPACE::tanf(x), 0.5);
54 }
55 }
56
TEST_F(LlvmLibcTanfTest,SpecificBitPatterns)57 TEST_F(LlvmLibcTanfTest, SpecificBitPatterns) {
58 constexpr int N = 54;
59 constexpr uint32_t INPUTS[N] = {
60 0x3a7a'8d2fU, // x = 0x1.f51a5ep-11f
61 0x3f06'0a92U, // x = pi/6
62 0x3f3a'dc51U, // x = 0x1.75b8a2p-1f
63 0x3f49'0fdbU, // x = pi/4
64 0x3f86'0a92U, // x = pi/3
65 0x3f8a'1f62U, // x = 0x1.143ec4p+0f
66 0x3fa7'832aU, // x = 0x1.4f0654p+0f
67 0x3fc9'0fdbU, // x = pi/2
68 0x4017'1973U, // x = 0x1.2e32e6p+1f
69 0x4049'0fdbU, // x = pi
70 0x4096'cbe4U, // x = 0x1.2d97c8p+2f
71 0x40c9'0fdbU, // x = 2*pi
72 0x433b'7490U, // x = 0x1.76e92p+7f
73 0x437c'e5f1U, // x = 0x1.f9cbe2p+7f
74 0x4619'9998U, // x = 0x1.33333p+13f
75 0x474d'246fU, // x = 0x1.9a48dep+15f
76 0x4afd'ece4U, // x = 0x1.fbd9c8p+22f
77 0x4c23'32e9U, // x = 0x1.4665d2p+25f
78 0x4d56'd355U, // x = 0x1.ada6aap+27f
79 0x5043'1032U, // x = 0x1.862064p+33f
80 0x50a3'e87fU, // x = 0x1.47d0fep+34f
81 0x5239'47f6U, // x = 0x1.728fecp+37f
82 0x531d'744cU, // x = 0x1.3ae898p+39f
83 0x53b1'46a6U, // x = 0x1.628d4cp+40f
84 0x5532'5019U, // x = 0x1.64a032p+43f
85 0x55ca'fb2aU, // x = 0x1.95f654p+44f
86 0x57d7'b0edU, // x = 0x1.af61dap+48f
87 0x588e'f060U, // x = 0x1.1de0cp+50f
88 0x5922'aa80U, // x = 0x1.4555p+51f
89 0x5980'445eU, // x = 0x1.0088bcp+52f
90 0x5aa4'542cU, // x = 0x1.48a858p+54f
91 0x5c07'bcd0U, // x = 0x1.0f79ap+57f
92 0x5ebc'fddeU, // x = 0x1.79fbbcp+62f
93 0x5f18'b878U, // x = 0x1.3170fp+63f
94 0x5fa6'eba7U, // x = 0x1.4dd74ep+64f
95 0x6115'cb11U, // x = 0x1.2b9622p+67f
96 0x61a4'0b40U, // x = 0x1.48168p+68f
97 0x6386'134eU, // x = 0x1.0c269cp+72f
98 0x63fc'86feU, // x = 0x1.f90dfcp+72f
99 0x6589'8498U, // x = 0x1.13093p+76f
100 0x65ee'8695U, // x = 0x1.dd0d2ap+76f
101 0x6600'0001U, // x = 0x1.000002p+77f
102 0x664e'46e4U, // x = 0x1.9c8dc8p+77f
103 0x66b0'14aaU, // x = 0x1.602954p+78f
104 0x6798'fe4fU, // x = 0x1.31fc9ep+80f
105 0x67a9'242bU, // x = 0x1.524856p+80f
106 0x6a19'76f1U, // x = 0x1.32ede2p+85f
107 0x6ad3'6709U, // x = 0x1.a6ce12p+86f
108 0x6c55'da58U, // x = 0x1.abb4bp+89f
109 0x6f79'be45U, // x = 0x1.f37c8ap+95f
110 0x7276'69d4U, // x = 0x1.ecd3a8p+101f
111 0x72b5'05bbU, // x = 0x1.6a0b76p+102f
112 0x7758'4625U, // x = 0x1.b08c4ap+111f
113 0x7bee'f5efU, // x = 0x1.ddebdep+120f
114 };
115
116 for (int i = 0; i < N; ++i) {
117 float x = FPBits(INPUTS[i]).get_val();
118 EXPECT_MPFR_MATCH_ALL_ROUNDING(mpfr::Operation::Tan, x,
119 LIBC_NAMESPACE::tanf(x), 0.5);
120 EXPECT_MPFR_MATCH_ALL_ROUNDING(mpfr::Operation::Tan, -x,
121 LIBC_NAMESPACE::tanf(-x), 0.5);
122 }
123 }
124
125 // SDCOMP-26094: check tanf in the cases for which the range reducer
126 // returns values furthest beyond its nominal upper bound of pi/4.
TEST_F(LlvmLibcTanfTest,SDCOMP_26094)127 TEST_F(LlvmLibcTanfTest, SDCOMP_26094) {
128 for (uint32_t v : SDCOMP26094_VALUES) {
129 float x = FPBits(v).get_val();
130 ASSERT_MPFR_MATCH(mpfr::Operation::Tan, x, LIBC_NAMESPACE::tanf(x), 0.5);
131 }
132 }
133