1ca6b3542SSiva Chandra Reddy //===-- Unittests mod_2pi, mod_pi_over_4 and mod_pi_over_2 ----------------===//
2ca6b3542SSiva Chandra Reddy //
3ca6b3542SSiva Chandra Reddy // Part of the LLVM Project, under the Apache License v2.0 with LLVM Exceptions.
4ca6b3542SSiva Chandra Reddy // See https://llvm.org/LICENSE.txt for license information.
5ca6b3542SSiva Chandra Reddy // SPDX-License-Identifier: Apache-2.0 WITH LLVM-exception
6ca6b3542SSiva Chandra Reddy //
7ca6b3542SSiva Chandra Reddy //===----------------------------------------------------------------------===//
8ca6b3542SSiva Chandra Reddy 
9ca6b3542SSiva Chandra Reddy #include "src/math/generic/dp_trig.h"
10ca6b3542SSiva Chandra Reddy #include "utils/MPFRWrapper/MPFRUtils.h"
11*6c3f53c7SSiva Chandra Reddy #include "utils/UnitTest/FPMatcher.h"
12ca6b3542SSiva Chandra Reddy #include "utils/UnitTest/Test.h"
13ca6b3542SSiva Chandra Reddy 
14ca6b3542SSiva Chandra Reddy #include <math.h>
15ca6b3542SSiva Chandra Reddy 
16ca6b3542SSiva Chandra Reddy namespace mpfr = __llvm_libc::testing::mpfr;
17ca6b3542SSiva Chandra Reddy using FPBits = __llvm_libc::fputil::FPBits<double>;
18ca6b3542SSiva Chandra Reddy using UIntType = FPBits::UIntType;
19ca6b3542SSiva Chandra Reddy 
TEST(LlvmLibcMod2PITest,Range)20ca6b3542SSiva Chandra Reddy TEST(LlvmLibcMod2PITest, Range) {
21ca6b3542SSiva Chandra Reddy   constexpr UIntType count = 1000000000;
22ca6b3542SSiva Chandra Reddy   constexpr UIntType step = UIntType(-1) / count;
23ca6b3542SSiva Chandra Reddy   for (UIntType i = 0, v = 0; i <= count; ++i, v += step) {
24ca6b3542SSiva Chandra Reddy     double x = double(FPBits(v));
25ca6b3542SSiva Chandra Reddy     if (isnan(x) || isinf(x) || x <= 0.0)
26ca6b3542SSiva Chandra Reddy       continue;
27ca6b3542SSiva Chandra Reddy 
28ca6b3542SSiva Chandra Reddy     ASSERT_MPFR_MATCH(mpfr::Operation::Mod2PI, x, __llvm_libc::mod_2pi(x), 0);
29ca6b3542SSiva Chandra Reddy   }
30ca6b3542SSiva Chandra Reddy }
31ca6b3542SSiva Chandra Reddy 
TEST(LlvmLibcModPIOver2Test,Range)32ca6b3542SSiva Chandra Reddy TEST(LlvmLibcModPIOver2Test, Range) {
33ca6b3542SSiva Chandra Reddy   constexpr UIntType count = 1000000000;
34ca6b3542SSiva Chandra Reddy   constexpr UIntType step = UIntType(-1) / count;
35ca6b3542SSiva Chandra Reddy   for (UIntType i = 0, v = 0; i <= count; ++i, v += step) {
36ca6b3542SSiva Chandra Reddy     double x = double(FPBits(v));
37ca6b3542SSiva Chandra Reddy     if (isnan(x) || isinf(x) || x <= 0.0)
38ca6b3542SSiva Chandra Reddy       continue;
39ca6b3542SSiva Chandra Reddy 
40ca6b3542SSiva Chandra Reddy     ASSERT_MPFR_MATCH(mpfr::Operation::ModPIOver2, x,
41ca6b3542SSiva Chandra Reddy                       __llvm_libc::mod_pi_over_2(x), 0);
42ca6b3542SSiva Chandra Reddy   }
43ca6b3542SSiva Chandra Reddy }
44ca6b3542SSiva Chandra Reddy 
TEST(LlvmLibcModPIOver4Test,Range)45ca6b3542SSiva Chandra Reddy TEST(LlvmLibcModPIOver4Test, Range) {
46ca6b3542SSiva Chandra Reddy   constexpr UIntType count = 1000000000;
47ca6b3542SSiva Chandra Reddy   constexpr UIntType step = UIntType(-1) / count;
48ca6b3542SSiva Chandra Reddy   for (UIntType i = 0, v = 0; i <= count; ++i, v += step) {
49ca6b3542SSiva Chandra Reddy     double x = double(FPBits(v));
50ca6b3542SSiva Chandra Reddy     if (isnan(x) || isinf(x) || x <= 0.0)
51ca6b3542SSiva Chandra Reddy       continue;
52ca6b3542SSiva Chandra Reddy 
53ca6b3542SSiva Chandra Reddy     ASSERT_MPFR_MATCH(mpfr::Operation::ModPIOver4, x,
54ca6b3542SSiva Chandra Reddy                       __llvm_libc::mod_pi_over_4(x), 0);
55ca6b3542SSiva Chandra Reddy   }
56ca6b3542SSiva Chandra Reddy }
57