#include <stdlib.h>
#include <stdio.h>
#include <string.h>
#include <math.h>
#include <iostream>
#include <iomanip>
#include "ModelViewMatrix.h"
using namespace std;

struct Face{
	vertex vert[3];
	double ambient[3];
	double refractive;
	double reflective;
	double diffuse[3];
	double specular[3];
	double shininess;
	double emission[3];
	Mat4 matrix;
	Mat4 inverseMatrix;
	double dielectric;
	double reflectivity;
	double translucency;

};

class Obj{
public:
	vector<vertex> v;
	vector<vertex> vt;
	vector<vertex> vn;
	vector<Face> faces;
	double center[3];
	double radius;
	void setMatrices(Mat4 m, Mat4 mi){
		for(int i = 0; i < faces.size(); i++){
			faces[i].matrix = m;
			faces[i].inverseMatrix = mi;
		}
	}
	void setDielectric(double a){
		for(int i = 0; i < faces.size(); i++){
			faces[i].dielectric = a;
		}
	}
	void setReflectivity(double a){
		for(int i = 0; i < faces.size(); i++){
			faces[i].reflectivity = a;
		}
	}
	void setTranslucency(double a){
		for(int i = 0; i < faces.size(); i++){
			faces[i].translucency = a;
		}
	}
};
double quickDist(double x, double y, double z, double a, double b, double c){
	return (x-a)*(x-a) + (y-b)*(y-b) + (z-c)*(z-c);
}

//void printObj(Obj obj){
//	for(int i = 0; i < obj.faces.size(); i++){
//		cout
//	}
//}

Obj loadObj(string filename){
	char line[1000];
	char command[1000];
	Obj obj;
	for(int i = 0; i < 3; i++)
		obj.center[i] = 0;

	obj.radius = 0;
	double numverts = 0;
	FILE *fp ; 
	fp = fopen(filename.c_str(), "rt") ;
	while (!feof(fp)) {
		cout << setiosflags(ios::fixed) << setprecision(5);
		fgets(line,sizeof(line),fp) ;
		if (feof(fp)) break ;
		if (line[0] == '#') continue ; // Comment lines
	    
		int num = sscanf(line, "%s", command) ;
		if (num != 1) continue ; // Blank line etc.

		// Now, we simply parse the file by looking at the first line for the 
		// various commands
	    
		/**************   OBJ  **********/
	  
		if (!strcmp(command, "v")) {
			vertex v;
			double x, y, z;
			int num = sscanf(line, "%s %f %f %f", command, &v.x, &v.y, &v.z) ;
			//int num = sscanf(line, "%s %lf %lf %lf", command, &x, &y, &z) ;
			if (num != 4) {
				fprintf(stderr, "v x y z") ;
				exit(1) ;
			}
			assert(!strcmp(command,"v")) ;
			//v.init();
			//v.x = x;
			//v.y = y; 
			//v.z = z;
			obj.center[0] += v.x;
			obj.center[1] += v.y;
			obj.center[2] += v.z;
			numverts++;
			obj.v.push_back(v);
			cout << "v  " << v.x << " " << v.y<< " " << v.z << endl;
		}else if(!strcmp(command,"vt")){
			vertex vt;
			int num = sscanf(line, "%s %f %f", command, &vt.u, &vt.v);
			if (num != 3) {
				fprintf(stderr, "vt u v") ;
				exit(1) ;
			}
			assert(!strcmp(command,"vt")) ;
			obj.vt.push_back(vt);
			cout << "vt  " << vt.u << " " << vt.v << endl;
			
		}else if(!strcmp(command,"vn")){
			vertex vn;
			int num = sscanf(line, "%s %f %f %f", command, &vn.x, &vn.y, &vn.z);
			if (num != 4) {
				fprintf(stderr, "vn x y z") ;
				exit(1) ;
			}
			assert(!strcmp(command,"vn")) ;
			obj.vn.push_back(vn);
			cout << "vn  " << vn.x << " " << vn.y<< " " << vn.z << endl;
		}else if(!strcmp(command,"f")){
			Face f;

			int v0,t0,n0,v1,t1,n1,v2,t2,n2;
			int num = sscanf(line, "%s %d/%d/%d %d/%d/%d %d/%d/%d", 
				command, &v0,&t0,&n0,&v1,&t1,&n1,&v2,&t2,&n2);
			cout << "f " <<v0<<"/"<<t0<<"/"<<n0<<" "<<v1<<"/"<<t1<<"/"<<n1<<" "<<v2<<"/"<<t2<<"/"<<n2<< endl;
			//obj is base 1
			v0 -= 1;
			t0 -= 1;
			n0 -= 1;
			v1 -= 1;
			t1 -= 1;
			n1 -= 1;
			v2 -= 1;
			t2 -= 1;
			n2 -= 1;

			if (num != 10) {
				fprintf(stderr, "f v1 t1 n1 v2 t2 n2 v3 t3 n3") ;
				exit(1) ;
			}
			assert(!strcmp(command,"f")) ;

			f.vert[0].x = obj.v[v0].x;
			f.vert[0].y = obj.v[v0].y;
			f.vert[0].z = obj.v[v0].z;
			f.vert[1].x = obj.v[v1].x;
			f.vert[1].y = obj.v[v1].y;
			f.vert[1].z = obj.v[v1].z;
			f.vert[2].x = obj.v[v2].x;
			f.vert[2].y = obj.v[v2].y;
			f.vert[2].z = obj.v[v2].z;

			f.vert[0].u = obj.vt[t0].u;
			f.vert[0].v = obj.vt[t0].v;
			f.vert[1].u = obj.vt[t1].u;
			f.vert[1].v = obj.vt[t1].v;
			f.vert[2].u = obj.vt[t2].u;
			f.vert[2].v = obj.vt[t2].v;

			f.vert[0].nx = obj.vn[n0].x;
			f.vert[0].ny = obj.vn[n0].y;
			f.vert[0].nz = obj.vn[n0].z;
			f.vert[1].nx = obj.vn[n1].x;
			f.vert[1].ny = obj.vn[n1].y;
			f.vert[1].nz = obj.vn[n1].z;
			f.vert[2].nx = obj.vn[n2].x;
			f.vert[2].ny = obj.vn[n2].y;
			f.vert[2].nz = obj.vn[n2].z;
			obj.faces.push_back(f);
			
		}
	}
	for(int i = 0; i < 3; i++)
		obj.center[i] /= numverts;
	for(int i = 0; i < obj.v.size(); i++){
		double newDist = quickDist(obj.v[i].x, obj.v[i].y, obj.v[i].z, obj.center[0], obj.center[1], obj.center[2]);
		obj.radius = obj.radius < newDist ? newDist : obj.radius;
	}
	obj.radius = sqrt(obj.radius);
//	cout << '\EOF';
	//printObj(obj);
	return obj;
}