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