shylie/llvm-project
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

Relax muldc3 test to avoid precision issue (#151663)

Browse Source
This commit is contained in:
jinge90 2025-08-02 13:09:21 +08:00 committed by GitHub
parent ee67f78776
commit de7ee884a6
No known key found for this signature in database
GPG Key ID: B5690EEEBB952194
1 changed files with 20 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

22
compiler-rt/test/builtins/Unit/muldc3_test.c
View File

@ -7,6 +7,7 @@
#include <complex.h>
#include <stdio.h>
#define RELATIVE_TOLERANCE 1e-9
// Returns: the product of a + ib and c + id
@ -15,6 +16,19 @@ __muldc3(double __a, double __b, double __c, double __d);
enum {zero, non_zero, inf, NaN, non_zero_nan};
int check_complex_equal(double _Complex r1, double _Complex r2)
{
double max_magnitude = fmax(cabs(r1), cabs(r2));
double real_diff = fabs(creal(r1) - creal(r2));
double imag_diff = fabs(cimag(r1) - cimag(r2));
if (real_diff >= max_magnitude * RELATIVE_TOLERANCE)
return 0;
if (imag_diff >= max_magnitude * RELATIVE_TOLERANCE)
return 0;
return 1;
}
int
classify(double _Complex x)
{
@ -46,11 +60,15 @@ int test__muldc3(double a, double b, double c, double d)
// a, b, c, d, creal(r), cimag(r));
double _Complex dividend;
double _Complex divisor;
double _Complex temp;
__real__ dividend = a;
__imag__ dividend = b;
__real__ divisor = c;
__imag__ divisor = d;
__real__ temp = a * c - b * d;
__imag__ temp = a * d + b * c;
switch (classify(dividend))
{
@ -89,7 +107,7 @@ int test__muldc3(double a, double b, double c, double d)
case non_zero:
if (classify(r) != non_zero)
return 1;
if (r != a * c - b * d + _Complex_I*(a * d + b * c))
if (!check_complex_equal(r, temp))
return 1;
break;
case inf: