#include "glos.h"
#include <GL/gl.h>
#include <GL/glaux.h>
#include "windows.h"
#include "math.h"
#include "stdio.h"
#include <fstream>
#include <iostream>
#include <strstream>//#include "../require.h"void myinit(void);
using namespace std;
#define SZ 500 // Buffer size
void CALLBACK  display(void);
void CALLBACK  reshape(GLsizei w,GLsizei h);
void CALLBACK  left(void);//按 LEFT
void CALLBACK  right(void);//按 RIGHT
void CALLBACK  up(void);//按 UP
void CALLBACK  down(void);//按 DOWN
void CALLBACK  big(void);//按 b
void CALLBACK  small1(void);//按 s



//////////////////////////////////////
#define MAXFILE 6000
#define MINFILE 0
#define STEPFILE 25

#define XFACTOR .05
#define YFACTOR .05
#define ZFACTOR .05

#define CRYSTCONST 5.65315*XFACTOR
//////////////////////////////////////




static int z_motion=0,rotate=0,yinzi=1;
int maxfile,minfile,stepfile,totalpoint=200;
float enlargerate;
/*
GLfloat points[4][4][3]={
        {{-1.5,-1.5,2.0},{-0.5,-1.5,2.0},
        {0.5,-1.5,-1.0},{1.5,-1.5,2.0}},
        {{-1.5,-0.5,1.0},{-0.5,1.5,2.0},
        {0.5,0.5,1.0},{1.5,-0.5,-1.0}},
        {{-1.5,0.5,2.0},{-0.5,0.5,1.0},
        {0.5,0.5,3.0},{1.5,-1.5,1.5}},
        {{-1.5,1.5,-2.0},{-0.5,1.5,-2.0},
        {0.5,0.5,1.0},{1.5,1.5,-1.0}}};
*/
int datanew=0;
/////////////////////////////
//	MATRIX SIZE!	
//	SHOULD BE MODIFIED 
	GLfloat temp[3200][5];
/////////////////////////////
GLfloat color[4][3]={
        {1.0,0.0,0.0},{0.0,1.0,0.0},{0.0,0.0,1.0},{1.0,1.0,1.0}};


//convert!
void convert(int ii,char* filename)
{
	int j,k,i;
	i=ii;
	filename[0]='P';
	filename[1]='O';
	filename[2]='S';
	filename[8]='.';
	filename[9]='D';
	filename[10]='A';
	filename[11]='T';
	filename[12]='\0';
	filename += 7 ;
	for(j=0;j<5;j++)
	{
		k=fmod(i,10);
		i=(i-k)/10;
		*filename = k+'0'-0;
		filename--;
	}
	return;
}
		
void getdata(int snoffile)
{

	float x,y,z,sn,label;
	float *p;
	
	char buf[SZ],filename[13];

	convert(snoffile,filename);
	fstream in(filename);
	p=(float *)temp;
	while(in.getline(buf, SZ)) { // Removes \n
	istrstream s(buf);
	s>>sn>>label>>x>>y>>z;
	*p=(float)sn;p++;
	*p=label;	 p++;
	*p=(x-0)*XFACTOR;p++;
	*p=(y-0)*YFACTOR;p++;
	*p=(z-0)*ZFACTOR;p++;
	//noofdata++;
	}
	
	
	
	
	
	cout<<"file " <<filename<< " input successfully!"<<"\n";
//	cout<<noofdata<<"input successfully!"<<"\n";
	return;
}
		
//初始化光照、材质的过程


void CALLBACK left(void)
{
//每当按下LEFT，rotate数值 +1
//以下函数都类似
rotate++;
}
void CALLBACK big(void)
{
//每当按下LEFT，rotate数值 +1
//以下函数都类似
yinzi++;
}
void CALLBACK small1(void)
{
//每当按下LEFT，rotate数值 +1
//以下函数都类似
	yinzi--;
}
void CALLBACK right(void)
{
rotate--;
}
void CALLBACK up(void)
{
z_motion++;
}
void CALLBACK down(void)
{
z_motion--;
}

		void initlights(void)
{
        GLfloat ambient[]={0.4,0.6,0.2,1.0};
        GLfloat position[]={.0,.2,1.0,0.0};
        glEnable(GL_LIGHTING);
        glEnable(GL_LIGHT0);
        glLightfv(GL_LIGHT0,GL_AMBIENT,ambient);
        glLightfv(GL_LIGHT0,GL_POSITION,position);
}

void myinit(void)
{

        auxInitDisplayMode(AUX_SINGLE|AUX_RGBA);
        auxInitPosition(0,0,500,500);
        auxInitWindow("sample1");
        glClearColor(0.0,0.5,0.0,0.3);
        glClear(GL_COLOR_BUFFER_BIT);
        //glEnable(GL_MAP2_VERTEX_3);
        glEnable(GL_AUTO_NORMAL);

      glEnable(GL_NORMALIZE);
        //glMapGrid2f(50,0.0,1.0,50,0.0,1.0);

        glEnable(GL_DEPTH_TEST);
//初始化光照、材质
        initlights();
//      glShadeModel(GL_FLAT);



}

void CALLBACK reshape(GLsizei w,GLsizei h)
{

glViewport(0,0,w,h);
glMatrixMode(GL_PROJECTION);
glLoadIdentity();

if(w<=h)
 glOrtho(-5.0,5.0,-5.0*(GLfloat)h/(GLfloat)w,
         5.0*(GLfloat)h/(GLfloat)w,-5.0,5.0);
else
  glOrtho(-5.0*(GLfloat)h/(GLfloat)w,
         5.0*(GLfloat)h/(GLfloat)w,-5.0,5.0,-5.0,5.0);
glMatrixMode(GL_MODELVIEW);
glLoadIdentity();
}


void CALLBACK display(void)
{
glClear(GL_COLOR_BUFFER_BIT|GL_DEPTH_BUFFER_BIT);
glColor3f(0.0,0.0,1.0);
//glScalef(1+yinzi/10.,1+yinzi/10.,1+yinzi/10.);
glPushMatrix();
//glRotatef(40+10*rotate,1.0,1.0,.0);
glRotatef(90*fmod(yinzi,2),.0,1.,.0);
glTranslatef(-.0,-.0+z_motion,0.0);	//glScalef(2,2,1);
GLfloat mat_specular[]={1.0,1.0,1.0,1.0};
GLfloat mat_shininess[]={80.0};
glMaterialfv(GL_FRONT,GL_SPECULAR,mat_specular);
glMaterialfv(GL_FRONT,GL_SHININESS,mat_shininess);

GLfloat mat_diffusered[]={0.9,0.1,0.1,1.0};
GLfloat mat_diffuseblack[]={0.1,0.1,0.1,1.0};
GLfloat mat_diffuseblue[]={0.1,0.1,0.9,1.0};

glMaterialfv(GL_FRONT,GL_DIFFUSE,mat_diffusered);
glEnable(GL_MAP2_VERTEX_3);
glEnable(GL_AUTO_NORMAL);

glTranslatef(-.8,-.8,-.8);	//glScalef(2,2,1);

//===========================
glPushMatrix();
glScalef(enlargerate,enlargerate,enlargerate);


//getdata(1750);
glPushMatrix();	
glTranslatef(-.706644*sqrt(totalpoint/200),-.706644*sqrt(totalpoint/200),-.706644*sqrt(totalpoint/200));

//display();
int filei;
for(filei=minfile;filei<=maxfile;filei=filei+stepfile)
{
getdata(filei);
glClear(GL_COLOR_BUFFER_BIT|GL_DEPTH_BUFFER_BIT);

glPushMatrix();	
glTranslatef(.706644*sqrt(totalpoint/200),.706644*sqrt(totalpoint/200),.706644*sqrt(totalpoint/200));
auxWireCube(1.4132875*sqrt(totalpoint/200));
glPopMatrix();

//       draw the next layer
/*
for (int ii=0;ii<5;ii++)
for(int j=0;j<5;j++)
{
glPushMatrix();	
glTranslatef(CRYSTCONST*.25,CRYSTCONST*(ii+.25),CRYSTCONST*(j+.25));
glMaterialfv(GL_FRONT,GL_DIFFUSE,mat_diffuseblue);
auxSolidSphere(.04);
glPopMatrix();
 
glPushMatrix();	
glTranslatef(CRYSTCONST*.25,CRYSTCONST*(ii+.75),CRYSTCONST*(j+.75));
glMaterialfv(GL_FRONT,GL_DIFFUSE,mat_diffuseblue);
auxSolidSphere(.04);
glPopMatrix();

glPushMatrix();	
glTranslatef(CRYSTCONST*.75,CRYSTCONST*(ii+.25),CRYSTCONST*(j+.75));
glMaterialfv(GL_FRONT,GL_DIFFUSE,mat_diffuseblue);
auxSolidSphere(.02);
glPopMatrix();

glPushMatrix();	
glTranslatef(CRYSTCONST*.75,CRYSTCONST*(ii+.75),CRYSTCONST*(j+.25));
glMaterialfv(GL_FRONT,GL_DIFFUSE,mat_diffuseblue);
auxSolidSphere(.02);
glPopMatrix();

}
*/


for(int i=0;i<totalpoint;i++)
{ 
float x,y,z;
x=temp[i][2];y=temp[i][3];z=temp[i][4];
glPushMatrix();	
glTranslatef(x,y,z);
if(temp[i][1]==1) 
{
	glMaterialfv(GL_FRONT,GL_DIFFUSE,mat_diffusered);
	if((int)fmod(temp[i][0],8)<2.5)
	auxSolidSphere(.06);
	else auxSolidSphere(.04);
}
else if(temp[i][1]==2)
{
	glMaterialfv(GL_FRONT,GL_DIFFUSE,mat_diffuseblue);
	if((int)fmod(temp[i][0],8)<5.5)
	auxSolidSphere(.05);
	else auxSolidSphere(.03);
}

glPopMatrix();
}

for(int ij=0;ij<10000;ij++)
for(int ik=0;ik<1000;ik++);

}


glPopMatrix();


glPopMatrix();
glPopMatrix();

  glFlush();


}
void main(void)
{

	float x,y,z,sn,label;
	float *p;
	int snoffile=0;
	printf("Pls input max,min,step,totalpoints,enlargerate\n");
	scanf("%d %d %d %d %f",&maxfile,&minfile,&stepfile,&totalpoint,&enlargerate);
	printf("%d %d %d %d %f",maxfile,minfile,stepfile,totalpoint,enlargerate);
	/*		double alf,xx,yy,zz,rr,pi=3.14159265;
		for(int i=0;i<11;i++){
			xx = -2+i*.4;
			rr = sqrt(2*2-fabs(xx)*fabs(xx));
			for(int j=0;j<4;j++){
			alf = (j+1)*pi/2.0;
			temp[i][j][0] = xx;
			temp[i][j][1] = rr*sin(alf);
			temp[i][j][2] = -rr*cos(alf);}
			}
*/
	
/*	char buf[SZ],filename[13];
	snoffile=1750;
	convert(snoffile,filename);
	fstream in(filename);
	p=(float *)temp;
	while(in.getline(buf, SZ)) { // Removes \n
	istrstream s(buf);
	s>>sn>>label>>x>>y>>z;
	*p=(float)sn;p++;
	*p=label;	 p++;
	*p=(x-0)*.03;p++;
	*p=(y-0)*.03;p++;
	*p=(z-0)*.03;p++;
	}
*/
  //printf(" %f %f %f\n",filename,temp[2][0],temp[2][1],temp[2][2]);
    //cout<<sn<<label<<x<<y<<z;  
	getdata(1750);
	myinit();
    auxKeyFunc(AUX_LEFT,left);
    auxKeyFunc(AUX_RIGHT,right);
    auxKeyFunc(AUX_UP,up);
    auxKeyFunc(AUX_DOWN,down);
    auxKeyFunc(AUX_RETURN,big);
    auxKeyFunc(AUX_SPACE,small1);

    
    auxReshapeFunc(reshape);
    auxMainLoop(display);
}
//end of sample
////////////////////////////////////////////////////////////////