81 lines · cpp
1//===-- Single-precision sinh function ------------------------------------===//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 "src/math/sinhf.h"10#include "src/__support/FPUtil/FEnvImpl.h"11#include "src/__support/FPUtil/FPBits.h"12#include "src/__support/FPUtil/rounding_mode.h"13#include "src/__support/macros/config.h"14#include "src/__support/macros/optimization.h" // LIBC_UNLIKELY15#include "src/__support/math/sinhfcoshf_utils.h"16 17namespace LIBC_NAMESPACE_DECL {18 19LLVM_LIBC_FUNCTION(float, sinhf, (float x)) {20 using FPBits = typename fputil::FPBits<float>;21 FPBits xbits(x);22 uint32_t x_abs = xbits.abs().uintval();23 24 // When |x| >= 90, or x is inf or nan25 if (LIBC_UNLIKELY(x_abs >= 0x42b4'0000U || x_abs <= 0x3da0'0000U)) {26 // |x| <= 0.07812527 if (x_abs <= 0x3da0'0000U) {28#ifndef LIBC_MATH_HAS_SKIP_ACCURATE_PASS29 // |x| = 0.000558942498173564672470092773437530 if (LIBC_UNLIKELY(x_abs == 0x3a12'85ffU)) {31 if (fputil::fenv_is_round_to_nearest())32 return x;33 }34#endif // !LIBC_MATH_HAS_SKIP_ACCURATE_PASS35 36 // |x| <= 2^-2637 if (LIBC_UNLIKELY(x_abs <= 0x3280'0000U)) {38 return static_cast<float>(39 LIBC_UNLIKELY(x_abs == 0) ? x : (x + 0.25 * x * x * x));40 }41 42 double xdbl = x;43 double x2 = xdbl * xdbl;44 // Sollya: fpminimax(sinh(x),[|3,5,7|],[|D...|],[-1/16-1/64;1/16+1/64],x);45 // Sollya output: x * (0x1p0 + x^0x1p1 * (0x1.5555555556583p-3 + x^0x1p146 // * (0x1.111110d239f1fp-747 // + x^0x1p1 * 0x1.a02b5a284013cp-13)))48 // Therefore, output of Sollya = x * pe;49 double pe = fputil::polyeval(x2, 0.0, 0x1.5555555556583p-3,50 0x1.111110d239f1fp-7, 0x1.a02b5a284013cp-13);51 return static_cast<float>(fputil::multiply_add(xdbl, pe, xdbl));52 }53 54 if (xbits.is_nan())55 return x + 1.0f; // sNaN to qNaN + signal56 57 if (xbits.is_inf())58 return x;59 60 int rounding = fputil::quick_get_round();61 if (xbits.is_neg()) {62 if (LIBC_UNLIKELY(rounding == FE_UPWARD || rounding == FE_TOWARDZERO))63 return -FPBits::max_normal().get_val();64 } else {65 if (LIBC_UNLIKELY(rounding == FE_DOWNWARD || rounding == FE_TOWARDZERO))66 return FPBits::max_normal().get_val();67 }68 69 fputil::set_errno_if_required(ERANGE);70 fputil::raise_except_if_required(FE_OVERFLOW);71 72 return x + FPBits::inf(xbits.sign()).get_val();73 }74 75 // sinh(x) = (e^x - e^(-x)) / 2.76 return static_cast<float>(77 math::sinhfcoshf_internal::exp_pm_eval</*is_sinh*/ true>(x));78}79 80} // namespace LIBC_NAMESPACE_DECL81