Subversion Repositories scarts

/* 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