d2ce3e9621
GCC provides these functions (e.g. __addtf3, etc.) in libgcc on x86_64. Since Clang supports float128, we can also enable the existing code by using float128 for fp_t if either __FLOAT128__ or __SIZEOF_FLOAT128__ is defined instead of only supporting these builtins for platforms with 128-bit IEEE long doubles. This commit defines a new tf_float typedef that matches a float with attribute((mode(TF)) on each given architecture. There are more tests that could be enabled for x86, but to keep the diff smaller, I restricted test changes to ones that started failing as part of this refactoring. This change has been tested on x86 (natively) and aarch64,powerpc64,riscv64 and sparc64 via qemu-user. This supersedes https://reviews.llvm.org/D98261 and should also cover the changes from https://github.com/llvm/llvm-project/pull/68041.
86 lines
2.3 KiB
C
86 lines
2.3 KiB
C
// RUN: %clang_builtins %s %librt -o %t && %run %t
|
|
|
|
#define QUAD_PRECISION
|
|
#include "fp_lib.h"
|
|
#include "int_lib.h"
|
|
#include <math.h>
|
|
#include <stdio.h>
|
|
|
|
#if defined(CRT_HAS_TF_MODE)
|
|
|
|
int test__compiler_rt_logbl(fp_t x) {
|
|
# if defined(__ve__)
|
|
if (fpclassify(x) == FP_SUBNORMAL)
|
|
return 0;
|
|
# endif
|
|
fp_t crt_value = __compiler_rt_logbl(x);
|
|
fp_t libm_value = logbl(x);
|
|
// Compare the values, considering all NaNs equivalent, as the spec doesn't
|
|
// specify the NaN signedness/payload.
|
|
if (crt_value != libm_value &&
|
|
!(crt_isnan(crt_value) && crt_isnan(libm_value))) {
|
|
// Split expected values into two for printf
|
|
twords x_t, crt_value_t, libm_value_t;
|
|
x_t.all = toRep(x);
|
|
crt_value_t.all = toRep(crt_value);
|
|
libm_value_t.all = toRep(libm_value);
|
|
printf("error: in __compiler_rt_logbl([%llX %llX]) = [%llX %llX] != "
|
|
"[%llX %llX]\n",
|
|
x_t.s.high, x_t.s.low, crt_value_t.s.high, crt_value_t.s.low,
|
|
libm_value_t.s.high, libm_value_t.s.low);
|
|
return 1;
|
|
}
|
|
return 0;
|
|
}
|
|
|
|
fp_t cases[] = {
|
|
1.e-6, -1.e-6, NAN, -NAN, INFINITY, -INFINITY, -1,
|
|
-0.0, 0.0, 1, -2, 2, -0.5, 0.5,
|
|
};
|
|
|
|
int main() {
|
|
const unsigned N = sizeof(cases) / sizeof(cases[0]);
|
|
for (unsigned i = 0; i < N; ++i) {
|
|
if (test__compiler_rt_logbl(cases[i]))
|
|
return 1;
|
|
}
|
|
|
|
// Test a moving 1 bit, especially to handle denormal values.
|
|
// Test the negation as well.
|
|
// Since we are comparing the compiler-rt IEEE implementation against libc's
|
|
// long double implementation, this test can only succeed if long double
|
|
// is an IEEE 128-bit floating point number (otherwise we will see mismatches
|
|
// once we reach numbers that cannot be precisely represented in long double
|
|
// format).
|
|
# if defined(CRT_LDBL_IEEE_F128)
|
|
rep_t x = signBit;
|
|
int i = 0;
|
|
while (x) {
|
|
if (test__compiler_rt_logbl(fromRep(x)))
|
|
return 1;
|
|
if (test__compiler_rt_logbl(fromRep(signBit ^ x)))
|
|
return 1;
|
|
x >>= 1;
|
|
printf("l1: %d\n", i++);
|
|
}
|
|
// Also try a couple moving ones
|
|
x = signBit | (signBit >> 1) | (signBit >> 2);
|
|
while (x) {
|
|
if (test__compiler_rt_logbl(fromRep(x)))
|
|
return 1;
|
|
if (test__compiler_rt_logbl(fromRep(signBit ^ x)))
|
|
return 1;
|
|
x >>= 1;
|
|
printf("l1: %d\n", i++);
|
|
}
|
|
#endif
|
|
|
|
return 0;
|
|
}
|
|
#else
|
|
int main() {
|
|
printf("skipped\n");
|
|
return 0;
|
|
}
|
|
#endif
|