winsup/mingw/mingwex/math/nextafterf.c


1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29

#include <math.h>

float
nextafterf (float x, float y)
{
  union
  {
    float f;
    unsigned int i;
  } u;
  if (isnan (y) || isnan (x))
    return x + y;
  if (x == y )
     return x;
  u.f = x; 
  if (u.i == 0u)
    {
      if (y > 0.0F)
	u.i = 1;
      else
	u.i = 0x80000001;
      return u.f;
    }
   if (((x > 0.0F) ^ (y > x)) == 0)
    u.i++;
  else
    u.i--;
  return u.f;
}