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