ly-random

ly_random.c

#include <stdint.h>

static unsigned rnd_head, rnd_tail;
static uint32_t rnd_pool[64];

static __inline
uint32_t random_poly (unsigned n)
{
        /* LY: доступ к ряду полинома через циклический буфер. */
        return rnd_pool [(rnd_tail + n) & 63];
}

static __inline
uint32_t random_rol (uint32_t v, unsigned s)
{
        /* LY: циклический сдвиг влево. */
        return (v << s) | (v >> (32 - s));
}

void random_salt (unsigned salt)
{
        /* LY: подсаливаем генератор случайностью. */
        rnd_pool[rnd_head++ & 63] += salt;
}

uint32_t random_get (void)
{
        /* LY: Генератор случайных числел на основе неприводимого
                полинома X^64 + X^4 + X^3 + X и нелинейной функции обратной связи
                отдаленно подобно вихрю Мерсенна (Mersenne twister).
                (C) 2003 Leo Yuriev
        */

        uint32_t v =
                random_poly (0)
                ^ random_poly (1)
                ^ random_poly (62)
                ^ random_poly (61);

        uint32_t s = 0xb7e15163ul
                + random_rol (random_poly (11), 5)
                + random_rol (random_poly (23), 7)
                + random_rol (random_poly (37), 1);

        v += s;
        rnd_pool[++rnd_tail & 63] = v;
        return v;
}

void random_init ( unsigned seed) {
        /* LY: инициализируем и разогреваем генератор. */
        int i;
        
        rnd_head = rnd_tail = 0;
        for(i = 0; i < 64; ++i)
                rnd_pool[i] = 0;
        
        random_salt (seed);
        for(i = 0; i < 64 * 64; ++i)
                random_get ();
}