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-exception 6 // 7 //===----------------------------------------------------------------------===// 8 9 #include "include/errno.h" 10 #include "include/math.h" 11 #include "src/errno/llvmlibc_errno.h" 12 #include "src/math/sinf.h" 13 #include "test/src/math/sdcomp26094.h" 14 #include "utils/CPP/Array.h" 15 #include "utils/FPUtil/BitPatterns.h" 16 #include "utils/FPUtil/FloatOperations.h" 17 #include "utils/FPUtil/FloatProperties.h" 18 #include "utils/MPFRWrapper/MPFRUtils.h" 19 #include "utils/UnitTest/Test.h" 20 21 #include <stdint.h> 22 23 using __llvm_libc::fputil::isNegativeQuietNaN; 24 using __llvm_libc::fputil::isQuietNaN; 25 using __llvm_libc::fputil::valueAsBits; 26 using __llvm_libc::fputil::valueFromBits; 27 28 using BitPatterns = __llvm_libc::fputil::BitPatterns<float>; 29 30 using __llvm_libc::testing::sdcomp26094Values; 31 32 namespace mpfr = __llvm_libc::testing::mpfr; 33 34 // 12 additional bits of precision over the base precision of a |float| 35 // value. 36 static constexpr mpfr::Tolerance tolerance{mpfr::Tolerance::floatPrecision, 12, 37 3 * 0x1000 / 4}; 38 39 TEST(SinfTest, SpecialNumbers) { 40 llvmlibc_errno = 0; 41 42 EXPECT_TRUE( 43 isQuietNaN(__llvm_libc::sinf(valueFromBits(BitPatterns::aQuietNaN)))); 44 EXPECT_EQ(llvmlibc_errno, 0); 45 46 EXPECT_TRUE(isNegativeQuietNaN( 47 __llvm_libc::sinf(valueFromBits(BitPatterns::aNegativeQuietNaN)))); 48 EXPECT_EQ(llvmlibc_errno, 0); 49 50 EXPECT_TRUE(isQuietNaN( 51 __llvm_libc::sinf(valueFromBits(BitPatterns::aSignallingNaN)))); 52 EXPECT_EQ(llvmlibc_errno, 0); 53 54 EXPECT_TRUE(isNegativeQuietNaN( 55 __llvm_libc::sinf(valueFromBits(BitPatterns::aNegativeSignallingNaN)))); 56 EXPECT_EQ(llvmlibc_errno, 0); 57 58 EXPECT_EQ(BitPatterns::zero, 59 valueAsBits(__llvm_libc::sinf(valueFromBits(BitPatterns::zero)))); 60 EXPECT_EQ(llvmlibc_errno, 0); 61 62 EXPECT_EQ(BitPatterns::negZero, valueAsBits(__llvm_libc::sinf( 63 valueFromBits(BitPatterns::negZero)))); 64 EXPECT_EQ(llvmlibc_errno, 0); 65 66 llvmlibc_errno = 0; 67 EXPECT_TRUE(isQuietNaN(__llvm_libc::sinf(valueFromBits(BitPatterns::inf)))); 68 EXPECT_EQ(llvmlibc_errno, EDOM); 69 70 llvmlibc_errno = 0; 71 EXPECT_TRUE(isNegativeQuietNaN( 72 __llvm_libc::sinf(valueFromBits(BitPatterns::negInf)))); 73 EXPECT_EQ(llvmlibc_errno, EDOM); 74 } 75 76 TEST(SinfTest, InFloatRange) { 77 constexpr uint32_t count = 1000000; 78 constexpr uint32_t step = UINT32_MAX / count; 79 for (uint32_t i = 0, v = 0; i <= count; ++i, v += step) { 80 float x = valueFromBits(v); 81 if (isnan(x) || isinf(x)) 82 continue; 83 ASSERT_MPFR_MATCH(mpfr::OP_Sin, x, __llvm_libc::sinf(x), tolerance); 84 } 85 } 86 87 TEST(SinfTest, SpecificBitPatterns) { 88 float x = valueFromBits(0xc70d39a1); 89 EXPECT_MPFR_MATCH(mpfr::OP_Sin, x, __llvm_libc::sinf(x), tolerance); 90 } 91 92 // For small values, sin(x) is x. 93 TEST(SinfTest, SmallValues) { 94 uint32_t bits = 0x17800000; 95 float x = valueFromBits(bits); 96 float result = __llvm_libc::sinf(x); 97 EXPECT_MPFR_MATCH(mpfr::OP_Sin, x, result, tolerance); 98 EXPECT_EQ(bits, valueAsBits(result)); 99 100 bits = 0x00400000; 101 x = valueFromBits(bits); 102 result = __llvm_libc::sinf(x); 103 EXPECT_MPFR_MATCH(mpfr::OP_Sin, x, result, tolerance); 104 EXPECT_EQ(bits, valueAsBits(result)); 105 } 106 107 // SDCOMP-26094: check sinf in the cases for which the range reducer 108 // returns values furthest beyond its nominal upper bound of pi/4. 109 TEST(SinfTest, SDCOMP_26094) { 110 for (uint32_t v : sdcomp26094Values) { 111 float x = valueFromBits(v); 112 EXPECT_MPFR_MATCH(mpfr::OP_Sin, x, __llvm_libc::sinf(x), tolerance); 113 } 114 } 115