Beispiel #1
0
/* initialize all variables at the beginning of the program */
void init_peaks ()
{
	movrandseed = 1;
	recent_change = 1; /* indicates that a change has just ocurred */
	offline_performance = 0.0;
	offline_error = 0.0;
	avg_error=0;    /* average error so far */ 
	current_error=0;/* error of the currently best individual */
	evals = 0;         /* number of evaluations so far */



  int i,j;
  double dummy;
  
  shift = (double *) calloc(geno_size, sizeof(double));
  coordinates = (double *) calloc(geno_size, sizeof(double));
  covered_peaks = (int *) calloc(number_of_peaks, sizeof(int));
  peak = (double **) calloc(number_of_peaks, sizeof(double*));
  prev_movement = (double **) calloc(number_of_peaks, sizeof(double*));
  for (i=0; i< number_of_peaks; i++){
    peak[i]= (double *) calloc(geno_size+2, sizeof(double));
    prev_movement[i] = (double *) calloc(geno_size, sizeof(double));
  }
  for (i=0; i< number_of_peaks; i++)
    for (j=0; j< geno_size; j++){
      peak[i][j] = 100.0*movrand();
      prev_movement[i][j] = movrand()-0.5;
    }
  if (standardheight <= 0.0)
    for (i=0; i< number_of_peaks; i++)
      peak[i][geno_size+1]= (maxheight-minheight)*movrand()+minheight;
  else
    for (i=0; i< number_of_peaks; i++)
      peak[i][geno_size+1]= standardheight;
  if (standardwidth <= 0.0)
    for (i=0; i< number_of_peaks; i++)
      peak[i][geno_size]= (maxwidth-minwidth)*movrand()+minwidth;
  else
    for (i=0; i< number_of_peaks; i++)
      peak[i][geno_size]= standardwidth;
  if(calculate_average_error){
    global_max = -100000.0;
    for (i=0;i<number_of_peaks; i++){
      for (j=0; j<geno_size; j++)
        coordinates[j]=peak[i][j];
      dummy = dummy_eval(coordinates);
      if (dummy>global_max)
        global_max = dummy;
    }
  }
}
/* initialize all variables at the beginning of the program */
void init_peaks ()
{
  int i,j;
  double dummy;
  
  shift = (double *) calloc(geno_size, sizeof(double));
  coordinates = (double *) calloc(geno_size, sizeof(double));
  covered_peaks = (int *) calloc(number_of_peaks, sizeof(int));
  peak = (double **) calloc(number_of_peaks, sizeof(double*));
  prev_movement = (double **) calloc(number_of_peaks, sizeof(double*));


  for (i=0; i< number_of_peaks; i++){
    peak[i]= (double *) calloc(geno_size+2, sizeof(double));
    prev_movement[i] = (double *) calloc(geno_size, sizeof(double));
  }

  for (i=0; i< number_of_peaks; i++)
    for (j=0; j< geno_size; j++){
// old      peak[i][j] = 100.0*movrand();
	  peak[i][j] = (maxcoordinate - mincoordinate) * movrand() + mincoordinate;
      prev_movement[i][j] = movrand()-0.5;
    }

  if (standardheight <= 0.0)
    for (i=0; i< number_of_peaks; i++)
      peak[i][geno_size+1]= (maxheight-minheight)*movrand()+minheight;
  else
    for (i=0; i< number_of_peaks; i++)
      peak[i][geno_size+1]= standardheight;


  if (standardwidth <= 0.0)
    for (i=0; i< number_of_peaks; i++)
      peak[i][geno_size]= (maxwidth-minwidth)*movrand()+minwidth;
  else
    for (i=0; i< number_of_peaks; i++)
      peak[i][geno_size]= standardwidth;


  if(calculate_average_error){
    global_max = -100000.0;
    for (i=0;i<number_of_peaks; i++){
      for (j=0; j<geno_size; j++)
        coordinates[j]=peak[i][j];
      dummy = dummy_eval(coordinates);
	  if (dummy>global_max)
        global_max = dummy;
		}
	}
}
/* whenever this function is called, the peaks are changed */ 
void change_peaks() 
{
  int i,j;
  double sum, sum2, offset, dummy;

  for(i=0; i<number_of_peaks; i++)
    {
      /* shift peak locations */
    sum = 0.0;
    for (j=0; j<geno_size; j++){
      shift[j]=movrand()-0.5;
      sum += shift[j]*shift[j];
      }
    if(sum>0.0)                 
      sum = vlength/sqrt(sum);
    else                           /* only in case of rounding errors */
      sum = 0.0;
    sum2=0.0;
    for (j=0; j<geno_size; j++){
      shift[j]=sum*(1.0-lambda)*shift[j]+lambda*prev_movement[i][j];
      sum2 += shift[j]*shift[j];
    }
    if(sum2>0.0)                 
      sum2 = vlength/sqrt(sum2);
    else                           /* only in case of rounding errors */
      sum2 = 0.0;


    for(j=0; j<geno_size; j++){
      shift[j]*=sum2;
      prev_movement[i][j]= shift[j];
      if (((peak[i][j]+prev_movement[i][j]) < mincoordinate) || ((peak[i][j]+prev_movement[i][j]) > mincoordinate)){
				dummy = (peak[i][j]+prev_movement[i][j]- mincoordinate)/(maxcoordinate-mincoordinate);
				if (((int)dummy % 2) == 0) {
					dummy = fabs(dummy-floor(dummy));
				}else {
					dummy=1-fabs(dummy-floor(dummy));
				}
				peak[i][j] = mincoordinate + (maxcoordinate - mincoordinate)* dummy;
			}
      else
        peak[i][j] += prev_movement[i][j];    
		}

    /* change peak width */
    j = geno_size;
    offset = movnrand()*width_severity;
    if ((peak[i][j]+offset) < minwidth)
      peak[i][j] = 2.0*minwidth-peak[i][j]-offset;
    else if ((peak[i][j]+offset) > maxwidth)
      peak[i][j]= 2.0*maxwidth-peak[i][j]-offset;
    else
      peak[i][j] += offset;
   /* change peak height */
    j++;
    offset = height_severity*movnrand();
    if ((peak[i][j]+offset) < minheight)
      peak[i][j] = 2.0*minheight-peak[i][j]-offset;
    else if ((peak[i][j]+offset) > maxheight)
      peak[i][j]= 2.0*maxheight-peak[i][j]-offset;
    else
      peak[i][j] += offset;    
  }
  if(calculate_average_error){
    global_max = -100000.0;
    for (i=0;i<number_of_peaks; i++){
      for (j=0; j<geno_size; j++)
        coordinates[j]=peak[i][j];
      dummy = dummy_eval(coordinates);
      if (dummy>global_max){
        global_max = dummy;
        maximum_peak = i;
      }
    }
  }
  recent_change = 1;
}