Skip to content

V
vlc-2-2
Commit fec6ded4 authored by Michel Lespinasse's avatar Michel Lespinasse
Browse files
Options

reecriture du deco mpeg audio 

* support du layer 1 mono et stereo
* support du layer 2 stereo (j'ai pas encore fait le mono)
* support des low bitrates proposes par le mpeg2 (non teste)
* passage des tests de compliance mpeg :)
* moins de macros qu'avant, mais je dois encore propriser un peu le code.
parent abe1132d
Hide whitespace changes
Inline Side-by-side
Showing with 990 additions and 1927 deletions
src/audio_decoder/audio_constants.h deleted 100644 → 0
View file @ abe1132d
/*****************************************************************************
* audio_constants.h : defines the MPEG1 Layer I-II audio constants and tables
*****************************************************************************
* Copyright (C) 1999, 2000 VideoLAN
*
* Authors:
*
* This program is free software; you can redistribute it and/or modify
* it under the terms of the GNU General Public License as published by
* the Free Software Foundation; either version 2 of the License, or
* (at your option) any later version.
*
* This program is distributed in the hope that it will be useful,
* but WITHOUT ANY WARRANTY; without even the implied warranty of
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
* General Public License for more details.
*
* You should have received a copy of the GNU General Public
* License along with this program; if not, write to the
* Free Software Foundation, Inc., 59 Temple Place - Suite 330,
* Boston, MA 02111-1307, USA.
*****************************************************************************/
/*****************************************************************************
* 32 bits (4 bytes) audio frame header definitions
*****************************************************************************/
/*
* syncword == `1111 1111 1111'
*/
#define ADEC_HEADER_SYNCWORD_MASK 0xFFF00000
#define ADEC_HEADER_SYNCWORD_SHIFT 20
/* ID :
*
* `0' == reserved
* `1' == ISO/CEI 11172-3
*/
#define ADEC_HEADER_ID_MASK 0x00080000
#define ADEC_HEADER_ID_SHIFT 19
/*
* Layer :
*
* - `00' == reserved
* - `01' == Layer III
* - `10' == Layer II
* - `11' == Layer I
*/
#define ADEC_HEADER_LAYER_MASK 0x00060000
#define ADEC_HEADER_LAYER_SHIFT 17
#define ADEC_HEADER_LAYER_1 0x00060000
#define ADEC_HEADER_LAYER_2 0x00040000
#define ADEC_HEADER_LAYER_3 0x00020000
/* protection_bit */
#define ADEC_HEADER_PROTECTION_BIT_MASK 0x00010000
#define ADEC_HEADER_PROTECTION_BIT_SHIFT 16
/* bitrate_index */
#define ADEC_HEADER_BITRATE_INDEX_MASK 0x0000F000
#define ADEC_HEADER_BITRATE_INDEX_SHIFT 12
/*
* sampling_frequency :
*
* - `00' == 44100 Hz
* - `01' == 48000 Hz
* - `10' == 32000 Hz
* - `11' == reserved
*/
#define ADEC_HEADER_SAMPLING_FREQUENCY_MASK 0x00000C00
#define ADEC_HEADER_SAMPLING_FREQUENCY_SHIFT 10
/* padding_bit */
#define ADEC_HEADER_PADDING_BIT_MASK 0x00000200
#define ADEC_HEADER_PADDING_BIT_SHIFT 9
/* private_bit */
#define ADEC_HEADER_PRIVATE_BIT_MASK 0x00000100
#define ADEC_HEADER_PRIVATE_BIT_SHIFT 8
/*
* mode :
*
* - `00' == stereo (stereo mode)
* - `01' == combined stereo (stereo mode)
* - `10' == two channels (stereo mode)
* - `11' == one channel (mono mode)
*/
#define ADEC_HEADER_MODE_MASK 0x000000C0
#define ADEC_HEADER_MODE_SHIFT 6
/* mode_extension */
#define ADEC_HEADER_MODE_EXTENSION_MASK 0x00000030
#define ADEC_HEADER_MODE_EXTENSION_SHIFT 4
/* copyright */
#define ADEC_HEADER_COPYRIGHT_MASK 0x00000008
#define ADEC_HEADER_COPYRIGHT_SHIFT 3
/* original/copy */
#define ADEC_HEADER_ORIGINAL_COPY_MASK 0x00000004
#define ADEC_HEADER_ORIGINAL_COPY_SHIFT 2
/* emphasis */
#define ADEC_HEADER_EMPHASIS_MASK 0x00000003
#define ADEC_HEADER_EMPHASIS_SHIFT 0
/*****************************************************************************
* frame sizes = f( layer, padding_bit, sampling_frequency, bitrate_index )
*****************************************************************************
* ISO/IEC 11172-3 2.4.3.1 explains how to find out the number of bytes between
* two consecutive syncwords. In order to work out the body size of the frame,
* we just have to substract 4 bytes for the header size.
*
* = Layer I : (slot_size == 4 bytes)
* - padding_bit == 0 :
* frame_size = ( floor( 12 * bitrate(layer, bitrate_index) / sampling_frequency ) * 4 ) - 4
* - padding_bit == 1 :
* frame_size = ( ceil( 12 * bitrate(layer, bitrate_index) / sampling_frequency ) * 4 ) - 4
*
* = Layer II : (slot_size == 1 byte)
* - padding_bit == 0 :
* frame_size = ( floor( 144 * bitrate(layer, bitrate_index) / sampling_frequency ) * 1 ) - 4
* - padding_bit == 1 :
* frame_size = ( ceil( 144 * bitrate(layer, bitrate_index) / sampling_frequency ) * 1 ) - 4
*
* The frame sizes are stored in the following array :
* frame_size = ADEC_FRAME_SIZE[ 128*layer + 64*padding_bit + 16*sampling_frequency + bitrate_index ]
*****************************************************************************/
#if 0
#define ADEC_FRAME_SIZE \
{ \
/* Layer == `00' (reserved) */ \
\
/* padding_bit == `0' */ \
\
/* sampling_frequency == `00' (44100 Hz) */ \
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, \
/* sampling_frequency == `01' (48000 Hz) */ \
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, \
/* sampling_frequency == `10' (32000 Hz) */ \
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, \
/* sampling_frequency == `11' (reserved) */ \
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, \
\
/* padding_bit == `1' */ \
\
/* sampling_frequency == `00' (44100 Hz) */ \
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, \
/* sampling_frequency == `01' (48000 Hz) */ \
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, \
/* sampling_frequency == `10' (32000 Hz) */ \
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, \
/* sampling_frequency == `11' (reserved) */ \
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, \
\
\
/* Layer == `01' (III) */ \
\
/* padding_bit == `0' */ \
\
/* sampling_frequency == `00' (44100 Hz) */ \
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, \
/* sampling_frequency == `01' (48000 Hz) */ \
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, \
/* sampling_frequency == `10' (32000 Hz) */ \
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, \
/* sampling_frequency == `11' (reserved) */ \
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, \
\
/* padding_bit == `1' */ \
\
/* sampling_frequency == `00' (44100 Hz) */ \
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, \
/* sampling_frequency == `01' (48000 Hz) */ \
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, \
/* sampling_frequency == `10' (32000 Hz) */ \
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, \
/* sampling_frequency == `11' (reserved) */ \
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, \
\
\
/* Layer == `10' (II) */ \
\
/* padding_bit == `0' */ \
\
/* sampling_frequency == `00' (44100 Hz) */ \
0, 100, 152, 178, 204, 257, 309, 361, 413, 518, 622, 727, 831, 1040, 1249, 0, \
/* sampling_frequency == `01' (48000 Hz) */ \
0, 92, 140, 164, 188, 236, 284, 332, 380, 476, 572, 668, 764, 956, 1148, 0, \
/* sampling_frequency == `10' (32000 Hz) */ \
0, 140, 212, 248, 284, 356, 428, 500, 572, 716, 860, 1004, 1148, 1436, 1724, 0, \
/* sampling_frequency == `11' (reserved) */ \
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, \
\
/* padding_bit == `1' */ \
\
/* sampling_frequency == `00' (44100 Hz) */ \
0, 101, 153, 179, 205, 258, 310, 362, 414, 519, 623, 728, 832, 1041, 1250, 0, \
/* sampling_frequency == `01' (48000 Hz) */ \
0, 92, 140, 164, 188, 236, 284, 332, 380, 476, 572, 668, 764, 956, 1148, 0, \
/* sampling_frequency == `10' (32000 Hz) */ \
0, 140, 212, 248, 284, 356, 428, 500, 572, 716, 860, 1004, 1148, 1436, 1724, 0, \
/* sampling_frequency == `11' (reserved) */ \
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, \
\
\
/* Layer == `11' (I) */ \
\
/* padding_bit == `0' */ \
\
/* sampling_frequency == `00' (44100 Hz) */ \
0, 28, 64, 100, 132, 168, 204, 236, 272, 308, 344, 376, 412, 448, 480, 0, \
/* sampling_frequency == `01' (48000 Hz) */ \
0, 28, 60, 92, 124, 156, 188, 220, 252, 284, 316, 348, 380, 412, 444, 0, \
/* sampling_frequency == `10' (32000 Hz) */ \
0, 44, 92, 140, 188, 236, 284, 332, 380, 428, 476, 524, 572, 620, 668, 0, \
/* sampling_frequency == `11' (reserved) */ \
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, \
\
/* padding_bit == `1' */ \
\
/* sampling_frequency == `00' (44100 Hz) */ \
0, 32, 68, 104, 136, 172, 208, 240, 276, 312, 348, 380, 416, 452, 484, 0, \
/* sampling_frequency == `01' (48000 Hz) */ \
0, 28, 60, 92, 124, 156, 188, 220, 252, 284, 316, 348, 380, 412, 444, 0, \
/* sampling_frequency == `10' (32000 Hz) */ \
0, 44, 92, 140, 188, 236, 284, 332, 380, 428, 476, 524, 572, 620, 668, 0, \
/* sampling_frequency == `11' (reserved) */ \
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 \
}
#endif
/*****************************************************************************
* scale factors = f( scalefactor ) (Layer I & II, see ISO/IEC 11172-3 2.4.1)
*****************************************************************************
* Theses values are 2^(1 - index/3) (see ISO/IEC 11172-3 2.4.2.5)
*****************************************************************************/
#define ADEC_SCALE_FACTOR \
{ \
/* 0*/ 2.00000000000000, /* 1*/ 1.58740105196820, /* 2*/ 1.25992104989487,\
/* 3*/ 1.00000000000000, /* 4*/ 0.79370052598410, /* 5*/ 0.62996052494744,\
/* 6*/ 0.50000000000000, /* 7*/ 0.39685026299205, /* 8*/ 0.31498026247372,\
/* 9*/ 0.25000000000000, /*10*/ 0.19842513149602, /*11*/ 0.15749013123686,\
/*12*/ 0.12500000000000, /*13*/ 0.09921256574801, /*14*/ 0.07874506561843,\
/*15*/ 0.06250000000000, /*16*/ 0.04960628287401, /*17*/ 0.03937253280921,\
/*18*/ 0.03125000000000, /*19*/ 0.02480314143700, /*20*/ 0.01968626640461,\
/*21*/ 0.01562500000000, /*22*/ 0.01240157071850, /*23*/ 0.00984313320230,\
/*24*/ 0.00781250000000, /*25*/ 0.00620078535925, /*26*/ 0.00492156660115,\
/*27*/ 0.00390625000000, /*28*/ 0.00310039267963, /*29*/ 0.00246078330058,\
/*30*/ 0.00195312500000, /*31*/ 0.00155019633981, /*32*/ 0.00123039165029,\
/*33*/ 0.00097656250000, /*34*/ 0.00077509816991, /*35*/ 0.00061519582514,\
/*36*/ 0.00048828125000, /*37*/ 0.00038754908495, /*38*/ 0.00030759791257,\
/*39*/ 0.00024414062500, /*40*/ 0.00019377454248, /*41*/ 0.00015379895629,\
/*42*/ 0.00012207031250, /*43*/ 0.00009688727124, /*44*/ 0.00007689947814,\
/*45*/ 0.00006103515625, /*46*/ 0.00004844363562, /*47*/ 0.00003844973907,\
/*48*/ 0.00003051757813, /*49*/ 0.00002422181781, /*50*/ 0.00001922486954,\
/*51*/ 0.00001525878906, /*52*/ 0.00001211090890, /*53*/ 0.00000961243477,\
/*54*/ 0.00000762939453, /*55*/ 0.00000605545445, /*56*/ 0.00000480621738,\
/*57*/ 0.00000381469727, /*58*/ 0.00000302772723, /*59*/ 0.00000240310869,\
/*60*/ 0.00000190734863, /*61*/ 0.00000151386361, /*62*/ 0.00000120155435,\
/*63*/ 0.0 /* XXX?? invalid scale factor ? */ \
}
/*****************************************************************************
* Layer I definitions
*****************************************************************************/
/*
* slope table = f( allocation[ch][sb] (see ISO/IEC 11171-3 2.4.1) )
*/
#define ADEC_LAYER1_SLOPE \
{ \
/* 0*/ 0.0, /* no sample */ \
/* 1*/ 2.0/3, \
/* 2*/ 2.0/7, \
/* 3*/ 2.0/15, \
/* 4*/ 2.0/31, \
/* 5*/ 2.0/63, \
/* 6*/ 2.0/127, \
/* 7*/ 2.0/255, \
/* 8*/ 2.0/511, \
/* 9*/ 2.0/1023, \
/*10*/ 2.0/2047, \
/*11*/ 2.0/4095, \
/*12*/ 2.0/8191, \
/*13*/ 2.0/16383, \
/*14*/ 2.0/32767, \
/*15*/ 0.0 /* invalid bit allocation */ \
}
/*
* offset table = f( allocation[ch][sb] (see ISO/IEC 11172-3 2.4.1) )
*/
#define ADEC_LAYER1_OFFSET \
{ \
/* 0*/ 0.0, /* no sample */ \
/* 1*/ -2.0/3, \
/* 2*/ -6.0/7, \
/* 3*/ -14.0/15, \
/* 4*/ -30.0/31, \
/* 5*/ -62.0/63, \
/* 6*/ -126.0/127, \
/* 7*/ -254.0/255, \
/* 8*/ -510.0/511, \
/* 9*/ -1022.0/1023, \
/*10*/ -2046.0/2047, \
/*11*/ -4094.0/4095, \
/*12*/ -8190.0/8191, \
/*13*/ -16382.0/16383, \
/*14*/ -32766.0/32767, \
/*15*/ 0.0 /* invalid bit allocation */ \
}
/*****************************************************************************
* Layer II definitions
*****************************************************************************/
/*
* Bitrate PER CHANNEL index = f( mode, bitrate_index )
* (see ISO/IEC 11172-3 2.4.2.3)
*
* - This index is used in the ADEC_LAYER2_SBLIMIT and ADEC_LAYER2_NBAL tables.
*
* - 0 == forbidden mode/bitrate_index combination
* - 1 == 32 kbits/s per channel
* - 2 == 48 kbits/s per channel
* - 3 == 56 kbits/s per channel
* - 4 == 64 kbits/s per channel
* - 5 == 80 kbits/s per channel
* - 6 == 96 kbits/s per channel
* - 7 == 112 kbits/s per channel
* - 8 == 128 kbits/s per channel
* - 9 == 160 kbits/s per channel
* - 10 == 192 kbits/s per channel
*/
#define ADEC_LAYER2_BITRATE_PER_CHANNEL_INDEX \
{ \
/* mode == `00' (stereo) */ \
{ 0, 0, 0, 0, 1, 0, 2, 3, 4, 5, 6, 7, 8, 9, 10 }, \
/* mode == `01' (combined stereo) */ \
{ 0, 0, 0, 0, 1, 0, 2, 3, 4, 5, 6, 7, 8, 9, 10 }, \
/* mode == `10' (two channels) */ \
{ 0, 0, 0, 0, 1, 0, 2, 3, 4, 5, 6, 7, 8, 9, 10 }, \
/* mode == `11' (one channel) */ \
{ 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 0, 0, 0, 0 } \
}
/*
* Number of subbands = f( sampling_frequency, bitrate_per_channel_index )
* (see ISO/IEC 11172-3 Annex B.2)
*/
#define ADEC_LAYER2_SBLIMIT \
{ \
/* sampling_frequency == `00' (44100 Hz) */ \
{ 0, 8, 8, 27, 27, 27, 30, 30, 30, 30, 30}, \
/* sampling_frequency == `01' (48000 Hz) */ \
{ 0, 8, 8, 27, 27, 27, 27, 27, 27, 27, 27}, \
/* sampling_frequency == `10' (32000 Hz) */ \
{ 0, 12, 12, 27, 27, 27, 30, 30, 30, 30, 30} \
}
/*
* Number of bits allocated = f( bitrate_per_channel_index, subband )
* (see ISO/IEC 11172-3 Annex B.2)
*/
#define ADEC_LAYER2_NBAL \
{ \
/* bitrate_per_channel_index <= 2 */ \
{ 4, 4, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}, \
/* bitrate_per_channel_index > 2 */ \
{ 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 2, 2, 2, 2, 2, 2, 2, 0, 0} \
}
/*
* = When 3 samples are grouped in one codeword, we have to ungroup them in
* order to dequantize the samples (see ISO/IEC 11172-3 2.4.3.3.4) :
* 1) Ungrouping
* for ( i = 0; i < 3; i++ )
* {
* s[i] = c % nlevels;
* c = c / nlevels;
* }
* 2) Requantization
*
* = We pre-calculated all this, and stored the results in the following
* ungroup`nlevels' tables. ISO/IEC 11172-3 Annex B.4 tells us that the
* samples are grouped only when nlevels == 3, 5, or 9.
*
* = ADEC_LAYER2_UNGROUPn = f(3 * n*n*n)
*/
#define ADEC_LAYER2_UNGROUP3 \
{ \
-2.0/3, -2.0/3, -2.0/3, \
.0 , -2.0/3, -2.0/3, \
2.0/3, -2.0/3, -2.0/3, \
-2.0/3, .0 , -2.0/3, \
.0 , .0 , -2.0/3, \
2.0/3, .0 , -2.0/3, \
-2.0/3, 2.0/3, -2.0/3, \
.0 , 2.0/3, -2.0/3, \
2.0/3, 2.0/3, -2.0/3, \
-2.0/3, -2.0/3, .0 , \
.0 , -2.0/3, .0 , \
2.0/3, -2.0/3, .0 , \
-2.0/3, .0 , .0 , \
.0 , .0 , .0 , \
2.0/3, .0 , .0 , \
-2.0/3, 2.0/3, .0 , \
.0 , 2.0/3, .0 , \
2.0/3, 2.0/3, .0 , \
-2.0/3, -2.0/3, 2.0/3, \
.0 , -2.0/3, 2.0/3, \
2.0/3, -2.0/3, 2.0/3, \
-2.0/3, .0 , 2.0/3, \
.0 , .0 , 2.0/3, \
2.0/3, .0 , 2.0/3, \
-2.0/3, 2.0/3, 2.0/3, \
.0 , 2.0/3, 2.0/3, \
2.0/3, 2.0/3, 2.0/3 \
}
#define ADEC_LAYER2_UNGROUP5 \
{ \
-.8, -.8, -.8, \
-.4, -.8, -.8, \
.0, -.8, -.8, \
.4, -.8, -.8, \
.8, -.8, -.8, \
-.8, -.4, -.8, \
-.4, -.4, -.8, \
.0, -.4, -.8, \
.4, -.4, -.8, \
.8, -.4, -.8, \
-.8, .0, -.8, \
-.4, .0, -.8, \
.0, .0, -.8, \
.4, .0, -.8, \
.8, .0, -.8, \
-.8, .4, -.8, \
-.4, .4, -.8, \
.0, .4, -.8, \
.4, .4, -.8, \
.8, .4, -.8, \
-.8, .8, -.8, \
-.4, .8, -.8, \
.0, .8, -.8, \
.4, .8, -.8, \
.8, .8, -.8, \
-.8, -.8, -.4, \
-.4, -.8, -.4, \
.0, -.8, -.4, \
.4, -.8, -.4, \
.8, -.8, -.4, \
-.8, -.4, -.4, \
-.4, -.4, -.4, \
.0, -.4, -.4, \
.4, -.4, -.4, \
.8, -.4, -.4, \
-.8, .0, -.4, \
-.4, .0, -.4, \
.0, .0, -.4, \
.4, .0, -.4, \
.8, .0, -.4, \
-.8, .4, -.4, \
-.4, .4, -.4, \
.0, .4, -.4, \
.4, .4, -.4, \
.8, .4, -.4, \
-.8, .8, -.4, \
-.4, .8, -.4, \
.0, .8, -.4, \
.4, .8, -.4, \
.8, .8, -.4, \
-.8, -.8, .0, \
-.4, -.8, .0, \
.0, -.8, .0, \
.4, -.8, .0, \
.8, -.8, .0, \
-.8, -.4, .0, \
-.4, -.4, .0, \
.0, -.4, .0, \
.4, -.4, .0, \
.8, -.4, .0, \
-.8, .0, .0, \
-.4, .0, .0, \
.0, .0, .0, \
.4, .0, .0, \
.8, .0, .0, \
-.8, .4, .0, \
-.4, .4, .0, \
.0, .4, .0, \
.4, .4, .0, \
.8, .4, .0, \
-.8, .8, .0, \
-.4, .8, .0, \
.0, .8, .0, \
.4, .8, .0, \
.8, .8, .0, \
-.8, -.8, .4, \
-.4, -.8, .4, \
.0, -.8, .4, \
.4, -.8, .4, \
.8, -.8, .4, \
-.8, -.4, .4, \
-.4, -.4, .4, \
.0, -.4, .4, \
.4, -.4, .4, \
.8, -.4, .4, \
-.8, .0, .4, \
-.4, .0, .4, \
.0, .0, .4, \
.4, .0, .4, \
.8, .0, .4, \
-.8, .4, .4, \
-.4, .4, .4, \
.0, .4, .4, \
.4, .4, .4, \
.8, .4, .4, \
-.8, .8, .4, \
-.4, .8, .4, \
.0, .8, .4, \
.4, .8, .4, \
.8, .8, .4, \
-.8, -.8, .8, \
-.4, -.8, .8, \
.0, -.8, .8, \
.4, -.8, .8, \
.8, -.8, .8, \
-.8, -.4, .8, \
-.4, -.4, .8, \
.0, -.4, .8, \
.4, -.4, .8, \
.8, -.4, .8, \
-.8, .0, .8, \
-.4, .0, .8, \
.0, .0, .8, \
.4, .0, .8, \
.8, .0, .8, \
-.8, .4, .8, \
-.4, .4, .8, \
.0, .4, .8, \
.4, .4, .8, \
.8, .4, .8, \
-.8, .8, .8, \
-.4, .8, .8, \
.0, .8, .8, \
.4, .8, .8, \
.8, .8, .8 \
}
#define ADEC_LAYER2_UNGROUP9 \
{ \
-8.0/9, -8.0/9, -8.0/9, \
-6.0/9, -8.0/9, -8.0/9, \
-4.0/9, -8.0/9, -8.0/9, \
-2.0/9, -8.0/9, -8.0/9, \
.0 , -8.0/9, -8.0/9, \
2.0/9, -8.0/9, -8.0/9, \
4.0/9, -8.0/9, -8.0/9, \
6.0/9, -8.0/9, -8.0/9, \
8.0/9, -8.0/9, -8.0/9, \
-8.0/9, -6.0/9, -8.0/9, \
-6.0/9, -6.0/9, -8.0/9, \
-4.0/9, -6.0/9, -8.0/9, \
-2.0/9, -6.0/9, -8.0/9, \
.0 , -6.0/9, -8.0/9, \
2.0/9, -6.0/9, -8.0/9, \
4.0/9, -6.0/9, -8.0/9, \
6.0/9, -6.0/9, -8.0/9, \
8.0/9, -6.0/9, -8.0/9, \
-8.0/9, -4.0/9, -8.0/9, \
-6.0/9, -4.0/9, -8.0/9, \
-4.0/9, -4.0/9, -8.0/9, \
-2.0/9, -4.0/9, -8.0/9, \
.0 , -4.0/9, -8.0/9, \
2.0/9, -4.0/9, -8.0/9, \
4.0/9, -4.0/9, -8.0/9, \
6.0/9, -4.0/9, -8.0/9, \
8.0/9, -4.0/9, -8.0/9, \
-8.0/9, -2.0/9, -8.0/9, \
-6.0/9, -2.0/9, -8.0/9, \
-4.0/9, -2.0/9, -8.0/9, \
-2.0/9, -2.0/9, -8.0/9, \
.0 , -2.0/9, -8.0/9, \
2.0/9, -2.0/9, -8.0/9, \
4.0/9, -2.0/9, -8.0/9, \
6.0/9, -2.0/9, -8.0/9, \
8.0/9, -2.0/9, -8.0/9, \
-8.0/9, .0 , -8.0/9, \
-6.0/9, .0 , -8.0/9, \
-4.0/9, .0 , -8.0/9, \
-2.0/9, .0 , -8.0/9, \
.0 , .0 , -8.0/9, \
2.0/9, .0 , -8.0/9, \
4.0/9, .0 , -8.0/9, \
6.0/9, .0 , -8.0/9, \
8.0/9, .0 , -8.0/9, \
-8.0/9, 2.0/9, -8.0/9, \
-6.0/9, 2.0/9, -8.0/9, \
-4.0/9, 2.0/9, -8.0/9, \
-2.0/9, 2.0/9, -8.0/9, \
.0 , 2.0/9, -8.0/9, \
2.0/9, 2.0/9, -8.0/9, \
4.0/9, 2.0/9, -8.0/9, \
6.0/9, 2.0/9, -8.0/9, \
8.0/9, 2.0/9, -8.0/9, \
-8.0/9, 4.0/9, -8.0/9, \
-6.0/9, 4.0/9, -8.0/9, \
-4.0/9, 4.0/9, -8.0/9, \
-2.0/9, 4.0/9, -8.0/9, \
.0 , 4.0/9, -8.0/9, \
2.0/9, 4.0/9, -8.0/9, \
4.0/9, 4.0/9, -8.0/9, \
6.0/9, 4.0/9, -8.0/9, \
8.0/9, 4.0/9, -8.0/9, \
-8.0/9, 6.0/9, -8.0/9, \
-6.0/9, 6.0/9, -8.0/9, \
-4.0/9, 6.0/9, -8.0/9, \
-2.0/9, 6.0/9, -8.0/9, \
.0 , 6.0/9, -8.0/9, \
2.0/9, 6.0/9, -8.0/9, \
4.0/9, 6.0/9, -8.0/9, \
6.0/9, 6.0/9, -8.0/9, \
8.0/9, 6.0/9, -8.0/9, \
-8.0/9, 8.0/9, -8.0/9, \
-6.0/9, 8.0/9, -8.0/9, \
-4.0/9, 8.0/9, -8.0/9, \
-2.0/9, 8.0/9, -8.0/9, \
.0 , 8.0/9, -8.0/9, \
2.0/9, 8.0/9, -8.0/9, \
4.0/9, 8.0/9, -8.0/9, \
6.0/9, 8.0/9, -8.0/9, \
8.0/9, 8.0/9, -8.0/9, \
-8.0/9, -8.0/9, -6.0/9, \
-6.0/9, -8.0/9, -6.0/9, \
-4.0/9, -8.0/9, -6.0/9, \
-2.0/9, -8.0/9, -6.0/9, \
.0 , -8.0/9, -6.0/9, \
2.0/9, -8.0/9, -6.0/9, \
4.0/9, -8.0/9, -6.0/9, \
6.0/9, -8.0/9, -6.0/9, \
8.0/9, -8.0/9, -6.0/9, \
-8.0/9, -6.0/9, -6.0/9, \
-6.0/9, -6.0/9, -6.0/9, \
-4.0/9, -6.0/9, -6.0/9, \
-2.0/9, -6.0/9, -6.0/9, \
.0 , -6.0/9, -6.0/9, \
2.0/9, -6.0/9, -6.0/9, \
4.0/9, -6.0/9, -6.0/9, \
6.0/9, -6.0/9, -6.0/9, \
8.0/9, -6.0/9, -6.0/9, \
-8.0/9, -4.0/9, -6.0/9, \
-6.0/9, -4.0/9, -6.0/9, \
-4.0/9, -4.0/9, -6.0/9, \
-2.0/9, -4.0/9, -6.0/9, \
.0 , -4.0/9, -6.0/9, \
2.0/9, -4.0/9, -6.0/9, \
4.0/9, -4.0/9, -6.0/9, \
6.0/9, -4.0/9, -6.0/9, \
8.0/9, -4.0/9, -6.0/9, \
-8.0/9, -2.0/9, -6.0/9, \
-6.0/9, -2.0/9, -6.0/9, \
-4.0/9, -2.0/9, -6.0/9, \
-2.0/9, -2.0/9, -6.0/9, \
.0 , -2.0/9, -6.0/9, \
2.0/9, -2.0/9, -6.0/9, \
4.0/9, -2.0/9, -6.0/9, \
6.0/9, -2.0/9, -6.0/9, \
8.0/9, -2.0/9, -6.0/9, \
-8.0/9, .0 , -6.0/9, \
-6.0/9, .0 , -6.0/9, \
-4.0/9, .0 , -6.0/9, \
-2.0/9, .0 , -6.0/9, \
.0 , .0 , -6.0/9, \
2.0/9, .0 , -6.0/9, \
4.0/9, .0 , -6.0/9, \
6.0/9, .0 , -6.0/9, \
8.0/9, .0 , -6.0/9, \
-8.0/9, 2.0/9, -6.0/9, \
-6.0/9, 2.0/9, -6.0/9, \
-4.0/9, 2.0/9, -6.0/9, \
-2.0/9, 2.0/9, -6.0/9, \
.0 , 2.0/9, -6.0/9, \
2.0/9, 2.0/9, -6.0/9, \
4.0/9, 2.0/9, -6.0/9, \
6.0/9, 2.0/9, -6.0/9, \
8.0/9, 2.0/9, -6.0/9, \
-8.0/9, 4.0/9, -6.0/9, \
-6.0/9, 4.0/9, -6.0/9, \
-4.0/9, 4.0/9, -6.0/9, \
-2.0/9, 4.0/9, -6.0/9, \
.0 , 4.0/9, -6.0/9, \
2.0/9, 4.0/9, -6.0/9, \
4.0/9, 4.0/9, -6.0/9, \
6.0/9, 4.0/9, -6.0/9, \
8.0/9, 4.0/9, -6.0/9, \
-8.0/9, 6.0/9, -6.0/9, \
-6.0/9, 6.0/9, -6.0/9, \
-4.0/9, 6.0/9, -6.0/9, \
-2.0/9, 6.0/9, -6.0/9, \
.0 , 6.0/9, -6.0/9, \
2.0/9, 6.0/9, -6.0/9, \
4.0/9, 6.0/9, -6.0/9, \
6.0/9, 6.0/9, -6.0/9, \
8.0/9, 6.0/9, -6.0/9, \
-8.0/9, 8.0/9, -6.0/9, \
-6.0/9, 8.0/9, -6.0/9, \
-4.0/9, 8.0/9, -6.0/9, \
-2.0/9, 8.0/9, -6.0/9, \
.0 , 8.0/9, -6.0/9, \
2.0/9, 8.0/9, -6.0/9, \
4.0/9, 8.0/9, -6.0/9, \
6.0/9, 8.0/9, -6.0/9, \
8.0/9, 8.0/9, -6.0/9, \
-8.0/9, -8.0/9, -4.0/9, \
-6.0/9, -8.0/9, -4.0/9, \
-4.0/9, -8.0/9, -4.0/9, \
-2.0/9, -8.0/9, -4.0/9, \
.0 , -8.0/9, -4.0/9, \
2.0/9, -8.0/9, -4.0/9, \
4.0/9, -8.0/9, -4.0/9, \
6.0/9, -8.0/9, -4.0/9, \
8.0/9, -8.0/9, -4.0/9, \
-8.0/9, -6.0/9, -4.0/9, \
-6.0/9, -6.0/9, -4.0/9, \
-4.0/9, -6.0/9, -4.0/9, \
-2.0/9, -6.0/9, -4.0/9, \
.0 , -6.0/9, -4.0/9, \
2.0/9, -6.0/9, -4.0/9, \
4.0/9, -6.0/9, -4.0/9, \
6.0/9, -6.0/9, -4.0/9, \
8.0/9, -6.0/9, -4.0/9, \
-8.0/9, -4.0/9, -4.0/9, \
-6.0/9, -4.0/9, -4.0/9, \
-4.0/9, -4.0/9, -4.0/9, \
-2.0/9, -4.0/9, -4.0/9, \
.0 , -4.0/9, -4.0/9, \
2.0/9, -4.0/9, -4.0/9, \
4.0/9, -4.0/9, -4.0/9, \
6.0/9, -4.0/9, -4.0/9, \
8.0/9, -4.0/9, -4.0/9, \
-8.0/9, -2.0/9, -4.0/9, \
-6.0/9, -2.0/9, -4.0/9, \
-4.0/9, -2.0/9, -4.0/9, \
-2.0/9, -2.0/9, -4.0/9, \
.0 , -2.0/9, -4.0/9, \
2.0/9, -2.0/9, -4.0/9, \
4.0/9, -2.0/9, -4.0/9, \
6.0/9, -2.0/9, -4.0/9, \
8.0/9, -2.0/9, -4.0/9, \
-8.0/9, .0 , -4.0/9, \
-6.0/9, .0 , -4.0/9, \
-4.0/9, .0 , -4.0/9, \
-2.0/9, .0 , -4.0/9, \
.0 , .0 , -4.0/9, \
2.0/9, .0 , -4.0/9, \
4.0/9, .0 , -4.0/9, \
6.0/9, .0 , -4.0/9, \
8.0/9, .0 , -4.0/9, \
-8.0/9, 2.0/9, -4.0/9, \
-6.0/9, 2.0/9, -4.0/9, \
-4.0/9, 2.0/9, -4.0/9, \
-2.0/9, 2.0/9, -4.0/9, \
.0 , 2.0/9, -4.0/9, \
2.0/9, 2.0/9, -4.0/9, \
4.0/9, 2.0/9, -4.0/9, \
6.0/9, 2.0/9, -4.0/9, \
8.0/9, 2.0/9, -4.0/9, \
-8.0/9, 4.0/9, -4.0/9, \
-6.0/9, 4.0/9, -4.0/9, \
-4.0/9, 4.0/9, -4.0/9, \
-2.0/9, 4.0/9, -4.0/9, \
.0 , 4.0/9, -4.0/9, \
2.0/9, 4.0/9, -4.0/9, \
4.0/9, 4.0/9, -4.0/9, \
6.0/9, 4.0/9, -4.0/9, \
8.0/9, 4.0/9, -4.0/9, \
-8.0/9, 6.0/9, -4.0/9, \
-6.0/9, 6.0/9, -4.0/9, \
-4.0/9, 6.0/9, -4.0/9, \
-2.0/9, 6.0/9, -4.0/9, \
.0 , 6.0/9, -4.0/9, \
2.0/9, 6.0/9, -4.0/9, \
4.0/9, 6.0/9, -4.0/9, \
6.0/9, 6.0/9, -4.0/9, \
8.0/9, 6.0/9, -4.0/9, \
-8.0/9, 8.0/9, -4.0/9, \
-6.0/9, 8.0/9, -4.0/9, \
-4.0/9, 8.0/9, -4.0/9, \
-2.0/9, 8.0/9, -4.0/9, \
.0 , 8.0/9, -4.0/9, \
2.0/9, 8.0/9, -4.0/9, \
4.0/9, 8.0/9, -4.0/9, \
6.0/9, 8.0/9, -4.0/9, \
8.0/9, 8.0/9, -4.0/9, \
-8.0/9, -8.0/9, -2.0/9, \
-6.0/9, -8.0/9, -2.0/9, \
-4.0/9, -8.0/9, -2.0/9, \
-2.0/9, -8.0/9, -2.0/9, \
.0 , -8.0/9, -2.0/9, \
2.0/9, -8.0/9, -2.0/9, \
4.0/9, -8.0/9, -2.0/9, \
6.0/9, -8.0/9, -2.0/9, \
8.0/9, -8.0/9, -2.0/9, \
-8.0/9, -6.0/9, -2.0/9, \
-6.0/9, -6.0/9, -2.0/9, \
-4.0/9, -6.0/9, -2.0/9, \
-2.0/9, -6.0/9, -2.0/9, \
.0 , -6.0/9, -2.0/9, \
2.0/9, -6.0/9, -2.0/9, \
4.0/9, -6.0/9, -2.0/9, \
6.0/9, -6.0/9, -2.0/9, \
8.0/9, -6.0/9, -2.0/9, \
-8.0/9, -4.0/9, -2.0/9, \
-6.0/9, -4.0/9, -2.0/9, \
-4.0/9, -4.0/9, -2.0/9, \
-2.0/9, -4.0/9, -2.0/9, \
.0 , -4.0/9, -2.0/9, \
2.0/9, -4.0/9, -2.0/9, \
4.0/9, -4.0/9, -2.0/9, \
6.0/9, -4.0/9, -2.0/9, \
8.0/9, -4.0/9, -2.0/9, \
-8.0/9, -2.0/9, -2.0/9, \
-6.0/9, -2.0/9, -2.0/9, \
-4.0/9, -2.0/9, -2.0/9, \
-2.0/9, -2.0/9, -2.0/9, \
.0 , -2.0/9, -2.0/9, \
2.0/9, -2.0/9, -2.0/9, \
4.0/9, -2.0/9, -2.0/9, \
6.0/9, -2.0/9, -2.0/9, \
8.0/9, -2.0/9, -2.0/9, \
-8.0/9, .0 , -2.0/9, \
-6.0/9, .0 , -2.0/9, \
-4.0/9, .0 , -2.0/9, \
-2.0/9, .0 , -2.0/9, \
.0 , .0 , -2.0/9, \
2.0/9, .0 , -2.0/9, \
4.0/9, .0 , -2.0/9, \
6.0/9, .0 , -2.0/9, \
8.0/9, .0 , -2.0/9, \
-8.0/9, 2.0/9, -2.0/9, \
-6.0/9, 2.0/9, -2.0/9, \
-4.0/9, 2.0/9, -2.0/9, \
-2.0/9, 2.0/9, -2.0/9, \
.0 , 2.0/9, -2.0/9, \
2.0/9, 2.0/9, -2.0/9, \
4.0/9, 2.0/9, -2.0/9, \
6.0/9, 2.0/9, -2.0/9, \
8.0/9, 2.0/9, -2.0/9, \
-8.0/9, 4.0/9, -2.0/9, \
-6.0/9, 4.0/9, -2.0/9, \
-4.0/9, 4.0/9, -2.0/9, \
-2.0/9, 4.0/9, -2.0/9, \
.0 , 4.0/9, -2.0/9, \
2.0/9, 4.0/9, -2.0/9, \
4.0/9, 4.0/9, -2.0/9, \
6.0/9, 4.0/9, -2.0/9, \
8.0/9, 4.0/9, -2.0/9, \
-8.0/9, 6.0/9, -2.0/9, \
-6.0/9, 6.0/9, -2.0/9, \
-4.0/9, 6.0/9, -2.0/9, \
-2.0/9, 6.0/9, -2.0/9, \
.0 , 6.0/9, -2.0/9, \
2.0/9, 6.0/9, -2.0/9, \
4.0/9, 6.0/9, -2.0/9, \
6.0/9, 6.0/9, -2.0/9, \
8.0/9, 6.0/9, -2.0/9, \
-8.0/9, 8.0/9, -2.0/9, \
-6.0/9, 8.0/9, -2.0/9, \
-4.0/9, 8.0/9, -2.0/9, \
-2.0/9, 8.0/9, -2.0/9, \
.0 , 8.0/9, -2.0/9, \
2.0/9, 8.0/9, -2.0/9, \
4.0/9, 8.0/9, -2.0/9, \
6.0/9, 8.0/9, -2.0/9, \
8.0/9, 8.0/9, -2.0/9, \
-8.0/9, -8.0/9, .0 , \
-6.0/9, -8.0/9, .0 , \
-4.0/9, -8.0/9, .0 , \
-2.0/9, -8.0/9, .0 , \
.0 , -8.0/9, .0 , \
2.0/9, -8.0/9, .0 , \
4.0/9, -8.0/9, .0 , \
6.0/9, -8.0/9, .0 , \
8.0/9, -8.0/9, .0 , \
-8.0/9, -6.0/9, .0 , \
-6.0/9, -6.0/9, .0 , \
-4.0/9, -6.0/9, .0 , \
-2.0/9, -6.0/9, .0 , \
.0 , -6.0/9, .0 , \
2.0/9, -6.0/9, .0 , \
4.0/9, -6.0/9, .0 , \
6.0/9, -6.0/9, .0 , \
8.0/9, -6.0/9, .0 , \
-8.0/9, -4.0/9, .0 , \
-6.0/9, -4.0/9, .0 , \
-4.0/9, -4.0/9, .0 , \
-2.0/9, -4.0/9, .0 , \
.0 , -4.0/9, .0 , \
2.0/9, -4.0/9, .0 , \
4.0/9, -4.0/9, .0 , \
6.0/9, -4.0/9, .0 , \
8.0/9, -4.0/9, .0 , \
-8.0/9, -2.0/9, .0 , \
-6.0/9, -2.0/9, .0 , \
-4.0/9, -2.0/9, .0 , \
-2.0/9, -2.0/9, .0 , \
.0 , -2.0/9, .0 , \
2.0/9, -2.0/9, .0 , \
4.0/9, -2.0/9, .0 , \
6.0/9, -2.0/9, .0 , \
8.0/9, -2.0/9, .0 , \
-8.0/9, .0 , .0 , \
-6.0/9, .0 , .0 , \
-4.0/9, .0 , .0 , \
-2.0/9, .0 , .0 , \
.0 , .0 , .0 , \
2.0/9, .0 , .0 , \
4.0/9, .0 , .0 , \
6.0/9, .0 , .0 , \
8.0/9, .0 , .0 , \
-8.0/9, 2.0/9, .0 , \
-6.0/9, 2.0/9, .0 , \
-4.0/9, 2.0/9, .0 , \
-2.0/9, 2.0/9, .0 , \
.0 , 2.0/9, .0 , \
2.0/9, 2.0/9, .0 , \
4.0/9, 2.0/9, .0 , \
6.0/9, 2.0/9, .0 , \
8.0/9, 2.0/9, .0 , \
-8.0/9, 4.0/9, .0 , \
-6.0/9, 4.0/9, .0 , \
-4.0/9, 4.0/9, .0 , \
-2.0/9, 4.0/9, .0 , \
.0 , 4.0/9, .0 , \
2.0/9, 4.0/9, .0 , \
4.0/9, 4.0/9, .0 , \
6.0/9, 4.0/9, .0 , \
8.0/9, 4.0/9, .0 , \
-8.0/9, 6.0/9, .0 , \
-6.0/9, 6.0/9, .0 , \
-4.0/9, 6.0/9, .0 , \
-2.0/9, 6.0/9, .0 , \
.0 , 6.0/9, .0 , \
2.0/9, 6.0/9, .0 , \
4.0/9, 6.0/9, .0 , \
6.0/9, 6.0/9, .0 , \
8.0/9, 6.0/9, .0 , \
-8.0/9, 8.0/9, .0 , \
-6.0/9, 8.0/9, .0 , \
-4.0/9, 8.0/9, .0 , \
-2.0/9, 8.0/9, .0 , \
.0 , 8.0/9, .0 , \
2.0/9, 8.0/9, .0 , \
4.0/9, 8.0/9, .0 , \
6.0/9, 8.0/9, .0 , \
8.0/9, 8.0/9, .0 , \
-8.0/9, -8.0/9, 2.0/9, \
-6.0/9, -8.0/9, 2.0/9, \
-4.0/9, -8.0/9, 2.0/9, \
-2.0/9, -8.0/9, 2.0/9, \
.0 , -8.0/9, 2.0/9, \
2.0/9, -8.0/9, 2.0/9, \
4.0/9, -8.0/9, 2.0/9, \
6.0/9, -8.0/9, 2.0/9, \
8.0/9, -8.0/9, 2.0/9, \
-8.0/9, -6.0/9, 2.0/9, \
-6.0/9, -6.0/9, 2.0/9, \
-4.0/9, -6.0/9, 2.0/9, \
-2.0/9, -6.0/9, 2.0/9, \
.0 , -6.0/9, 2.0/9, \
2.0/9, -6.0/9, 2.0/9, \
4.0/9, -6.0/9, 2.0/9, \
6.0/9, -6.0/9, 2.0/9, \
8.0/9, -6.0/9, 2.0/9, \
-8.0/9, -4.0/9, 2.0/9, \
-6.0/9, -4.0/9, 2.0/9, \
-4.0/9, -4.0/9, 2.0/9, \
-2.0/9, -4.0/9, 2.0/9, \
.0 , -4.0/9, 2.0/9, \
2.0/9, -4.0/9, 2.0/9, \
4.0/9, -4.0/9, 2.0/9, \
6.0/9, -4.0/9, 2.0/9, \
8.0/9, -4.0/9, 2.0/9, \
-8.0/9, -2.0/9, 2.0/9, \
-6.0/9, -2.0/9, 2.0/9, \
-4.0/9, -2.0/9, 2.0/9, \
-2.0/9, -2.0/9, 2.0/9, \
.0 , -2.0/9, 2.0/9, \
2.0/9, -2.0/9, 2.0/9, \
4.0/9, -2.0/9, 2.0/9, \
6.0/9, -2.0/9, 2.0/9, \
8.0/9, -2.0/9, 2.0/9, \
-8.0/9, .0 , 2.0/9, \
-6.0/9, .0 , 2.0/9, \
-4.0/9, .0 , 2.0/9, \
-2.0/9, .0 , 2.0/9, \
.0 , .0 , 2.0/9, \
2.0/9, .0 , 2.0/9, \
4.0/9, .0 , 2.0/9, \
6.0/9, .0 , 2.0/9, \
8.0/9, .0 , 2.0/9, \
-8.0/9, 2.0/9, 2.0/9, \
-6.0/9, 2.0/9, 2.0/9, \
-4.0/9, 2.0/9, 2.0/9, \
-2.0/9, 2.0/9, 2.0/9, \
.0 , 2.0/9, 2.0/9, \
2.0/9, 2.0/9, 2.0/9, \
4.0/9, 2.0/9, 2.0/9, \
6.0/9, 2.0/9, 2.0/9, \
8.0/9, 2.0/9, 2.0/9, \
-8.0/9, 4.0/9, 2.0/9, \
-6.0/9, 4.0/9, 2.0/9, \
-4.0/9, 4.0/9, 2.0/9, \
-2.0/9, 4.0/9, 2.0/9, \
.0 , 4.0/9, 2.0/9, \
2.0/9, 4.0/9, 2.0/9, \
4.0/9, 4.0/9, 2.0/9, \
6.0/9, 4.0/9, 2.0/9, \
8.0/9, 4.0/9, 2.0/9, \
-8.0/9, 6.0/9, 2.0/9, \
-6.0/9, 6.0/9, 2.0/9, \
-4.0/9, 6.0/9, 2.0/9, \
-2.0/9, 6.0/9, 2.0/9, \
.0 , 6.0/9, 2.0/9, \
2.0/9, 6.0/9, 2.0/9, \
4.0/9, 6.0/9, 2.0/9, \
6.0/9, 6.0/9, 2.0/9, \
8.0/9, 6.0/9, 2.0/9, \
-8.0/9, 8.0/9, 2.0/9, \
-6.0/9, 8.0/9, 2.0/9, \
-4.0/9, 8.0/9, 2.0/9, \
-2.0/9, 8.0/9, 2.0/9, \
.0 , 8.0/9, 2.0/9, \
2.0/9, 8.0/9, 2.0/9, \
4.0/9, 8.0/9, 2.0/9, \
6.0/9, 8.0/9, 2.0/9, \
8.0/9, 8.0/9, 2.0/9, \
-8.0/9, -8.0/9, 4.0/9, \
-6.0/9, -8.0/9, 4.0/9, \
-4.0/9, -8.0/9, 4.0/9, \
-2.0/9, -8.0/9, 4.0/9, \
.0 , -8.0/9, 4.0/9, \
2.0/9, -8.0/9, 4.0/9, \
4.0/9, -8.0/9, 4.0/9, \
6.0/9, -8.0/9, 4.0/9, \
8.0/9, -8.0/9, 4.0/9, \
-8.0/9, -6.0/9, 4.0/9, \
-6.0/9, -6.0/9, 4.0/9, \
-4.0/9, -6.0/9, 4.0/9, \
-2.0/9, -6.0/9, 4.0/9, \
.0 , -6.0/9, 4.0/9, \
2.0/9, -6.0/9, 4.0/9, \
4.0/9, -6.0/9, 4.0/9, \
6.0/9, -6.0/9, 4.0/9, \
8.0/9, -6.0/9, 4.0/9, \
-8.0/9, -4.0/9, 4.0/9, \
-6.0/9, -4.0/9, 4.0/9, \
-4.0/9, -4.0/9, 4.0/9, \
-2.0/9, -4.0/9, 4.0/9, \
.0 , -4.0/9, 4.0/9, \
2.0/9, -4.0/9, 4.0/9, \
4.0/9, -4.0/9, 4.0/9, \
6.0/9, -4.0/9, 4.0/9, \
8.0/9, -4.0/9, 4.0/9, \
-8.0/9, -2.0/9, 4.0/9, \
-6.0/9, -2.0/9, 4.0/9, \
-4.0/9, -2.0/9, 4.0/9, \
-2.0/9, -2.0/9, 4.0/9, \
.0 , -2.0/9, 4.0/9, \
2.0/9, -2.0/9, 4.0/9, \
4.0/9, -2.0/9, 4.0/9, \
6.0/9, -2.0/9, 4.0/9, \
8.0/9, -2.0/9, 4.0/9, \
-8.0/9, .0 , 4.0/9, \
-6.0/9, .0 , 4.0/9, \
-4.0/9, .0 , 4.0/9, \
-2.0/9, .0 , 4.0/9, \
.0 , .0 , 4.0/9, \
2.0/9, .0 , 4.0/9, \
4.0/9, .0 , 4.0/9, \
6.0/9, .0 , 4.0/9, \
8.0/9, .0 , 4.0/9, \
-8.0/9, 2.0/9, 4.0/9, \
-6.0/9, 2.0/9, 4.0/9, \
-4.0/9, 2.0/9, 4.0/9, \
-2.0/9, 2.0/9, 4.0/9, \
.0 , 2.0/9, 4.0/9, \
2.0/9, 2.0/9, 4.0/9, \
4.0/9, 2.0/9, 4.0/9, \
6.0/9, 2.0/9, 4.0/9, \
8.0/9, 2.0/9, 4.0/9, \
-8.0/9, 4.0/9, 4.0/9, \
-6.0/9, 4.0/9, 4.0/9, \
-4.0/9, 4.0/9, 4.0/9, \
-2.0/9, 4.0/9, 4.0/9, \
.0 , 4.0/9, 4.0/9, \
2.0/9, 4.0/9, 4.0/9, \
4.0/9, 4.0/9, 4.0/9, \
6.0/9, 4.0/9, 4.0/9, \
8.0/9, 4.0/9, 4.0/9, \
-8.0/9, 6.0/9, 4.0/9, \
-6.0/9, 6.0/9, 4.0/9, \
-4.0/9, 6.0/9, 4.0/9, \
-2.0/9, 6.0/9, 4.0/9, \
.0 , 6.0/9, 4.0/9, \
2.0/9, 6.0/9, 4.0/9, \
4.0/9, 6.0/9, 4.0/9, \
6.0/9, 6.0/9, 4.0/9, \
8.0/9, 6.0/9, 4.0/9, \
-8.0/9, 8.0/9, 4.0/9, \
-6.0/9, 8.0/9, 4.0/9, \
-4.0/9, 8.0/9, 4.0/9, \
-2.0/9, 8.0/9, 4.0/9, \
.0 , 8.0/9, 4.0/9, \
2.0/9, 8.0/9, 4.0/9, \
4.0/9, 8.0/9, 4.0/9, \
6.0/9, 8.0/9, 4.0/9, \
8.0/9, 8.0/9, 4.0/9, \
-8.0/9, -8.0/9, 6.0/9, \
-6.0/9, -8.0/9, 6.0/9, \
-4.0/9, -8.0/9, 6.0/9, \
-2.0/9, -8.0/9, 6.0/9, \
.0 , -8.0/9, 6.0/9, \
2.0/9, -8.0/9, 6.0/9, \
4.0/9, -8.0/9, 6.0/9, \
6.0/9, -8.0/9, 6.0/9, \
8.0/9, -8.0/9, 6.0/9, \
-8.0/9, -6.0/9, 6.0/9, \
-6.0/9, -6.0/9, 6.0/9, \
-4.0/9, -6.0/9, 6.0/9, \
-2.0/9, -6.0/9, 6.0/9, \
.0 , -6.0/9, 6.0/9, \
2.0/9, -6.0/9, 6.0/9, \
4.0/9, -6.0/9, 6.0/9, \
6.0/9, -6.0/9, 6.0/9, \
8.0/9, -6.0/9, 6.0/9, \
-8.0/9, -4.0/9, 6.0/9, \
-6.0/9, -4.0/9, 6.0/9, \
-4.0/9, -4.0/9, 6.0/9, \
-2.0/9, -4.0/9, 6.0/9, \
.0 , -4.0/9, 6.0/9, \
2.0/9, -4.0/9, 6.0/9, \
4.0/9, -4.0/9, 6.0/9, \
6.0/9, -4.0/9, 6.0/9, \
8.0/9, -4.0/9, 6.0/9, \
-8.0/9, -2.0/9, 6.0/9, \
-6.0/9, -2.0/9, 6.0/9, \
-4.0/9, -2.0/9, 6.0/9, \
-2.0/9, -2.0/9, 6.0/9, \
.0 , -2.0/9, 6.0/9, \
2.0/9, -2.0/9, 6.0/9, \
4.0/9, -2.0/9, 6.0/9, \
6.0/9, -2.0/9, 6.0/9, \
8.0/9, -2.0/9, 6.0/9, \
-8.0/9, .0 , 6.0/9, \
-6.0/9, .0 , 6.0/9, \
-4.0/9, .0 , 6.0/9, \
-2.0/9, .0 , 6.0/9, \
.0 , .0 , 6.0/9, \
2.0/9, .0 , 6.0/9, \
4.0/9, .0 , 6.0/9, \
6.0/9, .0 , 6.0/9, \
8.0/9, .0 , 6.0/9, \
-8.0/9, 2.0/9, 6.0/9, \
-6.0/9, 2.0/9, 6.0/9, \
-4.0/9, 2.0/9, 6.0/9, \
-2.0/9, 2.0/9, 6.0/9, \
.0 , 2.0/9, 6.0/9, \
2.0/9, 2.0/9, 6.0/9, \
4.0/9, 2.0/9, 6.0/9, \
6.0/9, 2.0/9, 6.0/9, \
8.0/9, 2.0/9, 6.0/9, \
-8.0/9, 4.0/9, 6.0/9, \
-6.0/9, 4.0/9, 6.0/9, \
-4.0/9, 4.0/9, 6.0/9, \
-2.0/9, 4.0/9, 6.0/9, \
.0 , 4.0/9, 6.0/9, \
2.0/9, 4.0/9, 6.0/9, \
4.0/9, 4.0/9, 6.0/9, \
6.0/9, 4.0/9, 6.0/9, \
8.0/9, 4.0/9, 6.0/9, \
-8.0/9, 6.0/9, 6.0/9, \
-6.0/9, 6.0/9, 6.0/9, \
-4.0/9, 6.0/9, 6.0/9, \
-2.0/9, 6.0/9, 6.0/9, \
.0 , 6.0/9, 6.0/9, \
2.0/9, 6.0/9, 6.0/9, \
4.0/9, 6.0/9, 6.0/9, \
6.0/9, 6.0/9, 6.0/9, \
8.0/9, 6.0/9, 6.0/9, \
-8.0/9, 8.0/9, 6.0/9, \
-6.0/9, 8.0/9, 6.0/9, \
-4.0/9, 8.0/9, 6.0/9, \
-2.0/9, 8.0/9, 6.0/9, \
.0 , 8.0/9, 6.0/9, \
2.0/9, 8.0/9, 6.0/9, \
4.0/9, 8.0/9, 6.0/9, \
6.0/9, 8.0/9, 6.0/9, \
8.0/9, 8.0/9, 6.0/9, \
-8.0/9, -8.0/9, 8.0/9, \
-6.0/9, -8.0/9, 8.0/9, \
-4.0/9, -8.0/9, 8.0/9, \
-2.0/9, -8.0/9, 8.0/9, \
.0 , -8.0/9, 8.0/9, \
2.0/9, -8.0/9, 8.0/9, \
4.0/9, -8.0/9, 8.0/9, \
6.0/9, -8.0/9, 8.0/9, \
8.0/9, -8.0/9, 8.0/9, \
-8.0/9, -6.0/9, 8.0/9, \
-6.0/9, -6.0/9, 8.0/9, \
-4.0/9, -6.0/9, 8.0/9, \
-2.0/9, -6.0/9, 8.0/9, \
.0 , -6.0/9, 8.0/9, \
2.0/9, -6.0/9, 8.0/9, \
4.0/9, -6.0/9, 8.0/9, \
6.0/9, -6.0/9, 8.0/9, \
8.0/9, -6.0/9, 8.0/9, \
-8.0/9, -4.0/9, 8.0/9, \
-6.0/9, -4.0/9, 8.0/9, \
-4.0/9, -4.0/9, 8.0/9, \
-2.0/9, -4.0/9, 8.0/9, \
.0 , -4.0/9, 8.0/9, \
2.0/9, -4.0/9, 8.0/9, \
4.0/9, -4.0/9, 8.0/9, \
6.0/9, -4.0/9, 8.0/9, \
8.0/9, -4.0/9, 8.0/9, \
-8.0/9, -2.0/9, 8.0/9, \
-6.0/9, -2.0/9, 8.0/9, \
-4.0/9, -2.0/9, 8.0/9, \
-2.0/9, -2.0/9, 8.0/9, \
.0 , -2.0/9, 8.0/9, \
2.0/9, -2.0/9, 8.0/9, \
4.0/9, -2.0/9, 8.0/9, \
6.0/9, -2.0/9, 8.0/9, \
8.0/9, -2.0/9, 8.0/9, \
-8.0/9, .0 , 8.0/9, \
-6.0/9, .0 , 8.0/9, \
-4.0/9, .0 , 8.0/9, \
-2.0/9, .0 , 8.0/9, \
.0 , .0 , 8.0/9, \
2.0/9, .0 , 8.0/9, \
4.0/9, .0 , 8.0/9, \
6.0/9, .0 , 8.0/9, \
8.0/9, .0 , 8.0/9, \
-8.0/9, 2.0/9, 8.0/9, \
-6.0/9, 2.0/9, 8.0/9, \
-4.0/9, 2.0/9, 8.0/9, \
-2.0/9, 2.0/9, 8.0/9, \
.0 , 2.0/9, 8.0/9, \
2.0/9, 2.0/9, 8.0/9, \
4.0/9, 2.0/9, 8.0/9, \
6.0/9, 2.0/9, 8.0/9, \
8.0/9, 2.0/9, 8.0/9, \
-8.0/9, 4.0/9, 8.0/9, \
-6.0/9, 4.0/9, 8.0/9, \
-4.0/9, 4.0/9, 8.0/9, \
-2.0/9, 4.0/9, 8.0/9, \
.0 , 4.0/9, 8.0/9, \
2.0/9, 4.0/9, 8.0/9, \
4.0/9, 4.0/9, 8.0/9, \
6.0/9, 4.0/9, 8.0/9, \
8.0/9, 4.0/9, 8.0/9, \
-8.0/9, 6.0/9, 8.0/9, \
-6.0/9, 6.0/9, 8.0/9, \
-4.0/9, 6.0/9, 8.0/9, \
-2.0/9, 6.0/9, 8.0/9, \
.0 , 6.0/9, 8.0/9, \
2.0/9, 6.0/9, 8.0/9, \
4.0/9, 6.0/9, 8.0/9, \
6.0/9, 6.0/9, 8.0/9, \
8.0/9, 6.0/9, 8.0/9, \
-8.0/9, 8.0/9, 8.0/9, \
-6.0/9, 8.0/9, 8.0/9, \
-4.0/9, 8.0/9, 8.0/9, \
-2.0/9, 8.0/9, 8.0/9, \
.0 , 8.0/9, 8.0/9, \
2.0/9, 8.0/9, 8.0/9, \
4.0/9, 8.0/9, 8.0/9, \
6.0/9, 8.0/9, 8.0/9, \
8.0/9, 8.0/9, 8.0/9 \
}
/*
* Requantization tables : see ISO/IEC 11172-3 Annex B.2 and B.4
*
* = We store the requantization information in the following structures :
* typedef struct requantization_s
* {
* byte_t i_bits_per_codeword;
* const float * pf_ungroup;
* float f_slope;
* float f_offset;
* } requantization_t;
*
* = Theses values depend on the bitrate per channel, on the subband number
* and on the allocation[ch][sb] (see ISO/IEC 11172-3 2.4.2.6 and Annex B.2).
*
* = But, in order to avoid data redundancy, we use the following properties :
*
* - When bitrate_per_channel == 32 or 48 kbits/s (ie when
* bitrate_per_channel_index == 1 or 2), the requantization values depend
* only on allocation[ch][sb] (see ISO/IEC 11172-3 Annex B.2c and B.2d).
* That's why ADEC_LAYER2_REQUANTIZATION_CD = f(allocation).
*
* - In the other cases (see ISO/IEC 11172-3 Annex B.2a and B.2b), we can
* divide the tables in 4 subtables :
* + ADEC_LAYER2_REQUANTIZATION_AB1 for sb in [0..2]
* + ADEC_LAYER2_REQUANTIZATION_AB2 for sb in [3..10]
* + ADEC_LAYER2_REQUANTIZATION_AB3 for sb in [11..22]
* + ADEC_LAYER2_REQUANTIZATION_AB4 for sb in [23..29]
* That's why ADEC_LAYER2_REQUANTIZATION_AB = f(sb, allocation)
*
* = When these tables are used, pf_ungroup3, pf_ungroup5, pf_ungroup9,
* requantization_ab1, requantization_ab2, requantization_ab3 and
* requantization_ab4 must already have been defined.
*/
#define ADEC_LAYER2_REQUANTIZATION_CD \
{ \
{ 0, NULL, .0 , .0 }, /* allocation 0 */ \
{ 5, pf_ungroup3, .0 , .0 }, /* allocation 1 */ \
{ 7, pf_ungroup5, .0 , .0 }, /* allocation 2 */ \
{10, pf_ungroup9, .0 , .0 }, /* allocation 3 */ \
{ 4, NULL, .133333333332 , -.933333333328}, /* allocation 4 */ \
{ 5, NULL, .0645161290325 , -.967741935488}, /* allocation 5 */ \
{ 6, NULL, .0317460317459 , -.984126984124}, /* allocation 6 */ \
{ 7, NULL, .0157480314961 , -.992125984254}, /* allocation 7 */ \
{ 8, NULL, .00784313725492 , -.996078431375}, /* allocation 8 */ \
{ 9, NULL, .00391389432484 , -.998043052835}, /* allocation 9 */ \
{10, NULL, .00195503421311 , -.999022482897}, /* allocation 10 */ \
{11, NULL, .000977039570107 , -.99951148022 }, /* allocation 11 */ \
{12, NULL, .000488400488398 , -.999755799752}, /* allocation 12 */ \
{13, NULL, .000244170430962 , -.999877914789}, /* allocation 13 */ \
{14, NULL, .000122077763535 , -.999938961117}, /* allocation 14 */ \
{15, NULL, .000061037018952 , -.999969481491} /* allocation 15 */ \
}
#define ADEC_LAYER2_REQUANTIZATION_AB1 \
{ \
{ 0, NULL, .0 , .0 }, /* allocation 0 */ \
{ 5, pf_ungroup3, .0 , .0 }, /* allocation 1 */ \
{ 3, NULL, .285714285715 , -.857142857145}, /* allocation 2 */ \
{ 4, NULL, .133333333332 , -.933333333328}, /* allocation 3 */ \
{ 5, NULL, .0645161290325 , -.967741935488}, /* allocation 4 */ \
{ 6, NULL, .0317460317459 , -.984126984124}, /* allocation 5 */ \
{ 7, NULL, .0157480314961 , -.992125984254}, /* allocation 6 */ \
{ 8, NULL, .00784313725492 , -.996078431375}, /* allocation 7 */ \
{ 9, NULL, .00391389432484 , -.998043052835}, /* allocation 8 */ \
{10, NULL, .00195503421311 , -.999022482897}, /* allocation 9 */ \
{11, NULL, .000977039570107 , -.99951148022 }, /* allocation 10 */ \
{12, NULL, .000488400488398 , -.999755799752}, /* allocation 11 */ \
{13, NULL, .000244170430962 , -.999877914789}, /* allocation 12 */ \
{14, NULL, .000122077763535 , -.999938961117}, /* allocation 13 */ \
{15, NULL, .000061037018952 , -.999969481491}, /* allocation 14 */ \
{16, NULL, .0000305180437933, -.999984740976} /* allocation 15 */ \
}
#define ADEC_LAYER2_REQUANTIZATION_AB2 \
{ \
{ 0, NULL, .0 , .0 }, /* allocation 0 */ \
{ 5, pf_ungroup3, .0 , .0 }, /* allocation 1 */ \
{ 7, pf_ungroup5, .0 , .0 }, /* allocation 2 */ \
{ 3, NULL, .285714285715 , -.857142857145}, /* allocation 3 */ \
{10, pf_ungroup9, .0 , .0 }, /* allocation 4 */ \
{ 4, NULL, .133333333332 , -.933333333328}, /* allocation 5 */ \
{ 5, NULL, .0645161290325 , -.967741935488}, /* allocation 6 */ \
{ 6, NULL, .0317460317459 , -.984126984124}, /* allocation 7 */ \
{ 7, NULL, .0157480314961 , -.992125984254}, /* allocation 8 */ \
{ 8, NULL, .00784313725492 , -.996078431375}, /* allocation 9 */ \
{ 9, NULL, .00391389432484 , -.998043052835}, /* allocation 10 */ \
{10, NULL, .00195503421311 , -.999022482897}, /* allocation 11 */ \
{11, NULL, .000977039570107 , -.99951148022 }, /* allocation 12 */ \
{12, NULL, .000488400488398 , -.999755799752}, /* allocation 13 */ \
{13, NULL, .000244170430962 , -.999877914789}, /* allocation 14 */ \
{16, NULL, .0000305180437933, -.999984740976} /* allocation 15 */ \
}
#define ADEC_LAYER2_REQUANTIZATION_AB3 \
{ \
{ 0, NULL, .0 , .0 }, /* allocation 0 */ \
{ 5, pf_ungroup3, .0 , .0 }, /* allocation 1 */ \
{ 7, pf_ungroup5, .0 , .0 }, /* allocation 2 */ \
{ 3, NULL, .285714285715 , -.857142857145}, /* allocation 3 */ \
{10, pf_ungroup9, .0 , .0 }, /* allocation 4 */ \
{ 4, NULL, .133333333332 , -.933333333328}, /* allocation 5 */ \
{ 5, NULL, .0645161290325 , -.967741935488}, /* allocation 6 */ \
{16, NULL, .0000305180437933, -.999984740976}, /* allocation 7 */ \
{ 0, NULL, .0 , .0 }, /* allocation 8 */ \
{ 0, NULL, .0 , .0 }, /* allocation 9 */ \
{ 0, NULL, .0 , .0 }, /* allocation 10 */ \
{ 0, NULL, .0 , .0 }, /* allocation 11 */ \
{ 0, NULL, .0 , .0 }, /* allocation 12 */ \
{ 0, NULL, .0 , .0 }, /* allocation 13 */ \
{ 0, NULL, .0 , .0 }, /* allocation 14 */ \
{ 0, NULL, .0 , .0 } /* allocation 15 */ \
}
#define ADEC_LAYER2_REQUANTIZATION_AB4 \
{ \
{ 0, NULL, .0 , .0 }, /* allocation 0 */ \
{ 5, pf_ungroup3, .0 , .0 }, /* allocation 1 */ \
{ 7, pf_ungroup5, .0 , .0 }, /* allocation 2 */ \
{16, NULL, .0000305180437933, -.999984740976}, /* allocation 3 */ \
{ 0, NULL, .0 , .0 }, /* allocation 4 */ \
{ 0, NULL, .0 , .0 }, /* allocation 5 */ \
{ 0, NULL, .0 , .0 }, /* allocation 6 */ \
{ 0, NULL, .0 , .0 }, /* allocation 7 */ \
{ 0, NULL, .0 , .0 }, /* allocation 8 */ \
{ 0, NULL, .0 , .0 }, /* allocation 9 */ \
{ 0, NULL, .0 , .0 }, /* allocation 10 */ \
{ 0, NULL, .0 , .0 }, /* allocation 11 */ \
{ 0, NULL, .0 , .0 }, /* allocation 12 */ \
{ 0, NULL, .0 , .0 }, /* allocation 13 */ \
{ 0, NULL, .0 , .0 }, /* allocation 14 */ \
{ 0, NULL, .0 , .0 } /* allocation 15 */ \
}
#define ADEC_LAYER2_REQUANTIZATION_AB \
{ \
p_requantization_ab1, /* subband == 0 */ \
p_requantization_ab1, /* subband == 1 */ \
p_requantization_ab1, /* subband == 2 */ \
p_requantization_ab2, /* subband == 3 */ \
p_requantization_ab2, /* subband == 4 */ \
p_requantization_ab2, /* subband == 5 */ \
p_requantization_ab2, /* subband == 6 */ \
p_requantization_ab2, /* subband == 7 */ \
p_requantization_ab2, /* subband == 8 */ \
p_requantization_ab2, /* subband == 9 */ \
p_requantization_ab2, /* subband == 10 */ \
p_requantization_ab3, /* subband == 11 */ \
p_requantization_ab3, /* subband == 12 */ \
p_requantization_ab3, /* subband == 13 */ \
p_requantization_ab3, /* subband == 14 */ \
p_requantization_ab3, /* subband == 15 */ \
p_requantization_ab3, /* subband == 16 */ \
p_requantization_ab3, /* subband == 17 */ \
p_requantization_ab3, /* subband == 18 */ \
p_requantization_ab3, /* subband == 19 */ \
p_requantization_ab3, /* subband == 20 */ \
p_requantization_ab3, /* subband == 21 */ \
p_requantization_ab3, /* subband == 22 */ \
p_requantization_ab4, /* subband == 23 */ \
p_requantization_ab4, /* subband == 24 */ \
p_requantization_ab4, /* subband == 25 */ \
p_requantization_ab4, /* subband == 26 */ \
p_requantization_ab4, /* subband == 27 */ \
p_requantization_ab4, /* subband == 28 */ \
p_requantization_ab4 /* subband == 29 */ \
}
src/audio_decoder/audio_decoder.c
View file @ fec6ded4
......@@ -29,478 +29,855 @@
*/
#include "int_types.h"
#include "audio_constants.h"
//#include "audio_constants.h"
#include "audio_decoder.h"
#include "audio_math.h" /* DCT32(), PCM() */
#include "audio_bit_stream.h"
#define NULL ((void *)0)
/*****************************************************************************
* adec_Layer`L'_`M': decodes an mpeg 1, layer `L', mode `M', audio frame
*****************************************************************************
* These functions decode the audio frame which has already its header loaded
* in the i_header member of the audio decoder thread structure and its first
* byte of data described by the bit stream structure of the audio decoder
* thread (there is no bit available in the bit buffer yet)
*****************************************************************************/
/**** wkn ****/
/*****************************************************************************
* adec_Layer1_Mono
*****************************************************************************/
static __inline__ int adec_Layer1_Mono (audiodec_t * p_adec)
{
p_adec->bit_stream.buffer = 0;
p_adec->bit_stream.i_available = 0;
return (0);
}
static float adec_scalefactor_table[64] = { /* 2 ^ (1 - i/3) */
2.0000000000000000, 1.5874010519681994, 1.2599210498948732,
1.0000000000000000, 0.7937005259840998, 0.6299605249474366,
0.5000000000000000, 0.3968502629920499, 0.3149802624737183,
0.2500000000000000, 0.1984251314960249, 0.1574901312368591,
0.1250000000000000, 0.0992125657480125, 0.0787450656184296,
0.0625000000000000, 0.0496062828740062, 0.0393725328092148,
0.0312500000000000, 0.0248031414370031, 0.0196862664046074,
0.0156250000000000, 0.0124015707185016, 0.0098431332023037,
0.0078125000000000, 0.0062007853592508, 0.0049215666011518,
0.0039062500000000, 0.0031003926796254, 0.0024607833005759,
0.0019531250000000, 0.0015501963398127, 0.0012303916502880,
0.0009765625000000, 0.0007750981699063, 0.0006151958251440,
0.0004882812500000, 0.0003875490849532, 0.0003075979125720,
0.0002441406250000, 0.0001937745424766, 0.0001537989562860,
0.0001220703125000, 0.0000968872712383, 0.0000768994781430,
0.0000610351562500, 0.0000484436356191, 0.0000384497390715,
0.0000305175781250, 0.0000242218178096, 0.0000192248695357,
0.0000152587890625, 0.0000121109089048, 0.0000096124347679,
0.0000076293945312, 0.0000060554544524, 0.0000048062173839,
0.0000038146972656, 0.0000030277272262, 0.0000024031086920,
0.0000019073486328, 0.0000015138636131, 0.0000012015543460,
0.0000009536743164 /* last element is not in the standard... invalid ??? */
};
/*****************************************************************************
* adec_Layer1_Stereo
*****************************************************************************/
static __inline__ int adec_Layer1_Stereo (audiodec_t * p_adec)
{
p_adec->bit_stream.buffer = 0;
p_adec->bit_stream.i_available = 0;
return (0);
}
static float adec_slope_table[15] = {
0.6666666666666666, 0.2857142857142857, 0.1333333333333333,
0.0645161290322581, 0.0317460317460317, 0.0157480314960630,
0.0078431372549020, 0.0039138943248532, 0.0019550342130987,
0.0009770395701026, 0.0004884004884005, 0.0002441704309608,
0.0001220777635354, 0.0000610370189520, 0.0000305180437934
};
/*****************************************************************************
* adec_Layer2_Mono
*****************************************************************************/
static __inline__ int adec_Layer2_Mono (audiodec_t * p_adec)
{
p_adec->bit_stream.buffer = 0;
p_adec->bit_stream.i_available = 0;
return (0);
}
static float adec_offset_table[15] = {
-0.6666666666666666, -0.8571428571428571, -0.9333333333333333,
-0.9677419354838710, -0.9841269841269841, -0.9921259842519685,
-0.9960784313725490, -0.9980430528375733, -0.9990224828934506,
-0.9995114802149487, -0.9997557997557998, -0.9998779147845196,
-0.9999389611182323, -0.9999694814905240, -0.9999847409781033
};
/*****************************************************************************
* adec_Layer2_Stereo
*****************************************************************************/
static __inline__ int adec_Layer2_Stereo (audiodec_t * p_adec, s16 * buffer)
static u8 adec_layer1_allocation_table[15] = {
0, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15
};
static int adec_bound_table[4] = {4, 8, 12, 16};
static int adec_layer1_mono (audiodec_t * p_adec, s16 * buffer)
{
typedef struct requantization_s
{
u8 i_bits_per_codeword;
const float * pf_ungroup;
float f_slope;
float f_offset;
} requantization_t;
static const float pf_scalefactor[64] = ADEC_SCALE_FACTOR;
static u32 i_header;
static int i_sampling_frequency, i_mode, i_bound;
static int pi_allocation_0[32], pi_allocation_1[32]; /* see ISO/IEC 11172-3 2.4.1.6 */
int i_sb, i_nbal;
float f_scalefactor_0, f_scalefactor_1;
static const u8 ppi_bitrate_per_channel_index[4][15] = ADEC_LAYER2_BITRATE_PER_CHANNEL_INDEX;
static const u8 ppi_sblimit[3][11] = ADEC_LAYER2_SBLIMIT;
static const u8 ppi_nbal[2][32] = ADEC_LAYER2_NBAL;
static const float pf_ungroup3[3*3*3 * 3] = ADEC_LAYER2_UNGROUP3;
static const float pf_ungroup5[5*5*5 * 3] = ADEC_LAYER2_UNGROUP5;
static const float pf_ungroup9[9*9*9 * 3] = ADEC_LAYER2_UNGROUP9;
static const requantization_t p_requantization_cd[16] = ADEC_LAYER2_REQUANTIZATION_CD;
static const requantization_t p_requantization_ab1[16] = ADEC_LAYER2_REQUANTIZATION_AB1;
static const requantization_t p_requantization_ab2[16] = ADEC_LAYER2_REQUANTIZATION_AB2;
static const requantization_t p_requantization_ab3[16] = ADEC_LAYER2_REQUANTIZATION_AB3;
static const requantization_t p_requantization_ab4[16] = ADEC_LAYER2_REQUANTIZATION_AB4;
static const requantization_t * pp_requantization_ab[30] = ADEC_LAYER2_REQUANTIZATION_AB;
static int i_sblimit, i_bitrate_per_channel_index;
static int pi_scfsi_0[30], pi_scfsi_1[30];
static const u8 * pi_nbal;
static float ppf_sample_0[3][32], ppf_sample_1[3][32];
static const requantization_t * pp_requantization_0[30];
static const requantization_t * pp_requantization_1[30];
static requantization_t requantization;
static const float * pf_ungroup;
static float pf_scalefactor_0_0[30], pf_scalefactor_0_1[30], pf_scalefactor_0_2[30];
static float pf_scalefactor_1_0[30], pf_scalefactor_1_1[30], pf_scalefactor_1_2[30];
int i_2nbal, i_gr;
float f_dummy;
s16 * p_s16;
int i_need = 0, i_dump = 0;
#if 0
static const int pi_framesize[512] = ADEC_FRAME_SIZE;
#endif
u8 allocation[32];
float slope[32];
float offset[32];
float sample[32];
/* Read the audio frame header and flush the bit buffer */
i_header = p_adec->header;
/* Read the sampling frequency (see ISO/IEC 11172-3 2.4.2.3) */
i_sampling_frequency = (int)((i_header & ADEC_HEADER_SAMPLING_FREQUENCY_MASK)
>> ADEC_HEADER_SAMPLING_FREQUENCY_SHIFT);
/* Read the mode (see ISO/IEC 11172-3 2.4.2.3) */
i_mode = (int)((i_header & ADEC_HEADER_MODE_MASK) >> ADEC_HEADER_MODE_SHIFT);
/* If a CRC can be found in the frame, get rid of it */
if ((i_header & ADEC_HEADER_PROTECTION_BIT_MASK) == 0)
{
NeedBits (&p_adec->bit_stream, 16);
DumpBits (&p_adec->bit_stream, 16);
int sb;
int s;
/* parse allocation */
for (sb = 0; sb < 32; sb += 2) {
u8 tmp;
tmp = GetByte (&p_adec->bit_stream);
if ((tmp >> 4) > 14)
return 1;
allocation[sb] = adec_layer1_allocation_table [tmp >> 4];
if ((tmp & 15) > 14)
return 1;
allocation[sb+1] = adec_layer1_allocation_table [tmp & 15];
}
/* Find out the bitrate per channel index */
i_bitrate_per_channel_index = (int)ppi_bitrate_per_channel_index[i_mode]
[(i_header & ADEC_HEADER_BITRATE_INDEX_MASK) >> ADEC_HEADER_BITRATE_INDEX_SHIFT];
/* Find out the number of subbands */
i_sblimit = (int)ppi_sblimit[i_sampling_frequency][i_bitrate_per_channel_index];
/* Check if the frame is valid or not */
if (i_sblimit == 0)
{
return (0); /* the frame is invalid */
/* parse scalefactors */
for (sb = 0; sb < 32; sb++) {
if (allocation[sb]) {
int index;
float scalefactor;
NeedBits (&p_adec->bit_stream, 6);
index = p_adec->bit_stream.buffer >> (32 - 6);
DumpBits (&p_adec->bit_stream, 6);
scalefactor = adec_scalefactor_table[index];
slope[sb] = adec_slope_table[allocation[sb]-2] * scalefactor;
offset[sb] = adec_offset_table[allocation[sb]-2] * scalefactor;
}
}
/* Find out the number of bits allocated */
pi_nbal = ppi_nbal[ (i_bitrate_per_channel_index <= 2) ? 0 : 1 ];
/* Find out the `bound' subband (see ISO/IEC 11172-3 2.4.2.3) */
if (i_mode == 1)
{
i_bound = (int)(((i_header & ADEC_HEADER_MODE_EXTENSION_MASK) >> (ADEC_HEADER_MODE_EXTENSION_SHIFT - 2)) + 4);
if (i_bound > i_sblimit)
{
i_bound = i_sblimit;
}
/* parse samples */
for (s = 0; s < 12; s++) {
s16 * XXX_buf;
for (sb = 0; sb < 32; sb++) {
if (!allocation[sb]) {
sample[sb] = 0;
} else {
int code;
NeedBits (&p_adec->bit_stream, allocation[sb]);
code = p_adec->bit_stream.buffer >> (32 - allocation[sb]);
DumpBits (&p_adec->bit_stream, allocation[sb]);
sample[sb] = slope[sb] * code + offset[sb];
}
}
DCT32 (sample, &p_adec->bank_0);
XXX_buf = buffer;
PCM (&p_adec->bank_0, &XXX_buf, 1);
buffer += 32;
}
else
{
i_bound = i_sblimit;
return 0;
}
static int adec_layer1_stereo (audiodec_t * p_adec, s16 * buffer)
{
u8 allocation_0[32], allocation_1[32];
float slope_0[32], slope_1[32];
float offset_0[32], offset_1[32];
float sample_0[32], sample_1[32];
int bound;
int sb;
int s;
/* calculate bound */
bound = 32;
if ((p_adec->header & 0xc0) == 0x40) { /* intensity stereo */
int index;
index = (p_adec->header >> 4) & 3;
bound = adec_bound_table[index];
}
/* Read the allocation information (see ISO/IEC 11172-3 2.4.1.6) */
for (i_sb = 0; i_sb < i_bound; i_sb++)
{
i_2nbal = 2 * (i_nbal = (int)pi_nbal[ i_sb ]);
NeedBits (&p_adec->bit_stream, i_2nbal);
i_need += i_2nbal;
pi_allocation_0[ i_sb ] = (int)(p_adec->bit_stream.buffer >> (32 - i_nbal));
p_adec->bit_stream.buffer <<= i_nbal;
pi_allocation_1[ i_sb ] = (int)(p_adec->bit_stream.buffer >> (32 - i_nbal));
p_adec->bit_stream.buffer <<= i_nbal;
p_adec->bit_stream.i_available -= i_2nbal;
i_dump += i_2nbal;
/* parse allocation */
for (sb = 0; sb < bound; sb++) {
u8 tmp;
tmp = GetByte (&p_adec->bit_stream);
if ((tmp >> 4) > 14)
return 1;
allocation_0[sb] = adec_layer1_allocation_table [tmp >> 4];
if ((tmp & 15) > 14)
return 1;
allocation_1[sb] = adec_layer1_allocation_table [tmp & 15];
}
for (; i_sb < i_sblimit; i_sb++)
{
i_nbal = (int)pi_nbal[ i_sb ];
NeedBits (&p_adec->bit_stream, i_nbal);
i_need += i_nbal;
pi_allocation_0[ i_sb ] = (int)(p_adec->bit_stream.buffer >> (32 - i_nbal));
DumpBits (&p_adec->bit_stream, i_nbal);
i_dump += i_nbal;
for (; sb < 32; sb += 2) {
u8 tmp;
tmp = GetByte (&p_adec->bit_stream);
if ((tmp >> 4) > 14)
return 1;
allocation_0[sb] = allocation_1[sb] = adec_layer1_allocation_table [tmp >> 4];
if ((tmp & 15) > 14)
return 1;
allocation_0[sb+1] = allocation_1[sb+1] = adec_layer1_allocation_table [tmp & 15];
}
#define MACRO(p_requantization) \
for (i_sb = 0; i_sb < i_bound; i_sb++) \
{ \
if (pi_allocation_0[i_sb]) \
{ \
pp_requantization_0[i_sb] = &((p_requantization)[pi_allocation_0[i_sb]]); \
NeedBits (&p_adec->bit_stream, 2); \
i_need += 2; \
pi_scfsi_0[i_sb] = (int)(p_adec->bit_stream.buffer >> (32 - 2)); \
DumpBits (&p_adec->bit_stream, 2); \
i_dump += 2; \
} \
else \
{ \
ppf_sample_0[0][i_sb] = .0; \
ppf_sample_0[1][i_sb] = .0; \
ppf_sample_0[2][i_sb] = .0; \
} \
\
if (pi_allocation_1[i_sb]) \
{ \
pp_requantization_1[i_sb] = &((p_requantization)[pi_allocation_1[i_sb]]); \
NeedBits (&p_adec->bit_stream, 2); \
i_need += 2; \
pi_scfsi_1[i_sb] = (int)(p_adec->bit_stream.buffer >> (32 - 2)); \
DumpBits (&p_adec->bit_stream, 2); \
i_dump += 2; \
} \
else \
{ \
ppf_sample_1[0][i_sb] = .0; \
ppf_sample_1[1][i_sb] = .0; \
ppf_sample_1[2][i_sb] = .0; \
} \
} \
\
for (; i_sb < i_sblimit; i_sb++) \
{ \
if (pi_allocation_0[i_sb]) \
{ \
pp_requantization_0[i_sb] = &((p_requantization)[pi_allocation_0[i_sb]]); \
NeedBits (&p_adec->bit_stream, 4); \
i_need += 4; \
pi_scfsi_0[i_sb] = (int)(p_adec->bit_stream.buffer >> (32 - 2)); \
p_adec->bit_stream.buffer <<= 2; \
pi_scfsi_1[i_sb] = (int)(p_adec->bit_stream.buffer >> (32 - 2)); \
p_adec->bit_stream.buffer <<= 2; \
p_adec->bit_stream.i_available -= 4; \
i_dump += 4; \
} \
else \
{ \
ppf_sample_0[0][i_sb] = .0; \
ppf_sample_0[1][i_sb] = .0; \
ppf_sample_0[2][i_sb] = .0; \
ppf_sample_1[0][i_sb] = .0; \
ppf_sample_1[1][i_sb] = .0; \
ppf_sample_1[2][i_sb] = .0; \
} \
}
/* #define MACRO */
/* parse scalefactors */
for (sb = 0; sb < 32; sb++) {
if (allocation_0[sb]) {
int index;
float scalefactor;
NeedBits (&p_adec->bit_stream, 6);
index = p_adec->bit_stream.buffer >> (32 - 6);
DumpBits (&p_adec->bit_stream, 6);
if (i_bitrate_per_channel_index <= 2)
{
MACRO (p_requantization_cd)
scalefactor = adec_scalefactor_table[index];
slope_0[sb] = adec_slope_table[allocation_0[sb]-2] * scalefactor;
offset_0[sb] = adec_offset_table[allocation_0[sb]-2] * scalefactor;
}
if (allocation_1[sb]) {
int index;
float scalefactor;
NeedBits (&p_adec->bit_stream, 6);
index = p_adec->bit_stream.buffer >> (32 - 6);
DumpBits (&p_adec->bit_stream, 6);
scalefactor = adec_scalefactor_table[index];
slope_1[sb] = adec_slope_table[allocation_1[sb]-2] * scalefactor;
offset_1[sb] = adec_offset_table[allocation_1[sb]-2] * scalefactor;
}
}
else
{
MACRO (pp_requantization_ab[i_sb])
/* parse samples */
for (s = 0; s < 12; s++) {
s16 * XXX_buf;
for (sb = 0; sb < bound; sb++) {
if (!allocation_0[sb])
sample_0[sb] = 0;
else {
int code;
NeedBits (&p_adec->bit_stream, allocation_0[sb]);
code = p_adec->bit_stream.buffer >> (32 - allocation_0[sb]);
DumpBits (&p_adec->bit_stream, allocation_0[sb]);
sample_0[sb] = slope_0[sb] * code + offset_0[sb];
}
if (!allocation_1[sb])
sample_1[sb] = 0;
else {
int code;
NeedBits (&p_adec->bit_stream, allocation_1[sb]);
code = p_adec->bit_stream.buffer >> (32 - allocation_1[sb]);
DumpBits (&p_adec->bit_stream, allocation_1[sb]);
sample_1[sb] = slope_1[sb] * code + offset_1[sb];
}
}
for (; sb < 32; sb++) {
if (!allocation_0[sb]) {
sample_0[sb] = 0;
sample_1[sb] = 0;
} else {
int sample;
NeedBits (&p_adec->bit_stream, allocation_0[sb]);
sample = p_adec->bit_stream.buffer >> (32 - allocation_0[sb]);
DumpBits (&p_adec->bit_stream, allocation_0[sb]);
sample_0[sb] = slope_0[sb] * sample + offset_0[sb];
sample_1[sb] = slope_1[sb] * sample + offset_1[sb];
}
}
DCT32 (sample_0, &p_adec->bank_0);
XXX_buf = buffer;
PCM (&p_adec->bank_0, &XXX_buf, 2);
DCT32 (sample_1, &p_adec->bank_1);
XXX_buf = buffer+1;
PCM (&p_adec->bank_1, &XXX_buf, 2);
buffer += 64;
}
#define SWITCH(pi_scfsi,pf_scalefactor_0,pf_scalefactor_1,pf_scalefactor_2) \
switch ((pi_scfsi)[i_sb]) \
{ \
case 0: \
NeedBits (&p_adec->bit_stream, (3*6)); \
i_need += 18; \
(pf_scalefactor_0)[i_sb] = pf_scalefactor[p_adec->bit_stream.buffer >> (32 - 6)]; \
p_adec->bit_stream.buffer <<= 6; \
(pf_scalefactor_1)[i_sb] = pf_scalefactor[p_adec->bit_stream.buffer >> (32 - 6)]; \
p_adec->bit_stream.buffer <<= 6; \
(pf_scalefactor_2)[i_sb] = pf_scalefactor[p_adec->bit_stream.buffer >> (32 - 6)]; \
p_adec->bit_stream.buffer <<= 6; \
p_adec->bit_stream.i_available -= (3*6); \
i_dump += 18; \
break; \
\
case 1: \
NeedBits (&p_adec->bit_stream, (2*6)); \
i_need += 12; \
(pf_scalefactor_0)[i_sb] = \
(pf_scalefactor_1)[i_sb] = pf_scalefactor[p_adec->bit_stream.buffer >> (32 - 6)]; \
p_adec->bit_stream.buffer <<= 6; \
(pf_scalefactor_2)[i_sb] = pf_scalefactor[p_adec->bit_stream.buffer >> (32 - 6)]; \
p_adec->bit_stream.buffer <<= 6; \
p_adec->bit_stream.i_available -= (2*6); \
i_dump += 12; \
break; \
\
case 2: \
NeedBits (&p_adec->bit_stream, (1*6)); \
i_need += 6; \
(pf_scalefactor_0)[i_sb] = \
(pf_scalefactor_1)[i_sb] = \
(pf_scalefactor_2)[i_sb] = pf_scalefactor[p_adec->bit_stream.buffer >> (32 - 6)]; \
DumpBits (&p_adec->bit_stream, (1*6)); \
i_dump += 6; \
break; \
\
case 3: \
NeedBits (&p_adec->bit_stream, (2*6)); \
i_need += 12; \
(pf_scalefactor_0)[i_sb] = pf_scalefactor[p_adec->bit_stream.buffer >> (32 - 6)]; \
p_adec->bit_stream.buffer <<= 6; \
(pf_scalefactor_1)[i_sb] = \
(pf_scalefactor_2)[i_sb] = pf_scalefactor[p_adec->bit_stream.buffer >> (32 - 6)]; \
p_adec->bit_stream.buffer <<= 6; \
p_adec->bit_stream.i_available -= (2*6); \
i_dump += 12; \
break; \
return 0;
}
typedef struct {
s8 nbal[32];
u8 * alloc[32];
} alloc_table_t;
#define L3 -1
#define L5 -2
#define L9 -3
static void adec_layer2_get_table (u32 header, u8 freq_table[15],
alloc_table_t ** alloc, int * sblimit)
{
static s8 table_ab0[16] = {
0, L3, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16
};
static s8 table_ab3[16] = {
0, L3, L5, 3, L9, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 16
};
static s8 table_ab11[8] = {
0, L3, L5, 3, L9, 4, 5, 16
};
static s8 table_ab23[8] = {
0, L3, L5, 16
};
static alloc_table_t mpeg1_ab = {
{4,4,4,4,4,4,4,4,4,4,4,3,3,3,3,3,3,3,3,3,3,3,3,2,2,2,2,2,2,2,0,0},
{table_ab0, table_ab0, table_ab0, table_ab3,
table_ab3, table_ab3, table_ab3, table_ab3,
table_ab3, table_ab3, table_ab3, table_ab11,
table_ab11, table_ab11, table_ab11, table_ab11,
table_ab11, table_ab11, table_ab11, table_ab11,
table_ab11, table_ab11, table_ab11, table_ab23,
table_ab23, table_ab23, table_ab23, table_ab23,
table_ab23, table_ab23, NULL, NULL}
};
static s8 table_cd[16] = {
0, L3, L5, L9, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15
};
static alloc_table_t mpeg1_cd = {
{4,4,3,3,3,3,3,3,3,3,3,3,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0},
{table_cd, table_cd, table_cd, table_cd,
table_cd, table_cd, table_cd, table_cd,
table_cd, table_cd, table_cd, table_cd,
NULL, NULL, NULL, NULL,
NULL, NULL, NULL, NULL, NULL, NULL, NULL, NULL,
NULL, NULL, NULL, NULL, NULL, NULL, NULL, NULL}
};
static s8 table_0[16] = {
0, L3, L5, 3, L9, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14
};
static s8 table_4[8] = {
0, L3, L5, L9, 4, 5, 6, 7
};
static alloc_table_t mpeg2 = {
{4,4,4,4,3,3,3,3,3,3,3,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,0,0},
{table_0, table_0, table_0, table_0,
table_4, table_4, table_4, table_4,
table_4, table_4, table_4, table_4,
table_4, table_4, table_4, table_4,
table_4, table_4, table_4, table_4,
table_4, table_4, table_4, table_4,
table_4, table_4, table_4, table_4,
table_4, table_4, NULL, NULL}
};
static alloc_table_t * alloc_table [4] = {
&mpeg2, &mpeg1_cd, &mpeg1_ab, &mpeg1_ab
};
static int sblimit_table[12] = {
30, 8, 27, 30, 30, 8, 27, 27, 30, 12, 27, 30
};
int index;
if (!(header & 0x80000))
index = 0; /* mpeg2 */
else {
index = (header >> 12) & 15; /* mpeg1, bitrate */
index = freq_table [index];
}
/* #define SWITCH */
for (i_sb = 0; i_sb < i_bound; i_sb++)
{
if (pi_allocation_0[i_sb])
{
SWITCH (pi_scfsi_0, pf_scalefactor_0_0, pf_scalefactor_0_1, pf_scalefactor_0_2)
}
if (pi_allocation_1[i_sb])
{
SWITCH (pi_scfsi_1, pf_scalefactor_1_0, pf_scalefactor_1_1, pf_scalefactor_1_2)
}
*alloc = alloc_table[index];
index |= (header >> 8) & 12;
*sblimit = sblimit_table[index];
}
static int adec_layer2_mono (audiodec_t * p_adec, s16 * buffer)
{
return 1;
}
static int adec_layer2_stereo (audiodec_t * p_adec, s16 * buffer)
{
static u8 freq_table[15] = {3, 0, 0, 0, 1, 0, 1, 2, 2, 2, 3, 3, 3, 3, 3};
static float L3_table[3] = {-2/3.0, 0, 2/3.0};
static float L5_table[5] = {-4/5.0, -2/5.0, 0, 2/5.0, 4/5.0};
static float L9_table[9] = {-8/9.0, -6/9.0, -4/9.0, -2/9.0, 0,
2/9.0, 4/9.0, 6/9.0, 8/9.0};
s8 allocation_0[32], allocation_1[32];
u8 scfsi_0[32], scfsi_1[32];
float slope_0[3][32], slope_1[3][32];
float offset_0[3][32], offset_1[3][32];
float sample_0[3][32], sample_1[3][32];
alloc_table_t * alloc_table;
int sblimit;
int bound;
int sb;
int gr0, gr1;
int s;
/* get the right allocation table */
adec_layer2_get_table (p_adec->header, freq_table, &alloc_table, &sblimit);
/* calculate bound */
bound = sblimit;
if ((p_adec->header & 0xc0) == 0x40) { /* intensity stereo */
int index;
index = (p_adec->header >> 4) & 3;
if (adec_bound_table[index] < sblimit)
bound = adec_bound_table[index];
}
for (; i_sb < i_sblimit; i_sb++)
{
if (pi_allocation_0[i_sb])
{
SWITCH (pi_scfsi_0, pf_scalefactor_0_0, pf_scalefactor_0_1, pf_scalefactor_0_2)
SWITCH (pi_scfsi_1, pf_scalefactor_1_0, pf_scalefactor_1_1, pf_scalefactor_1_2)
}
/* parse allocation */
for (sb = 0; sb < bound; sb++) {
int index;
NeedBits (&p_adec->bit_stream, alloc_table->nbal[sb]);
index = p_adec->bit_stream.buffer >> (32 - alloc_table->nbal[sb]);
DumpBits (&p_adec->bit_stream, alloc_table->nbal[sb]);
allocation_0[sb] = alloc_table->alloc[sb][index];
NeedBits (&p_adec->bit_stream, alloc_table->nbal[sb]);
index = p_adec->bit_stream.buffer >> (32 - alloc_table->nbal[sb]);
DumpBits (&p_adec->bit_stream, alloc_table->nbal[sb]);
allocation_1[sb] = alloc_table->alloc[sb][index];
}
for (; i_sb < 32; i_sb++)
{
ppf_sample_0[0][i_sb] = .0;
ppf_sample_0[1][i_sb] = .0;
ppf_sample_0[2][i_sb] = .0;
ppf_sample_1[0][i_sb] = .0;
ppf_sample_1[1][i_sb] = .0;
ppf_sample_1[2][i_sb] = .0;
for (; sb < sblimit; sb++) {
int index;
NeedBits (&p_adec->bit_stream, alloc_table->nbal[sb]);
index = p_adec->bit_stream.buffer >> (32 - alloc_table->nbal[sb]);
DumpBits (&p_adec->bit_stream, alloc_table->nbal[sb]);
allocation_0[sb] = allocation_1[sb] = alloc_table->alloc[sb][index];
}
#define GROUPTEST(pp_requantization,ppf_sample,pf_sf) \
requantization = *((pp_requantization)[i_sb]); \
if (requantization.pf_ungroup == NULL) \
{ \
NeedBits (&p_adec->bit_stream, requantization.i_bits_per_codeword); \
i_need += requantization.i_bits_per_codeword; \
(ppf_sample)[0][i_sb] = (f_scalefactor_0 = (pf_sf)[i_sb]) * (requantization.f_slope * \
(p_adec->bit_stream.buffer >> (32 - requantization.i_bits_per_codeword)) + requantization.f_offset); \
DumpBits (&p_adec->bit_stream, requantization.i_bits_per_codeword); \
i_dump += requantization.i_bits_per_codeword; \
\
NeedBits (&p_adec->bit_stream, requantization.i_bits_per_codeword); \
i_need += requantization.i_bits_per_codeword; \
(ppf_sample)[1][i_sb] = f_scalefactor_0 * (requantization.f_slope * \
(p_adec->bit_stream.buffer >> (32 - requantization.i_bits_per_codeword)) + requantization.f_offset); \
DumpBits (&p_adec->bit_stream, requantization.i_bits_per_codeword); \
i_dump += requantization.i_bits_per_codeword; \
\
NeedBits (&p_adec->bit_stream, requantization.i_bits_per_codeword); \
i_need += requantization.i_bits_per_codeword; \
(ppf_sample)[2][i_sb] = f_scalefactor_0 * (requantization.f_slope * \
(p_adec->bit_stream.buffer >> (32 - requantization.i_bits_per_codeword)) + requantization.f_offset); \
DumpBits (&p_adec->bit_stream, requantization.i_bits_per_codeword); \
i_dump += requantization.i_bits_per_codeword; \
} \
else \
{ \
NeedBits (&p_adec->bit_stream, requantization.i_bits_per_codeword); \
i_need += requantization.i_bits_per_codeword; \
pf_ungroup = requantization.pf_ungroup + 3 * \
(p_adec->bit_stream.buffer >> (32 - requantization.i_bits_per_codeword)); \
DumpBits (&p_adec->bit_stream, requantization.i_bits_per_codeword); \
i_dump += requantization.i_bits_per_codeword; \
(ppf_sample)[0][i_sb] = (f_scalefactor_0 = (pf_sf)[i_sb]) * pf_ungroup[0]; \
(ppf_sample)[1][i_sb] = f_scalefactor_0 * pf_ungroup[1]; \
(ppf_sample)[2][i_sb] = f_scalefactor_0 * pf_ungroup[2]; \
/* parse scfsi */
for (sb = 0; sb < sblimit; sb++) {
if (allocation_0[sb]) {
NeedBits (&p_adec->bit_stream, 2);
scfsi_0[sb] = p_adec->bit_stream.buffer >> (32 - 2);
DumpBits (&p_adec->bit_stream, 2);
}
if (allocation_1[sb]) {
NeedBits (&p_adec->bit_stream, 2);
scfsi_1[sb] = p_adec->bit_stream.buffer >> (32 - 2);
DumpBits (&p_adec->bit_stream, 2);
}
}
/* #define GROUPTEST */
#define READ_SAMPLE_L2S(pf_scalefactor_0,pf_scalefactor_1,i_grlimit) \
for (; i_gr < (i_grlimit); i_gr++) \
{ \
for (i_sb = 0; i_sb < i_bound; i_sb++) \
{ \
if (pi_allocation_0[i_sb]) \
{ \
GROUPTEST (pp_requantization_0, ppf_sample_0, (pf_scalefactor_0)) \
} \
if (pi_allocation_1[i_sb]) \
{ \
GROUPTEST (pp_requantization_1, ppf_sample_1, (pf_scalefactor_1)) \
} \
} \
for (; i_sb < i_sblimit; i_sb++) \
{ \
if (pi_allocation_0[i_sb]) \
{ \
requantization = *(pp_requantization_0[i_sb]); \
if (requantization.pf_ungroup == NULL) \
{ \
NeedBits (&p_adec->bit_stream, requantization.i_bits_per_codeword); \
i_need += requantization.i_bits_per_codeword; \
ppf_sample_0[0][i_sb] = (f_scalefactor_0 = (pf_scalefactor_0)[i_sb]) * \
(requantization.f_slope * (f_dummy = \
(float)(p_adec->bit_stream.buffer >> (32 - requantization.i_bits_per_codeword))) + \
requantization.f_offset); \
DumpBits (&p_adec->bit_stream, requantization.i_bits_per_codeword); \
i_dump += requantization.i_bits_per_codeword; \
ppf_sample_1[0][i_sb] = (f_scalefactor_1 = (pf_scalefactor_1)[i_sb]) * \
(requantization.f_slope * f_dummy + requantization.f_offset); \
\
NeedBits (&p_adec->bit_stream, requantization.i_bits_per_codeword); \
i_need += requantization.i_bits_per_codeword; \
ppf_sample_0[1][i_sb] = f_scalefactor_0 * \
(requantization.f_slope * (f_dummy = \
(float)(p_adec->bit_stream.buffer >> (32 - requantization.i_bits_per_codeword))) + \
requantization.f_offset); \
DumpBits (&p_adec->bit_stream, requantization.i_bits_per_codeword); \
i_dump += requantization.i_bits_per_codeword; \
ppf_sample_1[1][i_sb] = f_scalefactor_1 * \
(requantization.f_slope * f_dummy + requantization.f_offset); \
\
NeedBits (&p_adec->bit_stream, requantization.i_bits_per_codeword); \
i_need += requantization.i_bits_per_codeword; \
ppf_sample_0[2][i_sb] = f_scalefactor_0 * \
(requantization.f_slope * (f_dummy = \
(float)(p_adec->bit_stream.buffer >> (32 - requantization.i_bits_per_codeword))) + \
requantization.f_offset); \
DumpBits (&p_adec->bit_stream, requantization.i_bits_per_codeword); \
i_dump += requantization.i_bits_per_codeword; \
ppf_sample_1[2][i_sb] = f_scalefactor_1 * \
(requantization.f_slope * f_dummy + requantization.f_offset); \
} \
else \
{ \
NeedBits (&p_adec->bit_stream, requantization.i_bits_per_codeword); \
i_need += requantization.i_bits_per_codeword; \
pf_ungroup = requantization.pf_ungroup + 3 * \
(p_adec->bit_stream.buffer >> (32 - requantization.i_bits_per_codeword)); \
DumpBits (&p_adec->bit_stream, requantization.i_bits_per_codeword); \
i_dump += requantization.i_bits_per_codeword; \
\
ppf_sample_0[0][i_sb] = (f_scalefactor_0 = (pf_scalefactor_0)[i_sb]) * pf_ungroup[0]; \
ppf_sample_0[1][i_sb] = f_scalefactor_0 * pf_ungroup[1]; \
ppf_sample_0[2][i_sb] = f_scalefactor_0 * pf_ungroup[2]; \
\
ppf_sample_1[0][i_sb] = (f_scalefactor_1 = (pf_scalefactor_1)[i_sb]) * pf_ungroup[0]; \
ppf_sample_1[1][i_sb] = f_scalefactor_1 * pf_ungroup[1]; \
ppf_sample_1[2][i_sb] = f_scalefactor_1 * pf_ungroup[2]; \
} \
} \
} \
\
DCT32 (ppf_sample_0[0], &p_adec->bank_0); \
PCM (&p_adec->bank_0, &p_s16, 2); \
p_s16 -= 63; \
\
DCT32 (ppf_sample_1[0], &p_adec->bank_1); \
PCM (&p_adec->bank_1, &p_s16, 2); \
p_s16 -= 1; \
\
DCT32 (ppf_sample_0[1], &p_adec->bank_0); \
PCM (&p_adec->bank_0, &p_s16, 2); \
p_s16 -= 63; \
\
DCT32 (ppf_sample_1[1], &p_adec->bank_1); \
PCM (&p_adec->bank_1, &p_s16, 2); \
p_s16 -= 1; \
\
DCT32 (ppf_sample_0[2], &p_adec->bank_0); \
PCM (&p_adec->bank_0, &p_s16, 2); \
p_s16 -= 63; \
\
DCT32 (ppf_sample_1[2], &p_adec->bank_1); \
PCM (&p_adec->bank_1, &p_s16, 2); \
p_s16 -= 1; \
/* parse scalefactors */
for (sb = 0; sb < sblimit; sb++) {
if (allocation_0[sb]) {
int index_0, index_1, index_2;
switch (scfsi_0[sb]) {
case 0:
NeedBits (&p_adec->bit_stream, 18);
index_0 = p_adec->bit_stream.buffer >> (32 - 6);
index_1 = (p_adec->bit_stream.buffer >> (32 - 12)) & 63;
index_2 = (p_adec->bit_stream.buffer >> (32 - 18)) & 63;
DumpBits (&p_adec->bit_stream, 18);
if (allocation_0[sb] < 0) {
slope_0[0][sb] = adec_scalefactor_table[index_0];
slope_0[1][sb] = adec_scalefactor_table[index_1];
slope_0[2][sb] = adec_scalefactor_table[index_2];
} else {
float scalefactor;
float slope, offset;
slope = adec_slope_table[allocation_0[sb]-2];
offset = adec_offset_table[allocation_0[sb]-2];
scalefactor = adec_scalefactor_table[index_0];
slope_0[0][sb] = slope * scalefactor;
offset_0[0][sb] = offset * scalefactor;
scalefactor = adec_scalefactor_table[index_1];
slope_0[1][sb] = slope * scalefactor;
offset_0[1][sb] = offset * scalefactor;
scalefactor = adec_scalefactor_table[index_2];
slope_0[2][sb] = slope * scalefactor;
offset_0[2][sb] = offset * scalefactor;
}
break;
case 1:
NeedBits (&p_adec->bit_stream, 12);
index_0 = p_adec->bit_stream.buffer >> (32 - 6);
index_1 = (p_adec->bit_stream.buffer >> (32 - 12)) & 63;
DumpBits (&p_adec->bit_stream, 12);
if (allocation_0[sb] < 0) {
slope_0[0][sb] = slope_0[1][sb] =
adec_scalefactor_table[index_0];
slope_0[2][sb] = adec_scalefactor_table[index_1];
} else {
float scalefactor;
float slope, offset;
slope = adec_slope_table[allocation_0[sb]-2];
offset = adec_offset_table[allocation_0[sb]-2];
scalefactor = adec_scalefactor_table[index_0];
slope_0[0][sb] = slope_0[1][sb] = slope * scalefactor;
offset_0[0][sb] = offset_0[1][sb] = offset * scalefactor;
scalefactor = adec_scalefactor_table[index_1];
slope_0[2][sb] = slope * scalefactor;
offset_0[2][sb] = offset * scalefactor;
}
break;
case 2:
NeedBits (&p_adec->bit_stream, 6);
index_0 = p_adec->bit_stream.buffer >> (32 - 6);
DumpBits (&p_adec->bit_stream, 6);
if (allocation_0[sb] < 0) {
slope_0[0][sb] = slope_0[1][sb] = slope_0[2][sb] =
adec_scalefactor_table[index_0];
} else {
float scalefactor;
float slope, offset;
slope = adec_slope_table[allocation_0[sb]-2];
offset = adec_offset_table[allocation_0[sb]-2];
scalefactor = adec_scalefactor_table[index_0];
slope_0[0][sb] = slope_0[1][sb] = slope_0[2][sb] =
slope * scalefactor;
offset_0[0][sb] = offset_0[1][sb] = offset_0[2][sb] =
offset * scalefactor;
}
break;
case 3:
NeedBits (&p_adec->bit_stream, 12);
index_0 = p_adec->bit_stream.buffer >> (32 - 6);
index_1 = (p_adec->bit_stream.buffer >> (32 - 12)) & 63;
DumpBits (&p_adec->bit_stream, 12);
if (allocation_0[sb] < 0) {
slope_0[0][sb] = adec_scalefactor_table[index_0];
slope_0[1][sb] = slope_0[2][sb] =
adec_scalefactor_table[index_1];
} else {
float scalefactor;
float slope, offset;
slope = adec_slope_table[allocation_0[sb]-2];
offset = adec_offset_table[allocation_0[sb]-2];
scalefactor = adec_scalefactor_table[index_0];
slope_0[0][sb] = slope * scalefactor;
offset_0[0][sb] = offset * scalefactor;
scalefactor = adec_scalefactor_table[index_1];
slope_0[1][sb] = slope_0[2][sb] = slope * scalefactor;
offset_0[1][sb] = offset_0[2][sb] = offset * scalefactor;
}
break;
}
}
if (allocation_1[sb]) {
int index_0, index_1, index_2;
switch (scfsi_1[sb]) {
case 0:
NeedBits (&p_adec->bit_stream, 18);
index_0 = p_adec->bit_stream.buffer >> (32 - 6);
index_1 = (p_adec->bit_stream.buffer >> (32 - 12)) & 63;
index_2 = (p_adec->bit_stream.buffer >> (32 - 18)) & 63;
DumpBits (&p_adec->bit_stream, 18);
if (allocation_1[sb] < 0) {
slope_1[0][sb] = adec_scalefactor_table[index_0];
slope_1[1][sb] = adec_scalefactor_table[index_1];
slope_1[2][sb] = adec_scalefactor_table[index_2];
} else {
float scalefactor;
float slope, offset;
slope = adec_slope_table[allocation_1[sb]-2];
offset = adec_offset_table[allocation_1[sb]-2];
scalefactor = adec_scalefactor_table[index_0];
slope_1[0][sb] = slope * scalefactor;
offset_1[0][sb] = offset * scalefactor;
scalefactor = adec_scalefactor_table[index_1];
slope_1[1][sb] = slope * scalefactor;
offset_1[1][sb] = offset * scalefactor;
scalefactor = adec_scalefactor_table[index_2];
slope_1[2][sb] = slope * scalefactor;
offset_1[2][sb] = offset * scalefactor;
}
break;
case 1:
NeedBits (&p_adec->bit_stream, 12);
index_0 = p_adec->bit_stream.buffer >> (32 - 6);
index_1 = (p_adec->bit_stream.buffer >> (32 - 12)) & 63;
DumpBits (&p_adec->bit_stream, 12);
if (allocation_1[sb] < 0) {
slope_1[0][sb] = slope_1[1][sb] =
adec_scalefactor_table[index_0];
slope_1[2][sb] = adec_scalefactor_table[index_1];
} else {
float scalefactor;
float slope, offset;
slope = adec_slope_table[allocation_1[sb]-2];
offset = adec_offset_table[allocation_1[sb]-2];
scalefactor = adec_scalefactor_table[index_0];
slope_1[0][sb] = slope_1[1][sb] = slope * scalefactor;
offset_1[0][sb] = offset_1[1][sb] = offset * scalefactor;
scalefactor = adec_scalefactor_table[index_1];
slope_1[2][sb] = slope * scalefactor;
offset_1[2][sb] = offset * scalefactor;
}
break;
case 2:
NeedBits (&p_adec->bit_stream, 6);
index_0 = p_adec->bit_stream.buffer >> (32 - 6);
DumpBits (&p_adec->bit_stream, 6);
if (allocation_1[sb] < 0) {
slope_1[0][sb] = slope_1[1][sb] = slope_1[2][sb] =
adec_scalefactor_table[index_0];
} else {
float scalefactor;
float slope, offset;
slope = adec_slope_table[allocation_1[sb]-2];
offset = adec_offset_table[allocation_1[sb]-2];
scalefactor = adec_scalefactor_table[index_0];
slope_1[0][sb] = slope_1[1][sb] = slope_1[2][sb] =
slope * scalefactor;
offset_1[0][sb] = offset_1[1][sb] = offset_1[2][sb] =
offset * scalefactor;
}
break;
case 3:
NeedBits (&p_adec->bit_stream, 12);
index_0 = p_adec->bit_stream.buffer >> (32 - 6);
index_1 = (p_adec->bit_stream.buffer >> (32 - 12)) & 63;
DumpBits (&p_adec->bit_stream, 12);
if (allocation_1[sb] < 0) {
slope_1[0][sb] = adec_scalefactor_table[index_0];
slope_1[1][sb] = slope_1[2][sb] =
adec_scalefactor_table[index_1];
} else {
float scalefactor;
float slope, offset;
slope = adec_slope_table[allocation_1[sb]-2];
offset = adec_offset_table[allocation_1[sb]-2];
scalefactor = adec_scalefactor_table[index_0];
slope_1[0][sb] = slope * scalefactor;
offset_1[0][sb] = offset * scalefactor;
scalefactor = adec_scalefactor_table[index_1];
slope_1[1][sb] = slope_1[2][sb] = slope * scalefactor;
offset_1[1][sb] = offset_1[2][sb] = offset * scalefactor;
}
break;
}
}
}
/* #define READ_SAMPLE_L2S */
i_gr = 0;
p_s16 = buffer;
/* parse samples */
READ_SAMPLE_L2S (pf_scalefactor_0_0, pf_scalefactor_1_0, 4)
READ_SAMPLE_L2S (pf_scalefactor_0_1, pf_scalefactor_1_1, 8)
READ_SAMPLE_L2S (pf_scalefactor_0_2, pf_scalefactor_1_2, 12)
for (gr0 = 0; gr0 < 3; gr0++)
for (gr1 = 0; gr1 < 4; gr1++) {
s16 * XXX_buf;
p_adec->bit_stream.buffer = 0;
p_adec->bit_stream.i_available = 0;
return (6);
}
for (sb = 0; sb < bound; sb++) {
int code;
/**** wkn ****/
switch (allocation_0[sb]) {
case 0:
sample_0[0][sb] = sample_0[1][sb] = sample_0[2][sb] = 0;
break;
case L3:
NeedBits (&p_adec->bit_stream, 5);
code = p_adec->bit_stream.buffer >> (32 - 5);
DumpBits (&p_adec->bit_stream, 5);
sample_0[0][sb] = slope_0[gr0][sb] * L3_table[code % 3];
code /= 3;
sample_0[1][sb] = slope_0[gr0][sb] * L3_table[code % 3];
code /= 3;
sample_0[2][sb] = slope_0[gr0][sb] * L3_table[code];
break;
case L5:
NeedBits (&p_adec->bit_stream, 7);
code = p_adec->bit_stream.buffer >> (32 - 7);
DumpBits (&p_adec->bit_stream, 7);
sample_0[0][sb] = slope_0[gr0][sb] * L5_table[code % 5];
code /= 5;
sample_0[1][sb] = slope_0[gr0][sb] * L5_table[code % 5];
code /= 5;
sample_0[2][sb] = slope_0[gr0][sb] * L5_table[code];
break;
case L9:
NeedBits (&p_adec->bit_stream, 10);
code = p_adec->bit_stream.buffer >> (32 - 10);
DumpBits (&p_adec->bit_stream, 10);
sample_0[0][sb] = slope_0[gr0][sb] * L9_table[code % 9];
code /= 9;
sample_0[1][sb] = slope_0[gr0][sb] * L9_table[code % 9];
code /= 9;
sample_0[2][sb] = slope_0[gr0][sb] * L9_table[code];
break;
default:
for (s = 0; s < 3; s++) {
NeedBits (&p_adec->bit_stream, allocation_0[sb]);
code = (p_adec->bit_stream.buffer >>
(32 - allocation_0[sb]));
DumpBits (&p_adec->bit_stream, allocation_0[sb]);
sample_0[s][sb] =
slope_0[gr0][sb] * code + offset_0[gr0][sb];
}
}
switch (allocation_1[sb]) {
case 0:
sample_1[0][sb] = sample_1[1][sb] = sample_1[2][sb] = 0;
break;
case L3:
NeedBits (&p_adec->bit_stream, 5);
code = p_adec->bit_stream.buffer >> (32 - 5);
DumpBits (&p_adec->bit_stream, 5);
sample_1[0][sb] = slope_1[gr0][sb] * L3_table[code % 3];
code /= 3;
sample_1[1][sb] = slope_1[gr0][sb] * L3_table[code % 3];
code /= 3;
sample_1[2][sb] = slope_1[gr0][sb] * L3_table[code];
break;
case L5:
NeedBits (&p_adec->bit_stream, 7);
code = p_adec->bit_stream.buffer >> (32 - 7);
DumpBits (&p_adec->bit_stream, 7);
sample_1[0][sb] = slope_1[gr0][sb] * L5_table[code % 5];
code /= 5;
sample_1[1][sb] = slope_1[gr0][sb] * L5_table[code % 5];
code /= 5;
sample_1[2][sb] = slope_1[gr0][sb] * L5_table[code];
break;
case L9:
NeedBits (&p_adec->bit_stream, 10);
code = p_adec->bit_stream.buffer >> (32 - 10);
DumpBits (&p_adec->bit_stream, 10);
sample_1[0][sb] = slope_1[gr0][sb] * L9_table[code % 9];
code /= 9;
sample_1[1][sb] = slope_1[gr0][sb] * L9_table[code % 9];
code /= 9;
sample_1[2][sb] = slope_1[gr0][sb] * L9_table[code];
break;
default:
for (s = 0; s < 3; s++) {
NeedBits (&p_adec->bit_stream, allocation_1[sb]);
code = (p_adec->bit_stream.buffer >>
(32 - allocation_1[sb]));
DumpBits (&p_adec->bit_stream, allocation_1[sb]);
sample_1[s][sb] =
slope_1[gr0][sb] * code + offset_1[gr0][sb];
}
}
}
for (; sb < sblimit; sb++) {
int code;
switch (allocation_0[sb]) {
case 0:
sample_0[0][sb] = sample_0[1][sb] = sample_0[2][sb] = 0;
sample_1[0][sb] = sample_1[1][sb] = sample_1[2][sb] = 0;
break;
case L3:
NeedBits (&p_adec->bit_stream, 5);
code = p_adec->bit_stream.buffer >> (32 - 5);
DumpBits (&p_adec->bit_stream, 5);
sample_0[0][sb] = slope_0[gr0][sb] * L3_table[code % 3];
sample_1[0][sb] = slope_1[gr0][sb] * L3_table[code % 3];
code /= 3;
sample_0[1][sb] = slope_0[gr0][sb] * L3_table[code % 3];
sample_1[1][sb] = slope_1[gr0][sb] * L3_table[code % 3];
code /= 3;
sample_0[2][sb] = slope_0[gr0][sb] * L3_table[code];
sample_1[2][sb] = slope_1[gr0][sb] * L3_table[code];
break;
case L5:
NeedBits (&p_adec->bit_stream, 7);
code = p_adec->bit_stream.buffer >> (32 - 7);
DumpBits (&p_adec->bit_stream, 7);
sample_0[0][sb] = slope_0[gr0][sb] * L5_table[code % 5];
sample_1[0][sb] = slope_1[gr0][sb] * L5_table[code % 5];
code /= 5;
sample_0[1][sb] = slope_0[gr0][sb] * L5_table[code % 5];
sample_1[1][sb] = slope_1[gr0][sb] * L5_table[code % 5];
code /= 5;
sample_0[2][sb] = slope_0[gr0][sb] * L5_table[code];
sample_1[2][sb] = slope_1[gr0][sb] * L5_table[code];
break;
case L9:
NeedBits (&p_adec->bit_stream, 10);
code = p_adec->bit_stream.buffer >> (32 - 10);
DumpBits (&p_adec->bit_stream, 10);
sample_0[0][sb] = slope_0[gr0][sb] * L9_table[code % 9];
sample_1[0][sb] = slope_1[gr0][sb] * L9_table[code % 9];
code /= 9;
sample_0[1][sb] = slope_0[gr0][sb] * L9_table[code % 9];
sample_1[1][sb] = slope_1[gr0][sb] * L9_table[code % 9];
code /= 9;
sample_0[2][sb] = slope_0[gr0][sb] * L9_table[code];
sample_1[2][sb] = slope_1[gr0][sb] * L9_table[code];
break;
default:
for (s = 0; s < 3; s++) {
NeedBits (&p_adec->bit_stream, allocation_0[sb]);
code = (p_adec->bit_stream.buffer >>
(32 - allocation_0[sb]));
DumpBits (&p_adec->bit_stream, allocation_0[sb]);
sample_0[s][sb] =
slope_0[gr0][sb] * code + offset_0[gr0][sb];
sample_1[s][sb] =
slope_1[gr0][sb] * code + offset_1[gr0][sb];
}
}
}
for (; sb < 32; sb++) {
sample_0[0][sb] = sample_0[1][sb] = sample_0[2][sb] = 0;
sample_1[0][sb] = sample_1[1][sb] = sample_1[2][sb] = 0;
}
for (s = 0; s < 3; s++) {
DCT32 (sample_0[s], &p_adec->bank_0);
XXX_buf = buffer;
PCM (&p_adec->bank_0, &XXX_buf, 2);
DCT32 (sample_1[s], &p_adec->bank_1);
XXX_buf = buffer+1;
PCM (&p_adec->bank_1, &XXX_buf, 2);
buffer += 64;
}
}
return 0;
}
int adec_init (audiodec_t * p_adec)
{
......@@ -526,56 +903,64 @@ int adec_sync_frame (audiodec_t * p_adec, adec_sync_info_t * p_sync_info)
static int mpeg2_layer2_bit_rate[15] = {
0, 8, 16, 24, 32, 40, 48, 56, 64, 80, 96, 112, 128, 144, 160
};
static int * bit_rate_table[8] = {
NULL, NULL, mpeg2_layer2_bit_rate, mpeg2_layer1_bit_rate,
NULL, NULL, mpeg1_layer2_bit_rate, mpeg1_layer1_bit_rate
};
u32 header;
int index;
int * bit_rate_table;
int * bit_rates;
int sample_rate;
int bit_rate;
int frame_size;
p_adec->bit_stream.total_bytes_read = 0;
header = GetByte (&p_adec->bit_stream) << 24;
header |= GetByte (&p_adec->bit_stream) << 16;
header |= GetByte (&p_adec->bit_stream) << 8;
header |= GetByte (&p_adec->bit_stream);
p_adec->header = header;
/* basic header check : sync word, no emphasis */
if ((header & 0xfff00003) != 0xfff00000)
return 1;
index = (header >> 10) & 3; /* sample rate index */
if (index > 2)
return 1;
sample_rate = mpeg1_sample_rate[index];
/* calculate bit rate */
switch ((header >> 17) & 7) {
case 2: /* mpeg 2, layer 2 */
sample_rate >>= 1; /* half sample rate for mpeg2 */
bit_rate_table = mpeg2_layer2_bit_rate;
break;
case 3: /* mpeg 2, layer 1 */
sample_rate >>= 1; /* half sample rate for mpeg2 */
bit_rate_table = mpeg2_layer1_bit_rate;
break;
case 6: /* mpeg1, layer 2 */
bit_rate_table = mpeg1_layer2_bit_rate;
break;
case 7: /* mpeg1, layer 1 */
bit_rate_table = mpeg1_layer1_bit_rate;
break;
default: /* invalid layer */
return 1;
}
index = (header >> 17) & 7; /* mpeg ID + layer */
bit_rates = bit_rate_table[index];
if (bit_rate_table == NULL)
return 1; /* invalid layer */
index = (header >> 12) & 15; /* bit rate index */
if (index > 14)
return 1;
bit_rate = bit_rate_table[index];
bit_rate = bit_rates[index];
p_sync_info->sample_rate = sample_rate;
p_sync_info->bit_rate = bit_rate;
/* mpeg 1 layer 2 : check that bitrate per channel is valid */
if (bit_rates == mpeg1_layer2_bit_rate) {
if ((header & 0xc0) == 0xc0) { /* mono */
if (index > 10)
return 1; /* invalid bitrate per channel */
} else { /* stereo */
if ((1 << index) & 0x2e)
return 1; /* invalid bitrate per channel */
}
}
/* calculate sample rate */
index = (header >> 10) & 3; /* sample rate index */
if (index > 2)
return 1;
sample_rate = mpeg1_sample_rate[index];
if (!(header & 0x80000))
sample_rate >>= 1; /* half sample rate for mpeg2 */
/* calculate frame length */
if ((header & 0x60000) == 0x60000) { /* layer 1 */
frame_size = 48000 * bit_rate / sample_rate;
......@@ -587,6 +972,8 @@ int adec_sync_frame (audiodec_t * p_adec, adec_sync_info_t * p_sync_info)
frame_size ++;
}
p_sync_info->sample_rate = sample_rate;
p_sync_info->bit_rate = bit_rate;
p_sync_info->frame_size = frame_size;
p_adec->frame_size = frame_size;
......@@ -595,15 +982,53 @@ int adec_sync_frame (audiodec_t * p_adec, adec_sync_info_t * p_sync_info)
int adec_decode_frame (audiodec_t * p_adec, s16 * buffer)
{
if (!(p_adec->header & 0x10000)) { /* error check, skip it */
GetByte (&p_adec->bit_stream);
GetByte (&p_adec->bit_stream);
}
/* parse audio data */
p_adec->bit_stream.i_available = 0;
adec_Layer2_Stereo (p_adec, buffer);
switch ((p_adec->header >> 17) & 3) {
case 2: /* layer 2 */
if ((p_adec->header & 0xc0) == 0xc0) {
if (adec_layer2_mono (p_adec, buffer))
return 1;
} else {
if (adec_layer2_stereo (p_adec, buffer))
return 1;
}
break;
case 3: /* layer 1 */
if ((p_adec->header & 0xc0) == 0xc0) {
if (adec_layer1_mono (p_adec, buffer))
return 1;
} else {
if (adec_layer1_stereo (p_adec, buffer))
return 1;
}
break;
}
/* skip ancillary data */
if ((p_adec->header & 0xf000) == 0) /* free bitrate format */
return 0;
/* XXX rewrite the byte counting system to reduce overhead */
#if 0
printf ("skip %d\n",
p_adec->frame_size - p_adec->bit_stream.total_bytes_read);
#endif
if (p_adec->bit_stream.total_bytes_read > p_adec->frame_size)
return 1;
return 1; /* overrun */
while (p_adec->bit_stream.total_bytes_read < p_adec->frame_size)
GetByte (&p_adec->bit_stream);
GetByte (&p_adec->bit_stream); /* skip ancillary data */
return 0;
}
src/audio_decoder/audio_decoder_thread.c
View file @ fec6ded4
/*****************************************************************************
* audio_decoder.c: MPEG1 Layer I-II audio decoder thread
* audio_decoder_thread.c: MPEG1 Layer I-II audio decoder thread
*****************************************************************************
* Copyright (C) 1999, 2000 VideoLAN
*
......@@ -56,7 +56,6 @@
#include "audio_output.h" /* aout_fifo_t (for audio_decoder.h) */
#include "audio_constants.h"
#include "audio_decoder.h"
#include "audio_decoder_thread.h"
#include "audio_math.h" /* DCT32(), PCM() */
......
src/audio_decoder/audio_math.c
View file @ fec6ded4
......@@ -313,7 +313,7 @@ void DCT32(float *x, adec_bank_t *b)
void PCM(adec_bank_t *b, s16 **pcm, int jump)
{
/* scale factor */
#define F 32768
#define F -32768
/* These values are not in the same order as in Annex 3-B.3 of the ISO/IEC
DIS 11172-3 */
static const float c[512] = {
......
src/audio_decoder/audio_test.c 0 → 100644
View file @ fec6ded4
/*****************************************************************************
* audio_test.c: MPEG1 Layer I-II audio decoder test program
*****************************************************************************
* Copyright (C) 1999, 2000 VideoLAN
*
* Authors:
*
* This program is free software; you can redistribute it and/or modify
* it under the terms of the GNU General Public License as published by
* the Free Software Foundation; either version 2 of the License, or
* (at your option) any later version.
*
* This program is distributed in the hope that it will be useful,
* but WITHOUT ANY WARRANTY; without even the implied warranty of
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
* General Public License for more details.
*
* You should have received a copy of the GNU General Public
* License along with this program; if not, write to the
* Free Software Foundation, Inc., 59 Temple Place - Suite 330,
* Boston, MA 02111-1307, USA.
*****************************************************************************/
/*
* TODO :
*
* - optimiser les NeedBits() et les GetBits() du code l o c'est possible ;
* - vlc_cond_signal() / vlc_cond_wait() ;
*
*/
/*****************************************************************************
* Preamble
*****************************************************************************/
#include <stdio.h>
#include <string.h>
#include "int_types.h"
#include "audio_decoder.h"
#define ADEC_FRAME_SIZE (2*1152)
int main (void)
{
audiodec_t decoder;
adec_sync_info_t sync_info;
adec_byte_stream_t * stream;
s16 buffer [ADEC_FRAME_SIZE];
int framenum;
memset (&decoder, 0, sizeof (decoder));
if (adec_init (&decoder))
return 1;
stream = adec_byte_stream (&decoder);
stream->p_byte = NULL;
stream->p_end = NULL;
stream->info = stdin;
framenum = 0;
while (1) {
int i;
if (adec_sync_frame (&decoder, &sync_info))
return 1;
if (adec_decode_frame (&decoder, buffer))
return 1;
#if 1
for (i = 0; i < (2*1152); i++)
printf ("%04X\n",(u16)buffer[i]);
#endif
}
return 0;
}
void adec_byte_stream_next (adec_byte_stream_t * p_byte_stream)
{
static u8 buffer [1024];
static u8 dummy = 0;
FILE * fd;
int size;
fd = p_byte_stream->info;
size = fread (buffer, 1, 1024, fd);
if (size) {
p_byte_stream->p_byte = buffer;
p_byte_stream->p_end = buffer + size;
} else { /* end of stream, read dummy zeroes */
p_byte_stream->p_byte = &dummy;
p_byte_stream->p_end = &dummy + 1;
}
}
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment