1  /* s_ceilf.c -- float version of s_ceil.c.
       2   */
       3  
       4  /*
       5   * ====================================================
       6   * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
       7   *
       8   * Developed at SunPro, a Sun Microsystems, Inc. business.
       9   * Permission to use, copy, modify, and distribute this
      10   * software is freely granted, provided that this notice
      11   * is preserved.
      12   * ====================================================
      13   */
      14  
      15  #define NO_MATH_REDIRECT
      16  #include <math.h>
      17  #include <math_private.h>
      18  #include <libm-alias-float.h>
      19  #include <math-use-builtins.h>
      20  
      21  float
      22  __ceilf (float x)
      23  {
      24  #if USE_CEILF_BUILTIN
      25    return __builtin_ceilf (x);
      26  #else
      27    /* Use generic implementation.  */
      28    int32_t i0, j0;
      29    uint32_t i;
      30  
      31    GET_FLOAT_WORD (i0, x);
      32    j0 = ((i0 >> 23) & 0xff) - 0x7f;
      33    if (j0 < 23)
      34      {
      35        if (j0 < 0)
      36  	{
      37  	  /* return 0 * sign (x) if |x| < 1  */
      38  	  if (i0 < 0)
      39  	    i0 = 0x80000000;
      40  	  else if (i0 != 0)
      41  	    i0 = 0x3f800000;
      42  	}
      43        else
      44  	{
      45  	  i = (0x007fffff) >> j0;
      46  	  if ((i0 & i) == 0)
      47  	    return x;		/* x is integral  */
      48  	  if (i0 > 0)
      49  	    i0 += (0x00800000) >> j0;
      50  	  i0 &= (~i);
      51  	}
      52      }
      53    else
      54      {
      55        if (__glibc_unlikely (j0 == 0x80))
      56  	return x + x;		/* inf or NaN  */
      57        else
      58  	return x;		/* x is integral  */
      59      }
      60    SET_FLOAT_WORD (x, i0);
      61    return x;
      62  #endif /* ! USE_CEILF_BUILTIN  */
      63  }
      64  #ifndef __ceilf
      65  libm_alias_float (__ceil, ceil)
      66  #endif