// Henon.cpp: implementation of the Henon class.


//


//////////////////////////////////////////////////////////////////////





//#include ""


#include "Henon.h"


#include "neural.h" //for myrand()


#include<iostream.h>	//for debugging purposes!


//////////////////////////////////////////////////////////////////////


// Construction/Destruction


//////////////////////////////////////////////////////////////////////











Henon::Henon()


{


	b=0.3;


	c=1.4;


	x=0.1;


	y=0.1;





}





Henon::Henon(long double in_b, long double	in_c){





	b=in_b;


	c=in_c;


	x=0.1;


	y=0.1;








}





Henon::Henon(long double x0,long double y0,long double in_b,long double in_c){





	x=x0;


	y=y0;


	b=in_b;


	c=in_c;





}





long double inline Henon::henonx(long double x, long double y, long double b, long double c){








		return (1-c*x*x+b*y);





}








long double inline Henon::henony(long double x){





		return x;


}








void Henon::advance(){





	tempx=henonx(x,y,b,c);


	tempy=henony(x);


	x=tempx;


	y=tempy;


}





//the obligatory overloaded operator





Henon& Henon::operator=(Henon& in){





	x=in.x;
	y=in.y;


	b=in.b;


	c=in.c;





	return *this;





}





void Henon::discard(){








	int discards=400;


	for(int i=0;i<discards;i++){


		advance();


	}


}





void Henon::randomize(){





	x=myrand();


	y=myrand();








}








Henon::~Henon()


{





}








DecoupledMap::DecoupledMap(long n){


	num_maps=n;


	long i,j,k;


	map=new long double*[num_maps*2];


	henonmaps=new Henon[num_maps];


	for(i=0;i<num_maps;i++){


		if(i==12){
			cout<<"\nfor some reason, we run into trouble with the henon map 12"<<endl;
		
		
		}
		henonmaps[i].randomize();

	
		//??? 4/17/1998
		//ONLY when i==12, discard somehow sets henonmaps[12].x= - INFINITY
		//I do NOT get it...

		henonmaps[i].discard();


	}


	j=0;

	for(i=0;i<num_maps;i++){


		map[j]=&(henonmaps[i].x); //we can directly modify whatever a given map contains by modifying *map[i]


		map[j+1]=&(henonmaps[i].y);

		j+=2;


	}


}





DecoupledMap::~DecoupledMap(){


	delete[]map;


	delete[]henonmaps;


}





void DecoupledMap::advance(){


	for(long i=0;i<num_maps;i++){


		henonmaps[i].advance();


	}


}



//assumes stuff is already allocated!!!

DecoupledMap& DecoupledMap::operator=(DecoupledMap & in){





	long i,j,k;


	num_maps=in.num_maps;


	if(!map){

		cout<<"\nDecupledMap::opererator=() failed due to attempted assignment to null variable!!!"<<endl;

		//return NULL;

	}

	//delete[]map;


	//map=new long double*[num_maps*2];


	//delete[] henonmaps;


	//henonmaps=new Henon[num_maps];


	for(i=0;i<num_maps;i++){


		henonmaps[i]=in.henonmaps[i];


	}


	j=0;

	for(i=0;i<num_maps;i++){




		map[j]=&(henonmaps[i].x); //we can directly modify whatever a given map contains by modifying *map[i]


		map[j+1]=&(henonmaps[i].y);

		j+=2;


	}





	return *this;





}