#include <stdio.h>
#include <stdlib.h>
#include <GL/glut.h>
#include <GL/gl.h>
#include <sys/time.h>


void OpenGLInit(void);
static void Animate(void );
static void ResizeWindow(int w, int h);
struct timespec ts;
struct timeval starttime,endtime;

double fps = 0;
double time1 = 0;
double time2 = 0;


int rec_center_x = 0;
int rec_center_y = 0;
int rec_width = 100;
int rec_length = 100;
int rec_rotation = 0;
int uneven = 50;
int K, J;
int X;
int Y;  
int RES;
int DEPTH = 3;


struct screen {

	unsigned char red;
	unsigned char green;
	unsigned char blue;
	unsigned char alpha;	
};


//struct screen *mainwindow=( (struct screen *) malloc(sizeof(struct screen)*5) );
struct screen mainwindow[9437184];

double x_coord = 0;

static void Animate(void)
{
	double i;
	int x1 = 0;
	int y1 = 0;
	int z1 = 10;
	float rcalc;
	float gcalc;
	float bcalc;
	float acalc;

	glClear(GL_COLOR_BUFFER_BIT);

	x1 = 0;
	y1 = 0;
	
	/*
	glBegin(GL_POINTS);
	glColor3f(255, 0, 0);
	glVertex3i(50, 50, 1);
        glVertex3i(51, 50, 1);
        glVertex3i(51, 51, 1);
        glVertex3i(50, 51, 1);
	glEnd();

        glBegin(GL_POINTS);
        glColor3f(255*.5, 0, 255);
        glVertex3i(50, 50, 0);
	//glColor4f(0, 0, 255, 1); 
        glVertex3i(51, 50, 0);
        glVertex3i(51, 51, 0);
        glVertex3i(50, 51, 0);
        glEnd();
	*/


	z1 = 2;
	double p1, p2, p3, a1, a2;
	double t1, t2;
	long int elem1, elem2;
	for(i=0; i < RES*DEPTH; i+=X) {
		glBegin(GL_POINTS);
				
		for(x1=0; x1<X; x1++) {
			
			//if(mainwindow[x1+(y1*X)+(z1*RES)].red != 0 && mainwindow[x1+(y1*X)+(z1*RES)].green != 0 && 
                          //mainwindow[x1+(y1*X)+(z1*RES)].blue != 0) {
			if(z1 != 3 ) {
				if(z1 == 0) {
					elem1 = x1+(y1*X); 
					elem2 = RES+x1+(y1*X);
				}else if(z1 == 1){
					elem1 = RES+x1+(y1*X);
					elem2 =(2*RES)+x1+(y1*X);
				}else{
					elem1 =(2*RES)+x1+(y1*X);
					elem2 =(2*RES)+x1+(y1*X);
				}

				t1 = mainwindow[elem1].red+mainwindow[elem1].green+
                                    mainwindow[elem1].blue;
//if(t1 != 0) {
				a2 =  (255 - mainwindow[elem1].alpha);
				
				t2 = mainwindow[elem2].red+mainwindow[elem2].green+
                                    mainwindow[elem2].blue;

				p1 = (mainwindow[elem2].red/(t2+1))*a2;
				p2 = (mainwindow[elem2].green/(t2+1))*a2;
				p3 = (mainwindow[elem2].blue/(t2+1))*a2;

				//p1 = p1 *(a2/255);
				//p2 = p2 *(a2/255);
				//p3 = p3 *(a2/255);

				rcalc = (mainwindow[elem1].red/(t1+1))*(mainwindow[elem1].alpha); 
				gcalc = (mainwindow[elem1].green/(t1+1))*(mainwindow[elem1].alpha);
				bcalc = (mainwindow[elem1].blue/(t1+1))*(mainwindow[elem1].alpha);

				rcalc = rcalc + p1;
				gcalc = gcalc + p2;
				bcalc = bcalc + p3;

				rcalc = rcalc / 255;
				gcalc = gcalc / 255;
				bcalc = bcalc / 255;


				/*
                                scalc = mainwindow[x1+(y1*X)+(z1*RES)].red+
                                        mainwindow[x1+(y1*X)+(z1*RES)].green+mainwindow[x1+(y1*X)+(z1*RES)].blue;
				scalc = scalc*(mainwindow[x1+(y1*X)+((z1)*RES)].alpha/255);
				

				//result = ( ALPHA * ( srcPixel - destPixel ) ) / 256 + destPixel
				rcalc = (mainwindow[x1+(y1*X)+((z1)*RES)].alpha*
                                    (mainwindow[x1+(y1*X)+(z1*RES)].red - mainwindow[x1+(y1*X)+((z1+1)*RES)].red) ) /
                                    256 + mainwindow[x1+(y1*X)+((z1+1)*RES)].red;
				gcalc = (mainwindow[x1+(y1*X)+((z1)*RES)].alpha*
                                    (mainwindow[x1+(y1*X)+(z1*RES)].green - mainwindow[x1+(y1*X)+((z1+1)*RES)].green) ) /
                                    256 + mainwindow[x1+(y1*X)+((z1+1)*RES)].green;
				bcalc = (mainwindow[x1+(y1*X)+((z1)*RES)].alpha*
                                    (mainwindow[x1+(y1*X)+(z1*RES)].blue - mainwindow[x1+(y1*X)+((z1+1)*RES)].blue) ) /
                                    256 + mainwindow[x1+(y1*X)+((z1+1)*RES)].blue;

				
				rcalc = (mainwindow[x1+(y1*X)+(z1*RES)].red*mainwindow[x1+(y1*X)+(z1*RES)].alpha/255)+
                                        ( (mainwindow[x1+(y1*X)+((z1+1)*RES)].red )*
                                        (mainwindow[x1+(y1*X)+((z1)*RES)].alpha/255) ) ;
				rcalc = rcalc / mainwindow[x1+(y1*X)+(z1*RES)].alpha;

                                gcalc = (mainwindow[x1+(y1*X)+(z1*RES)].green*mainwindow[x1+(y1*X)+(z1*RES)].alpha/255)+
                                        ( (mainwindow[x1+(y1*X)+((z1+1)*RES)].green )*
                                        (mainwindow[x1+(y1*X)+((z1)*RES)].alpha/255) ) ;
                                gcalc = gcalc / mainwindow[x1+(y1*X)+(z1*RES)].alpha;

                                bcalc = (mainwindow[x1+(y1*X)+(z1*RES)].blue*mainwindow[x1+(y1*X)+(z1*RES)].alpha/255)+
                                        ( (mainwindow[x1+(y1*X)+((z1+1)*RES)].blue )*
                                        (mainwindow[x1+(y1*X)+((z1)*RES)].alpha/255) ) ;
                                bcalc = bcalc / mainwindow[x1+(y1*X)+(z1*RES)].alpha;

				*/



				/*
				glColor3f( (mainwindow[x1+(y1*X)+((z1+1)*RES)].red * mainwindow[x1+(y1*X)+((z1+1)*RES)].alpha)
                                            + (mainwindow[x1+(y1*X)+(z1*RES)].red * mainwindow[x1+(y1*X)+(z1*RES)].alpha),
                                            (mainwindow[x1+(y1*X)+((z1+1)*RES)].green * mainwindow[x1+(y1*X)+((z1+1)*RES)].alpha)
                                            + (mainwindow[x1+(y1*X)+(z1*RES)].green * mainwindow[x1+(y1*X)+(z1*RES)].alpha), 
                                            (mainwindow[x1+(y1*X)+((z1+1)*RES)].blue * mainwindow[x1+(y1*X)+((z1+1)*RES)].alpha)
                                            + (mainwindow[x1+(y1*X)+(z1*RES)].blue * mainwindow[x1+(y1*X)+(z1*RES)].alpha)
                                          );
				*/
				//glColor3f(rcalc, gcalc, bcalc);
//}else{
//glColor3f(0, 0, 0);
//}

				//glVertex2f(x1, y1);
			//}else{
				//glColor4f(mainwindow[x1+(y1*X)+((z1)*RES)].red, mainwindow[x1+(y1*X)+((z1)*RES)].green, 
                                //mainwindow[x1+(y1*X)+((z1)*RES)].blue, mainwindow[x1+(y1*X)+(z1*RES)].alpha);
			}
			//glLoadIdentity();
t1 = rcalc + gcalc + bcalc;
if(t1 != 0) {

			glColor3f(rcalc, gcalc, bcalc);
			glVertex3f(x1+.5, y1+.5, z1);
//}else{
//glColor3f(0, 0, 0);
//glVertex3f(x1+.5, y1+.5, z1);
}		

		}

		if(y1 == Y && x1 == X ) {
			y1 = 0;
			if(z1 > 0) {
				z1--;
			}

			//printf("here\n");
		}else{
			//if(y1 < Y) {
				y1++;
			//}else{
			//	y1 = 0;
			//}
			//z=10;
		}	


		glEnd();


	}

	glFlush();
	glutSwapBuffers();
	glutPostRedisplay();

	time2++;

	gettimeofday(&endtime, NULL);
	time1=((double)(endtime.tv_sec*1000000-starttime.tv_sec*1000000+endtime.tv_usec-starttime.tv_usec))/1000000;

        if (time2 == 5)
        {
                fps = time2 / time1;
                printf("fps: %lf - %lf frames in %lf seconds\n", fps, time2, time1);
                gettimeofday(&starttime, NULL);
                time2 = 0;

	



/*
glBegin(GL_LINES);
for(x_coord=0; x_coord<=X; x_coord+=50) {
        //glLoadIdentity();
        glColor4f( 0, 255, 0, 0);
        //glTranslatef(x_coord,0.0f,0.0f);

        //glBegin(GL_LINES);
        glVertex3f(x_coord, 0.0, 0);
        glVertex3f(x_coord, Y, 0);
        //glEnd( );
}
glEnd( );  




        glFlush();
        glutSwapBuffers();
        glutPostRedisplay();
*/


		z1 = 1;
		K = rec_center_x;
		J = rec_center_y;

		while(J < (rec_center_y + rec_length ) ) {
			while(K < (rec_center_x + rec_width ) ) {
				elem1 =(RES)+K+(X*J);
				mainwindow[elem1].red = 0;				
				mainwindow[elem1].green = 0;
				mainwindow[elem1].blue = 0;
				mainwindow[elem1].alpha = 0;
				K++;
			}
			J++;
			K = rec_center_x;
		}


		z1 = 0;
                K = rec_center_x+uneven;
                J = rec_center_y+uneven;
                 
                while(J < (rec_center_y+uneven + rec_length ) ) {
                        while(K < (rec_center_x+uneven + rec_width ) ) {
				elem1 = K+(X*J);
                                mainwindow[elem1].red = 0;
                                mainwindow[elem1].green = 0;
                                mainwindow[elem1].blue = 0;
                                mainwindow[elem1].alpha = 0;
                                K++;
                        }
                        J++;
                        K = rec_center_x+uneven;
                }
                                








		if(rec_center_x < 500) {
			rec_center_x+=25;
		} else {
			rec_center_x=0;
		}

                if(rec_center_y < 500) {
                        rec_center_y+=25;
                } else {
                        rec_center_y=0;
                }  

		if(uneven == 100) {
			uneven = 50;
		}else{
			uneven = 100;
		}





		z1 = 1;
                K = rec_center_x;
                J = rec_center_y;
                while(J < (rec_center_y + rec_length) ) {
                        while(K < (rec_center_x + rec_width) ) {
				elem1 = (RES)+K+(X*J);
                                mainwindow[elem1].red = 255;
                                mainwindow[elem1].green = 0;
                                mainwindow[elem1].blue = 0;
				mainwindow[elem1].alpha = 15;
                                K++;
                        }
                        J++;
                        K = rec_center_x;
                }


		z1 = 0;
                K = rec_center_x +uneven;
                J = rec_center_y +uneven;   
                while(J < (rec_center_y+uneven + rec_length) ) {
                        while(K < (rec_center_x+uneven + rec_width) ) {
				elem1 = K+(X*J);
                                mainwindow[elem1].red = 0;
                                mainwindow[elem1].green = 0;
                                mainwindow[elem1].blue = 255;
                                mainwindow[elem1].alpha = 85;
                                K++;
                        }
                        J++;
                        K = rec_center_x+uneven;
                }


}





}

void OpenGLInit(void)
{
	glShadeModel( GL_FLAT );
	glClearColor( 0.0, 0.0, 0.0, 0.0 );
	glClear(GL_COLOR_BUFFER_BIT);
	glDisable( GL_DEPTH_TEST );
}       
        

static void ResizeWindow(int w, int h)
{  
	float aspectRatio;
        h = (h == 0) ? 1 : h;
        w = (w == 0) ? 1 : w;
        glViewport( 0, 0, w, h );       // View port uses whole window
        aspectRatio = (float)w/(float)h;

	glMatrixMode( GL_PROJECTION );
	glLoadIdentity();
	glOrtho (0, w, h, 0, 0, 1);
	glMatrixMode( GL_MODELVIEW );
}




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

        glutInit(&argc,argv);
        glutInitDisplayMode(GLUT_DOUBLE | GLUT_RGB);
        glutInitWindowPosition( 0, 0 );

	Y = glutGet(GLUT_SCREEN_HEIGHT); 
	X = glutGet(GLUT_SCREEN_WIDTH);
	RES = X*Y;
	

	printf("res: %d x: %d y: %d\n", RES, X, Y);

//	if( (mainwindow = (screen *)realloc(mainwindow, sizeof(screen)*RES*DEPTH) ) ==NULL )
//	{
//		printf("Reallocation failed\n");
//		exit(1);
//	} 










	long int i;
	for(i=0; i < RES*DEPTH; i++) {

                        mainwindow[i].red = 0;
                        mainwindow[i].green = 0;
                        mainwindow[i].blue = 0;
			mainwindow[i].alpha = 0;
	}


	printf("res*depth %d\n", RES*DEPTH);

	mainwindow[9437183].red = 255;
	mainwindow[9437183].green = 255;
	printf("red: %d\n", mainwindow[9437183].red);

//9437184

	glutInitWindowSize( X, Y );
	glutCreateWindow( "lsrender test" );


	OpenGLInit();
	glutReshapeFunc( ResizeWindow );
	glutDisplayFunc( Animate );
	glutMainLoop(  );



	return 0;
}