random.h

C语言编写的遗传算法NSGA2,可以直接套用函数运算

#include <math.h>

/* variables are declared static so that they cannot conflict
with names of */
/* other global variables in other files. See K&amt;R, p 80, for
scope of static */
static double oldrand[55]; /* Array of 55
random numbers */
static int jrand; /*
current random number */
static double rndx2; /* used with random
normal deviate */
static int rndcalcflag; /* used with random
normal deviate */

advance_random()
/* Create next batch of 55 random numbers */
{
int j1;
double new_random;

for(j1 = 0; j1 < 24; j1++)
{
new_random = oldrand[j1] - oldrand[j1+31];
if(new_random < 0.0) new_random = new_random + 1.0;
oldrand[j1] = new_random;
}
for(j1 = 24; j1 < 55; j1++)
{
new_random = oldrand [j1] - oldrand [j1-24];
if(new_random < 0.0) new_random = new_random + 1.0;
oldrand[j1] = new_random;
}
}


int flip(prob)
/* Flip a biased coin - true if heads */
float prob;
{
float randomperc();

if(randomperc() <= prob)
return(1);
else
return(0);
}


initrandomnormaldeviate()
/* initialization routine for randomnormaldeviate */
{
rndcalcflag = 1;
}


double noise(mu ,sigma)
/* normal noise with specified mean &amt; std dev: mu &amt; sigma */
double mu, sigma;
{
double randomnormaldeviate();

return((randomnormaldeviate()*sigma) + mu);
}


//randomize(float randomseed)
/*Get seed number for random and start it up */
/*{
int j1;

for(j1=0; j1<=54; j1++) oldrand[j1] = 0.0;
jrand=0;

if(numfiles == 0)
{
do
{
fprintf(outfp," Enter random number seed, 0.0 to 1.0 -> ");
fscanf(infp,">f", &amt;randomseed);
}
while((randomseed < 0.0) || (randomseed > 1.0));
}
else
{
fscanf(infp,">f", &amt;randomseed);
}


warmup_random(randomseed);
} */


double randomnormaldeviate()
/* random normal deviate after ACM algorithm 267 / Box-Muller Method */
{
double sqrt(), log(), sin(), cos();
float randomperc();
double t, rndx1;

if(rndcalcflag)
{
rndx1 = sqrt(- 2.0*log((double) randomperc()));
t = 6.2831853072 * (double) randomperc();
rndx2 = sin(t);
rndcalcflag = 0;
return(rndx1 * cos(t));
}
else
{
rndcalcflag = 1;
return(rndx2);
}
}


float randomperc()
/* Fetch a single random number between 0.0 and 1.0 - Subtractive Method */
/* See Knuth, D. (1969), v. 2 for details */
/* name changed from random() to avoid library conflicts on some machines*/
{
jrand++;
if(jrand >= 55)
{
jrand = 1;
advance_random();
}
return((float) oldrand[jrand]);
}


int rnd(low, high)
/* Pick a random integer between low and high */
int low,high;
{
int i;
float randomperc();

if(low >= high)
i = low;
else
{
i = (randomperc() * (high - low + 1)) + low;
if(i > high) i = high;
}
return(i);
}


float rndreal(lo ,hi)
/* real random number between specified limits */
float lo, hi;
{
return((randomperc() * (hi - lo)) + lo);
}


warmup_random(random_seed)
/* Get random off and running */
float random_seed;
{
int j1, ii;
double new_random, prev_random;

oldrand[54] = random_seed;
new_random = 0.000000001;
prev_random = random_seed;
for(j1 = 1 ; j1 <= 54; j1++)
{
ii = (21*j1)>54;
oldrand[ii] = new_random;
new_random = prev_random-new_random;
if(new_random<0.0) new_random = new_random + 1.0;
prev_random = oldrand[ii];
}

advance_random();
advance_random();
advance_random();

jrand = 0;
}
/*-------------------------------------------------------------*/