floor.c File Reference

#include "libm.h"
#include "libm_util.h"
#include "libm_errno.h"
#include "libm_inlines.h"

Include dependency graph for floor.c:

Go to the source code of this file.

Macros
#define	USE_HANDLE_ERROR

Functions
double FN_PROTOTYPE()	floor (double x)

Macro Definition Documentation

USE_HANDLE_ERROR

#define USE_HANDLE_ERROR

Definition at line 31 of file floor.c.

Function Documentation

floor()

double FN_PROTOTYPE() floor

(

double

x

)

Definition at line 37 of file floor.c.

{
  double r;
  long long rexp, xneg;
 
 
  unsigned long long ux, ax, ur, mask;
 
  GET_BITS_DP64(x, ux);
  ax = ux & (~SIGNBIT_DP64);
  xneg = (ux != ax);
 
  if (ax >= 0x4340000000000000)
    {
      /* abs(x) is either NaN, infinity, or >= 2^53 */
      if (ax > 0x7ff0000000000000)
        /* x is NaN */
        return _handle_error("floor", OP_FLOOR, ux|0x0008000000000000, _DOMAIN,
                            0, EDOM, x, 0.0, 1);
      else
        return x;
    }
  else if (ax < 0x3ff0000000000000) /* abs(x) < 1.0 */
    {
      if (ax == 0x0000000000000000)
        /* x is +zero or -zero; return the same zero */
        return x;
      else if (xneg) /* x < 0.0 */
        return -1.0;
      else
        return 0.0;
    }
  else
    {
      r = x;
      rexp = ((ux & EXPBITS_DP64) >> EXPSHIFTBITS_DP64) - EXPBIAS_DP64;
      /* Mask out the bits of r that we don't want */
      mask = 1;
      mask = (mask << (EXPSHIFTBITS_DP64 - rexp)) - 1;
      ur = (ux & ~mask);
      PUT_BITS_DP64(ur, r);
      if (xneg && (ur != ux))
        /* We threw some bits away and x was negative */
        return r - 1.0;
      else
        return r;
    }
 
}