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

源码开发语言/平台
当前位置：首页 > 源码/资料 > Windows编程 > 多媒体编程 > 查看源码
inv_mdct.cpp
资源名称：Visual C++视频音频开发实用工程案例精选.rar　[点击查看]
上传用户：tuheem
上传日期：2007-05-01
资源大小：21889k
文件大小：13k
源码类别：

多媒体编程

开发平台：

Visual C++

inv_mdct.cpp：源码内容
  1. /* inv_mdct.cpp
  2.  *
  3.  * Inverse Discrete Cosine Transform for hybrid
  4.  * synthesis in MPEG Audio Layer III. Based on the public_c source,
  5.  * but almost completely redone by Jeff Tsay and Francois-Raymond Boyer.
  6.  *
  7.  * Tomilla's 9 point IDCT algorithm has been replaced with Boyer's, which
  8.  * is faster. The final output multiplications required by Lee's FDCT
  9.  * have been folded into the multiplication by the window.
  10.  *
  11.  * Last modified : 08/02/97
  12.  */
  14. #include "all.h"
  15. #include "inv_mdct.h"
  17. /*
  18.  * This uses Byeong Gi Lee's Fast Cosine Transform algorithm to
  19.  * decompose the 36 point and 12 point IDCT's into 9 point and 3
  20.  * point IDCT's, respectively. Then the 9 point IDCT is computed
  21.  * by Francois-Raymond Boyer's algorithm. See his comments before
  22.  * the first 9 point IDCT. The 3 point IDCT is already efficient
  23.  * to implement. -- Jeff Tsay. 
  24.  */
  26. #ifdef COOLPRO
  27. #pragma optimize("",off) // If full optimizations are done,
  28. #pragma optimize("t",on) // VC 5.0 generates bad code
  29. #endif
  31. real win[4][36] =
  32. {
  33.  { -1.6141214951E-02f, -5.3603178919E-02f, -1.0070713296E-01f, -1.6280817573E-01f,
  34.    -4.9999999679E-01f, -3.8388735032E-01f, -6.2061144372E-01f, -1.1659756083E+00f,
  35.    -3.8720752656E+00f, -4.2256286556E+00f, -1.5195289984E+00f, -9.7416483388E-01f,
  36.    -7.3744074053E-01f, -1.2071067773E+00f, -5.1636156596E-01f, -4.5426052317E-01f,
  37.    -4.0715656898E-01f, -3.6969460527E-01f, -3.3876269197E-01f, -3.1242222492E-01f,
  38.    -2.8939587111E-01f, -2.6880081906E-01f, -5.0000000266E-01f, -2.3251417468E-01f,
  39.    -2.1596714708E-01f, -2.0004979098E-01f, -1.8449493497E-01f, -1.6905846094E-01f,
  40.    -1.5350360518E-01f, -1.3758624925E-01f, -1.2103922149E-01f, -2.0710679058E-01f,
  41.    -8.4752577594E-02f, -6.4157525656E-02f, -4.1131172614E-02f, -1.4790705759E-02f },
  43.  { -1.6141214951E-02f, -5.3603178919E-02f, -1.0070713296E-01f, -1.6280817573E-01f,
  44.    -4.9999999679E-01f, -3.8388735032E-01f, -6.2061144372E-01f, -1.1659756083E+00f,
  45.    -3.8720752656E+00f, -4.2256286556E+00f, -1.5195289984E+00f, -9.7416483388E-01f,
  46.    -7.3744074053E-01f, -1.2071067773E+00f, -5.1636156596E-01f, -4.5426052317E-01f,
  47.    -4.0715656898E-01f, -3.6969460527E-01f, -3.3908542600E-01f, -3.1511810350E-01f,
  48.    -2.9642226150E-01f, -2.8184548650E-01f, -5.4119610000E-01f, -2.6213228100E-01f, 
  49.    -2.5387916537E-01f, -2.3296291359E-01f, -1.9852728987E-01f, -1.5233534808E-01f, 
  50.    -9.6496400054E-02f, -3.3423828516E-02f, 0.0000000000E+00f, 0.0000000000E+00f,
  51.    0.0000000000E+00f, 0.0000000000E+00f, 0.0000000000E+00f, 0.0000000000E+00f },
  53.  { -4.8300800645E-02f, -1.5715656932E-01f, -2.8325045177E-01f, -4.2953747763E-01f, 
  54.    -1.2071067795E+00f, -8.2426483178E-01f, -1.1451749106E+00f, -1.7695290101E+00f,
  55.    -4.5470225061E+00f, -3.4890531002E+00f, -7.3296292804E-01f, -1.5076514758E-01f, 
  56.    0.0000000000E+00f, 0.0000000000E+00f, 0.0000000000E+00f, 0.0000000000E+00f,
  57.    0.0000000000E+00f, 0.0000000000E+00f, 0.0000000000E+00f, 0.0000000000E+00f,
  58.    0.0000000000E+00f, 0.0000000000E+00f, 0.0000000000E+00f, 0.0000000000E+00f,
  59.    0.0000000000E+00f, 0.0000000000E+00f, 0.0000000000E+00f, 0.0000000000E+00f, 
  60.    0.0000000000E+00f, 0.0000000000E+00f, 0.0000000000E+00f, 0.0000000000E+00f,
  61.    0.0000000000E+00f, 0.0000000000E+00f, 0.0000000000E+00f, 0.0000000000E+00f },
  63.  { 0.0000000000E+00f, 0.0000000000E+00f, 0.0000000000E+00f, 0.0000000000E+00f,
  64.    0.0000000000E+00f, 0.0000000000E+00f, -1.5076513660E-01f, -7.3296291107E-01f,
  65.    -3.4890530566E+00f, -4.5470224727E+00f, -1.7695290031E+00f, -1.1451749092E+00f,
  66.    -8.3137738100E-01f, -1.3065629650E+00f, -5.4142014250E-01f, -4.6528974900E-01f,
  67.    -4.1066990750E-01f, -3.7004680800E-01f, -3.3876269197E-01f, -3.1242222492E-01f,
  68.    -2.8939587111E-01f, -2.6880081906E-01f, -5.0000000266E-01f, -2.3251417468E-01f,
  69.    -2.1596714708E-01f, -2.0004979098E-01f, -1.8449493497E-01f, -1.6905846094E-01f,
  70.    -1.5350360518E-01f, -1.3758624925E-01f, -1.2103922149E-01f, -2.0710679058E-01f,
  71.    -8.4752577594E-02f, -6.4157525656E-02f, -4.1131172614E-02f, -1.4790705759E-02f }
  72. };
  74. void inv_mdct(real *in, real *out, int32 block_type)
  75. {
  76.  real    tmp[18];
  77.  real    *win_bt;
  78.     int32   i;
  80.  if(block_type == 2){
  82.    for(int32 p=0;p<36;p+=9) {
  83.     out[p]   = out[p+1] = out[p+2] = out[p+3] =
  84.       out[p+4] = out[p+5] = out[p+6] = out[p+7] =
  85.       out[p+8] = 0.0f;
  86.    }
  88.      int32 six_i = 0;
  90.    for(i=0;i<3;i++)
  91.     {
  93.        // 12 point IMDCT
  95.          // Begin 12 point IDCT
  97.     // Input aliasing for 12 pt IDCT
  98.     in[15+i] += in[12+i]; in[12+i] += in[9+i]; in[9+i]  +=  in[6+i];
  99.     in[6+i]  += in[3+i];  in[3+i]  += in[0+i];
  101.     // Input aliasing on odd indices (for 6 point IDCT)
  102.     in[15+i] += in[9+i];  in[9+i]  += in[3+i];
  104.     // 3 point IDCT on even indices
  106.          {
  107.       real  pp1, pp2, sum;
  108. pp2 = in[12+i] * 0.500000000f;
  109.     pp1 = in[ 6+i] * 0.866025403f;
  110.     sum = in[0+i] + pp2;
  111.     tmp[1] = in[0+i] - in[12+i];
  112.     tmp[0] = sum + pp1;
  113.     tmp[2] = sum - pp1;
  115.        // End 3 point IDCT on even indices
  117.     // 3 point IDCT on odd indices (for 6 point IDCT)
  119. pp2 = in[15+i] * 0.500000000f;
  120.     pp1 = in[ 9+i] * 0.866025403f;
  121.     sum = in[ 3+i] + pp2;
  122.     tmp[4] = in[3+i] - in[15+i];
  123.     tmp[5] = sum + pp1;
  124.     tmp[3] = sum - pp1;
  125.          }
  127.        // End 3 point IDCT on odd indices
  129.     // Twiddle factors on odd indices (for 6 point IDCT)
  131.     tmp[3] *= 1.931851653f;
  132.     tmp[4] *= 0.707106781f;
  133.     tmp[5] *= 0.517638090f;
  135.     // Output butterflies on 2 3 point IDCT's (for 6 point IDCT)
  137.          {
  138.     real save = tmp[0];
  139.     tmp[0] += tmp[5];
  140.     tmp[5] = save - tmp[5];
  141.     save = tmp[1];
  142.     tmp[1] += tmp[4];
  143.     tmp[4] = save - tmp[4];
  144.     save = tmp[2];
  145.     tmp[2] += tmp[3];
  146.     tmp[3] = save - tmp[3];
  147.          }
  149.     // End 6 point IDCT
  151.     // Twiddle factors on indices (for 12 point IDCT)
  153.     tmp[0]  *=  0.504314480f;
  154.     tmp[1]  *=  0.541196100f;
  155.     tmp[2]  *=  0.630236207f;
  156.     tmp[3]  *=  0.821339815f;
  157.     tmp[4]  *=  1.306562965f;
  158.     tmp[5]  *=  3.830648788f;
  160.        // End 12 point IDCT
  162.     // Shift to 12 point modified IDCT, multiply by window type 2
  163.     tmp[8]  = -tmp[0] * 0.793353340f;
  164.     tmp[9]  = -tmp[0] * 0.608761429f;
  165.     tmp[7]  = -tmp[1] * 0.923879532f;
  166.     tmp[10] = -tmp[1] * 0.382683432f;
  167.     tmp[6]  = -tmp[2] * 0.991444861f;
  168.     tmp[11] = -tmp[2] * 0.130526192f;
  170.     tmp[0]  =  tmp[3];
  171.     tmp[1]  =  tmp[4] * 0.382683432f;
  172.     tmp[2]  =  tmp[5] * 0.608761429f;
  174.     tmp[3]  = -tmp[5] * 0.793353340f;
  175.     tmp[4]  = -tmp[4] * 0.923879532f;
  176.     tmp[5]  = -tmp[0] * 0.991444861f;
  178.     tmp[0] *= 0.130526192f;
  180.     out[six_i + 6]  += tmp[0];
  181. out[six_i + 7]  += tmp[1];
  182.     out[six_i + 8]  += tmp[2];
  183. out[six_i + 9]  += tmp[3];
  184.     out[six_i + 10] += tmp[4];
  185. out[six_i + 11] += tmp[5];
  186.     out[six_i + 12] += tmp[6];
  187. out[six_i + 13] += tmp[7];
  188.     out[six_i + 14] += tmp[8];
  189. out[six_i + 15] += tmp[9];
  190.     out[six_i + 16] += tmp[10];
  191. out[six_i + 17] += tmp[11];
  193.     six_i += 6;
  194.     }
  196.  } else {
  199.    // 36 point IDCT
  201.    // input aliasing for 36 point IDCT
  202.    in[17]+=in[16]; in[16]+=in[15]; in[15]+=in[14]; in[14]+=in[13];
  203.    in[13]+=in[12]; in[12]+=in[11]; in[11]+=in[10]; in[10]+=in[9];
  204.    in[9] +=in[8];  in[8] +=in[7];  in[7] +=in[6];  in[6] +=in[5];
  205.    in[5] +=in[4];  in[4] +=in[3];  in[3] +=in[2];  in[2] +=in[1];
  206.    in[1] +=in[0];
  208.    // 18 point IDCT for odd indices
  210.    // input aliasing for 18 point IDCT
  211.    in[17]+=in[15]; in[15]+=in[13]; in[13]+=in[11]; in[11]+=in[9];
  212.    in[9] +=in[7];  in[7] +=in[5];  in[5] +=in[3];  in[3] +=in[1];
  214. {
  215.    real tmp0,tmp1,tmp2,tmp3,tmp4,tmp0_,tmp1_,tmp2_,tmp3_;
  216.    real tmp0o,tmp1o,tmp2o,tmp3o,tmp4o,tmp0_o,tmp1_o,tmp2_o,tmp3_o;
  218. // Fast 9 Point Inverse Discrete Cosine Transform
  219. //
  220. // By  Francois-Raymond Boyer
  221. //         mailto:boyerf@iro.umontreal.ca
  222. //         http://www.iro.umontreal.ca/~boyerf
  223. //
  224. // The code has been optimized for Intel processors
  225. //  (takes a lot of time to convert float to and from iternal FPU representation)
  226. //
  227. // It is a simple "factorization" of the IDCT matrix.
  229.    // 9 point IDCT on even indices
  230.    {
  231. // 5 points on odd indices (not realy an IDCT)
  232.    real i0 = in[0]+in[0];
  233.    real i0p12 = i0 + in[12];
  235.    tmp0 = i0p12 + in[4]*1.8793852415718f  + in[8]*1.532088886238f   + in[16]*0.34729635533386f;
  236.    tmp1 = i0    + in[4]                   - in[8] - in[12] - in[12] - in[16];
  237.    tmp2 = i0p12 - in[4]*0.34729635533386f - in[8]*1.8793852415718f  + in[16]*1.532088886238f;
  238.    tmp3 = i0p12 - in[4]*1.532088886238f   + in[8]*0.34729635533386f - in[16]*1.8793852415718f;
  239.    tmp4 = in[0] - in[4]                   + in[8] - in[12]          + in[16];
  241. // 4 points on even indices
  242.    real i6_ = in[6]*1.732050808f; // Sqrt[3]
  244.    tmp0_ = in[2]*1.9696155060244f  + i6_ + in[10]*1.2855752193731f  + in[14]*0.68404028665134f;
  245.    tmp1_ = (in[2]                        - in[10]                   - in[14])*1.732050808f;
  246.    tmp2_ = in[2]*1.2855752193731f  - i6_ - in[10]*0.68404028665134f + in[14]*1.9696155060244f;
  247.    tmp3_ = in[2]*0.68404028665134f - i6_ + in[10]*1.9696155060244f  - in[14]*1.2855752193731f;
  248.    }
  250.    // 9 point IDCT on odd indices
  251.    {
  252. // 5 points on odd indices (not realy an IDCT)
  253.    real i0 = in[0+1]+in[0+1];
  254.    real i0p12 = i0 + in[12+1];
  256.    tmp0o = i0p12   + in[4+1]*1.8793852415718f  + in[8+1]*1.532088886238f       + in[16+1]*0.34729635533386f;
  257.    tmp1o = i0      + in[4+1]                   - in[8+1] - in[12+1] - in[12+1] - in[16+1];
  258.    tmp2o = i0p12   - in[4+1]*0.34729635533386f - in[8+1]*1.8793852415718f      + in[16+1]*1.532088886238f;
  259.    tmp3o = i0p12   - in[4+1]*1.532088886238f   + in[8+1]*0.34729635533386f     - in[16+1]*1.8793852415718f;
  260.    tmp4o = (in[0+1] - in[4+1]                   + in[8+1] - in[12+1]            + in[16+1])*0.707106781f; // Twiddled
  262. // 4 points on even indices
  263.    real i6_ = in[6+1]*1.732050808f; // Sqrt[3]
  265.    tmp0_o = in[2+1]*1.9696155060244f  + i6_ + in[10+1]*1.2855752193731f  + in[14+1]*0.68404028665134f;
  266.    tmp1_o = (in[2+1]                        - in[10+1]                   - in[14+1])*1.732050808f;
  267.    tmp2_o = in[2+1]*1.2855752193731f  - i6_ - in[10+1]*0.68404028665134f + in[14+1]*1.9696155060244f;
  268.    tmp3_o = in[2+1]*0.68404028665134f - i6_ + in[10+1]*1.9696155060244f  - in[14+1]*1.2855752193731f;
  269.    }
  271.    // Twiddle factors on odd indices
  272.    // and
  273.    // Butterflies on 9 point IDCT's
  274.    // and
  275.    // twiddle factors for 36 point IDCT
  277.    real e, o;
  278.    e = tmp0 + tmp0_; o = (tmp0o + tmp0_o)*0.501909918f; tmp[0] = e + o;    tmp[17] = e - o;
  279.    e = tmp1 + tmp1_; o = (tmp1o + tmp1_o)*0.517638090f; tmp[1] = e + o;    tmp[16] = e - o;
  280.    e = tmp2 + tmp2_; o = (tmp2o + tmp2_o)*0.551688959f; tmp[2] = e + o;    tmp[15] = e - o;
  281.    e = tmp3 + tmp3_; o = (tmp3o + tmp3_o)*0.610387294f; tmp[3] = e + o;    tmp[14] = e - o;
  282.    tmp[4] = tmp4 + tmp4o; tmp[13] = tmp4 - tmp4o;
  283.    e = tmp3 - tmp3_; o = (tmp3o - tmp3_o)*0.871723397f; tmp[5] = e + o;    tmp[12] = e - o;
  284.    e = tmp2 - tmp2_; o = (tmp2o - tmp2_o)*1.183100792f; tmp[6] = e + o;    tmp[11] = e - o;
  285.    e = tmp1 - tmp1_; o = (tmp1o - tmp1_o)*1.931851653f; tmp[7] = e + o;    tmp[10] = e - o;
  286.    e = tmp0 - tmp0_; o = (tmp0o - tmp0_o)*5.736856623f; tmp[8] = e + o;    tmp[9] =  e - o;
  287.    }
  289.    // end 36 point IDCT */
  291. // shift to modified IDCT
  292.    win_bt = win[block_type];
  294. out[0] =-tmp[9]  * win_bt[0];
  295.    out[1] =-tmp[10] * win_bt[1];
  296. out[2] =-tmp[11] * win_bt[2];
  297.    out[3] =-tmp[12] * win_bt[3];
  298.    out[4] =-tmp[13] * win_bt[4];
  299. out[5] =-tmp[14] * win_bt[5];
  300. out[6] =-tmp[15] * win_bt[6];
  301. out[7] =-tmp[16] * win_bt[7];
  302. out[8] =-tmp[17] * win_bt[8];
  304.    out[9] = tmp[17] * win_bt[9];
  305.    out[10]= tmp[16] * win_bt[10];
  306. out[11]= tmp[15] * win_bt[11];
  307. out[12]= tmp[14] * win_bt[12];
  308. out[13]= tmp[13] * win_bt[13];
  309. out[14]= tmp[12] * win_bt[14];
  310.    out[15]= tmp[11] * win_bt[15];
  311. out[16]= tmp[10] * win_bt[16];
  312. out[17]= tmp[9]  * win_bt[17];
  313. out[18]= tmp[8]  * win_bt[18];
  314.    out[19]= tmp[7]  * win_bt[19];
  315. out[20]= tmp[6]  * win_bt[20];
  316.    out[21]= tmp[5]  * win_bt[21];
  317. out[22]= tmp[4]  * win_bt[22];
  318. out[23]= tmp[3]  * win_bt[23];
  319.   out[24]= tmp[2]  * win_bt[24];
  320.    out[25]= tmp[1]  * win_bt[25];
  321. out[26]= tmp[0]  * win_bt[26];
  323.    out[27]= tmp[0]  * win_bt[27];
  324. out[28]= tmp[1]  * win_bt[28];
  325. out[29]= tmp[2]  * win_bt[29];
  326. out[30]= tmp[3]  * win_bt[30];
  327. out[31]= tmp[4]  * win_bt[31];
  328.    out[32]= tmp[5]  * win_bt[32];
  329. out[33]= tmp[6]  * win_bt[33];
  330. out[34]= tmp[7]  * win_bt[34];
  331. out[35]= tmp[8]  * win_bt[35];
  332. }
  333. }
  335. #ifdef COOLPRO
  336. #pragma optimize("",on)  // Restore optimizations for VC 5.0
  337. #endif