| //===-- 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 "src/__support/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(); |
| } |