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
|
//===-- Template for diffing remquo results ---------------------*- C++ -*-===//
//
// 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 "utils/FPUtil/FPBits.h"
#include <math.h>
#include <stddef.h>
#include <stdint.h>
template <typename T> using RemQuoFunc = T (*)(T, T, int *);
template <typename T>
void RemQuoDiff(RemQuoFunc<T> func1, RemQuoFunc<T> func2, const uint8_t *data,
size_t size) {
constexpr size_t typeSize = sizeof(T);
if (size < 2 * typeSize)
return;
T x = *reinterpret_cast<const T *>(data);
T y = *reinterpret_cast<const T *>(data + typeSize);
int q1, q2;
T remainder1 = func1(x, y, &q1);
T remainder2 = func2(x, y, &q2);
if (isnan(remainder1)) {
if (!isnan(remainder2))
__builtin_trap();
return;
}
if (isinf(remainder2) != isinf(remainder1))
__builtin_trap();
// Compare only the 3 LS bits of the quotient.
if ((q1 & 0x7) != (q2 & 0x7))
__builtin_trap();
__llvm_libc::fputil::FPBits<T> bits1(remainder1);
__llvm_libc::fputil::FPBits<T> bits2(remainder2);
if (bits1.uintval() != bits2.uintval())
__builtin_trap();
}
|
