#include <stdio.h>
#include <math.h>

int fesetround(int);
#define FE_DOWNWARD     (1)
#define FE_TONEAREST    (0)
#define FE_TOWARDZERO   (3)
#define FE_UPWARD       (2)


int main()
{
  int rnd[4] = { FE_TONEAREST, FE_TOWARDZERO, FE_UPWARD, FE_DOWNWARD };
  char Rnd[4] = "NZUD";
  float x = 0x1.ff07fep+127f;
  float y;
  for (int i = 0; i < 4; i++)
  {
    fesetround(rnd[i]);
    y = sqrtf (x);
    printf ("RND%c: %a\n", Rnd[i], y);
  }
}