#include "lib.h"

void erromedio(IMAGE img1, IMAGE img, int N)
{
  int i,j;
  float S=0.0;

  /* Convolui imagem de entrada */

  for(i=N/2; i < img1->info->nr-N/2; i++)
  {
     for(j=N/2; j < img1->info->nc-N/2; j++)
     {
	S = S + pow((img1->data[i][j] - img->data[i][j]),2);
     }
  }
  S = S/((img1->info->nr-N+1)*(img1->info->nc-N+1));
  S = sqrt(S);
  printf("Erro medio quadratico: %f\n",S);
}


/*********************** main ********************/
int main( int argc, char** argv ){
  IMAGE img1=0,img2=0;

  double tempo; 
  
  /* Ler imagem de entrada */

  //ler_imagem_P2(argv[1],argv[2]);

  /* aplica erro medio */

  img1=Input_PBM(argv[1]);
  img2=Input_PBM(argv[2]);

  MarkTime(NULL);
   
  erromedio(img1,img2,atoi(argv[3]));

  MarkTime(&tempo);

  ShowTime(tempo);
  printf("\n");

}