class anillo
{
	int num;
	pto3D centroide;
	int vert[];
	  anillo ()
	{
		this.num = 0;
		centroide = null;
		vert = new int[15];
	}

	/*Resumen de operaciones
	   añadir elemento
	   giro de numeros
	   ordenar anticlockwise
	   cambio origen


	 */


	void addVert (int pton)
	{
		if (this.num < 14) {
			this.num++;
			this.vert[num] = pton;
		}
	}



	void setCentro (pto3D cide)
	{
		this.centroide = cide;
	}

	String poncentroide (MoleculaT mol)
	{
		String st = "";
		pto3D ncide = new pto3D ();
		for (int i = 1; i <= num; i++) {
			Atomo at = (Atomo) mol.susatomos.get (this.vert[i]);
			ncide = ncide.mas (at.vert);
			st = st + "\n" + ncide.aTexto ();
		}
		centroide = ncide.escala (1 / (double) num);
		return st;
	}


	void centracentroide (MoleculaT mol)
	{
		pto3D ncide = new pto3D ();
		for (int i = 1; i <= num; i++) {
			Atomo at = (Atomo) mol.susatomos.get (this.vert[i]);
			ncide = ncide.mas (at.vert);
		}
		ncide.escala (1 / num);
		if (ncide.dista (this.centroide) < 0.7)
			centroide = ncide;
	}


	void ordena (pto3D vecref, MoleculaT mol)	//vector de referencia y molecula a la que pertenecen (con toda  informacion 3D
	{
		//creamos un vector que nos indique donde esta el exterior del tubo, y desde donde miramos
		pto3D ciderel = centroide.menos (mol.TEST1);
		pto3D vtop = ciderel.proyeccplano (mol.TESTC);
		//Buscamos atomo mas cercano por la derecha
		int dcha = 0;
		double dis = 10;
		for (int i = 1; i <= num; i++) {
			Atomo at = (Atomo) mol.susatomos.get (vert[i]);
			pto3D p = at.vert;
			pto3D a = p.menos (centroide).menos (vecref);
			double di = vecref.dista (p.menos (centroide));
			double dh = a.dihedrog (vecref, vtop);
			if (di < dis && dh > 30 && dh < 150) {
				dis = di;
				dcha = i;
			}
		}

		this.ordenaccw (vert[dcha], mol);

	}

	void ordenaccw (int ini, MoleculaT mol)	//vector de referencia y molecula a la que pertenecen (con toda  informacion 3D
	{
		int newvert[] = new int[num + 1];
		//creamos un vector que nos indique donde esta el exterior del tubo, y desde donde miramos
		pto3D ciderel = centroide.menos (mol.TEST1);
		pto3D vtop = ciderel.proyeccplano (mol.TESTC);
		//Comenzamos por el atomo de inicio
		newvert[1] = ini;

		for (int i = 2; i <= num; i++) {
			pto3D pobj = ((Atomo) mol.susatomos.get (newvert[i - 1])).vert;
			int nizq = 0;
			double dihedrotemp = 0;
			double distan = 10;
			for (int j = 1; j <= num; j++) {
				Atomo at = (Atomo) mol.susatomos.get (vert[j]);
				pto3D c = at.vert;
				pto3D a = c.menos (pobj);
				pto3D b = pobj.menos (centroide);
				double dist = a.modulo ();
				double dhdr = a.dihedrog (b, vtop);
				if (dist < distan && dist > 0.01 && dhdr > -150 && dhdr < -30) {
					distan = dist;
					nizq = j;
					dihedrotemp = dhdr;
				}
			}
			newvert[i] = vert[nizq];
		}
		for (int i = 1; i <= num; i++)
			vert[i] = newvert[i];
	}
	void ordenacw (int ini, MoleculaT mol)	//vector de referencia y molecula a la que pertenecen (con toda  informacion 3D
	{
		int newvert[] = new int[num + 1];
		//creamos un vector que nos indique donde esta el exterior del tubo, y desde donde miramos
		pto3D ciderel = centroide.menos (mol.TEST1);
		pto3D vtop = ciderel.proyeccplano (mol.TESTC);
		//Comenzamos por el atomo de inicio
		newvert[1] = ini;

		for (int i = 2; i <= num; i++) {
			pto3D pobj = ((Atomo) mol.susatomos.get (newvert[i - 1])).vert;
			int nizq = 0;
			double dihedrotemp = 0;
			double distan = 10;
			for (int j = 1; j <= num; j++) {
				pto3D c = ((Atomo) mol.susatomos.get (vert[j])).vert;
				pto3D a = c.menos (pobj);
				pto3D b = pobj.menos (centroide);
				double dist = a.modulo ();
				double dhdr = a.dihedrog (b, vtop);
				if (dist < distan && dist > 0.01 && dhdr < 150 && dhdr > 30) {
					distan = dist;
					nizq = j;
					dihedrotemp = dhdr;
				}
			}
			newvert[i] = vert[nizq];
		}
		for (int i = 1; i <= num; i++)
			vert[i] = newvert[i];
	}




	void rota (int giro)	//pasa los indices hacia abajo(si giro>0), con lo que se cambia la orientación ccw)
	{
		int newvert[] = new int[num + 1];
		for (int i = 1; i <= num; i++) {
			int nuevapos = i - giro;
			if (nuevapos > num && nuevapos < 2 * num)
				nuevapos = nuevapos - num;
			else if (nuevapos <= 0 && nuevapos > -num)
				nuevapos = nuevapos + num;
			newvert[i] = vert[nuevapos];
		}
		for (int i = 1; i <= num; i++)
			vert[i] = newvert[i];

	}



	String aCadena ()
	{
		String cad = "Miembros: ";
		for (int i = 1; i <= num; i++)
			cad = cad + " " + vert[i];
		if (centroide != null)
			cad = cad + " centroide en " + centroide.aTexto ();

		return cad;
	}



}