1 /* s_ilogbl.c -- long double version of s_ilogb.c.
2 * Conversion to long double by Ulrich Drepper,
3 * Cygnus Support, drepper@cygnus.com.
7 * ====================================================
8 * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
10 * Developed at SunPro, a Sun Microsystems, Inc. business.
11 * Permission to use, copy, modify, and distribute this
12 * software is freely granted, provided that this notice
14 * ====================================================
17 #if defined(LIBM_SCCS) && !defined(lint)
18 static char rcsid
[] = "$NetBSD: $";
21 /* ilogbl(long double x)
22 * return the binary exponent of non-zero x
23 * ilogbl(0) = FP_ILOGB0
24 * ilogbl(NaN) = FP_ILOGBNAN (no signal is raised)
25 * ilogbl(+-Inf) = INT_MAX (no signal is raised)
30 #include "math_private.h"
33 int __ilogbl(long double x
)
41 GET_LDOUBLE_EXP(es
,x
);
44 GET_LDOUBLE_WORDS(es
,hx
,lx
,x
);
46 return FP_ILOGB0
; /* ilogbl(0) = FP_ILOGB0 */
47 else /* subnormal x */
49 for (ix
= -16415; lx
>0; lx
<<=1) ix
-=1;
51 for (ix
= -16383; hx
>0; hx
<<=1) ix
-=1;
55 else if (es
<0x7fff) return es
-0x3fff;
56 else if (FP_ILOGBNAN
!= INT_MAX
)
58 GET_LDOUBLE_WORDS(es
,hx
,lx
,x
);
59 if (((hx
& 0x7fffffff)|lx
) == 0)
60 /* ISO C99 requires ilogbl(+-Inf) == INT_MAX. */
65 weak_alias (__ilogbl
, ilogbl
)