aleas.cxx File Reference

#include <time.h>
#include <sys/types.h>
#include <unistd.h>
#include <math.h>
#include "aleas.hxx"

Include dependency graph for aleas.cxx:

Go to the source code of this file.

Macros
#define	NOALEAS

Functions
void	init_by_array (unsigned long[], int)
double	genrand_real1 (void)
double	genrand_real3 (void)
static void	initrand (void)
static double	randGauss (void)

Macro Definition Documentation

#define NOALEAS

Definition at line 35 of file aleas.cxx.

Function Documentation

double genrand_real1

(

void

)

Definition at line 122 of file mt19937ar.cxx.

References genrand_int32().

Referenced by Sphere::fill(), and Cube::tire().

{
    return genrand_int32()*(1.0/4294967295.0); 
    /* divided by 2^32-1 */ 
}

double genrand_real3

(

void

)

Definition at line 136 of file mt19937ar.cxx.

References genrand_int32().

Referenced by randGauss().

{
    return (((double)genrand_int32()) + 0.5)*(1.0/4294967296.0); 
    /* divided by 2^32 */
}

void init_by_array	(	unsigned	long[],
		int
	)

Definition at line 52 of file mt19937ar.cxx.

References CORBAEngineTest::i, init_genrand(), mt, and N.

Referenced by initrand().

{
    int i, j, k;
    init_genrand(19650218UL);
    i=1; j=0;
    k = (N>key_length ? N : key_length);
    for (; k; k--) {
        mt[i] = (mt[i] ^ ((mt[i-1] ^ (mt[i-1] >> 30)) * 1664525UL))
          + init_key[j] + j; /* non linear */
        mt[i] &= 0xffffffffUL; /* for WORDSIZE > 32 machines */
        i++; j++;
        if (i>=N) { mt[0] = mt[N-1]; i=1; }
        if (j>=key_length) j=0;
    }
    for (k=N-1; k; k--) {
        mt[i] = (mt[i] ^ ((mt[i-1] ^ (mt[i-1] >> 30)) * 1566083941UL))
          - i; /* non linear */
        mt[i] &= 0xffffffffUL; /* for WORDSIZE > 32 machines */
        i++;
        if (i>=N) { mt[0] = mt[N-1]; i=1; }
    }

    mt[0] = 0x80000000UL; /* MSB is 1; assuring non-zero initial array */ 
}

static void initrand ( void  )

static

Definition at line 37 of file aleas.cxx.

References init_by_array().

Referenced by Aleatoire::Aleatoire().

{
     static long done=0;
     unsigned long vec[]={
             31081996, 21012006, 17921963, 0,
             11101998, 2112003, 25111964, 0
         };

     if (! done) {
#ifdef NOALEAS
         vec[3] = vec[7] = 2082007;
#else
         vec[3] = (unsigned long) getpid();
         vec[7] = (unsigned long) time(NULL);
#endif
         init_by_array(vec, 8);
         done += 1;
     }
} 

static double randGauss ( void  )

static

Definition at line 73 of file aleas.cxx.

References genrand_real3(), and gui.logview::log.

Referenced by Normale::tire(), and NormalePositif::tire().

{
    static double   next;
    static int      flag=0;
    double          v2, d, fac;
    
    if (flag) {
        flag = 0;
        return next;
    }
    else {
        do {
            next = 2.0*genrand_real3() - 1.0;
            v2 = 2.0*genrand_real3() - 1.0;
            d = next*next + v2*v2;
        } while (d >= 1.0 || d == 0.0);
        fac = sqrt(-2.0*log(d)/d);
        next *= fac;
        flag = 1;
        return v2 * fac;
    }

}