[compiler-rt] de7ee88 - Relax muldc3 test to avoid precision issue (#151663)

[via llvm-commits]

Author: jinge90
Date: 2025-08-02T13:09:21+08:00
New Revision: de7ee884a6ff50afc6379ee9febb4780ddedff11

URL: https://github.com/llvm/llvm-project/commit/de7ee884a6ff50afc6379ee9febb4780ddedff11
DIFF: https://github.com/llvm/llvm-project/commit/de7ee884a6ff50afc6379ee9febb4780ddedff11.diff

LOG: Relax muldc3 test to avoid precision issue (#151663)

Added: 
    

Modified: 
    compiler-rt/test/builtins/Unit/muldc3_test.c

Removed: 
    


################################################################################
diff  --git a/compiler-rt/test/builtins/Unit/muldc3_test.c b/compiler-rt/test/builtins/Unit/muldc3_test.c
index 4941078e99b06..1787b83f771b3 100644
--- a/compiler-rt/test/builtins/Unit/muldc3_test.c
+++ b/compiler-rt/test/builtins/Unit/muldc3_test.c
@@ -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: