/* Check that certain subnormal numbers (formerly known as denormalized numbers) are rounded to within 0.5 ulp. PR other/14354. */ /* This test requires that float and unsigned int are the same size and that the sign-bit of the float is at MSB of the unsigned int. */ #if __INT_MAX__ != 2147483647L int main () { exit (0); } #else union uf { unsigned int u; float f; }; static float u2f (unsigned int v) { union uf u; u.u = v; return u.f; } static unsigned int f2u (float v) { union uf u; u.f = v; return u.u; } int ok = 1; static void tstmul (unsigned int ux, unsigned int uy, unsigned int ur) { float x = u2f (ux); float y = u2f (uy); if (f2u (x * y) != ur) /* Set a variable rather than aborting here, to simplify tracing when several computations are wrong. */ ok = 0; } /* We don't want to make this const and static, or else we risk inlining causing the test to fold as constants at compile-time. */ struct { unsigned int p1, p2, res; } expected[] = { {0xfff, 0x3f800400, 0xfff}, {0xf, 0x3fc88888, 0x17}, {0xf, 0x3f844444, 0xf} }; int main (int argc, char *argv[], char *envp[]) { unsigned int i; for (i = 0; i < sizeof (expected) / sizeof (expected[0]); i++) { tstmul (expected[i].p1, expected[i].p2, expected[i].res); tstmul (expected[i].p2, expected[i].p1, expected[i].res); } if (!ok) abort (); exit (0); } #endif
