/* 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;
}