135 lines · cpp
1//===-- Unittests for sinf ------------------------------------------------===//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/FPUtil/FPBits.h"13#include "src/math/sinf.h"14#include "test/UnitTest/FPMatcher.h"15#include "test/UnitTest/Test.h"16#include "test/src/math/sdcomp26094.h"17#include "utils/MPFRWrapper/MPFRUtils.h"18 19#ifdef LIBC_MATH_HAS_SKIP_ACCURATE_PASS20#define TOLERANCE 321#define FP_ASSERT ASSERT_MPFR_MATCH22#else23#define TOLERANCE 024#define FP_ASSERT ASSERT_MPFR_MATCH_ALL_ROUNDING25#endif // LIBC_MATH_HAS_SKIP_ACCURATE_PASS26 27using LlvmLibcSinfTest = LIBC_NAMESPACE::testing::FPTest<float>;28 29using LIBC_NAMESPACE::testing::SDCOMP26094_VALUES;30 31namespace mpfr = LIBC_NAMESPACE::testing::mpfr;32 33TEST_F(LlvmLibcSinfTest, SpecialNumbers) {34 EXPECT_FP_EQ(aNaN, LIBC_NAMESPACE::sinf(aNaN));35 EXPECT_MATH_ERRNO(0);36 37 EXPECT_FP_EQ(0.0f, LIBC_NAMESPACE::sinf(0.0f));38 EXPECT_MATH_ERRNO(0);39 40// When accuracy is not required, signed zeros might not be preserved.41#ifdef LIBC_MATH_HAS_SKIP_ACCURATE_PASS42 EXPECT_TRUE(0.0f == LIBC_NAMESPACE::sinf(-0.0f));43#else44 EXPECT_FP_EQ(-0.0f, LIBC_NAMESPACE::sinf(-0.0f));45 EXPECT_MATH_ERRNO(0);46#endif // LIBC_MATH_HAS_SKIP_ACCURATE_PASS47 48 EXPECT_FP_EQ(aNaN, LIBC_NAMESPACE::sinf(inf));49 EXPECT_MATH_ERRNO(EDOM);50 51 EXPECT_FP_EQ(aNaN, LIBC_NAMESPACE::sinf(neg_inf));52 EXPECT_MATH_ERRNO(EDOM);53}54 55TEST_F(LlvmLibcSinfTest, InFloatRange) {56 constexpr uint32_t COUNT = 100'000;57 constexpr uint32_t STEP = UINT32_MAX / COUNT;58 for (uint32_t i = 0, v = 0; i <= COUNT; ++i, v += STEP) {59 float x = FPBits(v).get_val();60 if (FPBits(v).is_nan() || FPBits(v).is_inf())61 continue;62 FP_ASSERT(mpfr::Operation::Sin, x, LIBC_NAMESPACE::sinf(x),63 TOLERANCE + 0.5);64 }65}66 67TEST_F(LlvmLibcSinfTest, SpecificBitPatterns) {68 constexpr int N = 36;69 constexpr uint32_t INPUTS[N] = {70 0x3f06'0a92U, // x = pi/671 0x3f3a'dc51U, // x = 0x1.75b8a2p-1f72 0x3f49'0fdbU, // x = pi/473 0x3f86'0a92U, // x = pi/374 0x3fa7'832aU, // x = 0x1.4f0654p+0f75 0x3fc9'0fdbU, // x = pi/276 0x4017'1973U, // x = 0x1.2e32e6p+1f77 0x4049'0fdbU, // x = pi78 0x4096'cbe4U, // x = 0x1.2d97c8p+2f79 0x40c9'0fdbU, // x = 2*pi80 0x433b'7490U, // x = 0x1.76e92p+7f81 0x437c'e5f1U, // x = 0x1.f9cbe2p+7f82 0x4619'9998U, // x = 0x1.33333p+13f83 0x474d'246fU, // x = 0x1.9a48dep+15f84 0x4afd'ece4U, // x = 0x1.fbd9c8p+22f85 0x4c23'32e9U, // x = 0x1.4665d2p+25f86 0x50a3'e87fU, // x = 0x1.47d0fep+34f87 0x5239'47f6U, // x = 0x1.728fecp+37f88 0x53b1'46a6U, // x = 0x1.628d4cp+40f89 0x55ca'fb2aU, // x = 0x1.95f654p+44f90 0x588e'f060U, // x = 0x1.1de0cp+50f91 0x5c07'bcd0U, // x = 0x1.0f79ap+57f92 0x5ebc'fddeU, // x = 0x1.79fbbcp+62f93 0x5fa6'eba7U, // x = 0x1.4dd74ep+64f94 0x61a4'0b40U, // x = 0x1.48168p+68f95 0x6386'134eU, // x = 0x1.0c269cp+72f96 0x6589'8498U, // x = 0x1.13093p+76f97 0x6600'0001U, // x = 0x1.000002p+77f98 0x664e'46e4U, // x = 0x1.9c8dc8p+77f99 0x66b0'14aaU, // x = 0x1.602954p+78f100 0x67a9'242bU, // x = 0x1.524856p+80f101 0x6a19'76f1U, // x = 0x1.32ede2p+85f102 0x6c55'da58U, // x = 0x1.abb4bp+89f103 0x6f79'be45U, // x = 0x1.f37c8ap+95f104 0x7276'69d4U, // x = 0x1.ecd3a8p+101f105 0x7758'4625U, // x = 0x1.b08c4ap+111f106 };107 108 for (int i = 0; i < N; ++i) {109 float x = FPBits(INPUTS[i]).get_val();110 FP_ASSERT(mpfr::Operation::Sin, x, LIBC_NAMESPACE::sinf(x),111 TOLERANCE + 0.5);112 FP_ASSERT(mpfr::Operation::Sin, -x, LIBC_NAMESPACE::sinf(-x),113 TOLERANCE + 0.5);114 }115}116 117// For small values, sin(x) is x.118TEST_F(LlvmLibcSinfTest, SmallValues) {119 float x = FPBits(0x1780'0000U).get_val();120 FP_ASSERT(mpfr::Operation::Sin, x, LIBC_NAMESPACE::sinf(x), 0.5);121 122 x = FPBits(0x0040'0000U).get_val();123 FP_ASSERT(mpfr::Operation::Sin, x, LIBC_NAMESPACE::sinf(x), 0.5);124}125 126// SDCOMP-26094: check sinf in the cases for which the range reducer127// returns values furthest beyond its nominal upper bound of pi/4.128TEST_F(LlvmLibcSinfTest, SDCOMP_26094) {129 for (uint32_t v : SDCOMP26094_VALUES) {130 float x = FPBits((v)).get_val();131 FP_ASSERT(mpfr::Operation::Sin, x, LIBC_NAMESPACE::sinf(x),132 TOLERANCE + 0.5);133 }134}135