181 lines · cpp
1//===-- Unittests for sincosf ---------------------------------------------===//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/__support/macros/optimization.h"14#include "src/math/sincosf.h"15#include "test/UnitTest/FPMatcher.h"16#include "test/UnitTest/Test.h"17#include "test/src/math/sdcomp26094.h"18#include "utils/MPFRWrapper/MPFRUtils.h"19 20#ifdef LIBC_MATH_HAS_SKIP_ACCURATE_PASS21#define TOLERANCE 122#else23#define TOLERANCE 024#endif // LIBC_MATH_HAS_SKIP_ACCURATE_PASS25 26using LlvmLibcSinCosfTest = LIBC_NAMESPACE::testing::FPTest<float>;27 28using LIBC_NAMESPACE::testing::SDCOMP26094_VALUES;29 30namespace mpfr = LIBC_NAMESPACE::testing::mpfr;31 32TEST_F(LlvmLibcSinCosfTest, SpecialNumbers) {33 float sin, cos;34 35 LIBC_NAMESPACE::sincosf(aNaN, &sin, &cos);36 EXPECT_FP_EQ(aNaN, cos);37 EXPECT_FP_EQ(aNaN, sin);38 EXPECT_MATH_ERRNO(0);39 40 LIBC_NAMESPACE::sincosf(0.0f, &sin, &cos);41 EXPECT_FP_EQ(1.0f, cos);42 EXPECT_FP_EQ(0.0f, sin);43 EXPECT_MATH_ERRNO(0);44 45 LIBC_NAMESPACE::sincosf(-0.0f, &sin, &cos);46 EXPECT_FP_EQ(1.0f, cos);47 EXPECT_FP_EQ(-0.0f, sin);48 EXPECT_MATH_ERRNO(0);49 50 LIBC_NAMESPACE::sincosf(inf, &sin, &cos);51 EXPECT_FP_EQ(aNaN, cos);52 EXPECT_FP_EQ(aNaN, sin);53 EXPECT_MATH_ERRNO(EDOM);54 55 LIBC_NAMESPACE::sincosf(neg_inf, &sin, &cos);56 EXPECT_FP_EQ(aNaN, cos);57 EXPECT_FP_EQ(aNaN, sin);58 EXPECT_MATH_ERRNO(EDOM);59}60 61#define EXPECT_SINCOS_MATCH_ALL_ROUNDING(input) \62 { \63 float sin, cos; \64 namespace mpfr = LIBC_NAMESPACE::testing::mpfr; \65 \66 mpfr::ForceRoundingMode __r1(mpfr::RoundingMode::Nearest); \67 if (__r1.success) { \68 LIBC_NAMESPACE::sincosf(input, &sin, &cos); \69 EXPECT_MPFR_MATCH(mpfr::Operation::Sin, input, sin, TOLERANCE + 0.5, \70 mpfr::RoundingMode::Nearest); \71 EXPECT_MPFR_MATCH(mpfr::Operation::Cos, input, cos, TOLERANCE + 0.5, \72 mpfr::RoundingMode::Nearest); \73 } \74 \75 mpfr::ForceRoundingMode __r2(mpfr::RoundingMode::Upward); \76 if (__r2.success) { \77 LIBC_NAMESPACE::sincosf(input, &sin, &cos); \78 EXPECT_MPFR_MATCH(mpfr::Operation::Sin, input, sin, TOLERANCE + 0.5, \79 mpfr::RoundingMode::Upward); \80 EXPECT_MPFR_MATCH(mpfr::Operation::Cos, input, cos, TOLERANCE + 0.5, \81 mpfr::RoundingMode::Upward); \82 } \83 \84 mpfr::ForceRoundingMode __r3(mpfr::RoundingMode::Downward); \85 if (__r3.success) { \86 LIBC_NAMESPACE::sincosf(input, &sin, &cos); \87 EXPECT_MPFR_MATCH(mpfr::Operation::Sin, input, sin, TOLERANCE + 0.5, \88 mpfr::RoundingMode::Downward); \89 EXPECT_MPFR_MATCH(mpfr::Operation::Cos, input, cos, TOLERANCE + 0.5, \90 mpfr::RoundingMode::Downward); \91 } \92 \93 mpfr::ForceRoundingMode __r4(mpfr::RoundingMode::TowardZero); \94 if (__r4.success) { \95 LIBC_NAMESPACE::sincosf(input, &sin, &cos); \96 EXPECT_MPFR_MATCH(mpfr::Operation::Sin, input, sin, TOLERANCE + 0.5, \97 mpfr::RoundingMode::TowardZero); \98 EXPECT_MPFR_MATCH(mpfr::Operation::Cos, input, cos, TOLERANCE + 0.5, \99 mpfr::RoundingMode::TowardZero); \100 } \101 }102 103TEST_F(LlvmLibcSinCosfTest, InFloatRange) {104 constexpr uint32_t COUNT = 1'001;105 constexpr uint32_t STEP = UINT32_MAX / COUNT;106 for (uint32_t i = 0, v = 0; i <= COUNT; ++i, v += STEP) {107 float x = FPBits(v).get_val();108 if (FPBits(v).is_nan() || FPBits(v).is_inf())109 continue;110 111 EXPECT_SINCOS_MATCH_ALL_ROUNDING(x);112 EXPECT_SINCOS_MATCH_ALL_ROUNDING(-x);113 }114}115 116// For hard to round inputs.117TEST_F(LlvmLibcSinCosfTest, SpecialValues) {118 constexpr int N = 43;119 constexpr uint32_t INPUTS[N] = {120 0x3b56'37f5U, // x = 0x1.ac6feap-9f121 0x3f06'0a92U, // x = pi/6122 0x3f3a'dc51U, // x = 0x1.75b8a2p-1f123 0x3f49'0fdbU, // x = pi/4124 0x3f86'0a92U, // x = pi/3125 0x3fa7'832aU, // x = 0x1.4f0654p+0f126 0x3fc9'0fdbU, // x = pi/2127 0x4017'1973U, // x = 0x1.2e32e6p+1f128 0x4049'0fdbU, // x = pi129 0x4096'cbe4U, // x = 0x1.2d97c8p+2f130 0x40c9'0fdbU, // x = 2*pi131 0x433b'7490U, // x = 0x1.76e92p+7f132 0x437c'e5f1U, // x = 0x1.f9cbe2p+7f133 0x4619'9998U, // x = 0x1.33333p+13f134 0x474d'246fU, // x = 0x1.9a48dep+15f135 0x4afd'ece4U, // x = 0x1.fbd9c8p+22f136 0x4c23'32e9U, // x = 0x1.4665d2p+25f137 0x50a3'e87fU, // x = 0x1.47d0fep+34f138 0x5239'47f6U, // x = 0x1.728fecp+37f139 0x53b1'46a6U, // x = 0x1.628d4cp+40f140 0x5532'5019U, // x = 0x1.64a032p+43f141 0x55ca'fb2aU, // x = 0x1.95f654p+44f142 0x588e'f060U, // x = 0x1.1de0cp+50f143 0x5922'aa80U, // x = 0x1.4555p+51f144 0x5aa4'542cU, // x = 0x1.48a858p+54f145 0x5c07'bcd0U, // x = 0x1.0f79ap+57f146 0x5ebc'fddeU, // x = 0x1.79fbbcp+62f147 0x5f18'b878U, // x = 0x1.3170fp+63f148 0x5fa6'eba7U, // x = 0x1.4dd74ep+64f149 0x6115'cb11U, // x = 0x1.2b9622p+67f150 0x61a4'0b40U, // x = 0x1.48168p+68f151 0x6386'134eU, // x = 0x1.0c269cp+72f152 0x6589'8498U, // x = 0x1.13093p+76f153 0x6600'0001U, // x = 0x1.000002p+77f154 0x664e'46e4U, // x = 0x1.9c8dc8p+77f155 0x66b0'14aaU, // x = 0x1.602954p+78f156 0x67a9'242bU, // x = 0x1.524856p+80f157 0x6a19'76f1U, // x = 0x1.32ede2p+85f158 0x6c55'da58U, // x = 0x1.abb4bp+89f159 0x6f79'be45U, // x = 0x1.f37c8ap+95f160 0x7276'69d4U, // x = 0x1.ecd3a8p+101f161 0x7758'4625U, // x = 0x1.b08c4ap+111f162 0x7bee'f5efU, // x = 0x1.ddebdep+120f163 };164 165 for (int i = 0; i < N; ++i) {166 float x = FPBits(INPUTS[i]).get_val();167 EXPECT_SINCOS_MATCH_ALL_ROUNDING(x);168 EXPECT_SINCOS_MATCH_ALL_ROUNDING(-x);169 }170}171 172// SDCOMP-26094: check sincosf in the cases for which the range reducer173// returns values furthest beyond its nominal upper bound of pi/4.174TEST_F(LlvmLibcSinCosfTest, SDCOMP_26094) {175 for (uint32_t v : SDCOMP26094_VALUES) {176 float x = FPBits(v).get_val();177 EXPECT_SINCOS_MATCH_ALL_ROUNDING(x);178 EXPECT_SINCOS_MATCH_ALL_ROUNDING(-x);179 }180}181