27 lines
583 B
C
27 lines
583 B
C
|
/* Copyright (C) 1994 DJ Delorie, see COPYING.DJ for details */
|
||
|
#include <math.h>
|
||
|
|
||
|
double
|
||
|
frexp(double x, int *exptr)
|
||
|
{
|
||
|
union {
|
||
|
double d;
|
||
|
unsigned char c[8];
|
||
|
} u;
|
||
|
|
||
|
u.d = x;
|
||
|
/*
|
||
|
* The format of the number is:
|
||
|
* Sign, 12 exponent bits, 51 mantissa bits
|
||
|
* The exponent is 1023 biased and there is an implicit zero.
|
||
|
* We get the exponent from the upper bits and set the exponent
|
||
|
* to 0x3fe (1022).
|
||
|
*/
|
||
|
*exptr = (int)(((u.c[7] & 0x7f) << 4) | (u.c[6] >> 4)) - 1022;
|
||
|
u.c[7] &= 0x80;
|
||
|
u.c[7] |= 0x3f;
|
||
|
u.c[6] &= 0x0f;
|
||
|
u.c[6] |= 0xe0;
|
||
|
return u.d;
|
||
|
}
|