Mercurial > libavcodec.hg
annotate ac3dec.c @ 7601:69d5f318275f libavcodec
use LFG instead of Mersenne Twister for AC-3 PRNG
author | jbr |
---|---|
date | Sun, 17 Aug 2008 17:41:48 +0000 |
parents | cd6217c9ce92 |
children | 20c91136fe71 |
rev | line source |
---|---|
5328 | 1 /* |
2 * AC-3 Audio Decoder | |
5316
a09b700fcfca
AC-3 decoder, soc revision 56, Aug 17 14:21:58 2006 UTC by cloud9
jbr
parents:
5315
diff
changeset
|
3 * This code is developed as part of Google Summer of Code 2006 Program. |
a09b700fcfca
AC-3 decoder, soc revision 56, Aug 17 14:21:58 2006 UTC by cloud9
jbr
parents:
5315
diff
changeset
|
4 * |
5328 | 5 * Copyright (c) 2006 Kartikey Mahendra BHATT (bhattkm at gmail dot com). |
7594
cd6217c9ce92
update my email address to one which does not depend on my service provider
jbr
parents:
7564
diff
changeset
|
6 * Copyright (c) 2007 Justin Ruggles <justin.ruggles@gmail.com> |
5316
a09b700fcfca
AC-3 decoder, soc revision 56, Aug 17 14:21:58 2006 UTC by cloud9
jbr
parents:
5315
diff
changeset
|
7 * |
5377
2020fab0e37a
update license header. add copyright for myself. correct and complete the statement about liba52.
jbr
parents:
5376
diff
changeset
|
8 * Portions of this code are derived from liba52 |
5316
a09b700fcfca
AC-3 decoder, soc revision 56, Aug 17 14:21:58 2006 UTC by cloud9
jbr
parents:
5315
diff
changeset
|
9 * http://liba52.sourceforge.net |
5377
2020fab0e37a
update license header. add copyright for myself. correct and complete the statement about liba52.
jbr
parents:
5376
diff
changeset
|
10 * Copyright (C) 2000-2003 Michel Lespinasse <walken@zoy.org> |
2020fab0e37a
update license header. add copyright for myself. correct and complete the statement about liba52.
jbr
parents:
5376
diff
changeset
|
11 * Copyright (C) 1999-2000 Aaron Holtzman <aholtzma@ess.engr.uvic.ca> |
5316
a09b700fcfca
AC-3 decoder, soc revision 56, Aug 17 14:21:58 2006 UTC by cloud9
jbr
parents:
5315
diff
changeset
|
12 * |
5328 | 13 * This file is part of FFmpeg. |
5303
b8821ff5c30d
AC-3 decoder, soc revision 9, Jun 14 03:52:02 2006 UTC by cloud9
jbr
parents:
diff
changeset
|
14 * |
5328 | 15 * FFmpeg is free software; you can redistribute it and/or |
5329 | 16 * modify it under the terms of the GNU General Public |
5303
b8821ff5c30d
AC-3 decoder, soc revision 9, Jun 14 03:52:02 2006 UTC by cloud9
jbr
parents:
diff
changeset
|
17 * License as published by the Free Software Foundation; either |
b8821ff5c30d
AC-3 decoder, soc revision 9, Jun 14 03:52:02 2006 UTC by cloud9
jbr
parents:
diff
changeset
|
18 * version 2 of the License, or (at your option) any later version. |
b8821ff5c30d
AC-3 decoder, soc revision 9, Jun 14 03:52:02 2006 UTC by cloud9
jbr
parents:
diff
changeset
|
19 * |
5328 | 20 * FFmpeg is distributed in the hope that it will be useful, |
5303
b8821ff5c30d
AC-3 decoder, soc revision 9, Jun 14 03:52:02 2006 UTC by cloud9
jbr
parents:
diff
changeset
|
21 * but WITHOUT ANY WARRANTY; without even the implied warranty of |
b8821ff5c30d
AC-3 decoder, soc revision 9, Jun 14 03:52:02 2006 UTC by cloud9
jbr
parents:
diff
changeset
|
22 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU |
5326 | 23 * General Public License for more details. |
5303
b8821ff5c30d
AC-3 decoder, soc revision 9, Jun 14 03:52:02 2006 UTC by cloud9
jbr
parents:
diff
changeset
|
24 * |
5326 | 25 * You should have received a copy of the GNU General Public |
5328 | 26 * License along with FFmpeg; if not, write to the Free Software |
5303
b8821ff5c30d
AC-3 decoder, soc revision 9, Jun 14 03:52:02 2006 UTC by cloud9
jbr
parents:
diff
changeset
|
27 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA |
b8821ff5c30d
AC-3 decoder, soc revision 9, Jun 14 03:52:02 2006 UTC by cloud9
jbr
parents:
diff
changeset
|
28 */ |
b8821ff5c30d
AC-3 decoder, soc revision 9, Jun 14 03:52:02 2006 UTC by cloud9
jbr
parents:
diff
changeset
|
29 |
b8821ff5c30d
AC-3 decoder, soc revision 9, Jun 14 03:52:02 2006 UTC by cloud9
jbr
parents:
diff
changeset
|
30 #include <stdio.h> |
b8821ff5c30d
AC-3 decoder, soc revision 9, Jun 14 03:52:02 2006 UTC by cloud9
jbr
parents:
diff
changeset
|
31 #include <stddef.h> |
b8821ff5c30d
AC-3 decoder, soc revision 9, Jun 14 03:52:02 2006 UTC by cloud9
jbr
parents:
diff
changeset
|
32 #include <math.h> |
b8821ff5c30d
AC-3 decoder, soc revision 9, Jun 14 03:52:02 2006 UTC by cloud9
jbr
parents:
diff
changeset
|
33 #include <string.h> |
b8821ff5c30d
AC-3 decoder, soc revision 9, Jun 14 03:52:02 2006 UTC by cloud9
jbr
parents:
diff
changeset
|
34 |
6763 | 35 #include "libavutil/crc.h" |
5315
31033309ebc9
AC-3 decoder, soc revision 55, Aug 17 12:25:24 2006 UTC by cloud9
jbr
parents:
5314
diff
changeset
|
36 #include "avcodec.h" |
5386
d6347c3cb329
Use shared AC-3 header parsing. Move some initializations to block parsing function.
jbr
parents:
5385
diff
changeset
|
37 #include "ac3_parser.h" |
5303
b8821ff5c30d
AC-3 decoder, soc revision 9, Jun 14 03:52:02 2006 UTC by cloud9
jbr
parents:
diff
changeset
|
38 #include "bitstream.h" |
b8821ff5c30d
AC-3 decoder, soc revision 9, Jun 14 03:52:02 2006 UTC by cloud9
jbr
parents:
diff
changeset
|
39 #include "dsputil.h" |
7020
40a2db8dbaa2
move the decode context and some macro constants to a new file, ac3dec.h
jbr
parents:
7019
diff
changeset
|
40 #include "ac3dec.h" |
7033 | 41 #include "ac3dec_data.h" |
5309
0662a270aab7
AC-3 decoder, soc revision 38, Aug 7 00:03:00 2006 UTC by cloud9
jbr
parents:
5308
diff
changeset
|
42 |
6536
8f2186d5daca
additional protection from segmentation faults and memory access errors by
jbr
parents:
6529
diff
changeset
|
43 /** Maximum possible frame size when the specification limit is ignored */ |
8f2186d5daca
additional protection from segmentation faults and memory access errors by
jbr
parents:
6529
diff
changeset
|
44 #define AC3_MAX_FRAME_SIZE 21695 |
8f2186d5daca
additional protection from segmentation faults and memory access errors by
jbr
parents:
6529
diff
changeset
|
45 |
7471
1788013cf23b
use the same table for ungrouping exponents and bap=2 mantissas
jbr
parents:
7470
diff
changeset
|
46 /** |
1788013cf23b
use the same table for ungrouping exponents and bap=2 mantissas
jbr
parents:
7470
diff
changeset
|
47 * table for ungrouping 3 values in 7 bits. |
1788013cf23b
use the same table for ungrouping exponents and bap=2 mantissas
jbr
parents:
7470
diff
changeset
|
48 * used for exponents and bap=2 mantissas |
1788013cf23b
use the same table for ungrouping exponents and bap=2 mantissas
jbr
parents:
7470
diff
changeset
|
49 */ |
1788013cf23b
use the same table for ungrouping exponents and bap=2 mantissas
jbr
parents:
7470
diff
changeset
|
50 static uint8_t ungroup_3_in_7_bits_tab[128][3]; |
5325
0d963101e0d7
AC-3 decoder, soc revision 402, Jul 14 13:45:14 2007 UTC by jbr
jbr
parents:
5324
diff
changeset
|
51 |
0d963101e0d7
AC-3 decoder, soc revision 402, Jul 14 13:45:14 2007 UTC by jbr
jbr
parents:
5324
diff
changeset
|
52 |
5456 | 53 /** tables for ungrouping mantissas */ |
6430
de7502093922
Use 24-bit fixed-point transform coefficients until just before MDCT. This
jbr
parents:
6142
diff
changeset
|
54 static int b1_mantissas[32][3]; |
de7502093922
Use 24-bit fixed-point transform coefficients until just before MDCT. This
jbr
parents:
6142
diff
changeset
|
55 static int b2_mantissas[128][3]; |
de7502093922
Use 24-bit fixed-point transform coefficients until just before MDCT. This
jbr
parents:
6142
diff
changeset
|
56 static int b3_mantissas[8]; |
de7502093922
Use 24-bit fixed-point transform coefficients until just before MDCT. This
jbr
parents:
6142
diff
changeset
|
57 static int b4_mantissas[128][2]; |
de7502093922
Use 24-bit fixed-point transform coefficients until just before MDCT. This
jbr
parents:
6142
diff
changeset
|
58 static int b5_mantissas[16]; |
5325
0d963101e0d7
AC-3 decoder, soc revision 402, Jul 14 13:45:14 2007 UTC by jbr
jbr
parents:
5324
diff
changeset
|
59 |
5456 | 60 /** |
61 * Quantization table: levels for symmetric. bits for asymmetric. | |
62 * reference: Table 7.18 Mapping of bap to Quantizer | |
63 */ | |
6006 | 64 static const uint8_t quantization_tab[16] = { |
5456 | 65 0, 3, 5, 7, 11, 15, |
66 5, 6, 7, 8, 9, 10, 11, 12, 14, 16 | |
67 }; | |
5325
0d963101e0d7
AC-3 decoder, soc revision 402, Jul 14 13:45:14 2007 UTC by jbr
jbr
parents:
5324
diff
changeset
|
68 |
5457
0a5849a4858b
use table for converting dynamic range codes to scale factors
jbr
parents:
5456
diff
changeset
|
69 /** dynamic range table. converts codes to scale factors. */ |
6006 | 70 static float dynamic_range_tab[256]; |
5457
0a5849a4858b
use table for converting dynamic range codes to scale factors
jbr
parents:
5456
diff
changeset
|
71 |
5515
b499866a3397
comments modification: add, clean up, Doxygenize, and reword
jbr
parents:
5499
diff
changeset
|
72 /** Adjustments in dB gain */ |
7018
885cad616170
add more gain levels and adjust mix level tables accordingly.
jbr
parents:
7017
diff
changeset
|
73 #define LEVEL_PLUS_3DB 1.4142135623730950 |
885cad616170
add more gain levels and adjust mix level tables accordingly.
jbr
parents:
7017
diff
changeset
|
74 #define LEVEL_PLUS_1POINT5DB 1.1892071150027209 |
885cad616170
add more gain levels and adjust mix level tables accordingly.
jbr
parents:
7017
diff
changeset
|
75 #define LEVEL_MINUS_1POINT5DB 0.8408964152537145 |
5325
0d963101e0d7
AC-3 decoder, soc revision 402, Jul 14 13:45:14 2007 UTC by jbr
jbr
parents:
5324
diff
changeset
|
76 #define LEVEL_MINUS_3DB 0.7071067811865476 |
0d963101e0d7
AC-3 decoder, soc revision 402, Jul 14 13:45:14 2007 UTC by jbr
jbr
parents:
5324
diff
changeset
|
77 #define LEVEL_MINUS_4POINT5DB 0.5946035575013605 |
0d963101e0d7
AC-3 decoder, soc revision 402, Jul 14 13:45:14 2007 UTC by jbr
jbr
parents:
5324
diff
changeset
|
78 #define LEVEL_MINUS_6DB 0.5000000000000000 |
5480 | 79 #define LEVEL_MINUS_9DB 0.3535533905932738 |
5325
0d963101e0d7
AC-3 decoder, soc revision 402, Jul 14 13:45:14 2007 UTC by jbr
jbr
parents:
5324
diff
changeset
|
80 #define LEVEL_ZERO 0.0000000000000000 |
5480 | 81 #define LEVEL_ONE 1.0000000000000000 |
82 | |
7018
885cad616170
add more gain levels and adjust mix level tables accordingly.
jbr
parents:
7017
diff
changeset
|
83 static const float gain_levels[9] = { |
885cad616170
add more gain levels and adjust mix level tables accordingly.
jbr
parents:
7017
diff
changeset
|
84 LEVEL_PLUS_3DB, |
885cad616170
add more gain levels and adjust mix level tables accordingly.
jbr
parents:
7017
diff
changeset
|
85 LEVEL_PLUS_1POINT5DB, |
5480 | 86 LEVEL_ONE, |
7018
885cad616170
add more gain levels and adjust mix level tables accordingly.
jbr
parents:
7017
diff
changeset
|
87 LEVEL_MINUS_1POINT5DB, |
5480 | 88 LEVEL_MINUS_3DB, |
89 LEVEL_MINUS_4POINT5DB, | |
90 LEVEL_MINUS_6DB, | |
7018
885cad616170
add more gain levels and adjust mix level tables accordingly.
jbr
parents:
7017
diff
changeset
|
91 LEVEL_ZERO, |
5480 | 92 LEVEL_MINUS_9DB |
93 }; | |
94 | |
95 /** | |
7019
81d5c68233e5
move mix level tables from parser to decoder. have parser read bitstream value instead of using an index to a table in the decoder.
jbr
parents:
7018
diff
changeset
|
96 * Table for center mix levels |
81d5c68233e5
move mix level tables from parser to decoder. have parser read bitstream value instead of using an index to a table in the decoder.
jbr
parents:
7018
diff
changeset
|
97 * reference: Section 5.4.2.4 cmixlev |
81d5c68233e5
move mix level tables from parser to decoder. have parser read bitstream value instead of using an index to a table in the decoder.
jbr
parents:
7018
diff
changeset
|
98 */ |
81d5c68233e5
move mix level tables from parser to decoder. have parser read bitstream value instead of using an index to a table in the decoder.
jbr
parents:
7018
diff
changeset
|
99 static const uint8_t center_levels[4] = { 4, 5, 6, 5 }; |
81d5c68233e5
move mix level tables from parser to decoder. have parser read bitstream value instead of using an index to a table in the decoder.
jbr
parents:
7018
diff
changeset
|
100 |
81d5c68233e5
move mix level tables from parser to decoder. have parser read bitstream value instead of using an index to a table in the decoder.
jbr
parents:
7018
diff
changeset
|
101 /** |
81d5c68233e5
move mix level tables from parser to decoder. have parser read bitstream value instead of using an index to a table in the decoder.
jbr
parents:
7018
diff
changeset
|
102 * Table for surround mix levels |
81d5c68233e5
move mix level tables from parser to decoder. have parser read bitstream value instead of using an index to a table in the decoder.
jbr
parents:
7018
diff
changeset
|
103 * reference: Section 5.4.2.5 surmixlev |
81d5c68233e5
move mix level tables from parser to decoder. have parser read bitstream value instead of using an index to a table in the decoder.
jbr
parents:
7018
diff
changeset
|
104 */ |
81d5c68233e5
move mix level tables from parser to decoder. have parser read bitstream value instead of using an index to a table in the decoder.
jbr
parents:
7018
diff
changeset
|
105 static const uint8_t surround_levels[4] = { 4, 6, 7, 6 }; |
81d5c68233e5
move mix level tables from parser to decoder. have parser read bitstream value instead of using an index to a table in the decoder.
jbr
parents:
7018
diff
changeset
|
106 |
81d5c68233e5
move mix level tables from parser to decoder. have parser read bitstream value instead of using an index to a table in the decoder.
jbr
parents:
7018
diff
changeset
|
107 /** |
5480 | 108 * Table for default stereo downmixing coefficients |
109 * reference: Section 7.8.2 Downmixing Into Two Channels | |
110 */ | |
111 static const uint8_t ac3_default_coeffs[8][5][2] = { | |
7018
885cad616170
add more gain levels and adjust mix level tables accordingly.
jbr
parents:
7017
diff
changeset
|
112 { { 2, 7 }, { 7, 2 }, }, |
885cad616170
add more gain levels and adjust mix level tables accordingly.
jbr
parents:
7017
diff
changeset
|
113 { { 4, 4 }, }, |
885cad616170
add more gain levels and adjust mix level tables accordingly.
jbr
parents:
7017
diff
changeset
|
114 { { 2, 7 }, { 7, 2 }, }, |
885cad616170
add more gain levels and adjust mix level tables accordingly.
jbr
parents:
7017
diff
changeset
|
115 { { 2, 7 }, { 5, 5 }, { 7, 2 }, }, |
885cad616170
add more gain levels and adjust mix level tables accordingly.
jbr
parents:
7017
diff
changeset
|
116 { { 2, 7 }, { 7, 2 }, { 6, 6 }, }, |
885cad616170
add more gain levels and adjust mix level tables accordingly.
jbr
parents:
7017
diff
changeset
|
117 { { 2, 7 }, { 5, 5 }, { 7, 2 }, { 8, 8 }, }, |
885cad616170
add more gain levels and adjust mix level tables accordingly.
jbr
parents:
7017
diff
changeset
|
118 { { 2, 7 }, { 7, 2 }, { 6, 7 }, { 7, 6 }, }, |
885cad616170
add more gain levels and adjust mix level tables accordingly.
jbr
parents:
7017
diff
changeset
|
119 { { 2, 7 }, { 5, 5 }, { 7, 2 }, { 6, 7 }, { 7, 6 }, }, |
5480 | 120 }; |
5325
0d963101e0d7
AC-3 decoder, soc revision 402, Jul 14 13:45:14 2007 UTC by jbr
jbr
parents:
5324
diff
changeset
|
121 |
5310
9aa9197034d7
AC-3 decoder, soc revision 40, Aug 9 00:10:14 2006 UTC by cloud9
jbr
parents:
5309
diff
changeset
|
122 /** |
5515
b499866a3397
comments modification: add, clean up, Doxygenize, and reword
jbr
parents:
5499
diff
changeset
|
123 * Symmetrical Dequantization |
b499866a3397
comments modification: add, clean up, Doxygenize, and reword
jbr
parents:
5499
diff
changeset
|
124 * reference: Section 7.3.3 Expansion of Mantissas for Symmetrical Quantization |
b499866a3397
comments modification: add, clean up, Doxygenize, and reword
jbr
parents:
5499
diff
changeset
|
125 * Tables 7.19 to 7.23 |
b499866a3397
comments modification: add, clean up, Doxygenize, and reword
jbr
parents:
5499
diff
changeset
|
126 */ |
6430
de7502093922
Use 24-bit fixed-point transform coefficients until just before MDCT. This
jbr
parents:
6142
diff
changeset
|
127 static inline int |
5456 | 128 symmetric_dequant(int code, int levels) |
5309
0662a270aab7
AC-3 decoder, soc revision 38, Aug 7 00:03:00 2006 UTC by cloud9
jbr
parents:
5308
diff
changeset
|
129 { |
6430
de7502093922
Use 24-bit fixed-point transform coefficients until just before MDCT. This
jbr
parents:
6142
diff
changeset
|
130 return ((code - (levels >> 1)) << 24) / levels; |
5309
0662a270aab7
AC-3 decoder, soc revision 38, Aug 7 00:03:00 2006 UTC by cloud9
jbr
parents:
5308
diff
changeset
|
131 } |
0662a270aab7
AC-3 decoder, soc revision 38, Aug 7 00:03:00 2006 UTC by cloud9
jbr
parents:
5308
diff
changeset
|
132 |
5321
fcb2a52c88fe
AC-3 decoder, soc revision 71, Sep 4 09:47:12 2006 UTC by cloud9
jbr
parents:
5320
diff
changeset
|
133 /* |
fcb2a52c88fe
AC-3 decoder, soc revision 71, Sep 4 09:47:12 2006 UTC by cloud9
jbr
parents:
5320
diff
changeset
|
134 * Initialize tables at runtime. |
fcb2a52c88fe
AC-3 decoder, soc revision 71, Sep 4 09:47:12 2006 UTC by cloud9
jbr
parents:
5320
diff
changeset
|
135 */ |
6517
48759bfbd073
Apply 'cold' attribute to init/uninit functions in libavcodec
zuxy
parents:
6487
diff
changeset
|
136 static av_cold void ac3_tables_init(void) |
5305
5892b4a6380b
AC-3 decoder, soc revision 31, Jul 14 23:53:28 2006 UTC by cloud9
jbr
parents:
5304
diff
changeset
|
137 { |
5375 | 138 int i; |
5309
0662a270aab7
AC-3 decoder, soc revision 38, Aug 7 00:03:00 2006 UTC by cloud9
jbr
parents:
5308
diff
changeset
|
139 |
7471
1788013cf23b
use the same table for ungrouping exponents and bap=2 mantissas
jbr
parents:
7470
diff
changeset
|
140 /* generate table for ungrouping 3 values in 7 bits |
1788013cf23b
use the same table for ungrouping exponents and bap=2 mantissas
jbr
parents:
7470
diff
changeset
|
141 reference: Section 7.1.3 Exponent Decoding */ |
1788013cf23b
use the same table for ungrouping exponents and bap=2 mantissas
jbr
parents:
7470
diff
changeset
|
142 for(i=0; i<128; i++) { |
1788013cf23b
use the same table for ungrouping exponents and bap=2 mantissas
jbr
parents:
7470
diff
changeset
|
143 ungroup_3_in_7_bits_tab[i][0] = i / 25; |
1788013cf23b
use the same table for ungrouping exponents and bap=2 mantissas
jbr
parents:
7470
diff
changeset
|
144 ungroup_3_in_7_bits_tab[i][1] = (i % 25) / 5; |
1788013cf23b
use the same table for ungrouping exponents and bap=2 mantissas
jbr
parents:
7470
diff
changeset
|
145 ungroup_3_in_7_bits_tab[i][2] = (i % 25) % 5; |
1788013cf23b
use the same table for ungrouping exponents and bap=2 mantissas
jbr
parents:
7470
diff
changeset
|
146 } |
1788013cf23b
use the same table for ungrouping exponents and bap=2 mantissas
jbr
parents:
7470
diff
changeset
|
147 |
5456 | 148 /* generate grouped mantissa tables |
149 reference: Section 7.3.5 Ungrouping of Mantissas */ | |
150 for(i=0; i<32; i++) { | |
151 /* bap=1 mantissas */ | |
7472
83c16f531a90
Hardcode the table for ungrouping 3 values in 5 bits. It will be reused by
jbr
parents:
7471
diff
changeset
|
152 b1_mantissas[i][0] = symmetric_dequant(ff_ac3_ungroup_3_in_5_bits_tab[i][0], 3); |
83c16f531a90
Hardcode the table for ungrouping 3 values in 5 bits. It will be reused by
jbr
parents:
7471
diff
changeset
|
153 b1_mantissas[i][1] = symmetric_dequant(ff_ac3_ungroup_3_in_5_bits_tab[i][1], 3); |
83c16f531a90
Hardcode the table for ungrouping 3 values in 5 bits. It will be reused by
jbr
parents:
7471
diff
changeset
|
154 b1_mantissas[i][2] = symmetric_dequant(ff_ac3_ungroup_3_in_5_bits_tab[i][2], 3); |
5456 | 155 } |
156 for(i=0; i<128; i++) { | |
157 /* bap=2 mantissas */ | |
7471
1788013cf23b
use the same table for ungrouping exponents and bap=2 mantissas
jbr
parents:
7470
diff
changeset
|
158 b2_mantissas[i][0] = symmetric_dequant(ungroup_3_in_7_bits_tab[i][0], 5); |
1788013cf23b
use the same table for ungrouping exponents and bap=2 mantissas
jbr
parents:
7470
diff
changeset
|
159 b2_mantissas[i][1] = symmetric_dequant(ungroup_3_in_7_bits_tab[i][1], 5); |
1788013cf23b
use the same table for ungrouping exponents and bap=2 mantissas
jbr
parents:
7470
diff
changeset
|
160 b2_mantissas[i][2] = symmetric_dequant(ungroup_3_in_7_bits_tab[i][2], 5); |
5309
0662a270aab7
AC-3 decoder, soc revision 38, Aug 7 00:03:00 2006 UTC by cloud9
jbr
parents:
5308
diff
changeset
|
161 |
5456 | 162 /* bap=4 mantissas */ |
163 b4_mantissas[i][0] = symmetric_dequant(i / 11, 11); | |
164 b4_mantissas[i][1] = symmetric_dequant(i % 11, 11); | |
165 } | |
166 /* generate ungrouped mantissa tables | |
167 reference: Tables 7.21 and 7.23 */ | |
168 for(i=0; i<7; i++) { | |
169 /* bap=3 mantissas */ | |
170 b3_mantissas[i] = symmetric_dequant(i, 7); | |
171 } | |
172 for(i=0; i<15; i++) { | |
173 /* bap=5 mantissas */ | |
174 b5_mantissas[i] = symmetric_dequant(i, 15); | |
175 } | |
5321
fcb2a52c88fe
AC-3 decoder, soc revision 71, Sep 4 09:47:12 2006 UTC by cloud9
jbr
parents:
5320
diff
changeset
|
176 |
5457
0a5849a4858b
use table for converting dynamic range codes to scale factors
jbr
parents:
5456
diff
changeset
|
177 /* generate dynamic range table |
0a5849a4858b
use table for converting dynamic range codes to scale factors
jbr
parents:
5456
diff
changeset
|
178 reference: Section 7.7.1 Dynamic Range Control */ |
0a5849a4858b
use table for converting dynamic range codes to scale factors
jbr
parents:
5456
diff
changeset
|
179 for(i=0; i<256; i++) { |
0a5849a4858b
use table for converting dynamic range codes to scale factors
jbr
parents:
5456
diff
changeset
|
180 int v = (i >> 5) - ((i >> 7) << 3) - 5; |
6006 | 181 dynamic_range_tab[i] = powf(2.0f, v) * ((i & 0x1F) | 0x20); |
5457
0a5849a4858b
use table for converting dynamic range codes to scale factors
jbr
parents:
5456
diff
changeset
|
182 } |
5305
5892b4a6380b
AC-3 decoder, soc revision 31, Jul 14 23:53:28 2006 UTC by cloud9
jbr
parents:
5304
diff
changeset
|
183 } |
5892b4a6380b
AC-3 decoder, soc revision 31, Jul 14 23:53:28 2006 UTC by cloud9
jbr
parents:
5304
diff
changeset
|
184 |
5304
eff51058fe13
AC-3 decoder, soc revision 26, Jul 5 04:55:15 2006 UTC by cloud9
jbr
parents:
5303
diff
changeset
|
185 |
5515
b499866a3397
comments modification: add, clean up, Doxygenize, and reword
jbr
parents:
5499
diff
changeset
|
186 /** |
b499866a3397
comments modification: add, clean up, Doxygenize, and reword
jbr
parents:
5499
diff
changeset
|
187 * AVCodec initialization |
b499866a3397
comments modification: add, clean up, Doxygenize, and reword
jbr
parents:
5499
diff
changeset
|
188 */ |
6517
48759bfbd073
Apply 'cold' attribute to init/uninit functions in libavcodec
zuxy
parents:
6487
diff
changeset
|
189 static av_cold int ac3_decode_init(AVCodecContext *avctx) |
5304
eff51058fe13
AC-3 decoder, soc revision 26, Jul 5 04:55:15 2006 UTC by cloud9
jbr
parents:
5303
diff
changeset
|
190 { |
6090
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
191 AC3DecodeContext *s = avctx->priv_data; |
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
192 s->avctx = avctx; |
5304
eff51058fe13
AC-3 decoder, soc revision 26, Jul 5 04:55:15 2006 UTC by cloud9
jbr
parents:
5303
diff
changeset
|
193 |
5331 | 194 ac3_common_init(); |
5309
0662a270aab7
AC-3 decoder, soc revision 38, Aug 7 00:03:00 2006 UTC by cloud9
jbr
parents:
5308
diff
changeset
|
195 ac3_tables_init(); |
6090
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
196 ff_mdct_init(&s->imdct_256, 8, 1); |
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
197 ff_mdct_init(&s->imdct_512, 9, 1); |
6142
a35b838ab955
Add variable alpha and size of half window for Kaiser-Bessel Derived window
superdump
parents:
6139
diff
changeset
|
198 ff_kbd_window_init(s->window, 5.0, 256); |
6090
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
199 dsputil_init(&s->dsp, avctx); |
7601 | 200 av_lfg_init(&s->dith_state, 0); |
5304
eff51058fe13
AC-3 decoder, soc revision 26, Jul 5 04:55:15 2006 UTC by cloud9
jbr
parents:
5303
diff
changeset
|
201 |
5515
b499866a3397
comments modification: add, clean up, Doxygenize, and reword
jbr
parents:
5499
diff
changeset
|
202 /* set bias values for float to int16 conversion */ |
7550 | 203 if(s->dsp.float_to_int16_interleave == ff_float_to_int16_interleave_c) { |
6090
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
204 s->add_bias = 385.0f; |
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
205 s->mul_bias = 1.0f; |
5479
943c732c905d
use dsputil for float to signed 16-bit sample conversion
jbr
parents:
5477
diff
changeset
|
206 } else { |
6090
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
207 s->add_bias = 0.0f; |
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
208 s->mul_bias = 32767.0f; |
5479
943c732c905d
use dsputil for float to signed 16-bit sample conversion
jbr
parents:
5477
diff
changeset
|
209 } |
943c732c905d
use dsputil for float to signed 16-bit sample conversion
jbr
parents:
5477
diff
changeset
|
210 |
6112 | 211 /* allow downmixing to stereo or mono */ |
212 if (avctx->channels > 0 && avctx->request_channels > 0 && | |
213 avctx->request_channels < avctx->channels && | |
214 avctx->request_channels <= 2) { | |
215 avctx->channels = avctx->request_channels; | |
216 } | |
6480
6f01a499e785
downmix before imdct unless different size transforms are used. about 20%
jbr
parents:
6432
diff
changeset
|
217 s->downmixed = 1; |
6112 | 218 |
6538
b8f886bce6f5
only allocate context input buffer if AVCodecContext.error_reslience is greater than 0.
jbr
parents:
6537
diff
changeset
|
219 /* allocate context input buffer */ |
b8f886bce6f5
only allocate context input buffer if AVCodecContext.error_reslience is greater than 0.
jbr
parents:
6537
diff
changeset
|
220 if (avctx->error_resilience >= FF_ER_CAREFUL) { |
b8f886bce6f5
only allocate context input buffer if AVCodecContext.error_reslience is greater than 0.
jbr
parents:
6537
diff
changeset
|
221 s->input_buffer = av_mallocz(AC3_MAX_FRAME_SIZE + FF_INPUT_BUFFER_PADDING_SIZE); |
b8f886bce6f5
only allocate context input buffer if AVCodecContext.error_reslience is greater than 0.
jbr
parents:
6537
diff
changeset
|
222 if (!s->input_buffer) |
b8f886bce6f5
only allocate context input buffer if AVCodecContext.error_reslience is greater than 0.
jbr
parents:
6537
diff
changeset
|
223 return AVERROR_NOMEM; |
b8f886bce6f5
only allocate context input buffer if AVCodecContext.error_reslience is greater than 0.
jbr
parents:
6537
diff
changeset
|
224 } |
b8f886bce6f5
only allocate context input buffer if AVCodecContext.error_reslience is greater than 0.
jbr
parents:
6537
diff
changeset
|
225 |
7451
85ab7655ad4d
Modify all codecs to report their supported input and output sample format(s).
pross
parents:
7328
diff
changeset
|
226 avctx->sample_fmt = SAMPLE_FMT_S16; |
5304
eff51058fe13
AC-3 decoder, soc revision 26, Jul 5 04:55:15 2006 UTC by cloud9
jbr
parents:
5303
diff
changeset
|
227 return 0; |
eff51058fe13
AC-3 decoder, soc revision 26, Jul 5 04:55:15 2006 UTC by cloud9
jbr
parents:
5303
diff
changeset
|
228 } |
eff51058fe13
AC-3 decoder, soc revision 26, Jul 5 04:55:15 2006 UTC by cloud9
jbr
parents:
5303
diff
changeset
|
229 |
5386
d6347c3cb329
Use shared AC-3 header parsing. Move some initializations to block parsing function.
jbr
parents:
5385
diff
changeset
|
230 /** |
5515
b499866a3397
comments modification: add, clean up, Doxygenize, and reword
jbr
parents:
5499
diff
changeset
|
231 * Parse the 'sync info' and 'bit stream info' from the AC-3 bitstream. |
5321
fcb2a52c88fe
AC-3 decoder, soc revision 71, Sep 4 09:47:12 2006 UTC by cloud9
jbr
parents:
5320
diff
changeset
|
232 * GetBitContext within AC3DecodeContext must point to |
7470
1a93d3bbe3ee
cosmetics: make all references to AC-3 capitalized and hyphenated
jbr
parents:
7451
diff
changeset
|
233 * the start of the synchronized AC-3 bitstream. |
5321
fcb2a52c88fe
AC-3 decoder, soc revision 71, Sep 4 09:47:12 2006 UTC by cloud9
jbr
parents:
5320
diff
changeset
|
234 */ |
6090
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
235 static int ac3_parse_header(AC3DecodeContext *s) |
5303
b8821ff5c30d
AC-3 decoder, soc revision 9, Jun 14 03:52:02 2006 UTC by cloud9
jbr
parents:
diff
changeset
|
236 { |
7028 | 237 GetBitContext *gbc = &s->gbc; |
238 int i; | |
239 | |
240 /* read the rest of the bsi. read twice for dual mono mode. */ | |
241 i = !(s->channel_mode); | |
242 do { | |
243 skip_bits(gbc, 5); // skip dialog normalization | |
244 if (get_bits1(gbc)) | |
245 skip_bits(gbc, 8); //skip compression | |
246 if (get_bits1(gbc)) | |
247 skip_bits(gbc, 8); //skip language code | |
248 if (get_bits1(gbc)) | |
249 skip_bits(gbc, 7); //skip audio production information | |
250 } while (i--); | |
251 | |
252 skip_bits(gbc, 2); //skip copyright bit and original bitstream bit | |
253 | |
254 /* skip the timecodes (or extra bitstream information for Alternate Syntax) | |
255 TODO: read & use the xbsi1 downmix levels */ | |
256 if (get_bits1(gbc)) | |
257 skip_bits(gbc, 14); //skip timecode1 / xbsi1 | |
258 if (get_bits1(gbc)) | |
259 skip_bits(gbc, 14); //skip timecode2 / xbsi2 | |
260 | |
261 /* skip additional bitstream info */ | |
262 if (get_bits1(gbc)) { | |
263 i = get_bits(gbc, 6); | |
264 do { | |
265 skip_bits(gbc, 8); | |
266 } while(i--); | |
267 } | |
268 | |
269 return 0; | |
270 } | |
271 | |
272 /** | |
7470
1a93d3bbe3ee
cosmetics: make all references to AC-3 capitalized and hyphenated
jbr
parents:
7451
diff
changeset
|
273 * Common function to parse AC-3 or E-AC-3 frame header |
7028 | 274 */ |
275 static int parse_frame_header(AC3DecodeContext *s) | |
276 { | |
5386
d6347c3cb329
Use shared AC-3 header parsing. Move some initializations to block parsing function.
jbr
parents:
5385
diff
changeset
|
277 AC3HeaderInfo hdr; |
7028 | 278 int err; |
5309
0662a270aab7
AC-3 decoder, soc revision 38, Aug 7 00:03:00 2006 UTC by cloud9
jbr
parents:
5308
diff
changeset
|
279 |
7489 | 280 err = ff_ac3_parse_header(&s->gbc, &hdr); |
5386
d6347c3cb329
Use shared AC-3 header parsing. Move some initializations to block parsing function.
jbr
parents:
5385
diff
changeset
|
281 if(err) |
d6347c3cb329
Use shared AC-3 header parsing. Move some initializations to block parsing function.
jbr
parents:
5385
diff
changeset
|
282 return err; |
5305
5892b4a6380b
AC-3 decoder, soc revision 31, Jul 14 23:53:28 2006 UTC by cloud9
jbr
parents:
5304
diff
changeset
|
283 |
5386
d6347c3cb329
Use shared AC-3 header parsing. Move some initializations to block parsing function.
jbr
parents:
5385
diff
changeset
|
284 /* get decoding parameters from header info */ |
6090
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
285 s->bit_alloc_params.sr_code = hdr.sr_code; |
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
286 s->channel_mode = hdr.channel_mode; |
6091 | 287 s->lfe_on = hdr.lfe_on; |
6090
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
288 s->bit_alloc_params.sr_shift = hdr.sr_shift; |
6093 | 289 s->sample_rate = hdr.sample_rate; |
6090
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
290 s->bit_rate = hdr.bit_rate; |
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
291 s->channels = hdr.channels; |
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
292 s->fbw_channels = s->channels - s->lfe_on; |
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
293 s->lfe_ch = s->fbw_channels + 1; |
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
294 s->frame_size = hdr.frame_size; |
6661
a409fbf1f42b
change ff_ac3_parse_header() to take a GetBitContext instead of const char*
bwolowiec
parents:
6540
diff
changeset
|
295 s->center_mix_level = hdr.center_mix_level; |
a409fbf1f42b
change ff_ac3_parse_header() to take a GetBitContext instead of const char*
bwolowiec
parents:
6540
diff
changeset
|
296 s->surround_mix_level = hdr.surround_mix_level; |
7011
77f607fb4e8b
get the number of blocks from the ac3 parser and use in the ac3 decoder.
jbr
parents:
6999
diff
changeset
|
297 s->num_blocks = hdr.num_blocks; |
7013 | 298 s->frame_type = hdr.frame_type; |
7012 | 299 s->substreamid = hdr.substreamid; |
5412
3f6a6dda6b2d
remove broken downmixing. will add new implementation later.
jbr
parents:
5393
diff
changeset
|
300 |
6958
1c523b5212cd
set lfe params at frame init instead of in every block
jbr
parents:
6957
diff
changeset
|
301 if(s->lfe_on) { |
1c523b5212cd
set lfe params at frame init instead of in every block
jbr
parents:
6957
diff
changeset
|
302 s->start_freq[s->lfe_ch] = 0; |
1c523b5212cd
set lfe params at frame init instead of in every block
jbr
parents:
6957
diff
changeset
|
303 s->end_freq[s->lfe_ch] = 7; |
1c523b5212cd
set lfe params at frame init instead of in every block
jbr
parents:
6957
diff
changeset
|
304 s->num_exp_groups[s->lfe_ch] = 2; |
1c523b5212cd
set lfe params at frame init instead of in every block
jbr
parents:
6957
diff
changeset
|
305 s->channel_in_cpl[s->lfe_ch] = 0; |
1c523b5212cd
set lfe params at frame init instead of in every block
jbr
parents:
6957
diff
changeset
|
306 } |
1c523b5212cd
set lfe params at frame init instead of in every block
jbr
parents:
6957
diff
changeset
|
307 |
7473
e2f3016867a0
Read bitstream parameters for E-AC-3 streams before returning an error.
jbr
parents:
7472
diff
changeset
|
308 if(hdr.bitstream_id > 10) |
e2f3016867a0
Read bitstream parameters for E-AC-3 streams before returning an error.
jbr
parents:
7472
diff
changeset
|
309 return AC3_PARSE_ERROR_BSID; |
e2f3016867a0
Read bitstream parameters for E-AC-3 streams before returning an error.
jbr
parents:
7472
diff
changeset
|
310 |
7028 | 311 return ac3_parse_header(s); |
6126
c24cd5282108
split out init of downmix coeffs into a separate function that can be called
jbr
parents:
6125
diff
changeset
|
312 } |
c24cd5282108
split out init of downmix coeffs into a separate function that can be called
jbr
parents:
6125
diff
changeset
|
313 |
c24cd5282108
split out init of downmix coeffs into a separate function that can be called
jbr
parents:
6125
diff
changeset
|
314 /** |
c24cd5282108
split out init of downmix coeffs into a separate function that can be called
jbr
parents:
6125
diff
changeset
|
315 * Set stereo downmixing coefficients based on frame header info. |
c24cd5282108
split out init of downmix coeffs into a separate function that can be called
jbr
parents:
6125
diff
changeset
|
316 * reference: Section 7.8.2 Downmixing Into Two Channels |
c24cd5282108
split out init of downmix coeffs into a separate function that can be called
jbr
parents:
6125
diff
changeset
|
317 */ |
c24cd5282108
split out init of downmix coeffs into a separate function that can be called
jbr
parents:
6125
diff
changeset
|
318 static void set_downmix_coeffs(AC3DecodeContext *s) |
c24cd5282108
split out init of downmix coeffs into a separate function that can be called
jbr
parents:
6125
diff
changeset
|
319 { |
c24cd5282108
split out init of downmix coeffs into a separate function that can be called
jbr
parents:
6125
diff
changeset
|
320 int i; |
7019
81d5c68233e5
move mix level tables from parser to decoder. have parser read bitstream value instead of using an index to a table in the decoder.
jbr
parents:
7018
diff
changeset
|
321 float cmix = gain_levels[center_levels[s->center_mix_level]]; |
81d5c68233e5
move mix level tables from parser to decoder. have parser read bitstream value instead of using an index to a table in the decoder.
jbr
parents:
7018
diff
changeset
|
322 float smix = gain_levels[surround_levels[s->surround_mix_level]]; |
7554 | 323 float norm0, norm1; |
6126
c24cd5282108
split out init of downmix coeffs into a separate function that can be called
jbr
parents:
6125
diff
changeset
|
324 |
6090
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
325 for(i=0; i<s->fbw_channels; i++) { |
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
326 s->downmix_coeffs[i][0] = gain_levels[ac3_default_coeffs[s->channel_mode][i][0]]; |
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
327 s->downmix_coeffs[i][1] = gain_levels[ac3_default_coeffs[s->channel_mode][i][1]]; |
5480 | 328 } |
6090
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
329 if(s->channel_mode > 1 && s->channel_mode & 1) { |
6126
c24cd5282108
split out init of downmix coeffs into a separate function that can be called
jbr
parents:
6125
diff
changeset
|
330 s->downmix_coeffs[1][0] = s->downmix_coeffs[1][1] = cmix; |
5480 | 331 } |
6090
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
332 if(s->channel_mode == AC3_CHMODE_2F1R || s->channel_mode == AC3_CHMODE_3F1R) { |
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
333 int nf = s->channel_mode - 2; |
6126
c24cd5282108
split out init of downmix coeffs into a separate function that can be called
jbr
parents:
6125
diff
changeset
|
334 s->downmix_coeffs[nf][0] = s->downmix_coeffs[nf][1] = smix * LEVEL_MINUS_3DB; |
5480 | 335 } |
6090
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
336 if(s->channel_mode == AC3_CHMODE_2F2R || s->channel_mode == AC3_CHMODE_3F2R) { |
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
337 int nf = s->channel_mode - 4; |
6126
c24cd5282108
split out init of downmix coeffs into a separate function that can be called
jbr
parents:
6125
diff
changeset
|
338 s->downmix_coeffs[nf][0] = s->downmix_coeffs[nf+1][1] = smix; |
5480 | 339 } |
6431
305c49259b59
pre-calculate downmix coefficient sums. 2% faster 5.1-to-stereo downmixing.
jbr
parents:
6430
diff
changeset
|
340 |
7554 | 341 /* renormalize */ |
342 norm0 = norm1 = 0.0; | |
343 for(i=0; i<s->fbw_channels; i++) { | |
344 norm0 += s->downmix_coeffs[i][0]; | |
345 norm1 += s->downmix_coeffs[i][1]; | |
346 } | |
347 norm0 = 1.0f / norm0; | |
348 norm1 = 1.0f / norm1; | |
6431
305c49259b59
pre-calculate downmix coefficient sums. 2% faster 5.1-to-stereo downmixing.
jbr
parents:
6430
diff
changeset
|
349 for(i=0; i<s->fbw_channels; i++) { |
7554 | 350 s->downmix_coeffs[i][0] *= norm0; |
351 s->downmix_coeffs[i][1] *= norm1; | |
6431
305c49259b59
pre-calculate downmix coefficient sums. 2% faster 5.1-to-stereo downmixing.
jbr
parents:
6430
diff
changeset
|
352 } |
7554 | 353 |
354 if(s->output_mode == AC3_CHMODE_MONO) { | |
355 for(i=0; i<s->fbw_channels; i++) | |
356 s->downmix_coeffs[i][0] = (s->downmix_coeffs[i][0] + s->downmix_coeffs[i][1]) * LEVEL_MINUS_3DB; | |
357 } | |
5303
b8821ff5c30d
AC-3 decoder, soc revision 9, Jun 14 03:52:02 2006 UTC by cloud9
jbr
parents:
diff
changeset
|
358 } |
b8821ff5c30d
AC-3 decoder, soc revision 9, Jun 14 03:52:02 2006 UTC by cloud9
jbr
parents:
diff
changeset
|
359 |
5376 | 360 /** |
5515
b499866a3397
comments modification: add, clean up, Doxygenize, and reword
jbr
parents:
5499
diff
changeset
|
361 * Decode the grouped exponents according to exponent strategy. |
b499866a3397
comments modification: add, clean up, Doxygenize, and reword
jbr
parents:
5499
diff
changeset
|
362 * reference: Section 7.1.3 Exponent Decoding |
5304
eff51058fe13
AC-3 decoder, soc revision 26, Jul 5 04:55:15 2006 UTC by cloud9
jbr
parents:
5303
diff
changeset
|
363 */ |
6089 | 364 static void decode_exponents(GetBitContext *gbc, int exp_strategy, int ngrps, |
5389
cf99cf2d3af9
use consistent data type for exponents. needs to be int8_t to pass to bit allocation.
jbr
parents:
5388
diff
changeset
|
365 uint8_t absexp, int8_t *dexps) |
5304
eff51058fe13
AC-3 decoder, soc revision 26, Jul 5 04:55:15 2006 UTC by cloud9
jbr
parents:
5303
diff
changeset
|
366 { |
6006 | 367 int i, j, grp, group_size; |
5375 | 368 int dexp[256]; |
369 int expacc, prevexp; | |
5305
5892b4a6380b
AC-3 decoder, soc revision 31, Jul 14 23:53:28 2006 UTC by cloud9
jbr
parents:
5304
diff
changeset
|
370 |
5375 | 371 /* unpack groups */ |
6006 | 372 group_size = exp_strategy + (exp_strategy == EXP_D45); |
5375 | 373 for(grp=0,i=0; grp<ngrps; grp++) { |
6089 | 374 expacc = get_bits(gbc, 7); |
7471
1788013cf23b
use the same table for ungrouping exponents and bap=2 mantissas
jbr
parents:
7470
diff
changeset
|
375 dexp[i++] = ungroup_3_in_7_bits_tab[expacc][0]; |
1788013cf23b
use the same table for ungrouping exponents and bap=2 mantissas
jbr
parents:
7470
diff
changeset
|
376 dexp[i++] = ungroup_3_in_7_bits_tab[expacc][1]; |
1788013cf23b
use the same table for ungrouping exponents and bap=2 mantissas
jbr
parents:
7470
diff
changeset
|
377 dexp[i++] = ungroup_3_in_7_bits_tab[expacc][2]; |
5375 | 378 } |
5303
b8821ff5c30d
AC-3 decoder, soc revision 9, Jun 14 03:52:02 2006 UTC by cloud9
jbr
parents:
diff
changeset
|
379 |
5375 | 380 /* convert to absolute exps and expand groups */ |
381 prevexp = absexp; | |
382 for(i=0; i<ngrps*3; i++) { | |
383 prevexp = av_clip(prevexp + dexp[i]-2, 0, 24); | |
6006 | 384 for(j=0; j<group_size; j++) { |
385 dexps[(i*group_size)+j] = prevexp; | |
5304
eff51058fe13
AC-3 decoder, soc revision 26, Jul 5 04:55:15 2006 UTC by cloud9
jbr
parents:
5303
diff
changeset
|
386 } |
eff51058fe13
AC-3 decoder, soc revision 26, Jul 5 04:55:15 2006 UTC by cloud9
jbr
parents:
5303
diff
changeset
|
387 } |
5303
b8821ff5c30d
AC-3 decoder, soc revision 9, Jun 14 03:52:02 2006 UTC by cloud9
jbr
parents:
diff
changeset
|
388 } |
b8821ff5c30d
AC-3 decoder, soc revision 9, Jun 14 03:52:02 2006 UTC by cloud9
jbr
parents:
diff
changeset
|
389 |
5451
19bb5e64a57c
split out coupling coefficient reconstruction into a separate function
jbr
parents:
5435
diff
changeset
|
390 /** |
5515
b499866a3397
comments modification: add, clean up, Doxygenize, and reword
jbr
parents:
5499
diff
changeset
|
391 * Generate transform coefficients for each coupled channel in the coupling |
5451
19bb5e64a57c
split out coupling coefficient reconstruction into a separate function
jbr
parents:
5435
diff
changeset
|
392 * range using the coupling coefficients and coupling coordinates. |
19bb5e64a57c
split out coupling coefficient reconstruction into a separate function
jbr
parents:
5435
diff
changeset
|
393 * reference: Section 7.4.3 Coupling Coordinate Format |
19bb5e64a57c
split out coupling coefficient reconstruction into a separate function
jbr
parents:
5435
diff
changeset
|
394 */ |
7491
0dbdc638270e
cosmetics: use a better function name than uncouple_channels()
jbr
parents:
7490
diff
changeset
|
395 static void calc_transform_coeffs_cpl(AC3DecodeContext *s) |
5451
19bb5e64a57c
split out coupling coefficient reconstruction into a separate function
jbr
parents:
5435
diff
changeset
|
396 { |
19bb5e64a57c
split out coupling coefficient reconstruction into a separate function
jbr
parents:
5435
diff
changeset
|
397 int i, j, ch, bnd, subbnd; |
19bb5e64a57c
split out coupling coefficient reconstruction into a separate function
jbr
parents:
5435
diff
changeset
|
398 |
19bb5e64a57c
split out coupling coefficient reconstruction into a separate function
jbr
parents:
5435
diff
changeset
|
399 subbnd = -1; |
6090
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
400 i = s->start_freq[CPL_CH]; |
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
401 for(bnd=0; bnd<s->num_cpl_bands; bnd++) { |
5451
19bb5e64a57c
split out coupling coefficient reconstruction into a separate function
jbr
parents:
5435
diff
changeset
|
402 do { |
19bb5e64a57c
split out coupling coefficient reconstruction into a separate function
jbr
parents:
5435
diff
changeset
|
403 subbnd++; |
19bb5e64a57c
split out coupling coefficient reconstruction into a separate function
jbr
parents:
5435
diff
changeset
|
404 for(j=0; j<12; j++) { |
6090
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
405 for(ch=1; ch<=s->fbw_channels; ch++) { |
6123
8d88831a2777
fix phase flag processing for the case when coupling coordinates are reused
jbr
parents:
6117
diff
changeset
|
406 if(s->channel_in_cpl[ch]) { |
6430
de7502093922
Use 24-bit fixed-point transform coefficients until just before MDCT. This
jbr
parents:
6142
diff
changeset
|
407 s->fixed_coeffs[ch][i] = ((int64_t)s->fixed_coeffs[CPL_CH][i] * (int64_t)s->cpl_coords[ch][bnd]) >> 23; |
6123
8d88831a2777
fix phase flag processing for the case when coupling coordinates are reused
jbr
parents:
6117
diff
changeset
|
408 if (ch == 2 && s->phase_flags[bnd]) |
6430
de7502093922
Use 24-bit fixed-point transform coefficients until just before MDCT. This
jbr
parents:
6142
diff
changeset
|
409 s->fixed_coeffs[ch][i] = -s->fixed_coeffs[ch][i]; |
6123
8d88831a2777
fix phase flag processing for the case when coupling coordinates are reused
jbr
parents:
6117
diff
changeset
|
410 } |
5451
19bb5e64a57c
split out coupling coefficient reconstruction into a separate function
jbr
parents:
5435
diff
changeset
|
411 } |
19bb5e64a57c
split out coupling coefficient reconstruction into a separate function
jbr
parents:
5435
diff
changeset
|
412 i++; |
19bb5e64a57c
split out coupling coefficient reconstruction into a separate function
jbr
parents:
5435
diff
changeset
|
413 } |
6090
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
414 } while(s->cpl_band_struct[subbnd]); |
5451
19bb5e64a57c
split out coupling coefficient reconstruction into a separate function
jbr
parents:
5435
diff
changeset
|
415 } |
19bb5e64a57c
split out coupling coefficient reconstruction into a separate function
jbr
parents:
5435
diff
changeset
|
416 } |
19bb5e64a57c
split out coupling coefficient reconstruction into a separate function
jbr
parents:
5435
diff
changeset
|
417 |
5515
b499866a3397
comments modification: add, clean up, Doxygenize, and reword
jbr
parents:
5499
diff
changeset
|
418 /** |
b499866a3397
comments modification: add, clean up, Doxygenize, and reword
jbr
parents:
5499
diff
changeset
|
419 * Grouped mantissas for 3-level 5-level and 11-level quantization |
b499866a3397
comments modification: add, clean up, Doxygenize, and reword
jbr
parents:
5499
diff
changeset
|
420 */ |
b499866a3397
comments modification: add, clean up, Doxygenize, and reword
jbr
parents:
5499
diff
changeset
|
421 typedef struct { |
6430
de7502093922
Use 24-bit fixed-point transform coefficients until just before MDCT. This
jbr
parents:
6142
diff
changeset
|
422 int b1_mant[3]; |
de7502093922
Use 24-bit fixed-point transform coefficients until just before MDCT. This
jbr
parents:
6142
diff
changeset
|
423 int b2_mant[3]; |
de7502093922
Use 24-bit fixed-point transform coefficients until just before MDCT. This
jbr
parents:
6142
diff
changeset
|
424 int b4_mant[2]; |
5456 | 425 int b1ptr; |
426 int b2ptr; | |
427 int b4ptr; | |
5306
abc5c130b448
AC-3 decoder, soc revision 32, Jul 17 09:37:32 2006 UTC by cloud9
jbr
parents:
5305
diff
changeset
|
428 } mant_groups; |
abc5c130b448
AC-3 decoder, soc revision 32, Jul 17 09:37:32 2006 UTC by cloud9
jbr
parents:
5305
diff
changeset
|
429 |
5515
b499866a3397
comments modification: add, clean up, Doxygenize, and reword
jbr
parents:
5499
diff
changeset
|
430 /** |
b499866a3397
comments modification: add, clean up, Doxygenize, and reword
jbr
parents:
5499
diff
changeset
|
431 * Get the transform coefficients for a particular channel |
b499866a3397
comments modification: add, clean up, Doxygenize, and reword
jbr
parents:
5499
diff
changeset
|
432 * reference: Section 7.3 Quantization and Decoding of Mantissas |
b499866a3397
comments modification: add, clean up, Doxygenize, and reword
jbr
parents:
5499
diff
changeset
|
433 */ |
6953
64504967b891
get_transform_coeffs() never returns an error, so make the function and its children return void.
jbr
parents:
6951
diff
changeset
|
434 static void get_transform_coeffs_ch(AC3DecodeContext *s, int ch_index, mant_groups *m) |
5303
b8821ff5c30d
AC-3 decoder, soc revision 9, Jun 14 03:52:02 2006 UTC by cloud9
jbr
parents:
diff
changeset
|
435 { |
6090
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
436 GetBitContext *gbc = &s->gbc; |
5454
5b497d971546
dither zero-bit mantissas by default. remove dithering only if it's explicitly turned off.
jbr
parents:
5453
diff
changeset
|
437 int i, gcode, tbap, start, end; |
5310
9aa9197034d7
AC-3 decoder, soc revision 40, Aug 9 00:10:14 2006 UTC by cloud9
jbr
parents:
5309
diff
changeset
|
438 uint8_t *exps; |
9aa9197034d7
AC-3 decoder, soc revision 40, Aug 9 00:10:14 2006 UTC by cloud9
jbr
parents:
5309
diff
changeset
|
439 uint8_t *bap; |
6430
de7502093922
Use 24-bit fixed-point transform coefficients until just before MDCT. This
jbr
parents:
6142
diff
changeset
|
440 int *coeffs; |
5309
0662a270aab7
AC-3 decoder, soc revision 38, Aug 7 00:03:00 2006 UTC by cloud9
jbr
parents:
5308
diff
changeset
|
441 |
6090
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
442 exps = s->dexps[ch_index]; |
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
443 bap = s->bap[ch_index]; |
6430
de7502093922
Use 24-bit fixed-point transform coefficients until just before MDCT. This
jbr
parents:
6142
diff
changeset
|
444 coeffs = s->fixed_coeffs[ch_index]; |
6090
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
445 start = s->start_freq[ch_index]; |
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
446 end = s->end_freq[ch_index]; |
5303
b8821ff5c30d
AC-3 decoder, soc revision 9, Jun 14 03:52:02 2006 UTC by cloud9
jbr
parents:
diff
changeset
|
447 |
5452
190d793a22e7
merge decoding of coupling channel coefficients and regular channel coefficients
jbr
parents:
5451
diff
changeset
|
448 for (i = start; i < end; i++) { |
5309
0662a270aab7
AC-3 decoder, soc revision 38, Aug 7 00:03:00 2006 UTC by cloud9
jbr
parents:
5308
diff
changeset
|
449 tbap = bap[i]; |
0662a270aab7
AC-3 decoder, soc revision 38, Aug 7 00:03:00 2006 UTC by cloud9
jbr
parents:
5308
diff
changeset
|
450 switch (tbap) { |
5304
eff51058fe13
AC-3 decoder, soc revision 26, Jul 5 04:55:15 2006 UTC by cloud9
jbr
parents:
5303
diff
changeset
|
451 case 0: |
7601 | 452 coeffs[i] = (av_lfg_get(&s->dith_state) & 0x7FFFFF) - 0x400000; |
5453 | 453 break; |
5303
b8821ff5c30d
AC-3 decoder, soc revision 9, Jun 14 03:52:02 2006 UTC by cloud9
jbr
parents:
diff
changeset
|
454 |
5304
eff51058fe13
AC-3 decoder, soc revision 26, Jul 5 04:55:15 2006 UTC by cloud9
jbr
parents:
5303
diff
changeset
|
455 case 1: |
5456 | 456 if(m->b1ptr > 2) { |
6089 | 457 gcode = get_bits(gbc, 5); |
5456 | 458 m->b1_mant[0] = b1_mantissas[gcode][0]; |
459 m->b1_mant[1] = b1_mantissas[gcode][1]; | |
460 m->b1_mant[2] = b1_mantissas[gcode][2]; | |
461 m->b1ptr = 0; | |
5304
eff51058fe13
AC-3 decoder, soc revision 26, Jul 5 04:55:15 2006 UTC by cloud9
jbr
parents:
5303
diff
changeset
|
462 } |
5456 | 463 coeffs[i] = m->b1_mant[m->b1ptr++]; |
5453 | 464 break; |
5303
b8821ff5c30d
AC-3 decoder, soc revision 9, Jun 14 03:52:02 2006 UTC by cloud9
jbr
parents:
diff
changeset
|
465 |
5304
eff51058fe13
AC-3 decoder, soc revision 26, Jul 5 04:55:15 2006 UTC by cloud9
jbr
parents:
5303
diff
changeset
|
466 case 2: |
5456 | 467 if(m->b2ptr > 2) { |
6089 | 468 gcode = get_bits(gbc, 7); |
5456 | 469 m->b2_mant[0] = b2_mantissas[gcode][0]; |
470 m->b2_mant[1] = b2_mantissas[gcode][1]; | |
471 m->b2_mant[2] = b2_mantissas[gcode][2]; | |
472 m->b2ptr = 0; | |
5304
eff51058fe13
AC-3 decoder, soc revision 26, Jul 5 04:55:15 2006 UTC by cloud9
jbr
parents:
5303
diff
changeset
|
473 } |
5456 | 474 coeffs[i] = m->b2_mant[m->b2ptr++]; |
5453 | 475 break; |
5303
b8821ff5c30d
AC-3 decoder, soc revision 9, Jun 14 03:52:02 2006 UTC by cloud9
jbr
parents:
diff
changeset
|
476 |
5304
eff51058fe13
AC-3 decoder, soc revision 26, Jul 5 04:55:15 2006 UTC by cloud9
jbr
parents:
5303
diff
changeset
|
477 case 3: |
6089 | 478 coeffs[i] = b3_mantissas[get_bits(gbc, 3)]; |
5453 | 479 break; |
5303
b8821ff5c30d
AC-3 decoder, soc revision 9, Jun 14 03:52:02 2006 UTC by cloud9
jbr
parents:
diff
changeset
|
480 |
5304
eff51058fe13
AC-3 decoder, soc revision 26, Jul 5 04:55:15 2006 UTC by cloud9
jbr
parents:
5303
diff
changeset
|
481 case 4: |
5456 | 482 if(m->b4ptr > 1) { |
6089 | 483 gcode = get_bits(gbc, 7); |
5456 | 484 m->b4_mant[0] = b4_mantissas[gcode][0]; |
485 m->b4_mant[1] = b4_mantissas[gcode][1]; | |
486 m->b4ptr = 0; | |
5304
eff51058fe13
AC-3 decoder, soc revision 26, Jul 5 04:55:15 2006 UTC by cloud9
jbr
parents:
5303
diff
changeset
|
487 } |
5456 | 488 coeffs[i] = m->b4_mant[m->b4ptr++]; |
5453 | 489 break; |
5303
b8821ff5c30d
AC-3 decoder, soc revision 9, Jun 14 03:52:02 2006 UTC by cloud9
jbr
parents:
diff
changeset
|
490 |
5304
eff51058fe13
AC-3 decoder, soc revision 26, Jul 5 04:55:15 2006 UTC by cloud9
jbr
parents:
5303
diff
changeset
|
491 case 5: |
6089 | 492 coeffs[i] = b5_mantissas[get_bits(gbc, 4)]; |
5453 | 493 break; |
5303
b8821ff5c30d
AC-3 decoder, soc revision 9, Jun 14 03:52:02 2006 UTC by cloud9
jbr
parents:
diff
changeset
|
494 |
6430
de7502093922
Use 24-bit fixed-point transform coefficients until just before MDCT. This
jbr
parents:
6142
diff
changeset
|
495 default: { |
5515
b499866a3397
comments modification: add, clean up, Doxygenize, and reword
jbr
parents:
5499
diff
changeset
|
496 /* asymmetric dequantization */ |
6430
de7502093922
Use 24-bit fixed-point transform coefficients until just before MDCT. This
jbr
parents:
6142
diff
changeset
|
497 int qlevel = quantization_tab[tbap]; |
de7502093922
Use 24-bit fixed-point transform coefficients until just before MDCT. This
jbr
parents:
6142
diff
changeset
|
498 coeffs[i] = get_sbits(gbc, qlevel) << (24 - qlevel); |
5453 | 499 break; |
6430
de7502093922
Use 24-bit fixed-point transform coefficients until just before MDCT. This
jbr
parents:
6142
diff
changeset
|
500 } |
5304
eff51058fe13
AC-3 decoder, soc revision 26, Jul 5 04:55:15 2006 UTC by cloud9
jbr
parents:
5303
diff
changeset
|
501 } |
6430
de7502093922
Use 24-bit fixed-point transform coefficients until just before MDCT. This
jbr
parents:
6142
diff
changeset
|
502 coeffs[i] >>= exps[i]; |
5303
b8821ff5c30d
AC-3 decoder, soc revision 9, Jun 14 03:52:02 2006 UTC by cloud9
jbr
parents:
diff
changeset
|
503 } |
b8821ff5c30d
AC-3 decoder, soc revision 9, Jun 14 03:52:02 2006 UTC by cloud9
jbr
parents:
diff
changeset
|
504 } |
b8821ff5c30d
AC-3 decoder, soc revision 9, Jun 14 03:52:02 2006 UTC by cloud9
jbr
parents:
diff
changeset
|
505 |
5454
5b497d971546
dither zero-bit mantissas by default. remove dithering only if it's explicitly turned off.
jbr
parents:
5453
diff
changeset
|
506 /** |
5515
b499866a3397
comments modification: add, clean up, Doxygenize, and reword
jbr
parents:
5499
diff
changeset
|
507 * Remove random dithering from coefficients with zero-bit mantissas |
5454
5b497d971546
dither zero-bit mantissas by default. remove dithering only if it's explicitly turned off.
jbr
parents:
5453
diff
changeset
|
508 * reference: Section 7.3.4 Dither for Zero Bit Mantissas (bap=0) |
5b497d971546
dither zero-bit mantissas by default. remove dithering only if it's explicitly turned off.
jbr
parents:
5453
diff
changeset
|
509 */ |
6090
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
510 static void remove_dithering(AC3DecodeContext *s) { |
5454
5b497d971546
dither zero-bit mantissas by default. remove dithering only if it's explicitly turned off.
jbr
parents:
5453
diff
changeset
|
511 int ch, i; |
5b497d971546
dither zero-bit mantissas by default. remove dithering only if it's explicitly turned off.
jbr
parents:
5453
diff
changeset
|
512 int end=0; |
6430
de7502093922
Use 24-bit fixed-point transform coefficients until just before MDCT. This
jbr
parents:
6142
diff
changeset
|
513 int *coeffs; |
5454
5b497d971546
dither zero-bit mantissas by default. remove dithering only if it's explicitly turned off.
jbr
parents:
5453
diff
changeset
|
514 uint8_t *bap; |
5b497d971546
dither zero-bit mantissas by default. remove dithering only if it's explicitly turned off.
jbr
parents:
5453
diff
changeset
|
515 |
6090
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
516 for(ch=1; ch<=s->fbw_channels; ch++) { |
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
517 if(!s->dither_flag[ch]) { |
6430
de7502093922
Use 24-bit fixed-point transform coefficients until just before MDCT. This
jbr
parents:
6142
diff
changeset
|
518 coeffs = s->fixed_coeffs[ch]; |
6090
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
519 bap = s->bap[ch]; |
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
520 if(s->channel_in_cpl[ch]) |
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
521 end = s->start_freq[CPL_CH]; |
5454
5b497d971546
dither zero-bit mantissas by default. remove dithering only if it's explicitly turned off.
jbr
parents:
5453
diff
changeset
|
522 else |
6090
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
523 end = s->end_freq[ch]; |
5454
5b497d971546
dither zero-bit mantissas by default. remove dithering only if it's explicitly turned off.
jbr
parents:
5453
diff
changeset
|
524 for(i=0; i<end; i++) { |
6095 | 525 if(!bap[i]) |
6430
de7502093922
Use 24-bit fixed-point transform coefficients until just before MDCT. This
jbr
parents:
6142
diff
changeset
|
526 coeffs[i] = 0; |
5454
5b497d971546
dither zero-bit mantissas by default. remove dithering only if it's explicitly turned off.
jbr
parents:
5453
diff
changeset
|
527 } |
6090
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
528 if(s->channel_in_cpl[ch]) { |
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
529 bap = s->bap[CPL_CH]; |
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
530 for(; i<s->end_freq[CPL_CH]; i++) { |
6095 | 531 if(!bap[i]) |
6430
de7502093922
Use 24-bit fixed-point transform coefficients until just before MDCT. This
jbr
parents:
6142
diff
changeset
|
532 coeffs[i] = 0; |
5454
5b497d971546
dither zero-bit mantissas by default. remove dithering only if it's explicitly turned off.
jbr
parents:
5453
diff
changeset
|
533 } |
5b497d971546
dither zero-bit mantissas by default. remove dithering only if it's explicitly turned off.
jbr
parents:
5453
diff
changeset
|
534 } |
5b497d971546
dither zero-bit mantissas by default. remove dithering only if it's explicitly turned off.
jbr
parents:
5453
diff
changeset
|
535 } |
5b497d971546
dither zero-bit mantissas by default. remove dithering only if it's explicitly turned off.
jbr
parents:
5453
diff
changeset
|
536 } |
5b497d971546
dither zero-bit mantissas by default. remove dithering only if it's explicitly turned off.
jbr
parents:
5453
diff
changeset
|
537 } |
5b497d971546
dither zero-bit mantissas by default. remove dithering only if it's explicitly turned off.
jbr
parents:
5453
diff
changeset
|
538 |
5515
b499866a3397
comments modification: add, clean up, Doxygenize, and reword
jbr
parents:
5499
diff
changeset
|
539 /** |
b499866a3397
comments modification: add, clean up, Doxygenize, and reword
jbr
parents:
5499
diff
changeset
|
540 * Get the transform coefficients. |
5321
fcb2a52c88fe
AC-3 decoder, soc revision 71, Sep 4 09:47:12 2006 UTC by cloud9
jbr
parents:
5320
diff
changeset
|
541 */ |
6953
64504967b891
get_transform_coeffs() never returns an error, so make the function and its children return void.
jbr
parents:
6951
diff
changeset
|
542 static void get_transform_coeffs(AC3DecodeContext *s) |
5303
b8821ff5c30d
AC-3 decoder, soc revision 9, Jun 14 03:52:02 2006 UTC by cloud9
jbr
parents:
diff
changeset
|
543 { |
5476
e50f0d583a36
consolidate decoding of lfe and coupling channels with decoding of full-bandwidth channels.
jbr
parents:
5475
diff
changeset
|
544 int ch, end; |
5304
eff51058fe13
AC-3 decoder, soc revision 26, Jul 5 04:55:15 2006 UTC by cloud9
jbr
parents:
5303
diff
changeset
|
545 int got_cplchan = 0; |
5306
abc5c130b448
AC-3 decoder, soc revision 32, Jul 17 09:37:32 2006 UTC by cloud9
jbr
parents:
5305
diff
changeset
|
546 mant_groups m; |
abc5c130b448
AC-3 decoder, soc revision 32, Jul 17 09:37:32 2006 UTC by cloud9
jbr
parents:
5305
diff
changeset
|
547 |
5456 | 548 m.b1ptr = m.b2ptr = m.b4ptr = 3; |
5303
b8821ff5c30d
AC-3 decoder, soc revision 9, Jun 14 03:52:02 2006 UTC by cloud9
jbr
parents:
diff
changeset
|
549 |
6090
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
550 for (ch = 1; ch <= s->channels; ch++) { |
5515
b499866a3397
comments modification: add, clean up, Doxygenize, and reword
jbr
parents:
5499
diff
changeset
|
551 /* transform coefficients for full-bandwidth channel */ |
6953
64504967b891
get_transform_coeffs() never returns an error, so make the function and its children return void.
jbr
parents:
6951
diff
changeset
|
552 get_transform_coeffs_ch(s, ch, &m); |
5515
b499866a3397
comments modification: add, clean up, Doxygenize, and reword
jbr
parents:
5499
diff
changeset
|
553 /* tranform coefficients for coupling channel come right after the |
b499866a3397
comments modification: add, clean up, Doxygenize, and reword
jbr
parents:
5499
diff
changeset
|
554 coefficients for the first coupled channel*/ |
6090
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
555 if (s->channel_in_cpl[ch]) { |
5306
abc5c130b448
AC-3 decoder, soc revision 32, Jul 17 09:37:32 2006 UTC by cloud9
jbr
parents:
5305
diff
changeset
|
556 if (!got_cplchan) { |
6953
64504967b891
get_transform_coeffs() never returns an error, so make the function and its children return void.
jbr
parents:
6951
diff
changeset
|
557 get_transform_coeffs_ch(s, CPL_CH, &m); |
7491
0dbdc638270e
cosmetics: use a better function name than uncouple_channels()
jbr
parents:
7490
diff
changeset
|
558 calc_transform_coeffs_cpl(s); |
5306
abc5c130b448
AC-3 decoder, soc revision 32, Jul 17 09:37:32 2006 UTC by cloud9
jbr
parents:
5305
diff
changeset
|
559 got_cplchan = 1; |
abc5c130b448
AC-3 decoder, soc revision 32, Jul 17 09:37:32 2006 UTC by cloud9
jbr
parents:
5305
diff
changeset
|
560 } |
6090
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
561 end = s->end_freq[CPL_CH]; |
5461
c427d8950769
cosmetics: use braces when it makes the code easier to read.
jbr
parents:
5460
diff
changeset
|
562 } else { |
6090
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
563 end = s->end_freq[ch]; |
5461
c427d8950769
cosmetics: use braces when it makes the code easier to read.
jbr
parents:
5460
diff
changeset
|
564 } |
5306
abc5c130b448
AC-3 decoder, soc revision 32, Jul 17 09:37:32 2006 UTC by cloud9
jbr
parents:
5305
diff
changeset
|
565 do |
6906
86e3dce7f723
zero the upper frequencies of the correct coefficients
jbr
parents:
6889
diff
changeset
|
566 s->fixed_coeffs[ch][end] = 0; |
5306
abc5c130b448
AC-3 decoder, soc revision 32, Jul 17 09:37:32 2006 UTC by cloud9
jbr
parents:
5305
diff
changeset
|
567 while(++end < 256); |
abc5c130b448
AC-3 decoder, soc revision 32, Jul 17 09:37:32 2006 UTC by cloud9
jbr
parents:
5305
diff
changeset
|
568 } |
5303
b8821ff5c30d
AC-3 decoder, soc revision 9, Jun 14 03:52:02 2006 UTC by cloud9
jbr
parents:
diff
changeset
|
569 |
5454
5b497d971546
dither zero-bit mantissas by default. remove dithering only if it's explicitly turned off.
jbr
parents:
5453
diff
changeset
|
570 /* if any channel doesn't use dithering, zero appropriate coefficients */ |
6090
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
571 if(!s->dither_all) |
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
572 remove_dithering(s); |
5303
b8821ff5c30d
AC-3 decoder, soc revision 9, Jun 14 03:52:02 2006 UTC by cloud9
jbr
parents:
diff
changeset
|
573 } |
b8821ff5c30d
AC-3 decoder, soc revision 9, Jun 14 03:52:02 2006 UTC by cloud9
jbr
parents:
diff
changeset
|
574 |
5415 | 575 /** |
5515
b499866a3397
comments modification: add, clean up, Doxygenize, and reword
jbr
parents:
5499
diff
changeset
|
576 * Stereo rematrixing. |
5415 | 577 * reference: Section 7.5.4 Rematrixing : Decoding Technique |
578 */ | |
6090
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
579 static void do_rematrixing(AC3DecodeContext *s) |
5304
eff51058fe13
AC-3 decoder, soc revision 26, Jul 5 04:55:15 2006 UTC by cloud9
jbr
parents:
5303
diff
changeset
|
580 { |
5415 | 581 int bnd, i; |
5310
9aa9197034d7
AC-3 decoder, soc revision 40, Aug 9 00:10:14 2006 UTC by cloud9
jbr
parents:
5309
diff
changeset
|
582 int end, bndend; |
6430
de7502093922
Use 24-bit fixed-point transform coefficients until just before MDCT. This
jbr
parents:
6142
diff
changeset
|
583 int tmp0, tmp1; |
5304
eff51058fe13
AC-3 decoder, soc revision 26, Jul 5 04:55:15 2006 UTC by cloud9
jbr
parents:
5303
diff
changeset
|
584 |
6090
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
585 end = FFMIN(s->end_freq[1], s->end_freq[2]); |
5310
9aa9197034d7
AC-3 decoder, soc revision 40, Aug 9 00:10:14 2006 UTC by cloud9
jbr
parents:
5309
diff
changeset
|
586 |
6090
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
587 for(bnd=0; bnd<s->num_rematrixing_bands; bnd++) { |
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
588 if(s->rematrixing_flags[bnd]) { |
7033 | 589 bndend = FFMIN(end, ff_ac3_rematrix_band_tab[bnd+1]); |
590 for(i=ff_ac3_rematrix_band_tab[bnd]; i<bndend; i++) { | |
6430
de7502093922
Use 24-bit fixed-point transform coefficients until just before MDCT. This
jbr
parents:
6142
diff
changeset
|
591 tmp0 = s->fixed_coeffs[1][i]; |
de7502093922
Use 24-bit fixed-point transform coefficients until just before MDCT. This
jbr
parents:
6142
diff
changeset
|
592 tmp1 = s->fixed_coeffs[2][i]; |
de7502093922
Use 24-bit fixed-point transform coefficients until just before MDCT. This
jbr
parents:
6142
diff
changeset
|
593 s->fixed_coeffs[1][i] = tmp0 + tmp1; |
de7502093922
Use 24-bit fixed-point transform coefficients until just before MDCT. This
jbr
parents:
6142
diff
changeset
|
594 s->fixed_coeffs[2][i] = tmp0 - tmp1; |
5415 | 595 } |
596 } | |
5304
eff51058fe13
AC-3 decoder, soc revision 26, Jul 5 04:55:15 2006 UTC by cloud9
jbr
parents:
5303
diff
changeset
|
597 } |
eff51058fe13
AC-3 decoder, soc revision 26, Jul 5 04:55:15 2006 UTC by cloud9
jbr
parents:
5303
diff
changeset
|
598 } |
eff51058fe13
AC-3 decoder, soc revision 26, Jul 5 04:55:15 2006 UTC by cloud9
jbr
parents:
5303
diff
changeset
|
599 |
5515
b499866a3397
comments modification: add, clean up, Doxygenize, and reword
jbr
parents:
5499
diff
changeset
|
600 /** |
b499866a3397
comments modification: add, clean up, Doxygenize, and reword
jbr
parents:
5499
diff
changeset
|
601 * Inverse MDCT Transform. |
b499866a3397
comments modification: add, clean up, Doxygenize, and reword
jbr
parents:
5499
diff
changeset
|
602 * Convert frequency domain coefficients to time-domain audio samples. |
b499866a3397
comments modification: add, clean up, Doxygenize, and reword
jbr
parents:
5499
diff
changeset
|
603 * reference: Section 7.9.4 Transformation Equations |
b499866a3397
comments modification: add, clean up, Doxygenize, and reword
jbr
parents:
5499
diff
changeset
|
604 */ |
6480
6f01a499e785
downmix before imdct unless different size transforms are used. about 20%
jbr
parents:
6432
diff
changeset
|
605 static inline void do_imdct(AC3DecodeContext *s, int channels) |
5306
abc5c130b448
AC-3 decoder, soc revision 32, Jul 17 09:37:32 2006 UTC by cloud9
jbr
parents:
5305
diff
changeset
|
606 { |
5336 | 607 int ch; |
7561
be6d31643128
oops, non-normalized downmix to mono contaminated add_bias for non-simd float-to-int.
lorenm
parents:
7554
diff
changeset
|
608 float add_bias = s->add_bias; |
be6d31643128
oops, non-normalized downmix to mono contaminated add_bias for non-simd float-to-int.
lorenm
parents:
7554
diff
changeset
|
609 if(s->out_channels==1 && channels>1) |
be6d31643128
oops, non-normalized downmix to mono contaminated add_bias for non-simd float-to-int.
lorenm
parents:
7554
diff
changeset
|
610 add_bias *= LEVEL_MINUS_3DB; // compensate for the gain in downmix |
5476
e50f0d583a36
consolidate decoding of lfe and coupling channels with decoding of full-bandwidth channels.
jbr
parents:
5475
diff
changeset
|
611 |
6006 | 612 for (ch=1; ch<=channels; ch++) { |
6090
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
613 if (s->block_switch[ch]) { |
7549 | 614 int i; |
615 float *x = s->tmp_output+128; | |
616 for(i=0; i<128; i++) | |
617 x[i] = s->transform_coeffs[ch][2*i]; | |
618 ff_imdct_half(&s->imdct_256, s->tmp_output, x); | |
7561
be6d31643128
oops, non-normalized downmix to mono contaminated add_bias for non-simd float-to-int.
lorenm
parents:
7554
diff
changeset
|
619 s->dsp.vector_fmul_window(s->output[ch-1], s->delay[ch-1], s->tmp_output, s->window, add_bias, 128); |
7549 | 620 for(i=0; i<128; i++) |
621 x[i] = s->transform_coeffs[ch][2*i+1]; | |
622 ff_imdct_half(&s->imdct_256, s->delay[ch-1], x); | |
5461
c427d8950769
cosmetics: use braces when it makes the code easier to read.
jbr
parents:
5460
diff
changeset
|
623 } else { |
7549 | 624 ff_imdct_half(&s->imdct_512, s->tmp_output, s->transform_coeffs[ch]); |
7561
be6d31643128
oops, non-normalized downmix to mono contaminated add_bias for non-simd float-to-int.
lorenm
parents:
7554
diff
changeset
|
625 s->dsp.vector_fmul_window(s->output[ch-1], s->delay[ch-1], s->tmp_output, s->window, add_bias, 128); |
7549 | 626 memcpy(s->delay[ch-1], s->tmp_output+128, 128*sizeof(float)); |
5461
c427d8950769
cosmetics: use braces when it makes the code easier to read.
jbr
parents:
5460
diff
changeset
|
627 } |
5306
abc5c130b448
AC-3 decoder, soc revision 32, Jul 17 09:37:32 2006 UTC by cloud9
jbr
parents:
5305
diff
changeset
|
628 } |
abc5c130b448
AC-3 decoder, soc revision 32, Jul 17 09:37:32 2006 UTC by cloud9
jbr
parents:
5305
diff
changeset
|
629 } |
abc5c130b448
AC-3 decoder, soc revision 32, Jul 17 09:37:32 2006 UTC by cloud9
jbr
parents:
5305
diff
changeset
|
630 |
5480 | 631 /** |
5515
b499866a3397
comments modification: add, clean up, Doxygenize, and reword
jbr
parents:
5499
diff
changeset
|
632 * Downmix the output to mono or stereo. |
5480 | 633 */ |
7563 | 634 void ff_ac3_downmix_c(float (*samples)[256], float (*matrix)[2], int out_ch, int in_ch, int len) |
5480 | 635 { |
636 int i, j; | |
6431
305c49259b59
pre-calculate downmix coefficient sums. 2% faster 5.1-to-stereo downmixing.
jbr
parents:
6430
diff
changeset
|
637 float v0, v1; |
7563 | 638 if(out_ch == 2) { |
639 for(i=0; i<len; i++) { | |
7554 | 640 v0 = v1 = 0.0f; |
7563 | 641 for(j=0; j<in_ch; j++) { |
642 v0 += samples[j][i] * matrix[j][0]; | |
643 v1 += samples[j][i] * matrix[j][1]; | |
7554 | 644 } |
645 samples[0][i] = v0; | |
646 samples[1][i] = v1; | |
5480 | 647 } |
7563 | 648 } else if(out_ch == 1) { |
649 for(i=0; i<len; i++) { | |
7554 | 650 v0 = 0.0f; |
7563 | 651 for(j=0; j<in_ch; j++) |
652 v0 += samples[j][i] * matrix[j][0]; | |
7554 | 653 samples[0][i] = v0; |
5480 | 654 } |
655 } | |
656 } | |
657 | |
5515
b499866a3397
comments modification: add, clean up, Doxygenize, and reword
jbr
parents:
5499
diff
changeset
|
658 /** |
6480
6f01a499e785
downmix before imdct unless different size transforms are used. about 20%
jbr
parents:
6432
diff
changeset
|
659 * Upmix delay samples from stereo to original channel layout. |
6f01a499e785
downmix before imdct unless different size transforms are used. about 20%
jbr
parents:
6432
diff
changeset
|
660 */ |
6f01a499e785
downmix before imdct unless different size transforms are used. about 20%
jbr
parents:
6432
diff
changeset
|
661 static void ac3_upmix_delay(AC3DecodeContext *s) |
6f01a499e785
downmix before imdct unless different size transforms are used. about 20%
jbr
parents:
6432
diff
changeset
|
662 { |
7549 | 663 int channel_data_size = 128*sizeof(float); |
6480
6f01a499e785
downmix before imdct unless different size transforms are used. about 20%
jbr
parents:
6432
diff
changeset
|
664 switch(s->channel_mode) { |
6f01a499e785
downmix before imdct unless different size transforms are used. about 20%
jbr
parents:
6432
diff
changeset
|
665 case AC3_CHMODE_DUALMONO: |
6f01a499e785
downmix before imdct unless different size transforms are used. about 20%
jbr
parents:
6432
diff
changeset
|
666 case AC3_CHMODE_STEREO: |
6f01a499e785
downmix before imdct unless different size transforms are used. about 20%
jbr
parents:
6432
diff
changeset
|
667 /* upmix mono to stereo */ |
6f01a499e785
downmix before imdct unless different size transforms are used. about 20%
jbr
parents:
6432
diff
changeset
|
668 memcpy(s->delay[1], s->delay[0], channel_data_size); |
6f01a499e785
downmix before imdct unless different size transforms are used. about 20%
jbr
parents:
6432
diff
changeset
|
669 break; |
6f01a499e785
downmix before imdct unless different size transforms are used. about 20%
jbr
parents:
6432
diff
changeset
|
670 case AC3_CHMODE_2F2R: |
6f01a499e785
downmix before imdct unless different size transforms are used. about 20%
jbr
parents:
6432
diff
changeset
|
671 memset(s->delay[3], 0, channel_data_size); |
6f01a499e785
downmix before imdct unless different size transforms are used. about 20%
jbr
parents:
6432
diff
changeset
|
672 case AC3_CHMODE_2F1R: |
6f01a499e785
downmix before imdct unless different size transforms are used. about 20%
jbr
parents:
6432
diff
changeset
|
673 memset(s->delay[2], 0, channel_data_size); |
6f01a499e785
downmix before imdct unless different size transforms are used. about 20%
jbr
parents:
6432
diff
changeset
|
674 break; |
6f01a499e785
downmix before imdct unless different size transforms are used. about 20%
jbr
parents:
6432
diff
changeset
|
675 case AC3_CHMODE_3F2R: |
6f01a499e785
downmix before imdct unless different size transforms are used. about 20%
jbr
parents:
6432
diff
changeset
|
676 memset(s->delay[4], 0, channel_data_size); |
6f01a499e785
downmix before imdct unless different size transforms are used. about 20%
jbr
parents:
6432
diff
changeset
|
677 case AC3_CHMODE_3F1R: |
6f01a499e785
downmix before imdct unless different size transforms are used. about 20%
jbr
parents:
6432
diff
changeset
|
678 memset(s->delay[3], 0, channel_data_size); |
6f01a499e785
downmix before imdct unless different size transforms are used. about 20%
jbr
parents:
6432
diff
changeset
|
679 case AC3_CHMODE_3F: |
6f01a499e785
downmix before imdct unless different size transforms are used. about 20%
jbr
parents:
6432
diff
changeset
|
680 memcpy(s->delay[2], s->delay[1], channel_data_size); |
6f01a499e785
downmix before imdct unless different size transforms are used. about 20%
jbr
parents:
6432
diff
changeset
|
681 memset(s->delay[1], 0, channel_data_size); |
6f01a499e785
downmix before imdct unless different size transforms are used. about 20%
jbr
parents:
6432
diff
changeset
|
682 break; |
6f01a499e785
downmix before imdct unless different size transforms are used. about 20%
jbr
parents:
6432
diff
changeset
|
683 } |
6f01a499e785
downmix before imdct unless different size transforms are used. about 20%
jbr
parents:
6432
diff
changeset
|
684 } |
6f01a499e785
downmix before imdct unless different size transforms are used. about 20%
jbr
parents:
6432
diff
changeset
|
685 |
6f01a499e785
downmix before imdct unless different size transforms are used. about 20%
jbr
parents:
6432
diff
changeset
|
686 /** |
7490
bddd15a1d167
cosmetics: change function name and comments to refer to block decoding
jbr
parents:
7489
diff
changeset
|
687 * Decode a single audio block from the AC-3 bitstream. |
5321
fcb2a52c88fe
AC-3 decoder, soc revision 71, Sep 4 09:47:12 2006 UTC by cloud9
jbr
parents:
5320
diff
changeset
|
688 */ |
7490
bddd15a1d167
cosmetics: change function name and comments to refer to block decoding
jbr
parents:
7489
diff
changeset
|
689 static int decode_audio_block(AC3DecodeContext *s, int blk) |
5304
eff51058fe13
AC-3 decoder, soc revision 26, Jul 5 04:55:15 2006 UTC by cloud9
jbr
parents:
5303
diff
changeset
|
690 { |
6090
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
691 int fbw_channels = s->fbw_channels; |
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
692 int channel_mode = s->channel_mode; |
5471
37c1944ca4ee
make exponent decoding consistent (prepare for merging coupling and lfe code with normal channel code).
jbr
parents:
5461
diff
changeset
|
693 int i, bnd, seg, ch; |
6480
6f01a499e785
downmix before imdct unless different size transforms are used. about 20%
jbr
parents:
6432
diff
changeset
|
694 int different_transforms; |
6f01a499e785
downmix before imdct unless different size transforms are used. about 20%
jbr
parents:
6432
diff
changeset
|
695 int downmix_output; |
7026 | 696 int cpl_in_use; |
6090
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
697 GetBitContext *gbc = &s->gbc; |
5476
e50f0d583a36
consolidate decoding of lfe and coupling channels with decoding of full-bandwidth channels.
jbr
parents:
5475
diff
changeset
|
698 uint8_t bit_alloc_stages[AC3_MAX_CHANNELS]; |
5304
eff51058fe13
AC-3 decoder, soc revision 26, Jul 5 04:55:15 2006 UTC by cloud9
jbr
parents:
5303
diff
changeset
|
699 |
5476
e50f0d583a36
consolidate decoding of lfe and coupling channels with decoding of full-bandwidth channels.
jbr
parents:
5475
diff
changeset
|
700 memset(bit_alloc_stages, 0, AC3_MAX_CHANNELS); |
e50f0d583a36
consolidate decoding of lfe and coupling channels with decoding of full-bandwidth channels.
jbr
parents:
5475
diff
changeset
|
701 |
5515
b499866a3397
comments modification: add, clean up, Doxygenize, and reword
jbr
parents:
5499
diff
changeset
|
702 /* block switch flags */ |
6480
6f01a499e785
downmix before imdct unless different size transforms are used. about 20%
jbr
parents:
6432
diff
changeset
|
703 different_transforms = 0; |
6f01a499e785
downmix before imdct unless different size transforms are used. about 20%
jbr
parents:
6432
diff
changeset
|
704 for (ch = 1; ch <= fbw_channels; ch++) { |
6090
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
705 s->block_switch[ch] = get_bits1(gbc); |
6480
6f01a499e785
downmix before imdct unless different size transforms are used. about 20%
jbr
parents:
6432
diff
changeset
|
706 if(ch > 1 && s->block_switch[ch] != s->block_switch[1]) |
6f01a499e785
downmix before imdct unless different size transforms are used. about 20%
jbr
parents:
6432
diff
changeset
|
707 different_transforms = 1; |
6f01a499e785
downmix before imdct unless different size transforms are used. about 20%
jbr
parents:
6432
diff
changeset
|
708 } |
5309
0662a270aab7
AC-3 decoder, soc revision 38, Aug 7 00:03:00 2006 UTC by cloud9
jbr
parents:
5308
diff
changeset
|
709 |
5515
b499866a3397
comments modification: add, clean up, Doxygenize, and reword
jbr
parents:
5499
diff
changeset
|
710 /* dithering flags */ |
6090
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
711 s->dither_all = 1; |
6006 | 712 for (ch = 1; ch <= fbw_channels; ch++) { |
6090
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
713 s->dither_flag[ch] = get_bits1(gbc); |
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
714 if(!s->dither_flag[ch]) |
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
715 s->dither_all = 0; |
5454
5b497d971546
dither zero-bit mantissas by default. remove dithering only if it's explicitly turned off.
jbr
parents:
5453
diff
changeset
|
716 } |
5309
0662a270aab7
AC-3 decoder, soc revision 38, Aug 7 00:03:00 2006 UTC by cloud9
jbr
parents:
5308
diff
changeset
|
717 |
5497 | 718 /* dynamic range */ |
6090
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
719 i = !(s->channel_mode); |
5497 | 720 do { |
6089 | 721 if(get_bits1(gbc)) { |
6090
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
722 s->dynamic_range[i] = ((dynamic_range_tab[get_bits(gbc, 8)]-1.0) * |
6091 | 723 s->avctx->drc_scale)+1.0; |
5386
d6347c3cb329
Use shared AC-3 header parsing. Move some initializations to block parsing function.
jbr
parents:
5385
diff
changeset
|
724 } else if(blk == 0) { |
6090
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
725 s->dynamic_range[i] = 1.0f; |
5386
d6347c3cb329
Use shared AC-3 header parsing. Move some initializations to block parsing function.
jbr
parents:
5385
diff
changeset
|
726 } |
5497 | 727 } while(i--); |
5309
0662a270aab7
AC-3 decoder, soc revision 38, Aug 7 00:03:00 2006 UTC by cloud9
jbr
parents:
5308
diff
changeset
|
728 |
5515
b499866a3397
comments modification: add, clean up, Doxygenize, and reword
jbr
parents:
5499
diff
changeset
|
729 /* coupling strategy */ |
6089 | 730 if (get_bits1(gbc)) { |
5476
e50f0d583a36
consolidate decoding of lfe and coupling channels with decoding of full-bandwidth channels.
jbr
parents:
5475
diff
changeset
|
731 memset(bit_alloc_stages, 3, AC3_MAX_CHANNELS); |
7328
9309e19c2b98
read 'coupling in use' parameter from bitstream into decoder context first,
jbr
parents:
7311
diff
changeset
|
732 s->cpl_in_use[blk] = get_bits1(gbc); |
9309e19c2b98
read 'coupling in use' parameter from bitstream into decoder context first,
jbr
parents:
7311
diff
changeset
|
733 if (s->cpl_in_use[blk]) { |
5515
b499866a3397
comments modification: add, clean up, Doxygenize, and reword
jbr
parents:
5499
diff
changeset
|
734 /* coupling in use */ |
6006 | 735 int cpl_begin_freq, cpl_end_freq; |
5414 | 736 |
6868
5387096bf92a
add check for invalid channel mode when coupling is used
jbr
parents:
6763
diff
changeset
|
737 if (channel_mode < AC3_CHMODE_STEREO) { |
5387096bf92a
add check for invalid channel mode when coupling is used
jbr
parents:
6763
diff
changeset
|
738 av_log(s->avctx, AV_LOG_ERROR, "coupling not allowed in mono or dual-mono\n"); |
5387096bf92a
add check for invalid channel mode when coupling is used
jbr
parents:
6763
diff
changeset
|
739 return -1; |
5387096bf92a
add check for invalid channel mode when coupling is used
jbr
parents:
6763
diff
changeset
|
740 } |
5387096bf92a
add check for invalid channel mode when coupling is used
jbr
parents:
6763
diff
changeset
|
741 |
5515
b499866a3397
comments modification: add, clean up, Doxygenize, and reword
jbr
parents:
5499
diff
changeset
|
742 /* determine which channels are coupled */ |
6006 | 743 for (ch = 1; ch <= fbw_channels; ch++) |
6090
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
744 s->channel_in_cpl[ch] = get_bits1(gbc); |
5309
0662a270aab7
AC-3 decoder, soc revision 38, Aug 7 00:03:00 2006 UTC by cloud9
jbr
parents:
5308
diff
changeset
|
745 |
5515
b499866a3397
comments modification: add, clean up, Doxygenize, and reword
jbr
parents:
5499
diff
changeset
|
746 /* phase flags in use */ |
6005 | 747 if (channel_mode == AC3_CHMODE_STEREO) |
6090
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
748 s->phase_flags_in_use = get_bits1(gbc); |
5309
0662a270aab7
AC-3 decoder, soc revision 38, Aug 7 00:03:00 2006 UTC by cloud9
jbr
parents:
5308
diff
changeset
|
749 |
5515
b499866a3397
comments modification: add, clean up, Doxygenize, and reword
jbr
parents:
5499
diff
changeset
|
750 /* coupling frequency range and band structure */ |
6089 | 751 cpl_begin_freq = get_bits(gbc, 4); |
752 cpl_end_freq = get_bits(gbc, 4); | |
6006 | 753 if (3 + cpl_end_freq - cpl_begin_freq < 0) { |
6090
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
754 av_log(s->avctx, AV_LOG_ERROR, "3+cplendf = %d < cplbegf = %d\n", 3+cpl_end_freq, cpl_begin_freq); |
5309
0662a270aab7
AC-3 decoder, soc revision 38, Aug 7 00:03:00 2006 UTC by cloud9
jbr
parents:
5308
diff
changeset
|
755 return -1; |
0662a270aab7
AC-3 decoder, soc revision 38, Aug 7 00:03:00 2006 UTC by cloud9
jbr
parents:
5308
diff
changeset
|
756 } |
6090
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
757 s->num_cpl_bands = s->num_cpl_subbands = 3 + cpl_end_freq - cpl_begin_freq; |
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
758 s->start_freq[CPL_CH] = cpl_begin_freq * 12 + 37; |
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
759 s->end_freq[CPL_CH] = cpl_end_freq * 12 + 73; |
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
760 for (bnd = 0; bnd < s->num_cpl_subbands - 1; bnd++) { |
6089 | 761 if (get_bits1(gbc)) { |
6090
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
762 s->cpl_band_struct[bnd] = 1; |
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
763 s->num_cpl_bands--; |
5304
eff51058fe13
AC-3 decoder, soc revision 26, Jul 5 04:55:15 2006 UTC by cloud9
jbr
parents:
5303
diff
changeset
|
764 } |
5461
c427d8950769
cosmetics: use braces when it makes the code easier to read.
jbr
parents:
5460
diff
changeset
|
765 } |
6124
d06c62429572
make sure that the last coupling band stops at the end of the coupling range
jbr
parents:
6123
diff
changeset
|
766 s->cpl_band_struct[s->num_cpl_subbands-1] = 0; |
5388 | 767 } else { |
5515
b499866a3397
comments modification: add, clean up, Doxygenize, and reword
jbr
parents:
5499
diff
changeset
|
768 /* coupling not in use */ |
6006 | 769 for (ch = 1; ch <= fbw_channels; ch++) |
6090
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
770 s->channel_in_cpl[ch] = 0; |
5304
eff51058fe13
AC-3 decoder, soc revision 26, Jul 5 04:55:15 2006 UTC by cloud9
jbr
parents:
5303
diff
changeset
|
771 } |
6869 | 772 } else if (!blk) { |
773 av_log(s->avctx, AV_LOG_ERROR, "new coupling strategy must be present in block 0\n"); | |
774 return -1; | |
7026 | 775 } else { |
7328
9309e19c2b98
read 'coupling in use' parameter from bitstream into decoder context first,
jbr
parents:
7311
diff
changeset
|
776 s->cpl_in_use[blk] = s->cpl_in_use[blk-1]; |
5303
b8821ff5c30d
AC-3 decoder, soc revision 9, Jun 14 03:52:02 2006 UTC by cloud9
jbr
parents:
diff
changeset
|
777 } |
7328
9309e19c2b98
read 'coupling in use' parameter from bitstream into decoder context first,
jbr
parents:
7311
diff
changeset
|
778 cpl_in_use = s->cpl_in_use[blk]; |
5309
0662a270aab7
AC-3 decoder, soc revision 38, Aug 7 00:03:00 2006 UTC by cloud9
jbr
parents:
5308
diff
changeset
|
779 |
5515
b499866a3397
comments modification: add, clean up, Doxygenize, and reword
jbr
parents:
5499
diff
changeset
|
780 /* coupling coordinates */ |
7026 | 781 if (cpl_in_use) { |
6006 | 782 int cpl_coords_exist = 0; |
5309
0662a270aab7
AC-3 decoder, soc revision 38, Aug 7 00:03:00 2006 UTC by cloud9
jbr
parents:
5308
diff
changeset
|
783 |
6006 | 784 for (ch = 1; ch <= fbw_channels; ch++) { |
6090
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
785 if (s->channel_in_cpl[ch]) { |
6089 | 786 if (get_bits1(gbc)) { |
6006 | 787 int master_cpl_coord, cpl_coord_exp, cpl_coord_mant; |
788 cpl_coords_exist = 1; | |
6089 | 789 master_cpl_coord = 3 * get_bits(gbc, 2); |
6090
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
790 for (bnd = 0; bnd < s->num_cpl_bands; bnd++) { |
6089 | 791 cpl_coord_exp = get_bits(gbc, 4); |
792 cpl_coord_mant = get_bits(gbc, 4); | |
6006 | 793 if (cpl_coord_exp == 15) |
6430
de7502093922
Use 24-bit fixed-point transform coefficients until just before MDCT. This
jbr
parents:
6142
diff
changeset
|
794 s->cpl_coords[ch][bnd] = cpl_coord_mant << 22; |
5306
abc5c130b448
AC-3 decoder, soc revision 32, Jul 17 09:37:32 2006 UTC by cloud9
jbr
parents:
5305
diff
changeset
|
795 else |
6430
de7502093922
Use 24-bit fixed-point transform coefficients until just before MDCT. This
jbr
parents:
6142
diff
changeset
|
796 s->cpl_coords[ch][bnd] = (cpl_coord_mant + 16) << 21; |
de7502093922
Use 24-bit fixed-point transform coefficients until just before MDCT. This
jbr
parents:
6142
diff
changeset
|
797 s->cpl_coords[ch][bnd] >>= (cpl_coord_exp + master_cpl_coord); |
5304
eff51058fe13
AC-3 decoder, soc revision 26, Jul 5 04:55:15 2006 UTC by cloud9
jbr
parents:
5303
diff
changeset
|
798 } |
6869 | 799 } else if (!blk) { |
800 av_log(s->avctx, AV_LOG_ERROR, "new coupling coordinates must be present in block 0\n"); | |
801 return -1; | |
5304
eff51058fe13
AC-3 decoder, soc revision 26, Jul 5 04:55:15 2006 UTC by cloud9
jbr
parents:
5303
diff
changeset
|
802 } |
5461
c427d8950769
cosmetics: use braces when it makes the code easier to read.
jbr
parents:
5460
diff
changeset
|
803 } |
c427d8950769
cosmetics: use braces when it makes the code easier to read.
jbr
parents:
5460
diff
changeset
|
804 } |
5515
b499866a3397
comments modification: add, clean up, Doxygenize, and reword
jbr
parents:
5499
diff
changeset
|
805 /* phase flags */ |
6123
8d88831a2777
fix phase flag processing for the case when coupling coordinates are reused
jbr
parents:
6117
diff
changeset
|
806 if (channel_mode == AC3_CHMODE_STEREO && cpl_coords_exist) { |
6090
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
807 for (bnd = 0; bnd < s->num_cpl_bands; bnd++) { |
6123
8d88831a2777
fix phase flag processing for the case when coupling coordinates are reused
jbr
parents:
6117
diff
changeset
|
808 s->phase_flags[bnd] = s->phase_flags_in_use? get_bits1(gbc) : 0; |
5461
c427d8950769
cosmetics: use braces when it makes the code easier to read.
jbr
parents:
5460
diff
changeset
|
809 } |
c427d8950769
cosmetics: use braces when it makes the code easier to read.
jbr
parents:
5460
diff
changeset
|
810 } |
5309
0662a270aab7
AC-3 decoder, soc revision 38, Aug 7 00:03:00 2006 UTC by cloud9
jbr
parents:
5308
diff
changeset
|
811 } |
0662a270aab7
AC-3 decoder, soc revision 38, Aug 7 00:03:00 2006 UTC by cloud9
jbr
parents:
5308
diff
changeset
|
812 |
5515
b499866a3397
comments modification: add, clean up, Doxygenize, and reword
jbr
parents:
5499
diff
changeset
|
813 /* stereo rematrixing strategy and band structure */ |
6005 | 814 if (channel_mode == AC3_CHMODE_STEREO) { |
6125 | 815 if (get_bits1(gbc)) { |
6090
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
816 s->num_rematrixing_bands = 4; |
7026 | 817 if(cpl_in_use && s->start_freq[CPL_CH] <= 61) |
6090
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
818 s->num_rematrixing_bands -= 1 + (s->start_freq[CPL_CH] == 37); |
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
819 for(bnd=0; bnd<s->num_rematrixing_bands; bnd++) |
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
820 s->rematrixing_flags[bnd] = get_bits1(gbc); |
6869 | 821 } else if (!blk) { |
822 av_log(s->avctx, AV_LOG_ERROR, "new rematrixing strategy must be present in block 0\n"); | |
823 return -1; | |
5305
5892b4a6380b
AC-3 decoder, soc revision 31, Jul 14 23:53:28 2006 UTC by cloud9
jbr
parents:
5304
diff
changeset
|
824 } |
5303
b8821ff5c30d
AC-3 decoder, soc revision 9, Jun 14 03:52:02 2006 UTC by cloud9
jbr
parents:
diff
changeset
|
825 } |
5309
0662a270aab7
AC-3 decoder, soc revision 38, Aug 7 00:03:00 2006 UTC by cloud9
jbr
parents:
5308
diff
changeset
|
826 |
5515
b499866a3397
comments modification: add, clean up, Doxygenize, and reword
jbr
parents:
5499
diff
changeset
|
827 /* exponent strategies for each channel */ |
7027 | 828 s->exp_strategy[blk][CPL_CH] = EXP_REUSE; |
829 s->exp_strategy[blk][s->lfe_ch] = EXP_REUSE; | |
7026 | 830 for (ch = !cpl_in_use; ch <= s->channels; ch++) { |
7027 | 831 s->exp_strategy[blk][ch] = get_bits(gbc, 2 - (ch == s->lfe_ch)); |
832 if(s->exp_strategy[blk][ch] != EXP_REUSE) | |
5476
e50f0d583a36
consolidate decoding of lfe and coupling channels with decoding of full-bandwidth channels.
jbr
parents:
5475
diff
changeset
|
833 bit_alloc_stages[ch] = 3; |
e50f0d583a36
consolidate decoding of lfe and coupling channels with decoding of full-bandwidth channels.
jbr
parents:
5475
diff
changeset
|
834 } |
5309
0662a270aab7
AC-3 decoder, soc revision 38, Aug 7 00:03:00 2006 UTC by cloud9
jbr
parents:
5308
diff
changeset
|
835 |
5515
b499866a3397
comments modification: add, clean up, Doxygenize, and reword
jbr
parents:
5499
diff
changeset
|
836 /* channel bandwidth */ |
6006 | 837 for (ch = 1; ch <= fbw_channels; ch++) { |
6090
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
838 s->start_freq[ch] = 0; |
7027 | 839 if (s->exp_strategy[blk][ch] != EXP_REUSE) { |
6957
6eb895971766
only calculate number of exponent groups when exponents are not reused.
jbr
parents:
6956
diff
changeset
|
840 int group_size; |
6090
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
841 int prev = s->end_freq[ch]; |
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
842 if (s->channel_in_cpl[ch]) |
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
843 s->end_freq[ch] = s->start_freq[CPL_CH]; |
5305
5892b4a6380b
AC-3 decoder, soc revision 31, Jul 14 23:53:28 2006 UTC by cloud9
jbr
parents:
5304
diff
changeset
|
844 else { |
6089 | 845 int bandwidth_code = get_bits(gbc, 6); |
6006 | 846 if (bandwidth_code > 60) { |
6090
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
847 av_log(s->avctx, AV_LOG_ERROR, "bandwidth code = %d > 60", bandwidth_code); |
5305
5892b4a6380b
AC-3 decoder, soc revision 31, Jul 14 23:53:28 2006 UTC by cloud9
jbr
parents:
5304
diff
changeset
|
848 return -1; |
5892b4a6380b
AC-3 decoder, soc revision 31, Jul 14 23:53:28 2006 UTC by cloud9
jbr
parents:
5304
diff
changeset
|
849 } |
6090
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
850 s->end_freq[ch] = bandwidth_code * 3 + 73; |
5304
eff51058fe13
AC-3 decoder, soc revision 26, Jul 5 04:55:15 2006 UTC by cloud9
jbr
parents:
5303
diff
changeset
|
851 } |
7027 | 852 group_size = 3 << (s->exp_strategy[blk][ch] - 1); |
6957
6eb895971766
only calculate number of exponent groups when exponents are not reused.
jbr
parents:
6956
diff
changeset
|
853 s->num_exp_groups[ch] = (s->end_freq[ch]+group_size-4) / group_size; |
6090
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
854 if(blk > 0 && s->end_freq[ch] != prev) |
5476
e50f0d583a36
consolidate decoding of lfe and coupling channels with decoding of full-bandwidth channels.
jbr
parents:
5475
diff
changeset
|
855 memset(bit_alloc_stages, 3, AC3_MAX_CHANNELS); |
5305
5892b4a6380b
AC-3 decoder, soc revision 31, Jul 14 23:53:28 2006 UTC by cloud9
jbr
parents:
5304
diff
changeset
|
856 } |
5461
c427d8950769
cosmetics: use braces when it makes the code easier to read.
jbr
parents:
5460
diff
changeset
|
857 } |
7027 | 858 if (cpl_in_use && s->exp_strategy[blk][CPL_CH] != EXP_REUSE) { |
6957
6eb895971766
only calculate number of exponent groups when exponents are not reused.
jbr
parents:
6956
diff
changeset
|
859 s->num_exp_groups[CPL_CH] = (s->end_freq[CPL_CH] - s->start_freq[CPL_CH]) / |
7027 | 860 (3 << (s->exp_strategy[blk][CPL_CH] - 1)); |
6957
6eb895971766
only calculate number of exponent groups when exponents are not reused.
jbr
parents:
6956
diff
changeset
|
861 } |
5309
0662a270aab7
AC-3 decoder, soc revision 38, Aug 7 00:03:00 2006 UTC by cloud9
jbr
parents:
5308
diff
changeset
|
862 |
5515
b499866a3397
comments modification: add, clean up, Doxygenize, and reword
jbr
parents:
5499
diff
changeset
|
863 /* decode exponents for each channel */ |
7026 | 864 for (ch = !cpl_in_use; ch <= s->channels; ch++) { |
7027 | 865 if (s->exp_strategy[blk][ch] != EXP_REUSE) { |
6090
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
866 s->dexps[ch][0] = get_bits(gbc, 4) << !ch; |
7027 | 867 decode_exponents(gbc, s->exp_strategy[blk][ch], |
6957
6eb895971766
only calculate number of exponent groups when exponents are not reused.
jbr
parents:
6956
diff
changeset
|
868 s->num_exp_groups[ch], s->dexps[ch][0], |
6090
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
869 &s->dexps[ch][s->start_freq[ch]+!!ch]); |
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
870 if(ch != CPL_CH && ch != s->lfe_ch) |
6089 | 871 skip_bits(gbc, 2); /* skip gainrng */ |
5304
eff51058fe13
AC-3 decoder, soc revision 26, Jul 5 04:55:15 2006 UTC by cloud9
jbr
parents:
5303
diff
changeset
|
872 } |
5461
c427d8950769
cosmetics: use braces when it makes the code easier to read.
jbr
parents:
5460
diff
changeset
|
873 } |
5309
0662a270aab7
AC-3 decoder, soc revision 38, Aug 7 00:03:00 2006 UTC by cloud9
jbr
parents:
5308
diff
changeset
|
874 |
5515
b499866a3397
comments modification: add, clean up, Doxygenize, and reword
jbr
parents:
5499
diff
changeset
|
875 /* bit allocation information */ |
6089 | 876 if (get_bits1(gbc)) { |
6090
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
877 s->bit_alloc_params.slow_decay = ff_ac3_slow_decay_tab[get_bits(gbc, 2)] >> s->bit_alloc_params.sr_shift; |
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
878 s->bit_alloc_params.fast_decay = ff_ac3_fast_decay_tab[get_bits(gbc, 2)] >> s->bit_alloc_params.sr_shift; |
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
879 s->bit_alloc_params.slow_gain = ff_ac3_slow_gain_tab[get_bits(gbc, 2)]; |
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
880 s->bit_alloc_params.db_per_bit = ff_ac3_db_per_bit_tab[get_bits(gbc, 2)]; |
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
881 s->bit_alloc_params.floor = ff_ac3_floor_tab[get_bits(gbc, 3)]; |
7026 | 882 for(ch=!cpl_in_use; ch<=s->channels; ch++) |
5476
e50f0d583a36
consolidate decoding of lfe and coupling channels with decoding of full-bandwidth channels.
jbr
parents:
5475
diff
changeset
|
883 bit_alloc_stages[ch] = FFMAX(bit_alloc_stages[ch], 2); |
6869 | 884 } else if (!blk) { |
885 av_log(s->avctx, AV_LOG_ERROR, "new bit allocation info must be present in block 0\n"); | |
886 return -1; | |
5304
eff51058fe13
AC-3 decoder, soc revision 26, Jul 5 04:55:15 2006 UTC by cloud9
jbr
parents:
5303
diff
changeset
|
887 } |
5309
0662a270aab7
AC-3 decoder, soc revision 38, Aug 7 00:03:00 2006 UTC by cloud9
jbr
parents:
5308
diff
changeset
|
888 |
5515
b499866a3397
comments modification: add, clean up, Doxygenize, and reword
jbr
parents:
5499
diff
changeset
|
889 /* signal-to-noise ratio offsets and fast gains (signal-to-mask ratios) */ |
6089 | 890 if (get_bits1(gbc)) { |
5393 | 891 int csnr; |
6089 | 892 csnr = (get_bits(gbc, 6) - 15) << 4; |
7026 | 893 for (ch = !cpl_in_use; ch <= s->channels; ch++) { /* snr offset and fast gain */ |
6090
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
894 s->snr_offset[ch] = (csnr + get_bits(gbc, 4)) << 2; |
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
895 s->fast_gain[ch] = ff_ac3_fast_gain_tab[get_bits(gbc, 3)]; |
5304
eff51058fe13
AC-3 decoder, soc revision 26, Jul 5 04:55:15 2006 UTC by cloud9
jbr
parents:
5303
diff
changeset
|
896 } |
5476
e50f0d583a36
consolidate decoding of lfe and coupling channels with decoding of full-bandwidth channels.
jbr
parents:
5475
diff
changeset
|
897 memset(bit_alloc_stages, 3, AC3_MAX_CHANNELS); |
6869 | 898 } else if (!blk) { |
899 av_log(s->avctx, AV_LOG_ERROR, "new snr offsets must be present in block 0\n"); | |
900 return -1; | |
5303
b8821ff5c30d
AC-3 decoder, soc revision 9, Jun 14 03:52:02 2006 UTC by cloud9
jbr
parents:
diff
changeset
|
901 } |
5309
0662a270aab7
AC-3 decoder, soc revision 38, Aug 7 00:03:00 2006 UTC by cloud9
jbr
parents:
5308
diff
changeset
|
902 |
5515
b499866a3397
comments modification: add, clean up, Doxygenize, and reword
jbr
parents:
5499
diff
changeset
|
903 /* coupling leak information */ |
7026 | 904 if (cpl_in_use) { |
6871 | 905 if (get_bits1(gbc)) { |
6872 | 906 s->bit_alloc_params.cpl_fast_leak = get_bits(gbc, 3); |
907 s->bit_alloc_params.cpl_slow_leak = get_bits(gbc, 3); | |
908 bit_alloc_stages[CPL_CH] = FFMAX(bit_alloc_stages[CPL_CH], 2); | |
909 } else if (!blk) { | |
910 av_log(s->avctx, AV_LOG_ERROR, "new coupling leak info must be present in block 0\n"); | |
911 return -1; | |
912 } | |
6871 | 913 } |
5309
0662a270aab7
AC-3 decoder, soc revision 38, Aug 7 00:03:00 2006 UTC by cloud9
jbr
parents:
5308
diff
changeset
|
914 |
5515
b499866a3397
comments modification: add, clean up, Doxygenize, and reword
jbr
parents:
5499
diff
changeset
|
915 /* delta bit allocation information */ |
6089 | 916 if (get_bits1(gbc)) { |
5515
b499866a3397
comments modification: add, clean up, Doxygenize, and reword
jbr
parents:
5499
diff
changeset
|
917 /* delta bit allocation exists (strategy) */ |
7026 | 918 for (ch = !cpl_in_use; ch <= fbw_channels; ch++) { |
6090
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
919 s->dba_mode[ch] = get_bits(gbc, 2); |
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
920 if (s->dba_mode[ch] == DBA_RESERVED) { |
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
921 av_log(s->avctx, AV_LOG_ERROR, "delta bit allocation strategy reserved\n"); |
5304
eff51058fe13
AC-3 decoder, soc revision 26, Jul 5 04:55:15 2006 UTC by cloud9
jbr
parents:
5303
diff
changeset
|
922 return -1; |
eff51058fe13
AC-3 decoder, soc revision 26, Jul 5 04:55:15 2006 UTC by cloud9
jbr
parents:
5303
diff
changeset
|
923 } |
5476
e50f0d583a36
consolidate decoding of lfe and coupling channels with decoding of full-bandwidth channels.
jbr
parents:
5475
diff
changeset
|
924 bit_alloc_stages[ch] = FFMAX(bit_alloc_stages[ch], 2); |
5304
eff51058fe13
AC-3 decoder, soc revision 26, Jul 5 04:55:15 2006 UTC by cloud9
jbr
parents:
5303
diff
changeset
|
925 } |
5515
b499866a3397
comments modification: add, clean up, Doxygenize, and reword
jbr
parents:
5499
diff
changeset
|
926 /* channel delta offset, len and bit allocation */ |
7026 | 927 for (ch = !cpl_in_use; ch <= fbw_channels; ch++) { |
6090
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
928 if (s->dba_mode[ch] == DBA_NEW) { |
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
929 s->dba_nsegs[ch] = get_bits(gbc, 3); |
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
930 for (seg = 0; seg <= s->dba_nsegs[ch]; seg++) { |
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
931 s->dba_offsets[ch][seg] = get_bits(gbc, 5); |
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
932 s->dba_lengths[ch][seg] = get_bits(gbc, 4); |
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
933 s->dba_values[ch][seg] = get_bits(gbc, 3); |
5304
eff51058fe13
AC-3 decoder, soc revision 26, Jul 5 04:55:15 2006 UTC by cloud9
jbr
parents:
5303
diff
changeset
|
934 } |
6954
57797b6c9d59
new dba values requires running last 2 bit allocation stages
jbr
parents:
6953
diff
changeset
|
935 /* run last 2 bit allocation stages if new dba values */ |
57797b6c9d59
new dba values requires running last 2 bit allocation stages
jbr
parents:
6953
diff
changeset
|
936 bit_alloc_stages[ch] = FFMAX(bit_alloc_stages[ch], 2); |
5304
eff51058fe13
AC-3 decoder, soc revision 26, Jul 5 04:55:15 2006 UTC by cloud9
jbr
parents:
5303
diff
changeset
|
937 } |
5461
c427d8950769
cosmetics: use braces when it makes the code easier to read.
jbr
parents:
5460
diff
changeset
|
938 } |
5386
d6347c3cb329
Use shared AC-3 header parsing. Move some initializations to block parsing function.
jbr
parents:
5385
diff
changeset
|
939 } else if(blk == 0) { |
6090
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
940 for(ch=0; ch<=s->channels; ch++) { |
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
941 s->dba_mode[ch] = DBA_NONE; |
5386
d6347c3cb329
Use shared AC-3 header parsing. Move some initializations to block parsing function.
jbr
parents:
5385
diff
changeset
|
942 } |
5303
b8821ff5c30d
AC-3 decoder, soc revision 9, Jun 14 03:52:02 2006 UTC by cloud9
jbr
parents:
diff
changeset
|
943 } |
5305
5892b4a6380b
AC-3 decoder, soc revision 31, Jul 14 23:53:28 2006 UTC by cloud9
jbr
parents:
5304
diff
changeset
|
944 |
5515
b499866a3397
comments modification: add, clean up, Doxygenize, and reword
jbr
parents:
5499
diff
changeset
|
945 /* Bit allocation */ |
7026 | 946 for(ch=!cpl_in_use; ch<=s->channels; ch++) { |
5476
e50f0d583a36
consolidate decoding of lfe and coupling channels with decoding of full-bandwidth channels.
jbr
parents:
5475
diff
changeset
|
947 if(bit_alloc_stages[ch] > 2) { |
e50f0d583a36
consolidate decoding of lfe and coupling channels with decoding of full-bandwidth channels.
jbr
parents:
5475
diff
changeset
|
948 /* Exponent mapping into PSD and PSD integration */ |
6090
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
949 ff_ac3_bit_alloc_calc_psd(s->dexps[ch], |
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
950 s->start_freq[ch], s->end_freq[ch], |
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
951 s->psd[ch], s->band_psd[ch]); |
5461
c427d8950769
cosmetics: use braces when it makes the code easier to read.
jbr
parents:
5460
diff
changeset
|
952 } |
5476
e50f0d583a36
consolidate decoding of lfe and coupling channels with decoding of full-bandwidth channels.
jbr
parents:
5475
diff
changeset
|
953 if(bit_alloc_stages[ch] > 1) { |
e50f0d583a36
consolidate decoding of lfe and coupling channels with decoding of full-bandwidth channels.
jbr
parents:
5475
diff
changeset
|
954 /* Compute excitation function, Compute masking curve, and |
e50f0d583a36
consolidate decoding of lfe and coupling channels with decoding of full-bandwidth channels.
jbr
parents:
5475
diff
changeset
|
955 Apply delta bit allocation */ |
6090
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
956 ff_ac3_bit_alloc_calc_mask(&s->bit_alloc_params, s->band_psd[ch], |
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
957 s->start_freq[ch], s->end_freq[ch], |
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
958 s->fast_gain[ch], (ch == s->lfe_ch), |
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
959 s->dba_mode[ch], s->dba_nsegs[ch], |
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
960 s->dba_offsets[ch], s->dba_lengths[ch], |
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
961 s->dba_values[ch], s->mask[ch]); |
5461
c427d8950769
cosmetics: use braces when it makes the code easier to read.
jbr
parents:
5460
diff
changeset
|
962 } |
5476
e50f0d583a36
consolidate decoding of lfe and coupling channels with decoding of full-bandwidth channels.
jbr
parents:
5475
diff
changeset
|
963 if(bit_alloc_stages[ch] > 0) { |
e50f0d583a36
consolidate decoding of lfe and coupling channels with decoding of full-bandwidth channels.
jbr
parents:
5475
diff
changeset
|
964 /* Compute bit allocation */ |
6090
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
965 ff_ac3_bit_alloc_calc_bap(s->mask[ch], s->psd[ch], |
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
966 s->start_freq[ch], s->end_freq[ch], |
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
967 s->snr_offset[ch], |
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
968 s->bit_alloc_params.floor, |
7017 | 969 ff_ac3_bap_tab, s->bap[ch]); |
5461
c427d8950769
cosmetics: use braces when it makes the code easier to read.
jbr
parents:
5460
diff
changeset
|
970 } |
5310
9aa9197034d7
AC-3 decoder, soc revision 40, Aug 9 00:10:14 2006 UTC by cloud9
jbr
parents:
5309
diff
changeset
|
971 } |
5309
0662a270aab7
AC-3 decoder, soc revision 38, Aug 7 00:03:00 2006 UTC by cloud9
jbr
parents:
5308
diff
changeset
|
972 |
5515
b499866a3397
comments modification: add, clean up, Doxygenize, and reword
jbr
parents:
5499
diff
changeset
|
973 /* unused dummy data */ |
6089 | 974 if (get_bits1(gbc)) { |
975 int skipl = get_bits(gbc, 9); | |
5309
0662a270aab7
AC-3 decoder, soc revision 38, Aug 7 00:03:00 2006 UTC by cloud9
jbr
parents:
5308
diff
changeset
|
976 while(skipl--) |
6089 | 977 skip_bits(gbc, 8); |
5304
eff51058fe13
AC-3 decoder, soc revision 26, Jul 5 04:55:15 2006 UTC by cloud9
jbr
parents:
5303
diff
changeset
|
978 } |
5499 | 979 |
5304
eff51058fe13
AC-3 decoder, soc revision 26, Jul 5 04:55:15 2006 UTC by cloud9
jbr
parents:
5303
diff
changeset
|
980 /* unpack the transform coefficients |
5515
b499866a3397
comments modification: add, clean up, Doxygenize, and reword
jbr
parents:
5499
diff
changeset
|
981 this also uncouples channels if coupling is in use. */ |
6953
64504967b891
get_transform_coeffs() never returns an error, so make the function and its children return void.
jbr
parents:
6951
diff
changeset
|
982 get_transform_coeffs(s); |
5306
abc5c130b448
AC-3 decoder, soc revision 32, Jul 17 09:37:32 2006 UTC by cloud9
jbr
parents:
5305
diff
changeset
|
983 |
5304
eff51058fe13
AC-3 decoder, soc revision 26, Jul 5 04:55:15 2006 UTC by cloud9
jbr
parents:
5303
diff
changeset
|
984 /* recover coefficients if rematrixing is in use */ |
6090
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
985 if(s->channel_mode == AC3_CHMODE_STEREO) |
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
986 do_rematrixing(s); |
5304
eff51058fe13
AC-3 decoder, soc revision 26, Jul 5 04:55:15 2006 UTC by cloud9
jbr
parents:
5303
diff
changeset
|
987 |
6032
f74202e7e896
Skip dialog normalization. It shouldn't be used by the decoder, and there is
jbr
parents:
6031
diff
changeset
|
988 /* apply scaling to coefficients (headroom, dynrng) */ |
6090
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
989 for(ch=1; ch<=s->channels; ch++) { |
6430
de7502093922
Use 24-bit fixed-point transform coefficients until just before MDCT. This
jbr
parents:
6142
diff
changeset
|
990 float gain = s->mul_bias / 4194304.0f; |
6090
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
991 if(s->channel_mode == AC3_CHMODE_DUALMONO) { |
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
992 gain *= s->dynamic_range[ch-1]; |
5412
3f6a6dda6b2d
remove broken downmixing. will add new implementation later.
jbr
parents:
5393
diff
changeset
|
993 } else { |
6090
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
994 gain *= s->dynamic_range[0]; |
5412
3f6a6dda6b2d
remove broken downmixing. will add new implementation later.
jbr
parents:
5393
diff
changeset
|
995 } |
7564 | 996 s->dsp.int32_to_float_fmul_scalar(s->transform_coeffs[ch], s->fixed_coeffs[ch], gain, 256); |
5412
3f6a6dda6b2d
remove broken downmixing. will add new implementation later.
jbr
parents:
5393
diff
changeset
|
997 } |
5318
ef4ef249ca72
AC-3 decoder, soc revision 58, Aug 19 14:20:45 2006 UTC by cloud9
jbr
parents:
5317
diff
changeset
|
998 |
6480
6f01a499e785
downmix before imdct unless different size transforms are used. about 20%
jbr
parents:
6432
diff
changeset
|
999 /* downmix and MDCT. order depends on whether block switching is used for |
6f01a499e785
downmix before imdct unless different size transforms are used. about 20%
jbr
parents:
6432
diff
changeset
|
1000 any channel in this block. this is because coefficients for the long |
6f01a499e785
downmix before imdct unless different size transforms are used. about 20%
jbr
parents:
6432
diff
changeset
|
1001 and short transforms cannot be mixed. */ |
6f01a499e785
downmix before imdct unless different size transforms are used. about 20%
jbr
parents:
6432
diff
changeset
|
1002 downmix_output = s->channels != s->out_channels && |
6f01a499e785
downmix before imdct unless different size transforms are used. about 20%
jbr
parents:
6432
diff
changeset
|
1003 !((s->output_mode & AC3_OUTPUT_LFEON) && |
6f01a499e785
downmix before imdct unless different size transforms are used. about 20%
jbr
parents:
6432
diff
changeset
|
1004 s->fbw_channels == s->out_channels); |
6f01a499e785
downmix before imdct unless different size transforms are used. about 20%
jbr
parents:
6432
diff
changeset
|
1005 if(different_transforms) { |
6f01a499e785
downmix before imdct unless different size transforms are used. about 20%
jbr
parents:
6432
diff
changeset
|
1006 /* the delay samples have already been downmixed, so we upmix the delay |
6f01a499e785
downmix before imdct unless different size transforms are used. about 20%
jbr
parents:
6432
diff
changeset
|
1007 samples in order to reconstruct all channels before downmixing. */ |
6f01a499e785
downmix before imdct unless different size transforms are used. about 20%
jbr
parents:
6432
diff
changeset
|
1008 if(s->downmixed) { |
6f01a499e785
downmix before imdct unless different size transforms are used. about 20%
jbr
parents:
6432
diff
changeset
|
1009 s->downmixed = 0; |
6f01a499e785
downmix before imdct unless different size transforms are used. about 20%
jbr
parents:
6432
diff
changeset
|
1010 ac3_upmix_delay(s); |
6f01a499e785
downmix before imdct unless different size transforms are used. about 20%
jbr
parents:
6432
diff
changeset
|
1011 } |
5306
abc5c130b448
AC-3 decoder, soc revision 32, Jul 17 09:37:32 2006 UTC by cloud9
jbr
parents:
5305
diff
changeset
|
1012 |
6480
6f01a499e785
downmix before imdct unless different size transforms are used. about 20%
jbr
parents:
6432
diff
changeset
|
1013 do_imdct(s, s->channels); |
6f01a499e785
downmix before imdct unless different size transforms are used. about 20%
jbr
parents:
6432
diff
changeset
|
1014 |
6f01a499e785
downmix before imdct unless different size transforms are used. about 20%
jbr
parents:
6432
diff
changeset
|
1015 if(downmix_output) { |
7563 | 1016 s->dsp.ac3_downmix(s->output, s->downmix_coeffs, s->out_channels, s->fbw_channels, 256); |
6480
6f01a499e785
downmix before imdct unless different size transforms are used. about 20%
jbr
parents:
6432
diff
changeset
|
1017 } |
6f01a499e785
downmix before imdct unless different size transforms are used. about 20%
jbr
parents:
6432
diff
changeset
|
1018 } else { |
6f01a499e785
downmix before imdct unless different size transforms are used. about 20%
jbr
parents:
6432
diff
changeset
|
1019 if(downmix_output) { |
7563 | 1020 s->dsp.ac3_downmix(s->transform_coeffs+1, s->downmix_coeffs, s->out_channels, s->fbw_channels, 256); |
6480
6f01a499e785
downmix before imdct unless different size transforms are used. about 20%
jbr
parents:
6432
diff
changeset
|
1021 } |
6f01a499e785
downmix before imdct unless different size transforms are used. about 20%
jbr
parents:
6432
diff
changeset
|
1022 |
7562
ef456ee01ea2
don't mark the delayed samples for upmixing if they haven't been downmixed
lorenm
parents:
7561
diff
changeset
|
1023 if(downmix_output && !s->downmixed) { |
6480
6f01a499e785
downmix before imdct unless different size transforms are used. about 20%
jbr
parents:
6432
diff
changeset
|
1024 s->downmixed = 1; |
7563 | 1025 s->dsp.ac3_downmix(s->delay, s->downmix_coeffs, s->out_channels, s->fbw_channels, 128); |
6480
6f01a499e785
downmix before imdct unless different size transforms are used. about 20%
jbr
parents:
6432
diff
changeset
|
1026 } |
6f01a499e785
downmix before imdct unless different size transforms are used. about 20%
jbr
parents:
6432
diff
changeset
|
1027 |
6f01a499e785
downmix before imdct unless different size transforms are used. about 20%
jbr
parents:
6432
diff
changeset
|
1028 do_imdct(s, s->out_channels); |
5480 | 1029 } |
1030 | |
5479
943c732c905d
use dsputil for float to signed 16-bit sample conversion
jbr
parents:
5477
diff
changeset
|
1031 return 0; |
5306
abc5c130b448
AC-3 decoder, soc revision 32, Jul 17 09:37:32 2006 UTC by cloud9
jbr
parents:
5305
diff
changeset
|
1032 } |
abc5c130b448
AC-3 decoder, soc revision 32, Jul 17 09:37:32 2006 UTC by cloud9
jbr
parents:
5305
diff
changeset
|
1033 |
5515
b499866a3397
comments modification: add, clean up, Doxygenize, and reword
jbr
parents:
5499
diff
changeset
|
1034 /** |
b499866a3397
comments modification: add, clean up, Doxygenize, and reword
jbr
parents:
5499
diff
changeset
|
1035 * Decode a single AC-3 frame. |
5321
fcb2a52c88fe
AC-3 decoder, soc revision 71, Sep 4 09:47:12 2006 UTC by cloud9
jbr
parents:
5320
diff
changeset
|
1036 */ |
6487
8f1b27d73a07
make input buffer const, as it should be. fixes a warning.
jbr
parents:
6480
diff
changeset
|
1037 static int ac3_decode_frame(AVCodecContext * avctx, void *data, int *data_size, |
8f1b27d73a07
make input buffer const, as it should be. fixes a warning.
jbr
parents:
6480
diff
changeset
|
1038 const uint8_t *buf, int buf_size) |
5304
eff51058fe13
AC-3 decoder, soc revision 26, Jul 5 04:55:15 2006 UTC by cloud9
jbr
parents:
5303
diff
changeset
|
1039 { |
6092 | 1040 AC3DecodeContext *s = avctx->priv_data; |
5305
5892b4a6380b
AC-3 decoder, soc revision 31, Jul 14 23:53:28 2006 UTC by cloud9
jbr
parents:
5304
diff
changeset
|
1041 int16_t *out_samples = (int16_t *)data; |
7550 | 1042 int blk, ch, err; |
5318
ef4ef249ca72
AC-3 decoder, soc revision 58, Aug 19 14:20:45 2006 UTC by cloud9
jbr
parents:
5317
diff
changeset
|
1043 |
5515
b499866a3397
comments modification: add, clean up, Doxygenize, and reword
jbr
parents:
5499
diff
changeset
|
1044 /* initialize the GetBitContext with the start of valid AC-3 Frame */ |
6538
b8f886bce6f5
only allocate context input buffer if AVCodecContext.error_reslience is greater than 0.
jbr
parents:
6537
diff
changeset
|
1045 if (s->input_buffer) { |
6536
8f2186d5daca
additional protection from segmentation faults and memory access errors by
jbr
parents:
6529
diff
changeset
|
1046 /* copy input buffer to decoder context to avoid reading past the end |
8f2186d5daca
additional protection from segmentation faults and memory access errors by
jbr
parents:
6529
diff
changeset
|
1047 of the buffer, which can be caused by a damaged input stream. */ |
8f2186d5daca
additional protection from segmentation faults and memory access errors by
jbr
parents:
6529
diff
changeset
|
1048 memcpy(s->input_buffer, buf, FFMIN(buf_size, AC3_MAX_FRAME_SIZE)); |
8f2186d5daca
additional protection from segmentation faults and memory access errors by
jbr
parents:
6529
diff
changeset
|
1049 init_get_bits(&s->gbc, s->input_buffer, buf_size * 8); |
8f2186d5daca
additional protection from segmentation faults and memory access errors by
jbr
parents:
6529
diff
changeset
|
1050 } else { |
6537 | 1051 init_get_bits(&s->gbc, buf, buf_size * 8); |
6536
8f2186d5daca
additional protection from segmentation faults and memory access errors by
jbr
parents:
6529
diff
changeset
|
1052 } |
5305
5892b4a6380b
AC-3 decoder, soc revision 31, Jul 14 23:53:28 2006 UTC by cloud9
jbr
parents:
5304
diff
changeset
|
1053 |
5515
b499866a3397
comments modification: add, clean up, Doxygenize, and reword
jbr
parents:
5499
diff
changeset
|
1054 /* parse the syncinfo */ |
6950 | 1055 *data_size = 0; |
7028 | 1056 err = parse_frame_header(s); |
6950 | 1057 |
1058 /* check that reported frame size fits in input buffer */ | |
1059 if(s->frame_size > buf_size) { | |
1060 av_log(avctx, AV_LOG_ERROR, "incomplete frame\n"); | |
1061 err = AC3_PARSE_ERROR_FRAME_SIZE; | |
1062 } | |
1063 | |
1064 /* check for crc mismatch */ | |
1065 if(err != AC3_PARSE_ERROR_FRAME_SIZE && avctx->error_resilience >= FF_ER_CAREFUL) { | |
1066 if(av_crc(av_crc_get_table(AV_CRC_16_ANSI), 0, &buf[2], s->frame_size-2)) { | |
1067 av_log(avctx, AV_LOG_ERROR, "frame CRC mismatch\n"); | |
6951 | 1068 err = AC3_PARSE_ERROR_CRC; |
6950 | 1069 } |
1070 } | |
1071 | |
6951 | 1072 if(err && err != AC3_PARSE_ERROR_CRC) { |
5680 | 1073 switch(err) { |
1074 case AC3_PARSE_ERROR_SYNC: | |
6950 | 1075 av_log(avctx, AV_LOG_ERROR, "frame sync error\n"); |
6999 | 1076 return -1; |
5680 | 1077 case AC3_PARSE_ERROR_BSID: |
1078 av_log(avctx, AV_LOG_ERROR, "invalid bitstream id\n"); | |
1079 break; | |
1080 case AC3_PARSE_ERROR_SAMPLE_RATE: | |
1081 av_log(avctx, AV_LOG_ERROR, "invalid sample rate\n"); | |
1082 break; | |
1083 case AC3_PARSE_ERROR_FRAME_SIZE: | |
1084 av_log(avctx, AV_LOG_ERROR, "invalid frame size\n"); | |
1085 break; | |
6540
b0d44aec1ec0
change name from stream type to frame type in AC3 code
bwolowiec
parents:
6538
diff
changeset
|
1086 case AC3_PARSE_ERROR_FRAME_TYPE: |
7014 | 1087 /* skip frame if CRC is ok. otherwise use error concealment. */ |
1088 /* TODO: add support for substreams and dependent frames */ | |
1089 if(s->frame_type == EAC3_FRAME_TYPE_DEPENDENT || s->substreamid) { | |
1090 av_log(avctx, AV_LOG_ERROR, "unsupported frame type : skipping frame\n"); | |
1091 return s->frame_size; | |
1092 } else { | |
7015 | 1093 av_log(avctx, AV_LOG_ERROR, "invalid frame type\n"); |
7014 | 1094 } |
6529 | 1095 break; |
5680 | 1096 default: |
1097 av_log(avctx, AV_LOG_ERROR, "invalid header\n"); | |
1098 break; | |
1099 } | |
5304
eff51058fe13
AC-3 decoder, soc revision 26, Jul 5 04:55:15 2006 UTC by cloud9
jbr
parents:
5303
diff
changeset
|
1100 } |
5303
b8821ff5c30d
AC-3 decoder, soc revision 9, Jun 14 03:52:02 2006 UTC by cloud9
jbr
parents:
diff
changeset
|
1101 |
6889 | 1102 /* if frame is ok, set audio parameters */ |
1103 if (!err) { | |
6885 | 1104 avctx->sample_rate = s->sample_rate; |
1105 avctx->bit_rate = s->bit_rate; | |
6098 | 1106 |
6885 | 1107 /* channel config */ |
1108 s->out_channels = s->channels; | |
6889 | 1109 s->output_mode = s->channel_mode; |
1110 if(s->lfe_on) | |
1111 s->output_mode |= AC3_OUTPUT_LFEON; | |
6885 | 1112 if (avctx->request_channels > 0 && avctx->request_channels <= 2 && |
1113 avctx->request_channels < s->channels) { | |
1114 s->out_channels = avctx->request_channels; | |
1115 s->output_mode = avctx->request_channels == 1 ? AC3_CHMODE_MONO : AC3_CHMODE_STEREO; | |
1116 } | |
1117 avctx->channels = s->out_channels; | |
5303
b8821ff5c30d
AC-3 decoder, soc revision 9, Jun 14 03:52:02 2006 UTC by cloud9
jbr
parents:
diff
changeset
|
1118 |
6885 | 1119 /* set downmixing coefficients if needed */ |
1120 if(s->channels != s->out_channels && !((s->output_mode & AC3_OUTPUT_LFEON) && | |
1121 s->fbw_channels == s->out_channels)) { | |
1122 set_downmix_coeffs(s); | |
1123 } | |
6889 | 1124 } else if (!s->out_channels) { |
1125 s->out_channels = avctx->channels; | |
1126 if(s->out_channels < s->channels) | |
1127 s->output_mode = s->out_channels == 1 ? AC3_CHMODE_MONO : AC3_CHMODE_STEREO; | |
6884 | 1128 } |
6126
c24cd5282108
split out init of downmix coeffs into a separate function that can be called
jbr
parents:
6125
diff
changeset
|
1129 |
7490
bddd15a1d167
cosmetics: change function name and comments to refer to block decoding
jbr
parents:
7489
diff
changeset
|
1130 /* decode the audio blocks */ |
7011
77f607fb4e8b
get the number of blocks from the ac3 parser and use in the ac3 decoder.
jbr
parents:
6999
diff
changeset
|
1131 for (blk = 0; blk < s->num_blocks; blk++) { |
7550 | 1132 const float *output[s->out_channels]; |
7490
bddd15a1d167
cosmetics: change function name and comments to refer to block decoding
jbr
parents:
7489
diff
changeset
|
1133 if (!err && decode_audio_block(s, blk)) { |
bddd15a1d167
cosmetics: change function name and comments to refer to block decoding
jbr
parents:
7489
diff
changeset
|
1134 av_log(avctx, AV_LOG_ERROR, "error decoding the audio block\n"); |
5304
eff51058fe13
AC-3 decoder, soc revision 26, Jul 5 04:55:15 2006 UTC by cloud9
jbr
parents:
5303
diff
changeset
|
1135 } |
7550 | 1136 for (ch = 0; ch < s->out_channels; ch++) |
1137 output[ch] = s->output[ch]; | |
1138 s->dsp.float_to_int16_interleave(out_samples, output, 256, s->out_channels); | |
1139 out_samples += 256 * s->out_channels; | |
5304
eff51058fe13
AC-3 decoder, soc revision 26, Jul 5 04:55:15 2006 UTC by cloud9
jbr
parents:
5303
diff
changeset
|
1140 } |
7011
77f607fb4e8b
get the number of blocks from the ac3 parser and use in the ac3 decoder.
jbr
parents:
6999
diff
changeset
|
1141 *data_size = s->num_blocks * 256 * avctx->channels * sizeof (int16_t); |
6090
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
1142 return s->frame_size; |
5303
b8821ff5c30d
AC-3 decoder, soc revision 9, Jun 14 03:52:02 2006 UTC by cloud9
jbr
parents:
diff
changeset
|
1143 } |
5304
eff51058fe13
AC-3 decoder, soc revision 26, Jul 5 04:55:15 2006 UTC by cloud9
jbr
parents:
5303
diff
changeset
|
1144 |
5515
b499866a3397
comments modification: add, clean up, Doxygenize, and reword
jbr
parents:
5499
diff
changeset
|
1145 /** |
b499866a3397
comments modification: add, clean up, Doxygenize, and reword
jbr
parents:
5499
diff
changeset
|
1146 * Uninitialize the AC-3 decoder. |
5321
fcb2a52c88fe
AC-3 decoder, soc revision 71, Sep 4 09:47:12 2006 UTC by cloud9
jbr
parents:
5320
diff
changeset
|
1147 */ |
6517
48759bfbd073
Apply 'cold' attribute to init/uninit functions in libavcodec
zuxy
parents:
6487
diff
changeset
|
1148 static av_cold int ac3_decode_end(AVCodecContext *avctx) |
5304
eff51058fe13
AC-3 decoder, soc revision 26, Jul 5 04:55:15 2006 UTC by cloud9
jbr
parents:
5303
diff
changeset
|
1149 { |
6092 | 1150 AC3DecodeContext *s = avctx->priv_data; |
6090
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
1151 ff_mdct_end(&s->imdct_512); |
b7b19fbc4746
cosmetics: rename all AC3DecodeContext variables from ctx to s
jbr
parents:
6089
diff
changeset
|
1152 ff_mdct_end(&s->imdct_256); |
5321
fcb2a52c88fe
AC-3 decoder, soc revision 71, Sep 4 09:47:12 2006 UTC by cloud9
jbr
parents:
5320
diff
changeset
|
1153 |
6538
b8f886bce6f5
only allocate context input buffer if AVCodecContext.error_reslience is greater than 0.
jbr
parents:
6537
diff
changeset
|
1154 av_freep(&s->input_buffer); |
b8f886bce6f5
only allocate context input buffer if AVCodecContext.error_reslience is greater than 0.
jbr
parents:
6537
diff
changeset
|
1155 |
5304
eff51058fe13
AC-3 decoder, soc revision 26, Jul 5 04:55:15 2006 UTC by cloud9
jbr
parents:
5303
diff
changeset
|
1156 return 0; |
eff51058fe13
AC-3 decoder, soc revision 26, Jul 5 04:55:15 2006 UTC by cloud9
jbr
parents:
5303
diff
changeset
|
1157 } |
eff51058fe13
AC-3 decoder, soc revision 26, Jul 5 04:55:15 2006 UTC by cloud9
jbr
parents:
5303
diff
changeset
|
1158 |
5334
d406baea3275
no need for the decoder to have the license as a prefix
jbr
parents:
5333
diff
changeset
|
1159 AVCodec ac3_decoder = { |
5315
31033309ebc9
AC-3 decoder, soc revision 55, Aug 17 12:25:24 2006 UTC by cloud9
jbr
parents:
5314
diff
changeset
|
1160 .name = "ac3", |
31033309ebc9
AC-3 decoder, soc revision 55, Aug 17 12:25:24 2006 UTC by cloud9
jbr
parents:
5314
diff
changeset
|
1161 .type = CODEC_TYPE_AUDIO, |
31033309ebc9
AC-3 decoder, soc revision 55, Aug 17 12:25:24 2006 UTC by cloud9
jbr
parents:
5314
diff
changeset
|
1162 .id = CODEC_ID_AC3, |
31033309ebc9
AC-3 decoder, soc revision 55, Aug 17 12:25:24 2006 UTC by cloud9
jbr
parents:
5314
diff
changeset
|
1163 .priv_data_size = sizeof (AC3DecodeContext), |
31033309ebc9
AC-3 decoder, soc revision 55, Aug 17 12:25:24 2006 UTC by cloud9
jbr
parents:
5314
diff
changeset
|
1164 .init = ac3_decode_init, |
31033309ebc9
AC-3 decoder, soc revision 55, Aug 17 12:25:24 2006 UTC by cloud9
jbr
parents:
5314
diff
changeset
|
1165 .close = ac3_decode_end, |
31033309ebc9
AC-3 decoder, soc revision 55, Aug 17 12:25:24 2006 UTC by cloud9
jbr
parents:
5314
diff
changeset
|
1166 .decode = ac3_decode_frame, |
7040
e943e1409077
Make AVCodec long_names definition conditional depending on CONFIG_SMALL.
stefano
parents:
7033
diff
changeset
|
1167 .long_name = NULL_IF_CONFIG_SMALL("ATSC A/52 / AC-3"), |
5304
eff51058fe13
AC-3 decoder, soc revision 26, Jul 5 04:55:15 2006 UTC by cloud9
jbr
parents:
5303
diff
changeset
|
1168 }; |