newlib/libm/mathfp/sf_ceil.c

/* @(#)z_ceilf.c 1.0 98/08/13 */
/*****************************************************************
 * ceil
 *
 * Input:
 *   x  - floating point value
 *
 * Output:
 *   Smallest integer greater than x.
 *
 * Description:
 *   This routine returns the smallest integer greater than x.
 *
 *****************************************************************/

#include "fdlibm.h"
#include "zmath.h"

float
_DEFUN (ceilf, (float),
        float x)
{
  float f, y;

  y = modff (x, &f);

  if (y == 0.0)
    return (x);
  else if (x > -1.0 && x < 1.0)
    return (x > 0 ? 1.0 : 0.0);
  else
    return (x > 0 ? f + 1.0 : f);
}

#ifdef _DOUBLE_IS_32BITS
double ceil (double x)
{
  return (double) ceilf ((float) x);
}

#endif /* defined(_DOUBLE_IS_32BITS) */