1/**
 2 * This file has no copyright assigned and is placed in the Public Domain.
 3 * This file is part of the mingw-w64 runtime package.
 4 * No warranty is given; refer to the file DISCLAIMER.PD within this package.
 5 */
 6/*
 7   nexttowardf.c
 8   Contributed by Danny Smith <dannysmith@users.sourceforge.net>
 9   No copyright claimed, absolutely no warranties.
10
11   2005-05-10
12*/
13
14#include <math.h>
15
16float
17nexttowardf (float x, long double y)
18{
19  union
20  {
21    float f;
22    unsigned int i;
23  } u;
24
25  long double xx = x;
26
27  if (isnan (y) || isnan (x))
28    return x + y;
29  if (xx == y )
30     /* nextafter (0.0, -O.0) should return -0.0.  */
31     return y;
32  u.f = x; 
33  if (x == 0.0F)
34    {
35      u.i = 1;
36      return y > 0.0L ? u.f : -u.f;
37    }
38  if (((x > 0.0F) ^ (y > xx)) == 0)
39    u.i++;
40  else
41    u.i--;
42  return u.f;
43}