URL
https://opencores.org/ocsvn/mb-jpeg/mb-jpeg/trunk
Subversion Repositories mb-jpeg
[/] [mb-jpeg/] [tags/] [STEP2_1b/] [encoder/] [dct.c] - Rev 66
Compare with Previous | Blame | View Log
#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"); } */