1 //===-- Unittests for tanhf16 ---------------------------------------------===//
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/fenv_macros.h"
10 #include "src/__support/FPUtil/cast.h"
11 #include "src/errno/libc_errno.h"
12 #include "src/math/tanhf16.h"
13 #include "test/UnitTest/FPMatcher.h"
14 #include "test/UnitTest/Test.h"
15
16 using LlvmLibcTanhf16Test = LIBC_NAMESPACE::testing::FPTest<float16>;
17
TEST_F(LlvmLibcTanhf16Test,SpecialNumbers)18 TEST_F(LlvmLibcTanhf16Test, SpecialNumbers) {
19 LIBC_NAMESPACE::libc_errno = 0;
20
21 EXPECT_FP_EQ_ALL_ROUNDING(aNaN, LIBC_NAMESPACE::tanhf16(aNaN));
22 EXPECT_MATH_ERRNO(0);
23
24 EXPECT_FP_EQ_WITH_EXCEPTION(aNaN, LIBC_NAMESPACE::tanhf16(sNaN), FE_INVALID);
25 EXPECT_MATH_ERRNO(0);
26
27 EXPECT_FP_EQ_ALL_ROUNDING(LIBC_NAMESPACE::fputil::cast<float16>(1.0),
28 LIBC_NAMESPACE::tanhf16(inf));
29 EXPECT_MATH_ERRNO(0);
30
31 EXPECT_FP_EQ_ALL_ROUNDING(LIBC_NAMESPACE::fputil::cast<float16>(-1.0),
32 LIBC_NAMESPACE::tanhf16(neg_inf));
33 EXPECT_MATH_ERRNO(0);
34
35 EXPECT_FP_EQ_ALL_ROUNDING(zero, LIBC_NAMESPACE::tanhf16(zero));
36 EXPECT_MATH_ERRNO(0);
37
38 EXPECT_FP_EQ_ALL_ROUNDING(neg_zero, LIBC_NAMESPACE::tanhf16(neg_zero));
39 EXPECT_MATH_ERRNO(0);
40 }
41
TEST_F(LlvmLibcTanhf16Test,ResultNearBounds)42 TEST_F(LlvmLibcTanhf16Test, ResultNearBounds) {
43 LIBC_NAMESPACE::libc_errno = 0;
44
45 EXPECT_FP_EQ_WITH_EXCEPTION(LIBC_NAMESPACE::fputil::cast<float16>(1.0),
46 LIBC_NAMESPACE::tanhf16(max_normal), FE_INEXACT);
47 EXPECT_MATH_ERRNO(0);
48
49 EXPECT_FP_EQ_WITH_EXCEPTION(LIBC_NAMESPACE::fputil::cast<float16>(-1.0),
50 LIBC_NAMESPACE::tanhf16(neg_max_normal),
51 FE_INEXACT);
52 EXPECT_MATH_ERRNO(0);
53
54 // round(atanh(1 - 2^-11), HP, RU);
55 float16 x = LIBC_NAMESPACE::fputil::cast<float16>(0x1.0a4p+2);
56
57 EXPECT_FP_EQ_WITH_EXCEPTION_ROUNDING_NEAREST(
58 LIBC_NAMESPACE::fputil::cast<float16>(0x1.ffcp-1),
59 LIBC_NAMESPACE::tanhf16(x), FE_INEXACT);
60 EXPECT_MATH_ERRNO(0);
61
62 EXPECT_FP_EQ_WITH_EXCEPTION_ROUNDING_UPWARD(
63 LIBC_NAMESPACE::fputil::cast<float16>(1.0), LIBC_NAMESPACE::tanhf16(x),
64 FE_INEXACT);
65 EXPECT_MATH_ERRNO(0);
66
67 EXPECT_FP_EQ_WITH_EXCEPTION_ROUNDING_DOWNWARD(
68 LIBC_NAMESPACE::fputil::cast<float16>(0x1.ffcp-1),
69 LIBC_NAMESPACE::tanhf16(x), FE_INEXACT);
70 EXPECT_MATH_ERRNO(0);
71
72 EXPECT_FP_EQ_WITH_EXCEPTION_ROUNDING_TOWARD_ZERO(
73 LIBC_NAMESPACE::fputil::cast<float16>(0x1.ffcp-1),
74 LIBC_NAMESPACE::tanhf16(x), FE_INEXACT);
75 EXPECT_MATH_ERRNO(0);
76
77 x = LIBC_NAMESPACE::fputil::cast<float16>(0x1.208p+2);
78
79 EXPECT_FP_EQ_WITH_EXCEPTION_ROUNDING_NEAREST(
80 LIBC_NAMESPACE::fputil::cast<float16>(1.0), LIBC_NAMESPACE::tanhf16(x),
81 FE_INEXACT);
82 EXPECT_MATH_ERRNO(0);
83
84 EXPECT_FP_EQ_WITH_EXCEPTION_ROUNDING_UPWARD(
85 LIBC_NAMESPACE::fputil::cast<float16>(1.0), LIBC_NAMESPACE::tanhf16(x),
86 FE_INEXACT);
87 EXPECT_MATH_ERRNO(0);
88
89 EXPECT_FP_EQ_WITH_EXCEPTION_ROUNDING_DOWNWARD(
90 LIBC_NAMESPACE::fputil::cast<float16>(0x1.ffcp-1),
91 LIBC_NAMESPACE::tanhf16(x), FE_INEXACT);
92 EXPECT_MATH_ERRNO(0);
93
94 EXPECT_FP_EQ_WITH_EXCEPTION_ROUNDING_TOWARD_ZERO(
95 LIBC_NAMESPACE::fputil::cast<float16>(0x1.ffcp-1),
96 LIBC_NAMESPACE::tanhf16(x), FE_INEXACT);
97 EXPECT_MATH_ERRNO(0);
98
99 // round(atanh(-1 + 2^-11), HP, RD);
100 x = LIBC_NAMESPACE::fputil::cast<float16>(-0x1.0a4p+2);
101
102 EXPECT_FP_EQ_WITH_EXCEPTION_ROUNDING_NEAREST(
103 LIBC_NAMESPACE::fputil::cast<float16>(-0x1.ffcp-1),
104 LIBC_NAMESPACE::tanhf16(x), FE_INEXACT);
105 EXPECT_MATH_ERRNO(0);
106
107 EXPECT_FP_EQ_WITH_EXCEPTION_ROUNDING_UPWARD(
108 LIBC_NAMESPACE::fputil::cast<float16>(-0x1.ffcp-1),
109 LIBC_NAMESPACE::tanhf16(x), FE_INEXACT);
110 EXPECT_MATH_ERRNO(0);
111
112 EXPECT_FP_EQ_WITH_EXCEPTION_ROUNDING_DOWNWARD(
113 LIBC_NAMESPACE::fputil::cast<float16>(-1.0), LIBC_NAMESPACE::tanhf16(x),
114 FE_INEXACT);
115 EXPECT_MATH_ERRNO(0);
116
117 EXPECT_FP_EQ_WITH_EXCEPTION_ROUNDING_TOWARD_ZERO(
118 LIBC_NAMESPACE::fputil::cast<float16>(-0x1.ffcp-1),
119 LIBC_NAMESPACE::tanhf16(x), FE_INEXACT);
120 EXPECT_MATH_ERRNO(0);
121
122 x = LIBC_NAMESPACE::fputil::cast<float16>(-0x1.208p+2);
123
124 EXPECT_FP_EQ_WITH_EXCEPTION_ROUNDING_NEAREST(
125 LIBC_NAMESPACE::fputil::cast<float16>(-1.0), LIBC_NAMESPACE::tanhf16(x),
126 FE_INEXACT);
127 EXPECT_MATH_ERRNO(0);
128
129 EXPECT_FP_EQ_WITH_EXCEPTION_ROUNDING_UPWARD(
130 LIBC_NAMESPACE::fputil::cast<float16>(-0x1.ffcp-1),
131 LIBC_NAMESPACE::tanhf16(x), FE_INEXACT);
132 EXPECT_MATH_ERRNO(0);
133
134 EXPECT_FP_EQ_WITH_EXCEPTION_ROUNDING_DOWNWARD(
135 LIBC_NAMESPACE::fputil::cast<float16>(-1.0), LIBC_NAMESPACE::tanhf16(x),
136 FE_INEXACT);
137 EXPECT_MATH_ERRNO(0);
138
139 EXPECT_FP_EQ_WITH_EXCEPTION_ROUNDING_TOWARD_ZERO(
140 LIBC_NAMESPACE::fputil::cast<float16>(-0x1.ffcp-1),
141 LIBC_NAMESPACE::tanhf16(x), FE_INEXACT);
142 EXPECT_MATH_ERRNO(0);
143 }
144