Forum automatisme

Code : Tout sélectionner

uint32_t rand32(void)
{
    uint32_t result = 0;
    result = rand();
    result <<= 16;
    result |= rand();
    return result;
}

// Init all engine parts with seed
void initCMWC(unsigned int seed)
{
    srand(seed);
    for (int i = 0; i < CMWC_CYCLE; i++) Q[i] = rand32();
    do c = rand32(); while (c >= CMWC_C_MAX);
}

// CMWC engine
uint32_t randCMWC(void)
{
    static uint32_t i = CMWC_CYCLE - 1;
    uint64_t t = 0;
    uint64_t a = 18782; // as Marsaglia recommends
    uint32_t r = 0xfffffffe; // as Marsaglia recommends
    uint32_t x = 0;

    i = (i + 1) & (CMWC_CYCLE - 1);
    t = a * Q[i] + c;
    c = t >> 32;
    x = t + c;
    if (x < c)
    {
        x++;
        c++;
    }

    return Q[i] = r - x;
}

int main()
{
    unsigned int seed = time(NULL);
    initCMWC(seed);
    printf("Random CMWC: %u\n", randCMWC());
}