1 | //===-- Unittests for sinhf -----------------------------------------------===// |
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/math_macros.h" |
10 | #include "src/__support/CPP/array.h" |
11 | #include "src/__support/FPUtil/FPBits.h" |
12 | #include "src/__support/libc_errno.h" |
13 | #include "src/math/sinhf.h" |
14 | #include "test/UnitTest/FPMatcher.h" |
15 | #include "test/UnitTest/Test.h" |
16 | #include "utils/MPFRWrapper/MPFRUtils.h" |
17 | |
18 | #include <stdint.h> |
19 | |
20 | using LlvmLibcSinhfTest = LIBC_NAMESPACE::testing::FPTest<float>; |
21 | |
22 | namespace mpfr = LIBC_NAMESPACE::testing::mpfr; |
23 | |
24 | TEST_F(LlvmLibcSinhfTest, SpecialNumbers) { |
25 | libc_errno = 0; |
26 | |
27 | EXPECT_FP_EQ(aNaN, LIBC_NAMESPACE::sinhf(aNaN)); |
28 | EXPECT_MATH_ERRNO(0); |
29 | |
30 | EXPECT_FP_EQ(0.0f, LIBC_NAMESPACE::sinhf(0.0f)); |
31 | EXPECT_MATH_ERRNO(0); |
32 | |
33 | EXPECT_FP_EQ(-0.0f, LIBC_NAMESPACE::sinhf(-0.0f)); |
34 | EXPECT_MATH_ERRNO(0); |
35 | |
36 | EXPECT_FP_EQ(inf, LIBC_NAMESPACE::sinhf(inf)); |
37 | EXPECT_MATH_ERRNO(0); |
38 | |
39 | EXPECT_FP_EQ(neg_inf, LIBC_NAMESPACE::sinhf(neg_inf)); |
40 | EXPECT_MATH_ERRNO(0); |
41 | } |
42 | |
43 | TEST_F(LlvmLibcSinhfTest, InFloatRange) { |
44 | constexpr uint32_t COUNT = 100'000; |
45 | constexpr uint32_t STEP = UINT32_MAX / COUNT; |
46 | for (uint32_t i = 0, v = 0; i <= COUNT; ++i, v += STEP) { |
47 | float x = FPBits(v).get_val(); |
48 | if (FPBits(v).is_nan() || FPBits(v).is_inf()) |
49 | continue; |
50 | ASSERT_MPFR_MATCH(mpfr::Operation::Sinh, x, LIBC_NAMESPACE::sinhf(x), 0.5); |
51 | } |
52 | } |
53 | |
54 | // For small values, sinh(x) is x. |
55 | TEST_F(LlvmLibcSinhfTest, SmallValues) { |
56 | float x = FPBits(uint32_t(0x17800000)).get_val(); |
57 | float result = LIBC_NAMESPACE::sinhf(x); |
58 | EXPECT_MPFR_MATCH(mpfr::Operation::Sinh, x, result, 0.5); |
59 | EXPECT_FP_EQ(x, result); |
60 | |
61 | x = FPBits(uint32_t(0x00400000)).get_val(); |
62 | result = LIBC_NAMESPACE::sinhf(x); |
63 | EXPECT_MPFR_MATCH(mpfr::Operation::Sinh, x, result, 0.5); |
64 | EXPECT_FP_EQ(x, result); |
65 | } |
66 | |
67 | TEST_F(LlvmLibcSinhfTest, Overflow) { |
68 | libc_errno = 0; |
69 | EXPECT_FP_EQ_WITH_EXCEPTION( |
70 | inf, LIBC_NAMESPACE::sinhf(FPBits(0x7f7fffffU).get_val()), FE_OVERFLOW); |
71 | EXPECT_MATH_ERRNO(ERANGE); |
72 | |
73 | EXPECT_FP_EQ_WITH_EXCEPTION( |
74 | inf, LIBC_NAMESPACE::sinhf(FPBits(0x42cffff8U).get_val()), FE_OVERFLOW); |
75 | EXPECT_MATH_ERRNO(ERANGE); |
76 | |
77 | EXPECT_FP_EQ_WITH_EXCEPTION( |
78 | inf, LIBC_NAMESPACE::sinhf(FPBits(0x42d00008U).get_val()), FE_OVERFLOW); |
79 | EXPECT_MATH_ERRNO(ERANGE); |
80 | } |
81 | |
82 | TEST_F(LlvmLibcSinhfTest, ExceptionalValues) { |
83 | float x = FPBits(uint32_t(0x3a12'85ffU)).get_val(); |
84 | EXPECT_MPFR_MATCH_ALL_ROUNDING(mpfr::Operation::Sinh, x, |
85 | LIBC_NAMESPACE::sinhf(x), 0.5); |
86 | |
87 | x = -FPBits(uint32_t(0x3a12'85ffU)).get_val(); |
88 | EXPECT_MPFR_MATCH_ALL_ROUNDING(mpfr::Operation::Sinh, x, |
89 | LIBC_NAMESPACE::sinhf(x), 0.5); |
90 | } |
91 | |