#include <inttypes.h>
#include <stdlib.h>

#include "rand48_impl.h"

double erand48(unsigned short s[3]) {
    union {
        uint64_t u;
        double f;
    } x = {0x3ff0000000000000ULL | __rand48_step(s, __seed48 + 3) << 4};
    return x.f - 1.0;
}

double drand48(void) {
    return erand48(__seed48);
}