comparison src/floatfns.c @ 6359:800c035273e9

(Flogb): Check for 0.0. Emulate logb if needed.

6358:250e69a75938

{
  Lisp_Object val;
  int value;
  double f = extract_float (arg);

#ifdef HAVE_LOGB
  IN_FLOAT (value = logb (f), "logb", arg);
  XSET (val, Lisp_Int, value);
#else
#ifdef HAVE_FREXP
  {
    int exp;  

    IN_FLOAT (frexp (f, &exp), "logb", arg);
    XSET (val, Lisp_Int, exp-1);
  }
#else
  /* Would someone like to write code to emulate logb?  */
  error ("`logb' not implemented on this operating system");
#endif
#endif

  return val;
}

/* the rounding functions  */

6359:800c035273e9

{
  Lisp_Object val;
  int value;
  double f = extract_float (arg);

  if (f == 0.0)
    value = -(VALMASK >> 1);
  else
    {
#ifdef HAVE_LOGB
      IN_FLOAT (value = logb (f), "logb", arg);
#else
#ifdef HAVE_FREXP
      IN_FLOAT (frexp (f, &value), "logb", arg);
      value--;
#else
      int i;
      double d;
      if (f < 0.0)
	f = -f;
      value = -1;
      while (f < 0.5)
	{
	  for (i = 1, d = 0.5; d * d >= f; i += i)
	    d *= d;
	  f /= d;
	  value -= i;
	}
      while (f >= 1.0)
	{
	  for (i = 1, d = 2.0; d * d <= f; i += i)
	    d *= d;
	  f /= d;
	  value += i;
	}
#endif
#endif
    }
  XSET (val, Lisp_Int, value);
  return val;
}

/* the rounding functions  */