logb.c File Reference

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

Include dependency graph for logb.c:

Go to the source code of this file.

Macros
#define	USE_INFINITY_WITH_FLAGS
#define	USE_HANDLE_ERROR

Functions
double	_logb (double x)

Macro Definition Documentation

USE_HANDLE_ERROR

#define USE_HANDLE_ERROR

Definition at line 31 of file logb.c.

USE_INFINITY_WITH_FLAGS

#define USE_INFINITY_WITH_FLAGS

Definition at line 30 of file logb.c.

Function Documentation

_logb()

double _logb

(

double

x

)

Definition at line 38 of file logb.c.

{
 
  unsigned long long ux;
  long long u;
  GET_BITS_DP64(x, ux);
  u = ((ux & EXPBITS_DP64) >> EXPSHIFTBITS_DP64) - EXPBIAS_DP64;
  if ((ux & ~SIGNBIT_DP64) == 0)
    /* x is +/-zero. Return -infinity with div-by-zero flag. */
    return _handle_error("_logb", OP_LOGB, NINFBITPATT_DP64, _SING,
                        AMD_F_DIVBYZERO, ERANGE, x, 0.0, 1);
  else if (EMIN_DP64 <= u && u <= EMAX_DP64)
    /* x is a normal number */
    return (double)u;
  else if (u > EMAX_DP64)
    {
      /* x is infinity or NaN */
      if ((ux & MANTBITS_DP64) == 0)
        /* x is +/-infinity. For VC++, return infinity of same sign. */
        return x;
      else
        /* x is NaN, result is NaN */
        return _handle_error("_logb", OP_LOGB, ux|0x0008000000000000, _DOMAIN,
                            0, EDOM, x, 0.0, 1);
    }
  else
    {
      /* x is denormalized. */
#ifdef FOLLOW_IEEE754_LOGB
      /* Return the value of the minimum exponent to ensure that
         the relationship between logb and scalb, defined in
         IEEE 754, holds. */
      return EMIN_DP64;
#else
      /* Follow the rule set by IEEE 854 for logb */
      ux &= MANTBITS_DP64;
      u = EMIN_DP64;
      while (ux < IMPBIT_DP64)
        {
          ux <<= 1;
          u--;
        }
      return (double)u;
#endif
    }
 
}