Index: test/builtins/Unit/divsc3_test.c =================================================================== --- test/builtins/Unit/divsc3_test.c +++ test/builtins/Unit/divsc3_test.c @@ -14,6 +14,7 @@ #include "int_lib.h" #include #include +#include #include // Returns: the quotient of (a + ib) / (c + id) @@ -47,6 +48,25 @@ return non_zero; } +// check for equality assuming that both real and imaginary parts +// can differ by exactly 1 representable value, in order to handle +// different levels of accuracy on 32-bit x86 +static bool approx_equal(float _Complex a, float _Complex b) { + if (a != b) { + float ra = __real__ a; + float ia = __imag__ a; + float rb = __real__ b; + float ib = __imag__ b; + + if (ra != rb && nextafterf(ra, rb) != rb) + return false; + if (ia != ib && nextafterf(ia, ib) != ib) + return false; + } + + return true; +} + int test__divsc3(float a, float b, float c, float d) { float _Complex r = __divsc3(a, b, c, d); @@ -100,7 +120,7 @@ { float _Complex z = (a * c + b * d) / (c * c + d * d) + (b * c - a * d) / (c * c + d * d) * _Complex_I; - if (r != z) + if (!approx_equal(r, z)) return 1; } break;