Function eulpf

test/math/lgamma.c:52–60

};

static int eulpf(float x)
{
  union { float f; uint32_t i; } u = { x };
  int e = u.i>>23 & 0xff;

  if (!e)
    e++;
  return e - 0x7f - 23;
}

static int checkulp(float d, int r)
{