开通博客赚积分
发布资源赚积分
分类

源码开发语言/平台
当前位置：首页 > 源码/资料 > 嵌入式/单片机编程 > 其他嵌入式/单片机内容 > 查看源码
idctref.c
资源名称：mpeg2.rar　[点击查看]
上传用户：ma_junhua
上传日期：2008-04-11
资源大小：2752k
文件大小：3k
源码类别：

其他嵌入式/单片机内容

开发平台：

C/C++

idctref.c：源码内容
  1. /* Reference_IDCT.c, Inverse Discrete Fourier Transform, double precision          */
  3. /* Copyright (C) 1996, MPEG Software Simulation Group. All Rights Reserved. */
  5. /*
  6.  * Disclaimer of Warranty
  7.  *
  8.  * These software programs are available to the user without any license fee or
  9.  * royalty on an "as is" basis.  The MPEG Software Simulation Group disclaims
  10.  * any and all warranties, whether express, implied, or statuary, including any
  11.  * implied warranties or merchantability or of fitness for a particular
  12.  * purpose.  In no event shall the copyright-holder be liable for any
  13.  * incidental, punitive, or consequential damages of any kind whatsoever
  14.  * arising from the use of these programs.
  15.  *
  16.  * This disclaimer of warranty extends to the user of these programs and user's
  17.  * customers, employees, agents, transferees, successors, and assigns.
  18.  *
  19.  * The MPEG Software Simulation Group does not represent or warrant that the
  20.  * programs furnished hereunder are free of infringement of any third-party
  21.  * patents.
  22.  *
  23.  * Commercial implementations of MPEG-1 and MPEG-2 video, including shareware,
  24.  * are subject to royalty fees to patent holders.  Many of these patents are
  25.  * general enough such that they are unavoidable regardless of implementation
  26.  * design.
  27.  *
  28.  */
  30. /*  Perform IEEE 1180 reference (64-bit floating point, separable 8x1
  31.  *  direct matrix multiply) Inverse Discrete Cosine Transform
  32. */
  35. /* Here we use math.h to generate constants.  Compiler results may
  36.    vary a little */
  38. #include <math.h>
  40. #include "config.h"
  42. #ifndef PI
  43. # ifdef M_PI
  44. #  define PI M_PI
  45. # else
  46. #  define PI 3.14159265358979323846
  47. # endif
  48. #endif
  50. /* global declarations */
  51. void Initialize_Fast_IDCTref _ANSI_ARGS_((void));
  52. void Reference_IDCT _ANSI_ARGS_((short *block));
  54. /* private data */
  56. /* cosine transform matrix for 8x1 IDCT */
  57. static double c[8][8];
  59. /* initialize DCT coefficient matrix */
  61. void Initialize_Reference_IDCT()
  62. {
  63.   int freq, time;
  64.   double scale;
  66.   for (freq=0; freq < 8; freq++)
  67.   {
  68.     scale = (freq == 0) ? sqrt(0.125) : 0.5;
  69.     for (time=0; time<8; time++)
  70.       c[freq][time] = scale*cos((PI/8.0)*freq*(time + 0.5));
  71.   }
  72. }
  74. /* perform IDCT matrix multiply for 8x8 coefficient block */
  76. void Reference_IDCT(block)
  77. short *block;
  78. {
  79.   int i, j, k, v;
  80.   double partial_product;
  81.   double tmp[64];
  83.   for (i=0; i<8; i++)
  84.     for (j=0; j<8; j++)
  85.     {
  86.       partial_product = 0.0;
  88.       for (k=0; k<8; k++)
  89.         partial_product+= c[k][j]*block[8*i+k];
  91.       tmp[8*i+j] = partial_product;
  92.     }
  94.   /* Transpose operation is integrated into address mapping by switching 
  95.      loop order of i and j */
  97.   for (j=0; j<8; j++)
  98.     for (i=0; i<8; i++)
  99.     {
  100.       partial_product = 0.0;
  102.       for (k=0; k<8; k++)
  103.         partial_product+= c[k][i]*tmp[8*k+j];
  105.       v = (int) floor(partial_product+0.5);
  106.       block[8*i+j] = (v<-256) ? -256 : ((v>255) ? 255 : v);
  107.     }
  108. }