12345678910111213141516171819202122232425262728293031323334353637383940414243444546 |
- #include "math.h"
- #include "math_private.h"
- static const double
- two54 = 1.80143985094819840000e+16;
- double frexp(double x, int *eptr)
- {
- int32_t hx, ix, lx;
- EXTRACT_WORDS(hx,lx,x);
- ix = 0x7fffffff&hx;
- *eptr = 0;
- if(ix>=0x7ff00000||((ix|lx)==0)) return x;
- if (ix<0x00100000) {
- x *= two54;
- GET_HIGH_WORD(hx,x);
- ix = hx&0x7fffffff;
- *eptr = -54;
- }
- *eptr += (ix>>20)-1022;
- hx = (hx&0x800fffff)|0x3fe00000;
- SET_HIGH_WORD(x,hx);
- return x;
- }
- libm_hidden_def(frexp)
|