/* Programa calcula as arestas e angulos dados 7 vertices visiveis
    na ordem correta de um paralelepipedo.
   O arquivo de entrada dos vertices deve estar no formato
   x y z    para cada vertice.  */

#include <stdio.h>
#include <math.h>

/* calcula vetores dados 2 ptos . retorna o modulo */
float calcula_arestas(v1, v2, aresta)
float *v1,     /* vertice 1  */
      *v2,     /* vertice 2  */
      *aresta; /* aresta v1v2 = v2 - v1*/
{
 int i;
 float dist = 0.0;
 
 aresta[0] = v2[0] - v1[0];
 aresta[1] = v2[1] - v1[1];
 aresta[2] = v2[2] - v1[2];
 
 for(i=0; i<3; i++)
   dist+= aresta[i]*aresta[i];
   
 dist = sqrt(dist);
 return(dist);
}

/* para visualizar uma matriz- forma diferente de alocacao*/
void ver1(m, nlinhas, ncol)
int nlinhas, ncol;
float **m;
{
  int i, j;
  printf(stderr,"\n");
  for(i = 0; i<nlinhas; i++)
    {
      printf(" (%d)",i);
      for(j=0; j<ncol; j++)
	{
	  printf("   %.2f",m[i][j]);
        }

      printf("\n\n");
    }
}

 
void main()
{
  float modulo[9];
  float **mat_aresta;
  float **mat_vert;
  FILE *entrada;
  char file_in[50];
  int nvertices;  /* num. de vertices obtidos dos dados */
  int i;
  char s[255];
  
  fprintf(stderr,"\n\n   Calculo de arestas dados os vertices em certa ordem.\n");
  fprintf(stderr,"\n\n       Entrada de dados.\n");
  fprintf(stderr,"\n Arquivo de vertices de entrada .vert = ");
  scanf("%s",file_in);
  
  printf("\n Arquivo de vertices de entrada: %s\n",file_in); 
  
  entrada = fopen(file_in, "r");
  nvertices = 0;

  while((fgets(s, 254, entrada)) != NULL)
    nvertices++;

  fclose(entrada);
  fprintf(stderr,"\nnvertices = %d\n\n",nvertices);
  
  if(nvertices != 7)
    fprintf(stderr," O PROGRAMA DEVE SER INTERROMPIDO. NVERTICES =! 7 \n");
  
 
  mat_vert = (float **)malloc(nvertices * sizeof(float *));
  for(i=0; i<nvertices; i++)
    mat_vert[i] = (float *)malloc(3 * sizeof(float));
    
  mat_aresta = (float **)malloc(9 * sizeof(float *));
  for(i=0; i<9; i++)
    mat_aresta[i] = (float *)malloc(3 * sizeof(float));    
  
  entrada = fopen(file_in, "r");
  for(i=0; i<nvertices; i++)
    fscanf(entrada,"%f %f %f\n",&mat_vert[i][0],&mat_vert[i][1],&mat_vert[i][2]);
  

  modulo[0] = calcula_arestas(mat_vert[1], mat_vert[0], mat_aresta[0]);
  modulo[1] = calcula_arestas(mat_vert[3], mat_vert[5], mat_aresta[1]);
  modulo[2] = calcula_arestas(mat_vert[4], mat_vert[6], mat_aresta[2]);
  modulo[3] = calcula_arestas(mat_vert[2], mat_vert[1], mat_aresta[3]);
  modulo[4] = calcula_arestas(mat_vert[4], mat_vert[3], mat_aresta[4]);
  modulo[5] = calcula_arestas(mat_vert[6], mat_vert[5], mat_aresta[5]);
  modulo[6] = calcula_arestas(mat_vert[4], mat_vert[2], mat_aresta[6]);
  modulo[7] = calcula_arestas(mat_vert[3], mat_vert[1], mat_aresta[7]);
  modulo[8] = calcula_arestas(mat_vert[5], mat_vert[0], mat_aresta[8]);
  
  
  printf("  Vertices calculados\n");
  ver1(mat_vert,nvertices,3);


  printf("aresta 1-0 : [%.2f  %.2f  %.2f] mod = %.2f\n",
          mat_aresta[0][0],mat_aresta[0][1],mat_aresta[0][2],modulo[0]);

  printf("\n");
	  
  printf("aresta 3-5 : [%.2f  %.2f  %.2f] mod = %.2f\n",
          mat_aresta[1][0],mat_aresta[1][1],mat_aresta[1][2],modulo[1]);

  printf("\n");
  
  printf("aresta 4-6 : [%.2f  %.2f  %.2f] mod = %.2f\n",
          mat_aresta[2][0],mat_aresta[2][1],mat_aresta[2][2],modulo[2]);

  printf("\n");
  
  printf("aresta 2-1 : [%.2f  %.2f  %.2f] mod = %.2f\n",
          mat_aresta[3][0],mat_aresta[3][1],mat_aresta[3][2],modulo[3]);

  printf("\n");
  
  printf("aresta 4-3 : [%.2f  %.2f  %.2f] mod = %.2f\n",
          mat_aresta[4][0],mat_aresta[4][1],mat_aresta[4][2],modulo[4]);

  printf("\n");
  
  printf("aresta 6-5 : [%.2f  %.2f  %.2f] mod = %.2f\n",
          mat_aresta[5][0],mat_aresta[5][1],mat_aresta[5][2],modulo[5]);
	
  printf("\n");
  
  printf("aresta 4-2 : [%.2f  %.2f  %.2f] mod = %.2f\n",
          mat_aresta[6][0],mat_aresta[6][1],mat_aresta[6][2],modulo[6]);

  printf("\n");
  
  printf("aresta 3-1 : [%.2f  %.2f  %.2f] mod = %.2f\n",
          mat_aresta[7][0],mat_aresta[7][1],mat_aresta[7][2],modulo[7]);

  printf("\n");
  
  printf("aresta 5-0 : [%.2f  %.2f  %.2f] mod = %.2f\n",
          mat_aresta[8][0],mat_aresta[8][1],mat_aresta[8][2],modulo[8]);
	  
}