diff options
Diffstat (limited to 's_ilogb.c')
-rw-r--r-- | s_ilogb.c | 46 |
1 files changed, 0 insertions, 46 deletions
diff --git a/s_ilogb.c b/s_ilogb.c deleted file mode 100644 index 240f9bf..0000000 --- a/s_ilogb.c +++ /dev/null @@ -1,46 +0,0 @@ - -/* @(#)s_ilogb.c 1.3 95/01/18 */ -/* - * ==================================================== - * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved. - * - * Developed at SunSoft, a Sun Microsystems, Inc. business. - * Permission to use, copy, modify, and distribute this - * software is freely granted, provided that this notice - * is preserved. - * ==================================================== - */ - -/* ieee_ilogb(double x) - * return the binary exponent of non-zero x - * ieee_ilogb(0) = 0x80000001 - * ieee_ilogb(inf/NaN) = 0x7fffffff (no signal is raised) - */ - -#include "fdlibm.h" - -#ifdef __STDC__ - int ieee_ilogb(double x) -#else - int ieee_ilogb(x) - double x; -#endif -{ - int hx,lx,ix; - - hx = (__HI(x))&0x7fffffff; /* high word of x */ - if(hx<0x00100000) { - lx = __LO(x); - if((hx|lx)==0) - return 0x80000001; /* ieee_ilogb(0) = 0x80000001 */ - else /* subnormal x */ - if(hx==0) { - for (ix = -1043; lx>0; lx<<=1) ix -=1; - } else { - for (ix = -1022,hx<<=11; hx>0; hx<<=1) ix -=1; - } - return ix; - } - else if (hx<0x7ff00000) return (hx>>20)-1023; - else return 0x7fffffff; -} |