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