//life.cpp
//(c) 1997, Oguz Yetkin
//physics 505
//
//Game of life, and other good stuff


//#include "graph.h"


#include <graphics.h>	//will have to use borland c++ graphics directly :(
					//for getpixel and putpixel
#include <conio.h>
#include <stdio.h>
#include <iostream.h>
#include <stdlib.h>		//for exit



const int GEN=2;
const int I=70;
const int J=70;
const int offset=20;	//offset to plot
//const int ITER=300;


inline int fate(int value,int num_neighbors);		//return 0 or 1

inline int count_neighbor(int array[GEN][I][J], int i, int j,int gen); //assume I and J there

int main(int argc, char** argv){

	int world[GEN][I][J];	//2 2-d arrays in a 3-d array
	int world2[GEN][I][J];	//parallel universe

	int gen=0;	//either 0 or 1
	int iter=0;
	int curr_generation=0;

	int maxcolor;
	long int i,j,k;	//for iterating up to millions
	char ch;
	cout<<"\nLife 1.0, by Oguz Yetkin, oguz@writeme.com";
	cout<<"\npress q to stop the program";
	cout<<"\nThis is a straighforward implementation of LIFE,";
	cout<<"\nThe code is probably on";
	cout<<"\nhttp://www.cs.wisc.edu/~yetkin/code\n"<<endl;
	cout<<"\n\nmax iterations: ";
	cin>>iter;


	int maxx,maxy;

	//initialize graphics
	//************************************
	int driver=DETECT;
	int mode;
	registerbgidriver(EGAVGA_driver);
	//int maxx, maxy;	//will be 640x480


	//**************INIT GRAPHICS**********************
	initgraph(&driver, &mode, "");	//driver linked in
	//mode=CGAC0;
	//initgraph(&driver, &mode,"h:\\borlandc\\bgi");

	maxx=getmaxx();
	maxy=getmaxy();
	maxcolor=getmaxcolor();
	//cout<<maxcolor<<endl;
	//getch();

	//**************END GRAPHICS SETUP*****************


	//zero out worlds

	for(gen=0;gen<GEN;gen++){
		for(i=0;i<I;i++){
			for(j=0;j<J;j++){
				world[gen][i][j]=0;
				world2[gen][i][j]=0;
			}
		}
	}

	randomize();
	//fill main world with random nums;
	for(i=0;i<I;i++){
		for(j=0;j<J;j++){
			world[0][i][j]=random(2);
		}
	}
	//copy world into parallel universe

	for(i=0;i<I;i++){
		for(j=0;j<J;j++){
			world2[0][i][j]=world[0][i][j];
		}
	}
	//now alter one point in the parallel universe
	world2[0][15][25]=!world2[0][15][25];


	//start LIFE
	for(curr_generation=0;curr_generation<iter;curr_generation++){

		gen=curr_generation%2;

		for(i=0;i<I;i++){
			for(j=0;j<J;j++){
				//put whatever we want in next generation's array
				world[((gen+1)%2)][i][j]=
					fate(world[gen][i][j],
						count_neighbor(world,i,j,gen));

				world2[((gen+1)%2)][i][j]=
					fate(world2[gen][i][j],
						count_neighbor(world2,i,j,gen));
			}//end inner for
		}//end outer for

		//cout<<"\ngen: "<<curr_generation<<endl;
		for(i=0;i<I;i++){
			for(j=0;j<J;j++){

				if(world[((gen+1)%2)][i][j]==1){
					putpixel(i,j,WHITE);

				}else{
					putpixel(i,j,BLACK);
				}//end if,then,else
			}//end inner for
		}//end outer for

		//paint the parallel universe
		for(i=0;i<I;i++){
			for(j=0;j<J;j++){

				if(world2[((gen+1)%2)][i][j]==1){
					putpixel(i+I+offset,j,WHITE);

				}else{
					putpixel(i+I+offset,j,BLACK);
				}//end if,then,else
			}//end inner for
		}//end outer for
		//clrscr();
	}//end huge for loop

	getch();
	closegraph();
	return 1;

}//end main



inline int fate(int value,int num_neighbors){



	if(value==1){

		if((num_neighbors==2)||(num_neighbors==3)){
			return 1;	//stays alive
		}else{
			return 0;	//dies of loneliness or overcrowding
		}
	}else{

		if(num_neighbors==3){
			return 1;		//we have a new baby!!!
		}else{
			return 0;		//stay dead
		}
	}//end outer else


	cout<<"\nwe should never get here in fate()!!!\n"<<endl;
	return 0;


}



inline int count_neighbor(int array[GEN][I][J], int i, int j,int gen){


	int current=array[gen][i][j];
	//don't want any negatives
	int num;
	num=array[gen][((I+i+1)%I)][j]+array[gen][((I+i-1)%I)][j]	//W,E
		+array[gen][i][((J+j+1)%J)]+array[gen][i][((J+j-1)%J)] //N,S
		+array[gen][((I+i+1)%I)][((J+j+1)%J)]+array[gen][((I+i+1)%I)][((J+j-1)%J)]	//NE,SE
		+array[gen][((I+i-1)%I)][((J+j+1)%J)]+array[gen][((I+i-1)%I)][((J+j-1)%J)];
	return num;

}