//////////////////////////////////////////////////////////////////////
//
// 8.cpp
// for model of leakage 
//                       aka pampers =)
// issue N8
// Started: 16.07.03 9:10:12
//
//                                                    XIII
//////////////////////////////////////////////////////////////////////
// * - is wet cell 
// a - will wet
#ifndef _8_CPP_
#define _8_CPP_ 1

#include<fstream.h>
#include<stdio.h> // for getchar()
#include<stdlib.h> // for rand()
char *prog_name="pampers v0.02";
// size of massive
#define P_SIZE 20
#define P_SIZEY 60
const double p_LeakCoef=0.013; // coefficient of leakage 
							   // try to change it
typedef unsigned int uint;
// variables
char *file_name="pampers.txt";
fstream fs; // for manipulations with file
char p_mas[P_SIZE+2][P_SIZEY+2]; //massive
int p_i,p_j; // counter's
uint p_time=0; // age of model
uint p_coef=0; // coefficient for check cell
// functions
int in_file(char*); // get data from file and set massive
void calculate(uint,uint); // calculate one cell 
int main()
{
	// print Hello :)
	cout<<"+---------------+\n";
	cout<<"|  Pampers :)   |\n";
	cout<<"|               |\n";
	cout<<"|          XIII |\n";
	cout<<"+---------------+"<<endl;
	cout<<"press Enter for next step; for Exit press q"<<endl;
	// set border
	for(p_i=0;p_i<P_SIZEY+2;p_i++)
		p_mas[0][p_i]=p_mas[P_SIZE+1][p_i]='-';
	for(p_j=0;p_j<P_SIZE+2;p_j++)
		p_mas[p_j][0]=p_mas[p_j][P_SIZEY+1]='|';
	p_mas[P_SIZE+1][0]=p_mas[0][P_SIZEY+1]='+';
	p_mas[0][0]=p_mas[P_SIZE+1][P_SIZEY+1]='+';
	// set elements of p_mas zero
	for(p_i=1;p_i<P_SIZE+1;p_i++)
	{
		for(p_j=1;p_j<P_SIZEY+1;p_j++)
			p_mas[p_i][p_j]='0';
	}
	// get data
	if(in_file(file_name))
		return 1; // error open file
	while(getchar()!='q')
	{
		cout<<"Time: "<<p_time++<<endl;
		// print massive
		for(p_i=0;p_i<P_SIZE+2;p_i++)
		{
			for(p_j=0;p_j<P_SIZEY+2;p_j++)
			{
				if(p_mas[p_i][p_j]=='a')
					p_mas[p_i][p_j]='*'; // wet
				if(p_mas[p_i][p_j]=='0')
					cout<<" ";
				else
					cout<<p_mas[p_i][p_j];
			}
			cout<<endl;
		}
		// calculate next step
		for(p_i=1;p_i<P_SIZE+1;p_i++)
		{
			for(p_j=1;p_j<P_SIZEY+1;p_j++)
				calculate(p_i,p_j);
		}
	}
	////////////
	return 0;
}
int in_file(char*file)
{// get data from file and set life-massive
	fs.open(file,ios::in);
	if(!fs)
	{// can't open file
		cout<<"Error open file!\t"<<file_name<<endl;
		return 1;
	}
	// get data
	for(p_i=1;p_i<P_SIZE+1;p_i++)
	{
		for(p_j=1;p_j<P_SIZEY+1;p_j++)
			fs>>p_mas[p_i][p_j];
	}
	fs.close();
	return 0;
}
void calculate(uint x,uint y)
{// calculate one cell with x,y
	switch(p_mas[x][y])
	{
	case 'a':
	case '+':
	case '|':
	case '-':
	case '*':
	break;
	case '0':// only for dry cell
		//-------------------------------------//
		p_coef=0;
		if(p_mas[x][y-1]=='*') // west
			p_coef++;
		if(p_mas[x-1][y-1]=='*')// north-west
			p_coef++;
		if(p_mas[x-1][y]=='*')// north
			p_coef++;
		if(p_mas[x-1][y+1]=='*')// north-east
			p_coef++;
		if(p_mas[x][y+1]=='*')// east
			p_coef++;
		if(p_mas[x+1][y+1]=='*')// south-east
			p_coef++;
		if(p_mas[x+1][y]=='*')// south
			p_coef++;
		if(p_mas[x+1][y-1]=='*')// south-west
			p_coef++;
		//
		// take pseudorandom number in the range 0 to 1.
		// yes here very simple case; in reality,  
		// here must be equation with exp() ; i think so =)
		//
		if((double)rand()/RAND_MAX<p_LeakCoef*p_coef)
				p_mas[x][y]='a';
		//-------------------------------------//
		break;
	default:
		break;
	}
}
#endif