apps/codecs/libfaad/sbr_huff.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$
  26 **/
  27
  28 #include "common.h"
  29 #include "structs.h"
  30
  31 #ifdef SBR_DEC
  32
  33 #include "sbr_syntax.h"
  34 #include "bits.h"
  35 #include "sbr_huff.h"
  36 #include "sbr_e_nf.h"
  37
  38
  39 typedef const int8_t (*sbr_huff_tab)[2];
  40
  41 static const int8_t t_huffman_env_1_5dB[120][2] = {
  42     {   1,   2 },    { -64, -65 },    {   3,   4 },    { -63, -66 },
  43     {   5,   6 },    { -62, -67 },    {   7,   8 },    { -61, -68 },
  44     {   9,  10 },    { -60, -69 },    {  11,  12 },    { -59, -70 },
  45     {  13,  14 },    { -58, -71 },    {  15,  16 },    { -57, -72 },
  46     {  17,  18 },    { -73, -56 },    {  19,  21 },    { -74,  20 },
  47     { -55, -75 },    {  22,  26 },    {  23,  24 },    { -54, -76 },
  48     { -77,  25 },    { -53, -78 },    {  27,  34 },    {  28,  29 },
  49     { -52, -79 },    {  30,  31 },    { -80, -51 },    {  32,  33 },
  50     { -83, -82 },    { -81, -50 },    {  35,  57 },    {  36,  40 },
  51     {  37,  38 },    { -88, -84 },    { -48,  39 },    { -90, -85 },
  52     {  41,  46 },    {  42,  43 },    { -49, -87 },    {  44,  45 },
  53     { -89, -86 },    {-124,-123 },    {  47,  50 },    {  48,  49 },
  54     {-122,-121 },    {-120,-119 },    {  51,  54 },    {  52,  53 },
  55     {-118,-117 },    {-116,-115 },    {  55,  56 },    {-114,-113 },
  56     {-112,-111 },    {  58,  89 },    {  59,  74 },    {  60,  67 },
  57     {  61,  64 },    {  62,  63 },    {-110,-109 },    {-108,-107 },
  58     {  65,  66 },    {-106,-105 },    {-104,-103 },    {  68,  71 },
  59     {  69,  70 },    {-102,-101 },    {-100, -99 },    {  72,  73 },
  60     { -98, -97 },    { -96, -95 },    {  75,  82 },    {  76,  79 },
  61     {  77,  78 },    { -94, -93 },    { -92, -91 },    {  80,  81 },
  62     { -47, -46 },    { -45, -44 },    {  83,  86 },    {  84,  85 },
  63     { -43, -42 },    { -41, -40 },    {  87,  88 },    { -39, -38 },
  64     { -37, -36 },    {  90, 105 },    {  91,  98 },    {  92,  95 },
  65     {  93,  94 },    { -35, -34 },    { -33, -32 },    {  96,  97 },
  66     { -31, -30 },    { -29, -28 },    {  99, 102 },    { 100, 101 },
  67     { -27, -26 },    { -25, -24 },    { 103, 104 },    { -23, -22 },
  68     { -21, -20 },    { 106, 113 },    { 107, 110 },    { 108, 109 },
  69     { -19, -18 },    { -17, -16 },    { 111, 112 },    { -15, -14 },
  70     { -13, -12 },    { 114, 117 },    { 115, 116 },    { -11, -10 },
  71     {  -9,  -8 },    { 118, 119 },    {  -7,  -6 },    {  -5,  -4 }
  72 };
  73
  74 static const int8_t f_huffman_env_1_5dB[120][2] = {
  75     {   1,   2 },    { -64, -65 },    {   3,   4 },    { -63, -66 },
  76     {   5,   6 },    { -67, -62 },    {   7,   8 },    { -68, -61 },
  77     {   9,  10 },    { -69, -60 },    {  11,  13 },    { -70,  12 },
  78     { -59, -71 },    {  14,  16 },    { -58,  15 },    { -72, -57 },
  79     {  17,  19 },    { -73,  18 },    { -56, -74 },    {  20,  23 },
  80     {  21,  22 },    { -55, -75 },    { -54, -53 },    {  24,  27 },
  81     {  25,  26 },    { -76, -52 },    { -77, -51 },    {  28,  31 },
  82     {  29,  30 },    { -50, -78 },    { -79, -49 },    {  32,  36 },
  83     {  33,  34 },    { -48, -47 },    { -80,  35 },    { -81, -82 },
  84     {  37,  47 },    {  38,  41 },    {  39,  40 },    { -83, -46 },
  85     { -45, -84 },    {  42,  44 },    { -85,  43 },    { -44, -43 },
  86     {  45,  46 },    { -88, -87 },    { -86, -90 },    {  48,  66 },
  87     {  49,  56 },    {  50,  53 },    {  51,  52 },    { -92, -42 },
  88     { -41, -39 },    {  54,  55 },    {-105, -89 },    { -38, -37 },
  89     {  57,  60 },    {  58,  59 },    { -94, -91 },    { -40, -36 },
  90     {  61,  63 },    { -20,  62 },    {-115,-110 },    {  64,  65 },
  91     {-108,-107 },    {-101, -97 },    {  67,  89 },    {  68,  75 },
  92     {  69,  72 },    {  70,  71 },    { -95, -93 },    { -34, -27 },
  93     {  73,  74 },    { -22, -17 },    { -16,-124 },    {  76,  82 },
  94     {  77,  79 },    {-123,  78 },    {-122,-121 },    {  80,  81 },
  95     {-120,-119 },    {-118,-117 },    {  83,  86 },    {  84,  85 },
  96     {-116,-114 },    {-113,-112 },    {  87,  88 },    {-111,-109 },
  97     {-106,-104 },    {  90, 105 },    {  91,  98 },    {  92,  95 },
  98     {  93,  94 },    {-103,-102 },    {-100, -99 },    {  96,  97 },
  99     { -98, -96 },    { -35, -33 },    {  99, 102 },    { 100, 101 },
 100     { -32, -31 },    { -30, -29 },    { 103, 104 },    { -28, -26 },
 101     { -25, -24 },    { 106, 113 },    { 107, 110 },    { 108, 109 },
 102     { -23, -21 },    { -19, -18 },    { 111, 112 },    { -15, -14 },
 103     { -13, -12 },    { 114, 117 },    { 115, 116 },    { -11, -10 },
 104     {  -9,  -8 },    { 118, 119 },    {  -7,  -6 },    {  -5,  -4 }
 105 };
 106
 107 static const int8_t t_huffman_env_bal_1_5dB[48][2] = {
 108     { -64,   1 },    { -63,   2 },    { -65,   3 },    { -62,   4 },
 109     { -66,   5 },    { -61,   6 },    { -67,   7 },    { -60,   8 },
 110     { -68,   9 },    {  10,  11 },    { -69, -59 },    {  12,  13 },
 111     { -70, -58 },    {  14,  28 },    {  15,  21 },    {  16,  18 },
 112     { -57,  17 },    { -71, -56 },    {  19,  20 },    { -88, -87 },
 113     { -86, -85 },    {  22,  25 },    {  23,  24 },    { -84, -83 },
 114     { -82, -81 },    {  26,  27 },    { -80, -79 },    { -78, -77 },
 115     {  29,  36 },    {  30,  33 },    {  31,  32 },    { -76, -75 },
 116     { -74, -73 },    {  34,  35 },    { -72, -55 },    { -54, -53 },
 117     {  37,  41 },    {  38,  39 },    { -52, -51 },    { -50,  40 },
 118     { -49, -48 },    {  42,  45 },    {  43,  44 },    { -47, -46 },
 119     { -45, -44 },    {  46,  47 },    { -43, -42 },    { -41, -40 }
 120 };
 121
 122 static const int8_t f_huffman_env_bal_1_5dB[48][2] = {
 123     { -64,   1 },    { -65,   2 },    { -63,   3 },    { -66,   4 },
 124     { -62,   5 },    { -61,   6 },    { -67,   7 },    { -68,   8 },
 125     { -60,   9 },    {  10,  11 },    { -69, -59 },    { -70,  12 },
 126     { -58,  13 },    {  14,  17 },    { -71,  15 },    { -57,  16 },
 127     { -56, -73 },    {  18,  32 },    {  19,  25 },    {  20,  22 },
 128     { -72,  21 },    { -88, -87 },    {  23,  24 },    { -86, -85 },
 129     { -84, -83 },    {  26,  29 },    {  27,  28 },    { -82, -81 },
 130     { -80, -79 },    {  30,  31 },    { -78, -77 },    { -76, -75 },
 131     {  33,  40 },    {  34,  37 },    {  35,  36 },    { -74, -55 },
 132     { -54, -53 },    {  38,  39 },    { -52, -51 },    { -50, -49 },
 133     {  41,  44 },    {  42,  43 },    { -48, -47 },    { -46, -45 },
 134     {  45,  46 },    { -44, -43 },    { -42,  47 },    { -41, -40 }
 135 };
 136
 137 static const int8_t t_huffman_env_3_0dB[62][2] = {
 138     { -64,   1 },    { -65,   2 },    { -63,   3 },    { -66,   4 },
 139     { -62,   5 },    { -67,   6 },    { -61,   7 },    { -68,   8 },
 140     { -60,   9 },    {  10,  11 },    { -69, -59 },    {  12,  14 },
 141     { -70,  13 },    { -71, -58 },    {  15,  18 },    {  16,  17 },
 142     { -72, -57 },    { -73, -74 },    {  19,  22 },    { -56,  20 },
 143     { -55,  21 },    { -54, -77 },    {  23,  31 },    {  24,  25 },
 144     { -75, -76 },    {  26,  27 },    { -78, -53 },    {  28,  29 },
 145     { -52, -95 },    { -94,  30 },    { -93, -92 },    {  32,  47 },
 146     {  33,  40 },    {  34,  37 },    {  35,  36 },    { -91, -90 },
 147     { -89, -88 },    {  38,  39 },    { -87, -86 },    { -85, -84 },
 148     {  41,  44 },    {  42,  43 },    { -83, -82 },    { -81, -80 },
 149     {  45,  46 },    { -79, -51 },    { -50, -49 },    {  48,  55 },
 150     {  49,  52 },    {  50,  51 },    { -48, -47 },    { -46, -45 },
 151     {  53,  54 },    { -44, -43 },    { -42, -41 },    {  56,  59 },
 152     {  57,  58 },    { -40, -39 },    { -38, -37 },    {  60,  61 },
 153     { -36, -35 },    { -34, -33 }
 154 };
 155
 156 static const int8_t f_huffman_env_3_0dB[62][2] = {
 157     { -64,   1 },    { -65,   2 },    { -63,   3 },    { -66,   4 },
 158     { -62,   5 },    { -67,   6 },    {   7,   8 },    { -61, -68 },
 159     {   9,  10 },    { -60, -69 },    {  11,  12 },    { -59, -70 },
 160     {  13,  14 },    { -58, -71 },    {  15,  16 },    { -57, -72 },
 161     {  17,  19 },    { -56,  18 },    { -55, -73 },    {  20,  24 },
 162     {  21,  22 },    { -74, -54 },    { -53,  23 },    { -75, -76 },
 163     {  25,  30 },    {  26,  27 },    { -52, -51 },    {  28,  29 },
 164     { -77, -79 },    { -50, -49 },    {  31,  39 },    {  32,  35 },
 165     {  33,  34 },    { -78, -46 },    { -82, -88 },    {  36,  37 },
 166     { -83, -48 },    { -47,  38 },    { -86, -85 },    {  40,  47 },
 167     {  41,  44 },    {  42,  43 },    { -80, -44 },    { -43, -42 },
 168     {  45,  46 },    { -39, -87 },    { -84, -40 },    {  48,  55 },
 169     {  49,  52 },    {  50,  51 },    { -95, -94 },    { -93, -92 },
 170     {  53,  54 },    { -91, -90 },    { -89, -81 },    {  56,  59 },
 171     {  57,  58 },    { -45, -41 },    { -38, -37 },    {  60,  61 },
 172     { -36, -35 },    { -34, -33 }
 173 };
 174
 175 static const int8_t t_huffman_env_bal_3_0dB[24][2] = {
 176     { -64,   1 },    { -63,   2 },    { -65,   3 },    { -66,   4 },
 177     { -62,   5 },    { -61,   6 },    { -67,   7 },    { -68,   8 },
 178     { -60,   9 },    {  10,  16 },    {  11,  13 },    { -69,  12 },
 179     { -76, -75 },    {  14,  15 },    { -74, -73 },    { -72, -71 },
 180     {  17,  20 },    {  18,  19 },    { -70, -59 },    { -58, -57 },
 181     {  21,  22 },    { -56, -55 },    { -54,  23 },    { -53, -52 }
 182 };
 183
 184 static const int8_t f_huffman_env_bal_3_0dB[24][2] = {
 185     { -64,   1 },    { -65,   2 },    { -63,   3 },    { -66,   4 },
 186     { -62,   5 },    { -61,   6 },    { -67,   7 },    { -68,   8 },
 187     { -60,   9 },    {  10,  13 },    { -69,  11 },    { -59,  12 },
 188     { -58, -76 },    {  14,  17 },    {  15,  16 },    { -75, -74 },
 189     { -73, -72 },    {  18,  21 },    {  19,  20 },    { -71, -70 },
 190     { -57, -56 },    {  22,  23 },    { -55, -54 },    { -53, -52 }
 191 };
 192
 193 static const int8_t t_huffman_noise_3_0dB[62][2] = {
 194     { -64,   1 },    { -63,   2 },    { -65,   3 },    { -66,   4 },
 195     { -62,   5 },    { -67,   6 },    {   7,   8 },    { -61, -68 },
 196     {   9,  30 },    {  10,  15 },    { -60,  11 },    { -69,  12 },
 197     {  13,  14 },    { -59, -53 },    { -95, -94 },    {  16,  23 },
 198     {  17,  20 },    {  18,  19 },    { -93, -92 },    { -91, -90 },
 199     {  21,  22 },    { -89, -88 },    { -87, -86 },    {  24,  27 },
 200     {  25,  26 },    { -85, -84 },    { -83, -82 },    {  28,  29 },
 201     { -81, -80 },    { -79, -78 },    {  31,  46 },    {  32,  39 },
 202     {  33,  36 },    {  34,  35 },    { -77, -76 },    { -75, -74 },
 203     {  37,  38 },    { -73, -72 },    { -71, -70 },    {  40,  43 },
 204     {  41,  42 },    { -58, -57 },    { -56, -55 },    {  44,  45 },
 205     { -54, -52 },    { -51, -50 },    {  47,  54 },    {  48,  51 },
 206     {  49,  50 },    { -49, -48 },    { -47, -46 },    {  52,  53 },
 207     { -45, -44 },    { -43, -42 },    {  55,  58 },    {  56,  57 },
 208     { -41, -40 },    { -39, -38 },    {  59,  60 },    { -37, -36 },
 209     { -35,  61 },    { -34, -33 }
 210 };
 211
 212 static const int8_t t_huffman_noise_bal_3_0dB[24][2] = {
 213     { -64,   1 },    { -65,   2 },    { -63,   3 },    {   4,   9 },
 214     { -66,   5 },    { -62,   6 },    {   7,   8 },    { -76, -75 },
 215     { -74, -73 },    {  10,  17 },    {  11,  14 },    {  12,  13 },
 216     { -72, -71 },    { -70, -69 },    {  15,  16 },    { -68, -67 },
 217     { -61, -60 },    {  18,  21 },    {  19,  20 },    { -59, -58 },
 218     { -57, -56 },    {  22,  23 },    { -55, -54 },    { -53, -52 }
 219 };
 220
 221
 222 static INLINE int16_t sbr_huff_dec(bitfile *ld, sbr_huff_tab t_huff)
 223 {
 224     uint8_t bit;
 225     int16_t index = 0;
 226
 227     while (index >= 0)
 228     {
 229         bit = (uint8_t)faad_get1bit(ld);
 230         index = t_huff[index][bit];
 231     }
 232
 233     return index + 64;
 234 }
 235
 236 /* table 10 */
 237 void sbr_envelope(bitfile *ld, sbr_info *sbr, uint8_t ch)
 238 {
 239     uint8_t env, band;
 240     int8_t delta = 0;
 241     sbr_huff_tab t_huff, f_huff;
 242
 243     if ((sbr->L_E[ch] == 1) && (sbr->bs_frame_class[ch] == FIXFIX))
 244         sbr->amp_res[ch] = 0;
 245     else
 246         sbr->amp_res[ch] = sbr->bs_amp_res;
 247
 248     if ((sbr->bs_coupling) && (ch == 1))
 249     {
 250         delta = 1;
 251         if (sbr->amp_res[ch])
 252         {
 253             t_huff = t_huffman_env_bal_3_0dB;
 254             f_huff = f_huffman_env_bal_3_0dB;
 255         } else {
 256             t_huff = t_huffman_env_bal_1_5dB;
 257             f_huff = f_huffman_env_bal_1_5dB;
 258         }
 259     } else {
 260         delta = 0;
 261         if (sbr->amp_res[ch])
 262         {
 263             t_huff = t_huffman_env_3_0dB;
 264             f_huff = f_huffman_env_3_0dB;
 265         } else {
 266             t_huff = t_huffman_env_1_5dB;
 267             f_huff = f_huffman_env_1_5dB;
 268         }
 269     }
 270
 271     for (env = 0; env < sbr->L_E[ch]; env++)
 272     {
 273         if (sbr->bs_df_env[ch][env] == 0)
 274         {
 275             if ((sbr->bs_coupling == 1) && (ch == 1))
 276             {
 277                 if (sbr->amp_res[ch])
 278                 {
 279                     sbr->E[ch][0][env] = (uint16_t)(faad_getbits(ld, 5
 280                         DEBUGVAR(1,272,"sbr_envelope(): bs_data_env")) << delta);
 281                 } else {
 282                     sbr->E[ch][0][env] = (uint16_t)(faad_getbits(ld, 6
 283                         DEBUGVAR(1,273,"sbr_envelope(): bs_data_env")) << delta);
 284                 }
 285             } else {
 286                 if (sbr->amp_res[ch])
 287                 {
 288                     sbr->E[ch][0][env] = (uint16_t)(faad_getbits(ld, 6
 289                         DEBUGVAR(1,274,"sbr_envelope(): bs_data_env")) << delta);
 290                 } else {
 291                     sbr->E[ch][0][env] = (uint16_t)(faad_getbits(ld, 7
 292                         DEBUGVAR(1,275,"sbr_envelope(): bs_data_env")) << delta);
 293                 }
 294             }
 295
 296             for (band = 1; band < sbr->n[sbr->f[ch][env]]; band++)
 297             {
 298                 sbr->E[ch][band][env] = (sbr_huff_dec(ld, f_huff) << delta);
 299             }
 300
 301         } else {
 302             for (band = 0; band < sbr->n[sbr->f[ch][env]]; band++)
 303             {
 304                 sbr->E[ch][band][env] = (sbr_huff_dec(ld, t_huff) << delta);
 305             }
 306         }
 307     }
 308
 309     extract_envelope_data(sbr, ch);
 310 }
 311
 312 /* table 11 */
 313 void sbr_noise(bitfile *ld, sbr_info *sbr, uint8_t ch)
 314 {
 315     uint8_t noise, band;
 316     int8_t delta = 0;
 317     sbr_huff_tab t_huff, f_huff;
 318
 319     if ((sbr->bs_coupling == 1) && (ch == 1))
 320     {
 321         delta = 1;
 322         t_huff = t_huffman_noise_bal_3_0dB;
 323         f_huff = f_huffman_env_bal_3_0dB;
 324     } else {
 325         delta = 0;
 326         t_huff = t_huffman_noise_3_0dB;
 327         f_huff = f_huffman_env_3_0dB;
 328     }
 329
 330     for (noise = 0; noise < sbr->L_Q[ch]; noise++)
 331     {
 332         if(sbr->bs_df_noise[ch][noise] == 0)
 333         {
 334             if ((sbr->bs_coupling == 1) && (ch == 1))
 335             {
 336                 sbr->Q[ch][0][noise] = (faad_getbits(ld, 5
 337                     DEBUGVAR(1,276,"sbr_noise(): bs_data_noise")) << delta);
 338             } else {
 339                 sbr->Q[ch][0][noise] = (faad_getbits(ld, 5
 340                     DEBUGVAR(1,277,"sbr_noise(): bs_data_noise")) << delta);
 341             }
 342             for (band = 1; band < sbr->N_Q; band++)
 343             {
 344                 sbr->Q[ch][band][noise] = (sbr_huff_dec(ld, f_huff) << delta);
 345             }
 346         } else {
 347             for (band = 0; band < sbr->N_Q; band++)
 348             {
 349                 sbr->Q[ch][band][noise] = (sbr_huff_dec(ld, t_huff) << delta);
 350             }
 351         }
 352     }
 353
 354     extract_noise_floor_data(sbr, ch);
 355 }
 356
 357 #endif