/* ----------------------------------------------------------------------------
  ex1.C
  mbwall 28jul94
  Copyright (c) 1995-1996  Massachusetts Institute of Technology

 DESCRIPTION:
   Example program for the SimpleGA class and 2DBinaryStringGenome class.
This program tries to fill the 2Dgenome with alternating 1s and 0s. 
  This example uses the default crossover (single point), default mutator
(uniform random bit flip), and default initializer (uniform random) for the
2D genome.
  Notice that one-point crossover is not necessarily the best kind of crossover
to use if you want to generate a 'good' genome with this kind of objective 
function.  But it does work.
---------------------------------------------------------------------------- */
#include <stdio.h>
#include <iostream.h>
#include <math.h>

#include <ga/GASimpleGA.h>	// we're going to use the simple GA
#include <ga/GA2DBinStrGenome.h> // and the 2D binary string genome

float Objective(GAGenome &);	// This is the declaration of our obj function.
				// The definition comes later in the file.

int
main(int argc, char **argv)
{
  cout << "Example 1\n\n";
  cout << "This program tries to fill a 2DBinaryStringGenome with\n";
  cout << "alternating 1s and 0s using a SimpleGA\n\n"; cout.flush();

// See if we've been given a seed to use (for testing purposes).  When you
// specify a random seed, the evolution will be exactly the same each time
// you use that seed number.

  for(int ii=1; ii<argc; ii++) {
    if(strcmp(argv[ii++],"seed") == 0) {
      GARandomSeed((unsigned int)atoi(argv[ii]));
    }
  }

// Declare variables for the GA parameters and set them to some default values.

  int width    = 10;
  int height   = 2;
  int popsize  = 30;
  int ngen     = 3000;
  float pmut   = 0.001;
  float pcross = 0.9;

  /*GAParameterList params;
  GASimpleGA::registerDefaultParameters(params);
  params.set(gaNscoreFilename, "bog.dat");
  params.set(gaNflushFrequency, 50);*/

// Now create the GA and run it.  First we create a genome of the type that
// we want to use in the GA.  The ga doesn't operate on this genome in the
// optimization - it just uses it to clone a population of genomes.

  GA2DBinaryStringGenome genome(width, height, Objective);

// Now that we have the genome, we create the genetic algorithm and set
// its parameters - number of generations, mutation probability, and crossover
// probability.  And finally we tell it to evolve itself.

  GASimpleGA ga(genome);
  //ga.parameters(params);
  ga.populationSize(popsize);
  ga.nGenerations(ngen);
  ga.pMutation(pmut);
  ga.pCrossover(pcross);
  ga.scoreFilename("test1.dat");
  ga.scoreFrequency(10);	// keep the scores of every 10th generation
  ga.flushFrequency(50);	// dump scores to disk every 50th generation
  ga.evolve();
// Now we print out the best genome that the GA found.
  cout << ga.statistics() << endl;                   

  cout << "The GA found:\n" << ga.statistics().bestIndividual() << "\n";
  
  genome = ga.statistics().bestIndividual();
  cout << "the ga generated:\n";
  for(int j=0; j<height; j++){
    for(int i=0; i<width; i++)
      cout << (genome.gene(i,j) == 1 ? '*' : ' ') << " ";
    cout << "\n";
  }

// That's it!
  return 0;
}
 

/*long DecodeChromeSome(char *string,int point,int length)
{
	long decimal=0L;
	char *pointer;
	int i;

	for (i=0,pointer=string+point;i<length;i++,pointer++)
	{
		decimal+=(*pointer-'0')<<(length-1-i);
	}
	return (decimal);
}*/

// This is the objective function.  All it does is check for alternating 0s and
// 1s.  If the gene is odd and contains a 1, the fitness is incremented by 1.
// If the gene is even and contains a 0, the fitness is incremented by 1.  No
// penalties are assigned. 
//   We have to do the cast because a plain, generic GAGenome doesn't have 
// the members that a GA2DBinaryStringGenome has.  And it's ok to cast it
// because we know that we will only get GA2DBinaryStringGenomes and
// nothing else.

float
Objective(GAGenome& g) {
  GA2DBinaryStringGenome & genome = (GA2DBinaryStringGenome &)g;
  float amin=-2.048,amax=2.048;
  float score=0.0;
  int count=0;
  float temp[2],x[2];
  int width=10,height=2;
  /*int BinaryGenome[10][2];
  for (int i=0; i<height; i++)
  {
	BinaryGenome[0][i]=genome.gene(0,i);
	for(int j=1;j<width;j++)
	{
		if(BinaryGenome[j-1][i]==genome.gene(j,i))
			BinaryGenome[j][i]=0;
		else
			BinaryGenome[j][i]=1;
	}
  }*/
  for (int i=0;i<height;i++)
  {
		temp[i]=0;
		for (int j=0;j<width;j++)
		//	temp[i]+=pow(2,width-j-1)*BinaryGenome[j][i];
			temp[i]+=pow(2,width-j-1)*genome.gene(j,i);
  }

	x[0]=(amax-amin)*temp[0]/(pow(2,width)-1)+amin;
	x[1]=(amax-amin)*temp[1]/(pow(2,width)-1)+amin;
	score=100*(x[0]*x[0]-x[1])*(x[0]*x[0]-x[1])+(1-x[0])*(1-x[0]);
	//score=100*(x[0]*x[0])+x[1]*x[1];
   

  /*for(int i=0; i<genome.width(); i++){
		//score+=genome.gene(i,0)+genome.gene(i,1);
		//score+=genome.gene(i,1);
    for(int j=0; j<genome.height(); j++){
      if(genome.gene(i,j) == 0 && count%2 == 0)
	score += 1.0;
      if(genome.gene(i,j) == 1 && count%2 != 0)
	score += 1.0;
      count++;
    }
  }*/;
  return score;
}

/*void GlaytoBinary(char *GlayString,char *BinaryString,int length)
{
   *BinaryString=*GlayString;
   for(int i=1;i<length;i++,BinaryString++,GlayString++)
	   if(*BinaryString==*(GlayString+1))
		   *(BinaryString+1)=0;
	   else
		   *(BinaryString+1)=1;
}

*/