Subversion Repositories mb-jpeg

 
#include <stdio.h>
#include "dct.h"
#include "weights.h"
 
 
/*
	Function Name: dct
 
	Operation: Find the 8x8 DCT of an array using separable DCT
	First, finds 1-d DCT along rows, storing the result in inter[][]
	Then, 1-d DCT along columns of inter[][] is found
 
	Input: pixels is the 8x8 input array
 
	Output: dct is the 8x8 output array
*/
signed int max_val = 0;
static int __count=0;
 
void dct(signed char pixels[8][8], signed short dct[8][8])
{
        FILE * file;
	int inr, inc; 		/* rows and columns of input image */
	int intr, intc;		/* rows and columns of intermediate image */
	int outr, outc;		/* rows and columns of dct */
	int f_val;		/* cumulative sum */
	int inter[8][8];	/* stores intermediate result */
	int i,j,k;
        k=0;
    //    file = fopen("weights.h","w+");
      //  fprintf(file,"double weights1[512] = {");
	/* find 1-d dct along rows */
#if 0	
 
#ifdef __MICROBLAZE
       if (__count<4) {
	   	for (i=0; i<8; i++)
			for (j=0; j<8; j++)
	xil_printf("%x ", pixels[i][j]);
       xil_printf("\r\n");
       	}
#else
       if (__count<4) {
	   	for (i=0; i<8; i++)
			for (j=0; j<8; j++)
	printf(" %x ", pixels[i][j]);
       printf("\r\n");
       	}
#endif
#endif
 
 	for (intr=0; intr<8; intr++)
		for (intc=0; intc<8; intc++) {
			for (i=0,f_val=0; i<8; i++) {
 
			      	f_val += (pixels[intr][i]* weights[k]);//cos((double)(2*i+1)*(double)intc*PI/16);
                                k++;
                          //     fprintf(file, "\n%.0f,",cos((double)(2*i+1)*(double)intc*PI/16)*16384);
			}
                        if (intc!=0)
                                inter[intr][intc] =  f_val>>15;
                        else
                                inter[intr][intc] =  (11585*(f_val>>14))>>15;
 
                }
   //     fprintf(file,"\n};");
   //     fclose(file);
         k=0;
	/* find 1-d dct along columns */
 	for (outc=0; outc<8; outc++)
		for (outr=0; outr<8; outr++) {
			for (i=0,f_val=0; i<8; i++) {
				f_val += (inter[i][outc] *weights[k]);
                                k++;
			}
                        if (outr!=0)
			        dct[outr][outc] = f_val>>15;
                        else
                                dct[outr][outc] = (11585*(f_val>>14)>>15);
		}
#if 0
#ifdef __MICROBLAZE
       if (__count<4) {
	   	for (i=0; i<8; i++)
			for (j=0; j<8; j++)
	xil_printf("%x ", dct[i][j]);
       xil_printf("\r\n");
	__count++;
       	}
#else
       if (__count<4) {
	   	for (i=0; i<8; i++)
			for (j=0; j<8; j++)
	printf(" %x ", dct[i][j]);
       printf("\r\n");
	__count++;
       	}
#endif
#endif
 
}
 
/*****************************************************************
    UNCOMMENT THIS SECTION TO TEST 2D DCT 
*****************************************************************/
 
/*
main()
{
 
  unsigned char inputmatrix[8][8];
  unsigned char outputmatrix[8][8];
  unsigned int i,j;
 
 
  printf("Input Matrix (8*8) :-\n");
  for (i=0; i<8; i++){
	  printf("\n");
	  for (j=0;j<8;j++){
		   inputmatrix[i][j] = i*8+j;
		   printf("%4d",inputmatrix[i][j]);
	  }
 
  }
 
 
 
  dct(inputmatrix,outputmatrix);
 
  printf("\n\nOutput Matrix (8*8) :-\n");
 
   for (i=0; i<8; i++){
	   printf("\n");
	  for (j=0;j<8;j++){
	   printf("%4d",outputmatrix[i][j]);
 
	  }
 
  } 
printf("\n");
 
}
*/