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

源码开发语言/平台
当前位置：首页 > 源码/资料 > 多媒体 > mpeg/mp3 > 查看源码
encode.c
资源名称：bladeenc-0.90.0-src.zip　[点击查看]
上传用户：bjsgzm
上传日期：2007-01-08
资源大小：256k
文件大小：9k
源码类别：

mpeg/mp3

开发平台：

Visual C++

encode.c：源码内容
  1. /*
  2. (c) Copyright 1998, 1999 - Tord Jansson
  3. =======================================
  5. This file is part of the BladeEnc MP3 Encoder, based on
  6. ISO's reference code for MPEG Layer 3 compression, and might
  7. contain smaller or larger sections that are directly taken
  8. from ISO's reference code.
  10. All changes to the ISO reference code herein are either
  11. copyrighted by Tord Jansson (tord.jansson@swipnet.se)
  12. or sublicensed to Tord Jansson by a third party.
  14. BladeEnc is free software; you can redistribute this file
  15. and/or modify it under the terms of the GNU Lesser General Public
  16. License as published by the Free Software Foundation; either
  17. version 2.1 of the License, or (at your option) any later version.
  19. */
  21. #include "common.h"
  22. #include "encoder.h"
  25. typedef double MM[SBLIMIT][64];
  26. typedef double XX[2][512];
  28. static int  off[2];
  29. static int  half[2];
  31. static MM  m;
  32. static XX  x;
  35. /*____ rebuffer_audio() __________________________________________________*/
  37. void rebuffer_audio( short buffer[2][1152], short * insamp, unsigned int samples_read, int stereo )
  38. {
  39.   unsigned int j;
  41.   if(stereo == 2)
  42. { 
  44. for(j=0;j<samples_read/2;j++)
  45. {
  46.     buffer[0][j] = insamp[2*j];
  47.     buffer[1][j] = insamp[2*j+1];
  48. }
  49.   }
  50.   else
  51. {
  52. for(j=0;j<samples_read/2;j++)
  53. {
  54. buffer[0][j] = insamp[j];
  55. buffer[1][j] = 0;
  56. }
  57.   }
  59. for( ; j < 1152 ; j++ )
  60. {
  61. buffer[0][j] = 0;
  62. buffer[1][j] = 0;
  63. }
  65.   return;
  66. }
  69. /************************************************************************
  70. *
  71. * create_ana_filter()
  72. *
  73. * PURPOSE:  Calculates the analysis filter bank coefficients
  74. *
  75. * SEMANTICS:
  76. * Calculates the analysis filterbank coefficients and rounds to the
  77. * 9th decimal place accuracy of the filterbank tables in the ISO
  78. * document.  The coefficients are stored in #filter#
  80. ************************************************************************/
  82. void create_ana_filter(double filter[SBLIMIT][64])
  83. {
  84.    register int i,k;
  86.    for (i=0; i<32; i++)
  87.       for (k=0; k<64; k++) {
  88.           if ((filter[i][k] = 1e9*cos((double)((2*i+1)*(16-k)*PI64))) >= 0)
  89.              modf(filter[i][k]+0.5, &filter[i][k]);
  90.           else
  91.              modf(filter[i][k]-0.5, &filter[i][k]);
  92.           filter[i][k] *= 1e-9;
  93.    }
  94. }
  97. /*____ initWindowFilterSubband() ____________________________________________*/
  99. void initWindowFilterSubband( void )
  100. {
  101. int i, j;
  103. off[0] = 0;
  104. off[1] = 0;
  105. half[0] = 0;
  106. half[1] = 0;
  107.     for (i=0;i<2;i++)
  108. for (j=0;j<512;j++)
  109.         x[i][j] = 0;
  111.     create_ana_filter(m);
  112. }
  115. /************************************************************************
  116. *
  117. * filter_subband()
  118. *
  119. * PURPOSE:  Calculates the analysis filter bank coefficients
  120. *
  121. * SEMANTICS:
  122. *      The windowed samples #z# is filtered by the digital filter matrix #m#
  123. * to produce the subband samples #s#. This done by first selectively
  124. * picking out values from the windowed samples, and then multiplying
  125. * them by the filter matrix, producing 32 subband samples.
  126. *
  127. ************************************************************************/
  130. /*____ windowFilterSubband() ____________________________________________*/
  132. void  oldwindowFilterSubband( short * pBuffer, int k, double s[SBLIMIT] )
  133. {
  135.  double y[64];
  136.  int    i,j;
  138.   double  t;
  139.   double *  rpMM;
  141.   int offcache;
  143. double *dp;
  146.   dp = &x[k][0];
  147. offcache = off[k];
  149.   /* replace 32 oldest samples with 32 new samples */
  150.   for ( i=0 ; i<32 ; i++ )
  151. dp[31-i+offcache] = (double) pBuffer[i]/SCALE;
  152.   
  153.   for( i = 0 ; i<64 ; i++ )
  154.   {
  155.     t =  dp[(i+64*0+offcache)&(512-1)] * enwindow[i+64*0];
  156.     t += dp[(i+64*1+offcache)&(512-1)] * enwindow[i+64*1];
  157.     t += dp[(i+64*2+offcache)&(512-1)] * enwindow[i+64*2];
  158.     t += dp[(i+64*3+offcache)&(512-1)] * enwindow[i+64*3];
  159.     t += dp[(i+64*4+offcache)&(512-1)] * enwindow[i+64*4];
  160.     t += dp[(i+64*5+offcache)&(512-1)] * enwindow[i+64*5];
  161.     t += dp[(i+64*6+offcache)&(512-1)] * enwindow[i+64*6];
  162.     t += dp[(i+64*7+offcache)&(512-1)] * enwindow[i+64*7];
  163.     y[i] = t;
  164.   }
  166.   off[k] += 480;              /*offset is modulo (HAN_SIZE-1)*/
  167.   off[k] &= 512-1;
  171.   rpMM = (double *) (m);
  173.   for ( i=0 ; i<SBLIMIT ; i++ )
  174.   {
  175. t = 0;
  176. for( j = 0 ; j < 64 ; j++ )
  177. t += rpMM[j] * y[j];
  178. rpMM += j;
  179. s[i] = t;
  180. }  
  181.  
  182. }
  184. /*____ new windowFilterSubband() ____________________________________________*/
  186. void  windowFilterSubband( short * pBuffer, int k, double s[SBLIMIT] )
  187. {
  188.   double y[64];
  189. int    i,j;
  190. int a;
  192.   double  t;    /* TJ */
  193.   double *  rpMM;    /* TJ */
  195. double *dp, *dp2;
  196. double * pEnw;
  198.    dp = x[k]+off[k]+half[k]*256;
  200. /* replace 32 oldest samples with 32 new samples */
  202. for ( i=0 ; i<32 ; i++ )
  203. dp[(31-i)*8] = (double) pBuffer[i]/SCALE;
  206. for( j = 0 ; j < 2 ; j++ )
  207. {
  209. if( j == 0 )
  210. {
  211. dp = x[k] +half[k]*256;
  212. a = off[k];
  213. }
  214. else
  215. {
  216. if( half[k] == 0 )
  217. {
  218. dp = x[k]+256;
  219. a = off[k];
  220. }
  221. else
  222. {
  223. dp = x[k];
  224. a = (off[k]+1) & 7;
  225. }
  226. }
  228. switch( a )
  229. {
  230. case 0:
  231.   for( i = 0 ; i<32 ; i++ )
  232.    {
  233. dp2 = dp + i*8;
  234. pEnw = enwindow + j*32 + i;
  235.      t =  dp2[0] * pEnw[64*0];
  236.     t += dp2[1] * pEnw[64*1];
  237.     t += dp2[2] * pEnw[64*2];
  238.     t += dp2[3] * pEnw[64*3];
  239.     t += dp2[4] * pEnw[64*4];
  240.     t += dp2[5] * pEnw[64*5];
  241.     t += dp2[6] * pEnw[64*6];
  242.     t += dp2[7] * pEnw[64*7];
  243.     y[i+j*32] = t;
  244.   }
  245. break;
  246. case 1:
  247.   for( i = 0 ; i<32 ; i++ )
  248.    {
  249. dp2 = dp + i*8;
  250. pEnw = enwindow + j*32 + i;
  251.      t =  dp2[1] * pEnw[64*0];
  252.     t += dp2[2] * pEnw[64*1];
  253.     t += dp2[3] * pEnw[64*2];
  254.     t += dp2[4] * pEnw[64*3];
  255.     t += dp2[5] * pEnw[64*4];
  256.     t += dp2[6] * pEnw[64*5];
  257.     t += dp2[7] * pEnw[64*6];
  258.     t += dp2[0] * pEnw[64*7];
  259.     y[i+j*32] = t;
  260.   }
  261. break;
  262. case 2:
  263.   for( i = 0 ; i<32 ; i++ )
  264.    {
  265. dp2 = dp + i*8;
  266. pEnw = enwindow + j*32 + i;
  267.      t =  dp2[2] * pEnw[64*0];
  268.     t += dp2[3] * pEnw[64*1];
  269.     t += dp2[4] * pEnw[64*2];
  270.     t += dp2[5] * pEnw[64*3];
  271.     t += dp2[6] * pEnw[64*4];
  272.     t += dp2[7] * pEnw[64*5];
  273.     t += dp2[0] * pEnw[64*6];
  274.     t += dp2[1] * pEnw[64*7];
  275.     y[i+j*32] = t;
  276.   }
  277. break;
  278. case 3:
  279.   for( i = 0 ; i<32 ; i++ )
  280.    {
  281. dp2 = dp + i*8;
  282. pEnw = enwindow + j*32 + i;
  283.      t =  dp2[3] * pEnw[64*0];
  284.     t += dp2[4] * pEnw[64*1];
  285.     t += dp2[5] * pEnw[64*2];
  286.     t += dp2[6] * pEnw[64*3];
  287.     t += dp2[7] * pEnw[64*4];
  288.     t += dp2[0] * pEnw[64*5];
  289.     t += dp2[1] * pEnw[64*6];
  290.     t += dp2[2] * pEnw[64*7];
  291.     y[i+j*32] = t;
  292.   }
  293. break;
  294. case 4:
  295.   for( i = 0 ; i<32 ; i++ )
  296.    {
  297. dp2 = dp + i*8;
  298. pEnw = enwindow + j*32 + i;
  299.      t =  dp2[4] * pEnw[64*0];
  300.     t += dp2[5] * pEnw[64*1];
  301.     t += dp2[6] * pEnw[64*2];
  302.     t += dp2[7] * pEnw[64*3];
  303.     t += dp2[0] * pEnw[64*4];
  304.     t += dp2[1] * pEnw[64*5];
  305.     t += dp2[2] * pEnw[64*6];
  306.     t += dp2[3] * pEnw[64*7];
  307.     y[i+j*32] = t;
  308.   }
  309. break;
  310. case 5:
  311.   for( i = 0 ; i<32 ; i++ )
  312.    {
  313. dp2 = dp + i*8;
  314. pEnw = enwindow + j*32 + i;
  315.      t =  dp2[5] * pEnw[64*0];
  316.     t += dp2[6] * pEnw[64*1];
  317.     t += dp2[7] * pEnw[64*2];
  318.     t += dp2[0] * pEnw[64*3];
  319.     t += dp2[1] * pEnw[64*4];
  320.     t += dp2[2] * pEnw[64*5];
  321.     t += dp2[3] * pEnw[64*6];
  322.     t += dp2[4] * pEnw[64*7];
  323.     y[i+j*32] = t;
  324.   }
  325. break;
  326. case 6:
  327.   for( i = 0 ; i<32 ; i++ )
  328.    {
  329. dp2 = dp + i*8;
  330. pEnw = enwindow + j*32 + i;
  331.      t =  dp2[6] * pEnw[64*0];
  332.     t += dp2[7] * pEnw[64*1];
  333.     t += dp2[0] * pEnw[64*2];
  334.     t += dp2[1] * pEnw[64*3];
  335.     t += dp2[2] * pEnw[64*4];
  336.     t += dp2[3] * pEnw[64*5];
  337.     t += dp2[4] * pEnw[64*6];
  338.     t += dp2[5] * pEnw[64*7];
  339.     y[i+j*32] = t;
  340.   }
  341. break;
  342. case 7:
  343.   for( i = 0 ; i<32 ; i++ )
  344.    {
  345. dp2 = dp + i*8;
  346. pEnw = enwindow + j*32 + i;
  347.      t =  dp2[7] * pEnw[64*0];
  348.     t += dp2[0] * pEnw[64*1];
  349.     t += dp2[1] * pEnw[64*2];
  350.     t += dp2[2] * pEnw[64*3];
  351.     t += dp2[3] * pEnw[64*4];
  352.     t += dp2[4] * pEnw[64*5];
  353.     t += dp2[5] * pEnw[64*6];
  354.     t += dp2[6] * pEnw[64*7];
  355.     y[i+j*32] = t;
  356.   }
  357. break;
  358. }
  359. }
  361. half[k] = (half[k]+1) & 1;
  362. if( half[k] == 1 )
  363.   off[k] = (off[k]+7) & 7;              /*offset is modulo (HAN_SIZE-1)*/
  365. /* ======================================= */
  367.   rpMM = (double *) (m);
  369.   for ( i=0 ; i<SBLIMIT ; i++ )
  370.   {
  371. t = 0;
  372. for( j = 0 ; j < 64 ; j++ )
  373. t += rpMM[j] * y[j];
  374. rpMM += j;
  375. s[i] = t;
  376. }  
  377. }