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/__support/libc_errno.h"
13#include "src/math/expm1f16.h"
14#include "test/UnitTest/FPMatcher.h"
15#include "test/UnitTest/Test.h"
16
17using LlvmLibcExpm1f16Test = LIBC_NAMESPACE::testing::FPTest<float16>;
18
19TEST_F(LlvmLibcExpm1f16Test, SpecialNumbers) {
20 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
42TEST_F(LlvmLibcExpm1f16Test, Overflow) {
43 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
69TEST_F(LlvmLibcExpm1f16Test, ResultNearNegOne) {
70 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

source code of libc/test/src/math/smoke/expm1f16_test.cpp