brintos

brintos / llvm-project-archived public Read only

0
0
Text · 2.7 KiB · 20702c5 Raw
83 lines · cpp
1//===-- Unittests for asinf2//------------------------------------------------===//3//4// Part of the LLVM Project, under the Apache License v2.0 with LLVM Exceptions.5// See https://llvm.org/LICENSE.txt for license information.6// SPDX-License-Identifier: Apache-2.0 WITH LLVM-exception7//8//===----------------------------------------------------------------------===//9 10#include "hdr/errno_macros.h"11#include "hdr/math_macros.h"12#include "hdr/stdint_proxy.h"13#include "src/__support/FPUtil/FPBits.h"14#include "src/__support/macros/optimization.h"15#include "src/math/asinf.h"16#include "test/UnitTest/FPMatcher.h"17#include "test/UnitTest/Test.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 LlvmLibcAsinfTest = LIBC_NAMESPACE::testing::FPTest<float>;27 28namespace mpfr = LIBC_NAMESPACE::testing::mpfr;29 30TEST_F(LlvmLibcAsinfTest, SpecialNumbers) {31  EXPECT_FP_EQ_ALL_ROUNDING(aNaN, LIBC_NAMESPACE::asinf(aNaN));32  EXPECT_MATH_ERRNO(0);33 34  EXPECT_FP_EQ_ALL_ROUNDING(0.0f, LIBC_NAMESPACE::asinf(0.0f));35  EXPECT_MATH_ERRNO(0);36 37  EXPECT_FP_EQ_ALL_ROUNDING(-0.0f, LIBC_NAMESPACE::asinf(-0.0f));38  EXPECT_MATH_ERRNO(0);39 40  EXPECT_FP_EQ_ALL_ROUNDING(aNaN, LIBC_NAMESPACE::asinf(inf));41  EXPECT_MATH_ERRNO(EDOM);42 43  EXPECT_FP_EQ_ALL_ROUNDING(aNaN, LIBC_NAMESPACE::asinf(neg_inf));44  EXPECT_MATH_ERRNO(EDOM);45}46 47TEST_F(LlvmLibcAsinfTest, InFloatRange) {48  constexpr uint32_t COUNT = 100'000;49  constexpr uint32_t STEP = UINT32_MAX / COUNT;50  for (uint32_t i = 0, v = 0; i <= COUNT; ++i, v += STEP) {51    float x = FPBits(v).get_val();52    if (FPBits(v).is_nan() || FPBits(v).is_inf())53      continue;54    ASSERT_MPFR_MATCH_ALL_ROUNDING(mpfr::Operation::Asin, x,55                                   LIBC_NAMESPACE::asinf(x), 0.5);56  }57}58 59TEST_F(LlvmLibcAsinfTest, SpecificBitPatterns) {60  constexpr int N = 11;61  constexpr uint32_t INPUTS[N] = {62      0x3f000000, // x = 0.5f63      0x3f3504f3, // x = sqrt(2)/2, FE_DOWNWARD64      0x3f3504f4, // x = sqrt(2)/2, FE_UPWARD65      0x3f5db3d7, // x = sqrt(3)/2, FE_DOWNWARD66      0x3f5db3d8, // x = sqrt(3)/2, FE_UPWARD67      0x3f800000, // x = 1.0f68      0x40000000, // x = 2.0f69      0x3d09bf86, // x = 0x1.137f0cp-5f70      0x3de5fa1e, // x = 0x1.cbf43cp-4f71      0x3f083a1a, // x = 0x1.107434p-1f72      0x3f7741b6, // x = 0x1.ee836cp-1f73  };74 75  for (int i = 0; i < N; ++i) {76    float x = FPBits(INPUTS[i]).get_val();77    EXPECT_MPFR_MATCH_ALL_ROUNDING(mpfr::Operation::Asin, x,78                                   LIBC_NAMESPACE::asinf(x), TOLERANCE + 0.5);79    EXPECT_MPFR_MATCH_ALL_ROUNDING(mpfr::Operation::Asin, -x,80                                   LIBC_NAMESPACE::asinf(-x), TOLERANCE + 0.5);81  }82}83