1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
|
//===-- Exhaustive test for hypotf16 --------------------------------------===//
//
// Part of the LLVM Project, under the Apache License v2.0 with LLVM Exceptions.
// See https://llvm.org/LICENSE.txt for license information.
// SPDX-License-Identifier: Apache-2.0 WITH LLVM-exception
//
//===----------------------------------------------------------------------===//
#include "exhaustive_test.h"
#include "src/__support/FPUtil/FPBits.h"
#include "src/__support/FPUtil/Hypot.h"
#include "src/math/hypotf16.h"
#include "test/UnitTest/FPMatcher.h"
#include "utils/MPFRWrapper/MPFRUtils.h"
namespace mpfr = LIBC_NAMESPACE::testing::mpfr;
struct Hypotf16Checker : public virtual LIBC_NAMESPACE::testing::Test {
using FloatType = float16;
using FPBits = LIBC_NAMESPACE::fputil::FPBits<float16>;
using StorageType = typename FPBits::StorageType;
uint64_t check(uint16_t x_start, uint16_t x_stop,
[[maybe_unused]] uint16_t y_start, uint16_t y_stop,
mpfr::RoundingMode rounding) {
mpfr::ForceRoundingMode r(rounding);
if (!r.success)
return true;
uint16_t xbits = x_start;
uint64_t failed = 0;
do {
float16 x = FPBits(xbits).get_val();
uint16_t ybits = xbits;
do {
float16 y = FPBits(ybits).get_val();
bool correct = TEST_FP_EQ(LIBC_NAMESPACE::fputil::hypot<float16>(x, y),
LIBC_NAMESPACE::hypotf16(x, y));
// Using MPFR will be much slower.
// mpfr::BinaryInput<float16> input{x, y};
// bool correct = TEST_MPFR_MATCH_ROUNDING_SILENTLY(
// mpfr::Operation::Hypot, input, LIBC_NAMESPACE::hypotf16(x, y),
// 0.5,
// rounding);
failed += (!correct);
} while (ybits++ < y_stop);
} while (xbits++ < x_stop);
return failed;
}
};
using LlvmLibcHypotf16ExhaustiveTest =
LlvmLibcExhaustiveMathTest<Hypotf16Checker, 1 << 2>;
// Range of both inputs: [0, inf]
static constexpr uint16_t POS_START = 0x0000U;
static constexpr uint16_t POS_STOP = 0x7C00U;
TEST_F(LlvmLibcHypotf16ExhaustiveTest, PositiveRange) {
test_full_range_all_roundings(POS_START, POS_STOP, POS_START, POS_STOP);
}
// Range of both inputs: [-0, -inf]
static constexpr uint16_t NEG_START = 0x8000U;
static constexpr uint16_t NEG_STOP = 0xFC00U;
TEST_F(LlvmLibcHypotf16ExhaustiveTest, NegativeRange) {
test_full_range_all_roundings(NEG_START, NEG_STOP, NEG_START, NEG_STOP);
}
|
