1 //===-- Unittests for acosf -----------------------------------------------===// 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 "hdr/math_macros.h" 10 #include "src/__support/FPUtil/FPBits.h" 11 #include "src/errno/libc_errno.h" 12 #include "src/math/acosf.h" 13 #include "test/UnitTest/FPMatcher.h" 14 #include "test/UnitTest/Test.h" 15 #include "utils/MPFRWrapper/MPFRUtils.h" 16 17 #include <stdint.h> 18 19 namespace mpfr = LIBC_NAMESPACE::testing::mpfr; 20 21 using LlvmLibcAcosfTest = LIBC_NAMESPACE::testing::FPTest<float>; 22 23 TEST_F(LlvmLibcAcosfTest, SpecialNumbers) { 24 LIBC_NAMESPACE::libc_errno = 0; 25 26 EXPECT_FP_EQ_ALL_ROUNDING(aNaN, LIBC_NAMESPACE::acosf(aNaN)); 27 EXPECT_MATH_ERRNO(0); 28 29 EXPECT_FP_EQ_ALL_ROUNDING(aNaN, LIBC_NAMESPACE::acosf(inf)); 30 EXPECT_MATH_ERRNO(EDOM); 31 32 EXPECT_FP_EQ_ALL_ROUNDING(aNaN, LIBC_NAMESPACE::acosf(neg_inf)); 33 EXPECT_MATH_ERRNO(EDOM); 34 35 EXPECT_FP_EQ_ALL_ROUNDING(zero, LIBC_NAMESPACE::acosf(1.0f)); 36 EXPECT_MATH_ERRNO(0); 37 38 EXPECT_FP_EQ_ALL_ROUNDING(aNaN, LIBC_NAMESPACE::acosf(2.0f)); 39 EXPECT_MATH_ERRNO(EDOM); 40 41 EXPECT_FP_EQ_ALL_ROUNDING(aNaN, LIBC_NAMESPACE::acosf(-2.0f)); 42 EXPECT_MATH_ERRNO(EDOM); 43 } 44 45 TEST_F(LlvmLibcAcosfTest, InFloatRange) { 46 constexpr uint32_t COUNT = 100'000; 47 constexpr uint32_t STEP = UINT32_MAX / COUNT; 48 for (uint32_t i = 0, v = 0; i <= COUNT; ++i, v += STEP) { 49 float x = FPBits(v).get_val(); 50 if (FPBits(v).is_nan() || FPBits(v).is_inf()) 51 continue; 52 ASSERT_MPFR_MATCH_ALL_ROUNDING(mpfr::Operation::Acos, x, 53 LIBC_NAMESPACE::acosf(x), 0.5); 54 } 55 } 56 57 TEST_F(LlvmLibcAcosfTest, SpecificBitPatterns) { 58 constexpr int N = 13; 59 constexpr uint32_t INPUTS[N] = { 60 0x3f000000, // x = 0.5f 61 0x3f3504f3, // x = sqrt(2)/2, FE_DOWNWARD 62 0x3f3504f4, // x = sqrt(2)/2, FE_UPWARD 63 0x3f5db3d7, // x = sqrt(3)/2, FE_DOWNWARD 64 0x3f5db3d8, // x = sqrt(3)/2, FE_UPWARD 65 0x3f800000, // x = 1.0f 66 0x40000000, // x = 2.0f 67 0x328885a3, // x = 0x1.110b46p-26 68 0x39826222, // x = 0x1.04c444p-12 69 0x3d09bf86, // x = 0x1.137f0cp-5f 70 0x3de5fa1e, // x = 0x1.cbf43cp-4f 71 0x3f083a1a, // x = 0x1.107434p-1f 72 0x3f7741b6, // x = 0x1.ee836cp-1f 73 }; 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::Acos, x, 78 LIBC_NAMESPACE::acosf(x), 0.5); 79 EXPECT_MPFR_MATCH_ALL_ROUNDING(mpfr::Operation::Acos, -x, 80 LIBC_NAMESPACE::acosf(-x), 0.5); 81 } 82 } 83