projects
/
oweals
/
musl.git
/ blob
commit
grep
author
committer
pickaxe
?
search:
re
summary
|
shortlog
|
log
|
commit
|
commitdiff
|
tree
history
|
raw
|
HEAD
assembly optimizations for fmod/remainder functions
[oweals/musl.git]
/
src
/
math
/
ilogb.c
1
#include <limits.h>
2
#include "libm.h"
3
4
int ilogb(double x)
5
{
6
union dshape u = {x};
7
int e = u.bits>>52 & 0x7ff;
8
9
if (!e) {
10
u.bits <<= 12;
11
if (u.bits == 0)
12
return FP_ILOGB0;
13
/* subnormal x */
14
for (e = -0x3ff; u.bits < (uint64_t)1<<63; e--, u.bits<<=1);
15
return e;
16
}
17
if (e == 0x7ff)
18
return u.bits<<12 ? FP_ILOGBNAN : INT_MAX;
19
return e - 0x3ff;
20
}