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");

}

*/