#include <stdlib.h>
#include <string.h>
#include <stdio.h>
#include "lib.h"

/*********************** main ********************/
int main( int argc, char** argv ){
  float result;
  int val,i;
  float *lista_modelo, *lista_elem1, *lista_elem2;
  char dado[30],dado1[30];
  FILE *fp,*fp1;
  double S1=0.0, S2=0.0, S11=0.0, S22=0.0, S12=0.0, COOR=0.0;
  char t[30];
  int aux;
  double parte1,parte2,parte21,parte22,parte3,parte4;
    
   //--------------------------

   fp1 = fopen(argv[1],"r");
   if (fp1==NULL) {
	 printf("Nao consigo abrir o arquivo %s",argv[1]);
	 exit(-1);
   }

   fscanf(fp1,"%s\n",dado1);

   lista_elem1 = (float *) calloc(atoi(dado1),sizeof(float));

   for(i=0; i < atoi(dado1); i++)
   {
	 fscanf(fp1,"%s ",t);
	 lista_elem1[i] =  atof(t);
	 //printf("t[%d]=%s le=%f ",i,t,lista_elem1[i]);
   }

   //---------------------------

   S1=0.0; S11=0.0;
   for(i=0;i<atoi(dado1);i++)
   {
      S1 = S1 + lista_elem1[i];
	S11 = S11 + lista_elem1[i]*lista_elem1[i];
        
   }

   //--------------------------

   fp = fopen(argv[2],"r");
   if (fp==NULL) {
	 printf("Nao consigo abrir o arquivo %s",argv[2]);
	 exit(-1);
   }

   aux=0;
   while(!feof(fp))
   {
      fscanf(fp,"%s\n",dado);
      if(atoi(dado)!=0)
	{
      lista_elem2 = (float *) calloc(atoi(dado),sizeof(float));
      for(i=0; i < atoi(dado); i++)
      {
	 fscanf(fp,"%s ",t);
	 lista_elem2[i] =  atof(t);
      }

      //--------------------------

      S2=0.0;S22=0.0;
	for(i=0;i<atoi(dado);i++)
   	{
        	S2 = S2 + lista_elem2[i];
		S22 = S22 + lista_elem2[i]*lista_elem2[i];
   	}	

	double a=0.0,b=0.0;
      S12=0.0;
 	for(i=0;i<atoi(dado);i++)
   	{
        	S12 = S12 + lista_elem1[i]*lista_elem2[i];
   	}
      parte1 = (atoi(dado)*S12) - (S1*S2);
      parte2 = atoi(dado)*S11 - S1*S1;
      parte21 = atoi(dado)*S22;
      parte22 = S2*S2;
      parte3 = parte21 - parte22;
      parte4 = sqrt(parte2*parte3);
   	COOR = parte1 / parte4;
      aux++;
      printf("coor=%.5f ",COOR);
	
      //--------------------------
	}
      aux++;
   }

   printf("\n");
}