#include "vtkCara.h"
#include "vtkObjectFactory.h"
#include "math.h"

vtkCxxRevisionMacro(vtkCara, "$Revision: 0.1 $");
vtkStandardNewMacro(vtkCara);
vtkCara::vtkCara()
//Constructor por defecto
{
	vertices = NULL;
	numvertices = 0;
};

vtkCara::~vtkCara(){
	delete[] vertices;
	//this->UnRegister(this);
}

void vtkCara::Add_vertice(int p)
//Esta función añade un nuevo vértice a la cara tratada
{
	int i;
	int *vloc;

	numvertices++;
	if (numvertices>1)
	{
	vloc = (int *)new int[numvertices];
	for (i=0;i<(numvertices-1);i++)
		vloc[i]=vertices[i];
	vloc[numvertices-1]=p;
	delete[] vertices;
	vertices = (int *)new int[numvertices];
	for (i=0;i<numvertices;i++)
		vertices[i]=vloc[i];
	delete[] vloc;
	}
	else
	{
		vertices = (int *)new int[numvertices];
		vertices[numvertices-1]=p;
	}
};

int vtkCara::operator [](int i) 
//Se obtiene el vertice número i de la cara
{
	if (i<numvertices)
		return vertices[i];
	else
		return -1;

};

Punto vtkCara::normal(Punto *puntost2)
//Devuelve la normal de la cara
{
	Punto seg1,seg2;
	if (numvertices > 2)
	{
		seg1 = puntost2[vertices[1]] - puntost2[vertices[0]];
		seg2 = puntost2[vertices[2]] - puntost2[vertices[1]];
		return (seg1^seg2).unitario();
	}
	else
		return Punto();
};

vtkCara vtkCara::operator =(const vtkCara &c)
//Asigna caras
{
	int i;
	if (numvertices>0)
	{
		delete vertices;
		numvertices = 0;
	}
	numvertices = c.numvertices;
	vertices = (int *)new int[numvertices];
	for (i=0;i<numvertices;i++)
		vertices[i]=c.vertices[i];

	return *this;
};


BOOL vtkCara::contiene(int p)
//Comprueba si la cara contiene el vértice de índice p
{
	int i;
	BOOL encontrado = FALSE;

	i=0;
	while ((!(encontrado))&&(i<numvertices))
	{
		if (vertices[i]==p) 
			encontrado = TRUE;
		i++;
	}
	return encontrado;
}


BOOL vtkCara::dentroproy(Punto p, Punto *puntost2)
//Comprueba si el punto p pertenece a la proyección ortográfica de la cara
{
	int i;
	Punto *vp;		//Vectores del punto a los vértices de la proyección
	float ang=0.0f;	//Angulo acumulado para comprobar interioridad
	
	vp = (Punto *)new Punto[numvertices];
	for (i=0; i<numvertices;i++)
	{
		vp[i]=Punto(puntost2[vertices[i]].x-p.x,puntost2[vertices[i]].y-p.y,0.0f);
		vp[i]=vp[i].unitario();
	}
	for (i=0; i<numvertices-1;i++)
		ang += (float)acos(vp[i]*vp[i+1]);
	ang = (float)grad(ang);
	delete[] vp;
	if ((ang >= 359.9f)&&(ang <= 360.1f)) 
		return TRUE;
	else 
		return FALSE;
}

int vtkCara::indice(Punto p, Punto *puntost2)
{
	int i;

	for (i=0; i<numvertices; i++)
		if (puntost2[vertices[i]]==p)
			return i;

	return -1;
}

Punto vtkCara::centro(Punto * puntost2)
{
	Punto acum;
	int i;

	acum = puntost2[vertices[0]];
	for (i=1; i<numvertices;i++)
		acum = acum + puntost2[vertices[i]];
	acum = acum*(1.0f/numvertices);
	return acum;
}