125 lines · cpp
1//===-- Unittests for sinpif ----------------------------------------------===//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 "src/math/sinpif.h"11#include "test/UnitTest/FPMatcher.h"12#include "test/src/math/sdcomp26094.h"13#include "utils/MPFRWrapper/MPFRUtils.h"14 15#include "hdr/stdint_proxy.h"16 17using LlvmLibcSinpifTest = LIBC_NAMESPACE::testing::FPTest<float>;18 19using LIBC_NAMESPACE::testing::SDCOMP26094_VALUES;20 21namespace mpfr = LIBC_NAMESPACE::testing::mpfr;22 23TEST_F(LlvmLibcSinpifTest, SpecialNumbers) {24 EXPECT_FP_EQ(aNaN, LIBC_NAMESPACE::sinpif(aNaN));25 EXPECT_MATH_ERRNO(0);26 27 EXPECT_FP_EQ(0.0f, LIBC_NAMESPACE::sinpif(0.0f));28 EXPECT_MATH_ERRNO(0);29 30 EXPECT_FP_EQ(-0.0f, LIBC_NAMESPACE::sinpif(-0.0f));31 EXPECT_MATH_ERRNO(0);32 33 EXPECT_FP_EQ(aNaN, LIBC_NAMESPACE::sinpif(inf));34 EXPECT_MATH_ERRNO(EDOM);35 36 EXPECT_FP_EQ(aNaN, LIBC_NAMESPACE::sinpif(neg_inf));37 EXPECT_MATH_ERRNO(EDOM);38}39 40TEST_F(LlvmLibcSinpifTest, SpecificBitPatterns) {41 constexpr int N = 36;42 constexpr uint32_t INPUTS[N] = {43 0x3f06'0a92U, // x = pi/644 0x3f3a'dc51U, // x = 0x1.75b8a2p-1f45 0x3f49'0fdbU, // x = pi/446 0x3f86'0a92U, // x = pi/347 0x3fa7'832aU, // x = 0x1.4f0654p+0f48 0x3fc9'0fdbU, // x = pi/249 0x4017'1973U, // x = 0x1.2e32e6p+1f50 0x4049'0fdbU, // x = pi51 0x4096'cbe4U, // x = 0x1.2d97c8p+2f52 0x40c9'0fdbU, // x = 2*pi53 0x433b'7490U, // x = 0x1.76e92p+7f54 0x437c'e5f1U, // x = 0x1.f9cbe2p+7f55 0x4619'9998U, // x = 0x1.33333p+13f56 0x474d'246fU, // x = 0x1.9a48dep+15f57 0x4afd'ece4U, // x = 0x1.fbd9c8p+22f58 0x4c23'32e9U, // x = 0x1.4665d2p+25f59 0x50a3'e87fU, // x = 0x1.47d0fep+34f60 0x5239'47f6U, // x = 0x1.728fecp+37f61 0x53b1'46a6U, // x = 0x1.628d4cp+40f62 0x55ca'fb2aU, // x = 0x1.95f654p+44f63 0x588e'f060U, // x = 0x1.1de0cp+50f64 0x5c07'bcd0U, // x = 0x1.0f79ap+57f65 0x5ebc'fddeU, // x = 0x1.79fbbcp+62f66 0x5fa6'eba7U, // x = 0x1.4dd74ep+64f67 0x61a4'0b40U, // x = 0x1.48168p+68f68 0x6386'134eU, // x = 0x1.0c269cp+72f69 0x6589'8498U, // x = 0x1.13093p+76f70 0x6600'0001U, // x = 0x1.000002p+77f71 0x664e'46e4U, // x = 0x1.9c8dc8p+77f72 0x66b0'14aaU, // x = 0x1.602954p+78f73 0x67a9'242bU, // x = 0x1.524856p+80f74 0x6a19'76f1U, // x = 0x1.32ede2p+85f75 0x6c55'da58U, // x = 0x1.abb4bp+89f76 0x6f79'be45U, // x = 0x1.f37c8ap+95f77 0x7276'69d4U, // x = 0x1.ecd3a8p+101f78 0x7758'4625U, // x = 0x1.b08c4ap+111f79 };80 81 for (int i = 0; i < N; ++i) {82 float x = FPBits(INPUTS[i]).get_val();83 EXPECT_MPFR_MATCH_ALL_ROUNDING(mpfr::Operation::Sinpi, x,84 LIBC_NAMESPACE::sinpif(x), 0.5);85 EXPECT_MPFR_MATCH_ALL_ROUNDING(mpfr::Operation::Sinpi, -x,86 LIBC_NAMESPACE::sinpif(-x), 0.5);87 }88}89 90// For small values, sinpi(x) is pi * x.91TEST_F(LlvmLibcSinpifTest, SmallValues) {92 float x = FPBits(0x1780'0000U).get_val();93 EXPECT_MPFR_MATCH_ALL_ROUNDING(mpfr::Operation::Sinpi, x,94 LIBC_NAMESPACE::sinpif(x), 0.5);95 96 x = FPBits(0x0040'0000U).get_val();97 EXPECT_MPFR_MATCH_ALL_ROUNDING(mpfr::Operation::Sinpi, x,98 LIBC_NAMESPACE::sinpif(x), 0.5);99}100 101// SDCOMP-26094: check sinpif in the cases for which the range reducer102// returns values furthest beyond its nominal upper bound of pi/4.103TEST_F(LlvmLibcSinpifTest, SDCOMP_26094) {104 for (uint32_t v : SDCOMP26094_VALUES) {105 float x = FPBits((v)).get_val();106 EXPECT_MPFR_MATCH_ALL_ROUNDING(mpfr::Operation::Sinpi, x,107 LIBC_NAMESPACE::sinpif(x), 0.5);108 }109}110 111// sinpi(-n) = -0.0112// sinpi(+n) = +0.0113TEST_F(LlvmLibcSinpifTest, SignedZeros) {114 EXPECT_FP_EQ(-0.0, LIBC_NAMESPACE::sinpif(-0x99));115 EXPECT_FP_EQ(-0.0, LIBC_NAMESPACE::sinpif(-0x1p+43));116 EXPECT_FP_EQ(-0.0, LIBC_NAMESPACE::sinpif(-0x1.4p+64));117 EXPECT_FP_EQ(-0.0, LIBC_NAMESPACE::sinpif(-0x1.cp+106));118 EXPECT_FP_EQ(-0.0, LIBC_NAMESPACE::sinpif(-0x1.cp+21));119 EXPECT_FP_EQ(-0.0, LIBC_NAMESPACE::sinpif(-0x9999));120 EXPECT_FP_EQ(0.0, LIBC_NAMESPACE::sinpif(0x1p+43));121 EXPECT_FP_EQ(0.0, LIBC_NAMESPACE::sinpif(0x1.4p+64));122 EXPECT_FP_EQ(0.0, LIBC_NAMESPACE::sinpif(0x1.cp+106));123 EXPECT_FP_EQ(0.0, LIBC_NAMESPACE::sinpif(0x1.cp+21));124}125