void Smooth()
{
int imember; /* current member */
int i, j, igene; /* counters */
int begin, end; /* limit the operator */
static double *Vector_Aux;
static int first=1;
if (first)
{
Vector_Aux =
(double *) calloc((unsigned) Genes, sizeof(double));
if (Vector_Aux == NULL)
Error("Allocation failed for Vector_Aux");
first = 0;
}
/*
For all the members of the current population
*/
if (verbose)
fprintf(stderr,"\nSubpopulation %d smoothing trial solutions\n", instance);
for (imember = 0; imember < Popsize; imember++)
{
Unpack(New[imember].Gene, Bitstring, Length);
FloatRep(Bitstring, Vector, Genes);
for (igene = 0, Vector_Aux[igene] = 0.;
igene < NSOURCES + NRECEIVERS;
igene++, Vector_Aux[igene] = 0.)
{
if (igene < NSOURCES)
{
/* no smooth of source statics */
Vector_Aux[igene] = Vector[igene];
}
else
{
i = to_filter[igene - NSOURCES][0];
j = to_filter[igene - NSOURCES][1];
if (to_filter[igene - NSOURCES][2] == 2)
Vector_Aux[igene] =
(Vector[i] + Vector[j]) / 2.;
else
Vector_Aux[igene] =
(Vector[i] + Vector[igene] +
Vector[j]) / 3.;
}
}
/*
Packing
*/
StringRep(Vector_Aux, Bitstring, Genes);
Pack(Bitstring, New[imember].Gene, Length);
New[imember].Needs_evaluation = 1;
}
}
void Initialize () {
FILE *fp;
register int i, j; /* loop control */
int status; /* indicates end of file in initfile */
Trace("Initialize entered");
if (Totalexperiments > 1)
sprintf(Bestfile, "%s.%d", Minfile, Experiment+1);
/* prepare for new experiment */
Doneflag = 0;
Curr_dump = 0;
Bestsize = 0;
/* set next mutation location */
if (M_rate < 1.0)
Mu_next = ceil (log(Rand()) / log(1.0 - M_rate));
else
Mu_next = 1;
Trials = Gen = 0;
Lost = Conv = 0;
Plateau = 0;
Spin = 0;
Onsum = Offsum = 0.0;
for (i=0; i<Windowsize; i++) Window[i] = 0.0;
/* set up initial population */
i = 0; /* current structure */
if (Initflag) { /* get some structures from Initfile */
if ((fp = fopen(Initfile, "r")) == NULL) {
char msg[40];
sprintf(msg, "Initialize: can't open %s", Initfile);
Error(msg);
}
status = 1;
if (Floatflag) {
for (j = 0; j < Genes && status != EOF; j++) {
status = fscanf(fp, "%lf", &Vector[j]);
}
}
else
status = fscanf(fp, "%s", Bitstring);
while (status != EOF && i < Popsize) {
if (Floatflag)
StringRep(Vector, Bitstring, Genes);
Pack(Bitstring, New[i].Gene, Length);
New[i].Needs_evaluation = 1;
i++;
/* get the next structure */
if (Floatflag)
for (j = 0; j < Genes && status != EOF; j++)
status = fscanf(fp, "%lf", &Vector[j]);
else
status = fscanf(fp, "%s", Bitstring);
}
fclose(fp);
}
/********************************************************/
/* The seed for the random number generator is saved */
/* after the initialization of the first population */
/* in each experiment. The saved value is used as the */
/* Seed in subsequent experiments. The reason is: */
/* often we will run several experiments with one set */
/* of parameters, and compare the results with several */
/* experiments run with a different set of parameters. */
/* This means that for all runs which use the */
/* same population size, the initial populations for */
/* corresponding experiments will be the same. */
/********************************************************/
if ( Experiment > 0 ) Seed = Initseed;
for (; i < Popsize; i++) { /* initialize remainder randomly */
for (j = 0; j < Length; j++) {
if (Randint(0,1))
Bitstring[j] = '1';
else
Bitstring[j] = '0';
}
Pack(Bitstring , New[i].Gene , Length);
New[i].Needs_evaluation = 1;
}
Initseed = Seed;
Trace("Initialize completed");
} /* end Initialize */
