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/errno/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
TEST_F(LlvmLibcSinhfTest,SpecialNumbers)24 TEST_F(LlvmLibcSinhfTest, SpecialNumbers) {
25 LIBC_NAMESPACE::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
TEST_F(LlvmLibcSinhfTest,InFloatRange)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.
TEST_F(LlvmLibcSinhfTest,SmallValues)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
TEST_F(LlvmLibcSinhfTest,Overflow)67 TEST_F(LlvmLibcSinhfTest, Overflow) {
68 LIBC_NAMESPACE::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
TEST_F(LlvmLibcSinhfTest,ExceptionalValues)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