drc.c
上传用户:wstnjxml
上传日期:2014-04-03
资源大小:7248k
文件大小:5k
源码类别:

Windows CE

开发平台:

C/C++

  1. /*
  2. ** FAAD2 - Freeware Advanced Audio (AAC) Decoder including SBR decoding
  3. ** Copyright (C) 2003-2004 M. Bakker, Ahead Software AG, http://www.nero.com
  4. **  
  5. ** This program is free software; you can redistribute it and/or modify
  6. ** it under the terms of the GNU General Public License as published by
  7. ** the Free Software Foundation; either version 2 of the License, or
  8. ** (at your option) any later version.
  9. ** 
  10. ** This program is distributed in the hope that it will be useful,
  11. ** but WITHOUT ANY WARRANTY; without even the implied warranty of
  12. ** MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  13. ** GNU General Public License for more details.
  14. ** 
  15. ** You should have received a copy of the GNU General Public License
  16. ** along with this program; if not, write to the Free Software 
  17. ** Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA.
  18. **
  19. ** Any non-GPL usage of this software or parts of this software is strictly
  20. ** forbidden.
  21. **
  22. ** Commercial non-GPL licensing of this software is possible.
  23. ** For more info contact Ahead Software through Mpeg4AAClicense@nero.com.
  24. **
  25. ** $Id: drc.c,v 1.24 2004/09/04 14:56:28 menno Exp $
  26. **/
  27. #include "common.h"
  28. #include "structs.h"
  29. #include <stdlib.h>
  30. #include <string.h>
  31. #include "syntax.h"
  32. #include "drc.h"
  33. drc_info *drc_init(real_t cut, real_t boost)
  34. {
  35.     drc_info *drc = (drc_info*)faad_malloc(sizeof(drc_info));
  36.     memset(drc, 0, sizeof(drc_info));
  37.     drc->ctrl1 = cut;
  38.     drc->ctrl2 = boost;
  39.     drc->num_bands = 1;
  40.     drc->band_top[0] = 1024/4 - 1;
  41.     drc->dyn_rng_sgn[0] = 1;
  42.     drc->dyn_rng_ctl[0] = 0;
  43.     return drc;
  44. }
  45. void drc_end(drc_info *drc)
  46. {
  47.     if (drc) faad_free(drc);
  48. }
  49. #ifdef FIXED_POINT
  50. static const real_t drc_pow2_table[] =
  51. {
  52.     COEF_CONST(0.5146511183),
  53.     COEF_CONST(0.5297315472),
  54.     COEF_CONST(0.5452538663),
  55.     COEF_CONST(0.5612310242),
  56.     COEF_CONST(0.5776763484),
  57.     COEF_CONST(0.5946035575),
  58.     COEF_CONST(0.6120267717),
  59.     COEF_CONST(0.6299605249),
  60.     COEF_CONST(0.6484197773),
  61.     COEF_CONST(0.6674199271),
  62.     COEF_CONST(0.6869768237),
  63.     COEF_CONST(0.7071067812),
  64.     COEF_CONST(0.7278265914),
  65.     COEF_CONST(0.7491535384),
  66.     COEF_CONST(0.7711054127),
  67.     COEF_CONST(0.7937005260),
  68.     COEF_CONST(0.8169577266),
  69.     COEF_CONST(0.8408964153),
  70.     COEF_CONST(0.8655365610),
  71.     COEF_CONST(0.8908987181),
  72.     COEF_CONST(0.9170040432),
  73.     COEF_CONST(0.9438743127),
  74.     COEF_CONST(0.9715319412),
  75.     COEF_CONST(1.0000000000),
  76.     COEF_CONST(1.0293022366),
  77.     COEF_CONST(1.0594630944),
  78.     COEF_CONST(1.0905077327),
  79.     COEF_CONST(1.1224620483),
  80.     COEF_CONST(1.1553526969),
  81.     COEF_CONST(1.1892071150),
  82.     COEF_CONST(1.2240535433),
  83.     COEF_CONST(1.2599210499),
  84.     COEF_CONST(1.2968395547),
  85.     COEF_CONST(1.3348398542),
  86.     COEF_CONST(1.3739536475),
  87.     COEF_CONST(1.4142135624),
  88.     COEF_CONST(1.4556531828),
  89.     COEF_CONST(1.4983070769),
  90.     COEF_CONST(1.5422108254),
  91.     COEF_CONST(1.5874010520),
  92.     COEF_CONST(1.6339154532),
  93.     COEF_CONST(1.6817928305),
  94.     COEF_CONST(1.7310731220),
  95.     COEF_CONST(1.7817974363),
  96.     COEF_CONST(1.8340080864),
  97.     COEF_CONST(1.8877486254),
  98.     COEF_CONST(1.9430638823)
  99. };
  100. #endif
  101. void drc_decode(drc_info *drc, real_t *spec)
  102. {
  103.     uint16_t i, bd, top;
  104. #ifdef FIXED_POINT
  105.     int32_t exp, frac;
  106. #else
  107.     real_t factor, exp;
  108. #endif
  109.     uint16_t bottom = 0;
  110.     if (drc->num_bands == 1)
  111.         drc->band_top[0] = 1024/4 - 1;
  112.     for (bd = 0; bd < drc->num_bands; bd++)
  113.     {
  114.         top = 4 * (drc->band_top[bd] + 1);
  115. #ifndef FIXED_POINT
  116.         /* Decode DRC gain factor */
  117.         if (drc->dyn_rng_sgn[bd])  /* compress */
  118.             exp = -drc->ctrl1 * (drc->dyn_rng_ctl[bd] - (DRC_REF_LEVEL - drc->prog_ref_level))/REAL_CONST(24.0);
  119.         else /* boost */
  120.             exp = drc->ctrl2 * (drc->dyn_rng_ctl[bd] - (DRC_REF_LEVEL - drc->prog_ref_level))/REAL_CONST(24.0);
  121.         factor = (real_t)pow(2.0, exp);
  122.         /* Apply gain factor */
  123.         for (i = bottom; i < top; i++)
  124.             spec[i] *= factor;
  125. #else
  126.         /* Decode DRC gain factor */
  127.         if (drc->dyn_rng_sgn[bd])  /* compress */
  128.         {
  129.             exp = -1 * (drc->dyn_rng_ctl[bd] - (DRC_REF_LEVEL - drc->prog_ref_level))/ 24;
  130.             frac = -1 * (drc->dyn_rng_ctl[bd] - (DRC_REF_LEVEL - drc->prog_ref_level)) % 24;
  131.         } else { /* boost */
  132.             exp = (drc->dyn_rng_ctl[bd] - (DRC_REF_LEVEL - drc->prog_ref_level))/ 24;
  133.             frac = (drc->dyn_rng_ctl[bd] - (DRC_REF_LEVEL - drc->prog_ref_level)) % 24;
  134.         }
  135.         /* Apply gain factor */
  136.         if (exp < 0)
  137.         {
  138.             for (i = bottom; i < top; i++)
  139.             {
  140.                 spec[i] >>= -exp;
  141.                 if (frac)
  142.                     spec[i] = MUL_R(spec[i],drc_pow2_table[frac+23]);
  143.             }
  144.         } else {
  145.             for (i = bottom; i < top; i++)
  146.             {
  147.                 spec[i] <<= exp;
  148.                 if (frac)
  149.                     spec[i] = MUL_R(spec[i],drc_pow2_table[frac+23]);
  150.             }
  151.         }
  152. #endif
  153.         bottom = top;
  154.     }
  155. }