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