#include "prj.h"
#include <time.h>
#include <stdlib.h>

main()
{
  short frame,q,p;
  short j,i;
  short table[N][N];
  int dc;
  short qn;
  clock_t start,end;
  
  qn=1;
  set_eigen();
  set_zigzag();

  for(j=0;j<N;j++)
    for(i=0;i<N;i++){
      picel[0][j][i]=40+rand()%120;
      picel[1][j][i]=picel[0][j][i]+rand()%10-5;
    }

  printf("\npicel:\n");
  for(j=0;j<N;j++){
    for(i=0;i<N;i++)
      printf("%4d ",picel[0][j][i]);
    printf("\n");
  }
    
  printf("\npicel:\n");
  for(j=0;j<N;j++){
    for(i=0;i<N;i++)
      printf("%4d ",picel[1][j][i]);
    printf("\n");
  }
    
  p=0; q=0; frame=1; mv_x[1][0][0]=0; mv_y[1][0][0]=0;

	residue(1,0,0,0,0,table);

        dct(table,dct_value[frame][q*LY+p]);
        quant(qn,dct_value[frame][q*LY+p]);
  printf("\nresidue:\n");
  for(j=0;j<N;j++){
    for(i=0;i<N;i++)
      printf("%4d ",dct_value[frame][q*LY+p][j*8+i]);
    printf("\n");
  }


        bitlen[frame][q*LY+p]=tcoeff(dct_value[frame][q*LY+p],N*N);

        iquant(qn,dct_value[frame][q*LY+p]);
        idct(table,dct_value[frame][q*LY+p]);

	iresidue(frame,q*LY,p*LX,mv_y[frame][q/2][p/2],mv_x[frame][q/2][p/2],table);

  printf("\npicel:\n");
  for(j=0;j<N;j++){
    for(i=0;i<N;i++)
      printf("%4d ",picel[1][j][i]);
    printf("\n");
  }
    
  printf("\n%d\n",bitlen[1][q*LY+p]);
}