Генерация случайных вещественных чисел в диапазоне 0..1

// ///////////////////////////////////////////////////////////////////////// //
//                                                                           //
// L'Ecuyer algorithm for uniform random generator with practically endless  // 
 
// period. Combining 2 sequences.                                            //
//                                                                           //
// ///////////////////////////////////////////////////////////////////////// //
 
#include <stdio.h>
#include <stdlib.h>
#include <time.h>
 
// ///////////////////////////////////////////////////////////////////////// //
 
long dummy;
 
// ///////////////////////////////////////////////////////////////////////// //
 
void seed(long d)
 
{
 
    dummy = d;
 
}
 
// ///////////////////////////////////////////////////////////////////////// //
 
float rnd()
 
{
 
    const long IM1 = 2147483563;
 
    const long IM2 = 2147483399;
 
    const double AM = (1./IM1);
 
    const long IMM1 = (IM1 - 1);
 
    const int IA1 = 40014;
 
    const int IA2 = 40692;
 
    const long IQ1 = 53668;
 
    const long IQ2 = 52774;
 
    const int IR1 = 12211;
 
    const int IR2 = 3791;
 
    const int NTAB = 32;
 
    const int NWUP = 8;
 
    const double NDIV = (1 + IMM1 / NTAB);
 
    const double EPS = 1.2e-7;
 
    const double RNMX = (1.0 - EPS);
 
    
 
    static long dummy2 = 123456789;
 
    static long iy = 0;
 
    static long iv[32];
 
 
 
    long k = 0;
 
 
 
    // initialize the random sequence (first set of coefficients, the 
 
    // routine close to that in the function above
 
    if (dummy <= 0 || !iy)
 
    {
 
        // avoid negative or zero seed
 
        if (dummy < 0)
 
            dummy =- dummy;
 
        else
 
            if (!dummy)
 
                dummy = 1;
 
        dummy2 = dummy;
 
        // after NWUP warmups, initialize shuffle table
 
        for (int j = NTAB + NWUP - 1; j >= 0; --j)
 
        {
 
            k = dummy / IQ1;
 
            if ((dummy = IA1 * (dummy - k * IQ1) - IR1 * k) < 0) dummy += IM1;
 
            if (j < NTAB) iv[j] = dummy;
 
        }
 
        // first specimen from the table
 
        iy = iv[0];
 
    }
 
    // regular work: generate 2 sequences.
 
    k = dummy / IQ1;
 
    if ((dummy = IA1 * (dummy - k * IQ1) - IR1 * k) < 0) dummy+=IM1;
 
    k = dummy2 / IQ2;
 
    if ((dummy2 = IA2 * (dummy2 - k * IQ2) - IR2 * k) < 0) dummy2+=IM2;
 
    // shuffle output combining 2 sequences
 
    int j = 0;
 
    iy = iv[j = iy / NDIV] - dummy2;
 
    iv[j] = dummy;
 
    if(iy < 1) iy += IMM1;
 
    // return the result, as in the previous function
 
    float temp;
 
    return ((temp = AM * iy) > RNMX ? RNMX : temp);
 
}
 
// ///////////////////////////////////////////////////////////////////////// //
 
 
int main()
{
    seed(time(NULL));
    for (int i = 0; i < 100; ++i)
        printf("%f\n", rnd());
    return 0;
}
 
// ///////////////////////////////////////////////////////////////////////// //