1 //===-- Unittests for cospif ----------------------------------------------===// 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 "src/errno/libc_errno.h" 10 #include "src/math/cospif.h" 11 #include "test/UnitTest/FPMatcher.h" 12 #include "test/src/math/sdcomp26094.h" 13 #include "utils/MPFRWrapper/MPFRUtils.h" 14 15 using LlvmLibcCospifTest = LIBC_NAMESPACE::testing::FPTest<float>; 16 17 using LIBC_NAMESPACE::testing::SDCOMP26094_VALUES; 18 19 namespace mpfr = LIBC_NAMESPACE::testing::mpfr; 20 21 TEST_F(LlvmLibcCospifTest, SpecialNumbers) { 22 LIBC_NAMESPACE::libc_errno = 0; 23 24 EXPECT_FP_EQ(aNaN, LIBC_NAMESPACE::cospif(aNaN)); 25 EXPECT_MATH_ERRNO(0); 26 27 EXPECT_FP_EQ(1.0f, LIBC_NAMESPACE::cospif(0.0f)); 28 EXPECT_MATH_ERRNO(0); 29 30 EXPECT_FP_EQ(1.0f, LIBC_NAMESPACE::cospif(-0.0f)); 31 EXPECT_MATH_ERRNO(0); 32 33 EXPECT_FP_EQ(aNaN, LIBC_NAMESPACE::cospif(inf)); 34 EXPECT_MATH_ERRNO(EDOM); 35 36 EXPECT_FP_EQ(aNaN, LIBC_NAMESPACE::cospif(neg_inf)); 37 EXPECT_MATH_ERRNO(EDOM); 38 } 39 40 TEST_F(LlvmLibcCospifTest, SpecificBitPatterns) { 41 constexpr int N = 38; 42 constexpr uint32_t INPUTS[N] = { 43 0x3f00'0000U, // x = 0.5 44 0x461d'd600U, // x = 10101.5 45 0x3f06'0a92U, // x = pi/6 46 0x3f3a'dc51U, // x = 0x1.75b8a2p-1f 47 0x3f49'0fdbU, // x = pi/4 48 0x3f86'0a92U, // x = pi/3 49 0x3fa7'832aU, // x = 0x1.4f0654p+0f 50 0x3fc9'0fdbU, // x = pi/2 51 0x4017'1973U, // x = 0x1.2e32e6p+1f 52 0x4049'0fdbU, // x = pi 53 0x4096'cbe4U, // x = 0x1.2d97c8p+2f 54 0x40c9'0fdbU, // x = 2*pi 55 0x433b'7490U, // x = 0x1.76e92p+7f 56 0x437c'e5f1U, // x = 0x1.f9cbe2p+7f 57 0x4619'9998U, // x = 0x1.33333p+13f 58 0x474d'246fU, // x = 0x1.9a48dep+15f 59 0x4afd'ece4U, // x = 0x1.fbd9c8p+22f 60 0x4c23'32e9U, // x = 0x1.4665d2p+25f 61 0x50a3'e87fU, // x = 0x1.47d0fep+34f 62 0x5239'47f6U, // x = 0x1.728fecp+37f 63 0x53b1'46a6U, // x = 0x1.628d4cp+40f 64 0x55ca'fb2aU, // x = 0x1.95f654p+44f 65 0x588e'f060U, // x = 0x1.1de0cp+50f 66 0x5c07'bcd0U, // x = 0x1.0f79ap+57f 67 0x5ebc'fddeU, // x = 0x1.79fbbcp+62f 68 0x5fa6'eba7U, // x = 0x1.4dd74ep+64f 69 0x61a4'0b40U, // x = 0x1.48168p+68f 70 0x6386'134eU, // x = 0x1.0c269cp+72f 71 0x6589'8498U, // x = 0x1.13093p+76f 72 0x6600'0001U, // x = 0x1.000002p+77f 73 0x664e'46e4U, // x = 0x1.9c8dc8p+77f 74 0x66b0'14aaU, // x = 0x1.602954p+78f 75 0x67a9'242bU, // x = 0x1.524856p+80f 76 0x6a19'76f1U, // x = 0x1.32ede2p+85f 77 0x6c55'da58U, // x = 0x1.abb4bp+89f 78 0x6f79'be45U, // x = 0x1.f37c8ap+95f 79 0x7276'69d4U, // x = 0x1.ecd3a8p+101f 80 0x7758'4625U, // x = 0x1.b08c4ap+111f 81 }; 82 83 for (int i = 0; i < N; ++i) { 84 float x = FPBits(INPUTS[i]).get_val(); 85 EXPECT_MPFR_MATCH_ALL_ROUNDING(mpfr::Operation::Cospi, x, 86 LIBC_NAMESPACE::cospif(x), 0.5); 87 EXPECT_MPFR_MATCH_ALL_ROUNDING(mpfr::Operation::Cospi, -x, 88 LIBC_NAMESPACE::cospif(-x), 0.5); 89 } 90 } 91 92 // For small values, sinpi(x) is pi * x. 93 TEST_F(LlvmLibcCospifTest, SmallValues) { 94 float x = FPBits(0x1780'0000U).get_val(); 95 EXPECT_MPFR_MATCH_ALL_ROUNDING(mpfr::Operation::Cospi, x, 96 LIBC_NAMESPACE::cospif(x), 0.5); 97 98 x = FPBits(0x0040'0000U).get_val(); 99 EXPECT_MPFR_MATCH_ALL_ROUNDING(mpfr::Operation::Cospi, x, 100 LIBC_NAMESPACE::cospif(x), 0.5); 101 } 102 103 // SDCOMP-26094: check sinfpi in the cases for which the range reducer 104 // returns values furthest beyond its nominal upper bound of pi/4. 105 TEST_F(LlvmLibcCospifTest, SDCOMP_26094) { 106 for (uint32_t v : SDCOMP26094_VALUES) { 107 float x = FPBits((v)).get_val(); 108 EXPECT_MPFR_MATCH_ALL_ROUNDING(mpfr::Operation::Cospi, x, 109 LIBC_NAMESPACE::cospif(x), 0.5); 110 } 111 } 112 113 // sinpi(+n + 0.5) = sinpi(-n + 0.5) = +0.0 114 TEST_F(LlvmLibcCospifTest, SignedZeros) { 115 EXPECT_FP_EQ(0.0, LIBC_NAMESPACE::cospif(100.5f)); 116 EXPECT_FP_EQ(0.0, LIBC_NAMESPACE::cospif(-100.5f)); 117 EXPECT_FP_EQ(0.0, LIBC_NAMESPACE::cospif(45678.5f)); 118 EXPECT_FP_EQ(0.0, LIBC_NAMESPACE::cospif(-45678.5f)); 119 EXPECT_FP_EQ(0.0, LIBC_NAMESPACE::cospif(8000000.5f)); 120 EXPECT_FP_EQ(0.0, LIBC_NAMESPACE::cospif(-8000000.5f)); 121 } 122