#include "prj.h"

dct(short y[N][N],int Y[N*N])
{
  register short i,j,m,n;
  float t;

  for(j=0;j<N;j++)
    for(i=0;i<N;i++){
      t=0;
      for(n=0;n<N;n++)
	for(m=0;m<N;m++)
	  t+=eigen[n][j]*y[n][m]*eigen[m][i];
      Y[zz[j*N+i]]=(int) t;
    }

}

idct(short y[N][N],int Y[N*N])
{
  register short i,j,m,n;
  float t;

  for(n=0;n<N;n++)
    for(m=0;m<N;m++){
      t=0;
      for(j=0;j<N;j++)
	for(i=0;i<N;i++)
	  t+=eigen[n][j]*Y[zz[j*N+i]]*eigen[m][i];
      y[n][m]=(short) t;
    }

}

/*

main()
{
  int i,j,dc;
  short table[N][N];
  int Y[N*N];

  set_eigen();
  set_zigzag();

  table[0][0]=-1; table[0][1]=4; table[0][2]=6; table[0][3]=8;
  table[0][4]=4; table[0][5]=-5; table[0][6]=-4; table[0][7]=2;
  table[1][0]=-2; table[1][1]=-2; table[1][2]=1; table[1][3]=7;
  table[1][4]=9; table[1][5]=6; table[1][6]=3; table[1][7]=2;

  table[2][0]=-2; table[2][1]=-9;  table[2][2]=-6;  table[2][3]=-3;
  table[2][4]=4;  table[2][5]=20; table[2][6]=29; table[2][7]=22;
  table[3][0]=-1; table[3][1]=6;  table[3][2]=5;  table[3][3]=-7;
  table[3][4]=-13;  table[3][5]=12; table[3][6]=34; table[3][7]=26;

  table[4][0]=2; table[4][1]=-8;  table[4][2]=-3;  table[4][3]=2;
  table[4][4]=11;  table[4][5]=17; table[4][6]=-11; table[4][7]=-33;
  table[5][0]=4; table[5][1]=-8;  table[5][2]=-5;  table[5][3]=6;
  table[5][4]=13;  table[5][5]=17; table[5][6]=-14; table[5][7]=-36;

  table[6][0]=2; table[6][1]=-5;  table[6][2]=-7;  table[6][3]=6;
  table[6][4]=10;  table[6][5]=3; table[6][6]=-12; table[6][7]=11;
  table[7][0]=-3; table[7][1]=-4;  table[7][2]=-3;  table[7][3]=11;
  table[7][4]=12;  table[7][5]=-9; table[7][6]=-7; table[7][7]=29;


  for(j=0;j<N;j++){
    for(i=0;i<N;i++)
      printf(" %3d ",table[j][i]);
    printf("\n");
  }
  printf("\n");

  dct(table,Y);

  for(j=0;j<N;j++){
    for(i=0;i<N;i++)
      printf(" %4d ",Y[j*N+i]);
    printf("\n");
  }
  printf("\n");

  quant(1,Y);
  iquant(1,Y);

  for(j=0;j<N;j++){
    for(i=0;i<N;i++)
      printf(" %4d ",Y[j*N+i]);
    printf("\n");
  }
  printf("\n");

  for(j=0;j<N;j++)
    for(i=0;i<N;i++)
      table[i][j]=0;

  idct(table,Y);

  for(j=0;j<N;j++){
    for(i=0;i<N;i++)
      printf(" %3d ",table[j][i]);
    printf("\n");
  }
  printf("\n");
}


*/