1 //===-- Unittests for expm1f16 --------------------------------------------===//
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/errno_macros.h"
10 #include "hdr/fenv_macros.h"
11 #include "src/__support/FPUtil/cast.h"
12 #include "src/errno/libc_errno.h"
13 #include "src/math/expm1f16.h"
14 #include "test/UnitTest/FPMatcher.h"
15 #include "test/UnitTest/Test.h"
16
17 using LlvmLibcExpm1f16Test = LIBC_NAMESPACE::testing::FPTest<float16>;
18
TEST_F(LlvmLibcExpm1f16Test,SpecialNumbers)19 TEST_F(LlvmLibcExpm1f16Test, SpecialNumbers) {
20 LIBC_NAMESPACE::libc_errno = 0;
21
22 EXPECT_FP_EQ_ALL_ROUNDING(aNaN, LIBC_NAMESPACE::expm1f16(aNaN));
23 EXPECT_MATH_ERRNO(0);
24
25 EXPECT_FP_EQ_WITH_EXCEPTION(aNaN, LIBC_NAMESPACE::expm1f16(sNaN), FE_INVALID);
26 EXPECT_MATH_ERRNO(0);
27
28 EXPECT_FP_EQ_ALL_ROUNDING(inf, LIBC_NAMESPACE::expm1f16(inf));
29 EXPECT_MATH_ERRNO(0);
30
31 EXPECT_FP_EQ_ALL_ROUNDING(LIBC_NAMESPACE::fputil::cast<float16>(-1.0),
32 LIBC_NAMESPACE::expm1f16(neg_inf));
33 EXPECT_MATH_ERRNO(0);
34
35 EXPECT_FP_EQ_ALL_ROUNDING(zero, LIBC_NAMESPACE::expm1f16(zero));
36 EXPECT_MATH_ERRNO(0);
37
38 EXPECT_FP_EQ_ALL_ROUNDING(neg_zero, LIBC_NAMESPACE::expm1f16(neg_zero));
39 EXPECT_MATH_ERRNO(0);
40 }
41
TEST_F(LlvmLibcExpm1f16Test,Overflow)42 TEST_F(LlvmLibcExpm1f16Test, Overflow) {
43 LIBC_NAMESPACE::libc_errno = 0;
44
45 EXPECT_FP_EQ_WITH_EXCEPTION(inf, LIBC_NAMESPACE::expm1f16(max_normal),
46 FE_OVERFLOW | FE_INEXACT);
47 EXPECT_MATH_ERRNO(ERANGE);
48
49 // round(16 * log(2), HP, RN);
50 float16 x = LIBC_NAMESPACE::fputil::cast<float16>(0x1.63p+3);
51
52 EXPECT_FP_EQ_WITH_EXCEPTION_ROUNDING_NEAREST(inf, LIBC_NAMESPACE::expm1f16(x),
53 FE_OVERFLOW | FE_INEXACT);
54 EXPECT_MATH_ERRNO(ERANGE);
55
56 EXPECT_FP_EQ_WITH_EXCEPTION_ROUNDING_UPWARD(inf, LIBC_NAMESPACE::expm1f16(x),
57 FE_OVERFLOW | FE_INEXACT);
58 EXPECT_MATH_ERRNO(ERANGE);
59
60 EXPECT_FP_EQ_WITH_EXCEPTION_ROUNDING_DOWNWARD(
61 max_normal, LIBC_NAMESPACE::expm1f16(x), FE_INEXACT);
62 EXPECT_MATH_ERRNO(0);
63
64 EXPECT_FP_EQ_WITH_EXCEPTION_ROUNDING_TOWARD_ZERO(
65 max_normal, LIBC_NAMESPACE::expm1f16(x), FE_INEXACT);
66 EXPECT_MATH_ERRNO(0);
67 }
68
TEST_F(LlvmLibcExpm1f16Test,ResultNearNegOne)69 TEST_F(LlvmLibcExpm1f16Test, ResultNearNegOne) {
70 LIBC_NAMESPACE::libc_errno = 0;
71
72 EXPECT_FP_EQ_WITH_EXCEPTION(LIBC_NAMESPACE::fputil::cast<float16>(-1.0),
73 LIBC_NAMESPACE::expm1f16(neg_max_normal),
74 FE_INEXACT);
75
76 // round(-11 * log(2), HP, RN);
77 float16 x = LIBC_NAMESPACE::fputil::cast<float16>(-0x1.e8p+2);
78
79 EXPECT_FP_EQ_WITH_EXCEPTION_ROUNDING_NEAREST(
80 LIBC_NAMESPACE::fputil::cast<float16>(-0x1.ffcp-1),
81 LIBC_NAMESPACE::expm1f16(x), FE_INEXACT);
82
83 EXPECT_FP_EQ_WITH_EXCEPTION_ROUNDING_UPWARD(
84 LIBC_NAMESPACE::fputil::cast<float16>(-0x1.ffcp-1),
85 LIBC_NAMESPACE::expm1f16(x), FE_INEXACT);
86
87 EXPECT_FP_EQ_WITH_EXCEPTION_ROUNDING_DOWNWARD(
88 LIBC_NAMESPACE::fputil::cast<float16>(-1.0), LIBC_NAMESPACE::expm1f16(x),
89 FE_INEXACT);
90
91 EXPECT_FP_EQ_WITH_EXCEPTION_ROUNDING_TOWARD_ZERO(
92 LIBC_NAMESPACE::fputil::cast<float16>(-0x1.ffcp-1),
93 LIBC_NAMESPACE::expm1f16(x), FE_INEXACT);
94
95 x = LIBC_NAMESPACE::fputil::cast<float16>(-0x1.0a4p+3);
96
97 EXPECT_FP_EQ_WITH_EXCEPTION_ROUNDING_NEAREST(
98 LIBC_NAMESPACE::fputil::cast<float16>(-1.0), LIBC_NAMESPACE::expm1f16(x),
99 FE_INEXACT);
100
101 EXPECT_FP_EQ_WITH_EXCEPTION_ROUNDING_UPWARD(
102 LIBC_NAMESPACE::fputil::cast<float16>(-0x1.ffcp-1),
103 LIBC_NAMESPACE::expm1f16(x), FE_INEXACT);
104
105 EXPECT_FP_EQ_WITH_EXCEPTION_ROUNDING_DOWNWARD(
106 LIBC_NAMESPACE::fputil::cast<float16>(-1.0), LIBC_NAMESPACE::expm1f16(x),
107 FE_INEXACT);
108
109 EXPECT_FP_EQ_WITH_EXCEPTION_ROUNDING_TOWARD_ZERO(
110 LIBC_NAMESPACE::fputil::cast<float16>(-0x1.ffcp-1),
111 LIBC_NAMESPACE::expm1f16(x), FE_INEXACT);
112 }
113