Mercurial > libavcodec.hg
annotate qcelpdata.h @ 11352:6e0af2cfdcfe libavcodec
Do MC and IDCT in coding (hilbert) order
This increases the slice size to 64 pixels, due to having to decode an
entire chroma superblock row per slice.
This can be up to 6% slower depending on clip and CPU, but is necessary
for future optimizations that gain significantly more than was lost.
author | conrad |
---|---|
date | Wed, 03 Mar 2010 23:27:40 +0000 |
parents | 3c98f3e1b719 |
children | 7dd2a45249a9 |
rev | line source |
---|---|
8081 | 1 /* |
8147
3b9eefe120e6
Please keep original copyright heading and attribution lines...
reynaldo
parents:
8145
diff
changeset
|
2 * QCELP decoder |
8081 | 3 * Copyright (c) 2007 Reynaldo H. Verdejo Pinochet |
4 * | |
5 * This file is part of FFmpeg. | |
6 * | |
7 * FFmpeg is free software; you can redistribute it and/or | |
8 * modify it under the terms of the GNU Lesser General Public | |
9 * License as published by the Free Software Foundation; either | |
10 * version 2.1 of the License, or (at your option) any later version. | |
11 * | |
12 * FFmpeg is distributed in the hope that it will be useful, | |
13 * but WITHOUT ANY WARRANTY; without even the implied warranty of | |
14 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU | |
15 * Lesser General Public License for more details. | |
16 * | |
17 * You should have received a copy of the GNU Lesser General Public | |
18 * License along with FFmpeg; if not, write to the Free Software | |
19 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA | |
20 */ | |
21 | |
22 #ifndef AVCODEC_QCELPDATA_H | |
23 #define AVCODEC_QCELPDATA_H | |
24 | |
8147
3b9eefe120e6
Please keep original copyright heading and attribution lines...
reynaldo
parents:
8145
diff
changeset
|
25 /** |
8718
e9d9d946f213
Use full internal pathname in doxygen @file directives.
diego
parents:
8258
diff
changeset
|
26 * @file libavcodec/qcelpdata.h |
8149 | 27 * Data tables for the QCELP decoder |
8147
3b9eefe120e6
Please keep original copyright heading and attribution lines...
reynaldo
parents:
8145
diff
changeset
|
28 * @author Reynaldo H. Verdejo Pinochet |
8151
e2c068cb210a
Credit Kenan Gillet for his contributions towards merging the SoC QCELP decoder.
reynaldo
parents:
8149
diff
changeset
|
29 * @remark FFmpeg merging spearheaded by Kenan Gillet |
8258
b41482ad0ef5
COSMETICS, add missing remarks crediting Ben and Kenan
reynaldo
parents:
8237
diff
changeset
|
30 * @remark Development mentored by Benjamin Larson |
8147
3b9eefe120e6
Please keep original copyright heading and attribution lines...
reynaldo
parents:
8145
diff
changeset
|
31 */ |
3b9eefe120e6
Please keep original copyright heading and attribution lines...
reynaldo
parents:
8145
diff
changeset
|
32 |
8212
c245082d70f8
Add missing headers to get 'make checkheaders' closer to passing.
diego
parents:
8191
diff
changeset
|
33 #include <stddef.h> |
8081 | 34 #include <stdint.h> |
8212
c245082d70f8
Add missing headers to get 'make checkheaders' closer to passing.
diego
parents:
8191
diff
changeset
|
35 #include "libavutil/common.h" |
8081 | 36 |
8089 | 37 /** |
8230 | 38 * QCELP unpacked data frame |
39 */ | |
40 typedef struct { | |
41 /// @defgroup qcelp_codebook_parameters QCELP excitation codebook parameters | |
42 /// @{ | |
43 uint8_t cbsign[16]; ///!< sign of the codebook gain for each codebook subframe | |
44 uint8_t cbgain[16]; ///!< unsigned codebook gain for each codebook subframe | |
45 uint8_t cindex[16]; ///!< codebook index for each codebook subframe | |
46 /// @} | |
47 | |
48 /// @defgroup qcelp_pitch_parameters QCELP pitch prediction parameters | |
49 /// @{ | |
50 uint8_t plag[4]; ///!< pitch lag for each pitch subframe | |
51 uint8_t pfrac[4]; ///!< fractional pitch lag for each pitch subframe | |
52 uint8_t pgain[4]; ///!< pitch gain for each pitch subframe | |
53 /// @} | |
54 | |
55 /** | |
56 * line spectral pair frequencies (LSP) for RATE_OCTAVE, | |
57 * line spectral pair frequencies grouped into five vectors | |
58 * of dimension two (LSPV) for other rates | |
59 */ | |
60 uint8_t lspv[10]; | |
61 | |
8236 | 62 /** |
63 * reserved bits only present in bitrate 1, 1/4 and 1/8 packets | |
64 */ | |
65 uint8_t reserved; | |
8230 | 66 } QCELPFrame; |
67 | |
68 /** | |
8089 | 69 * pre-calculated table for hammsinc function |
8092 | 70 * Only half of the table is needed because of symmetry. |
8089 | 71 * |
72 * TIA/EIA/IS-733 2.4.5.2-2/3 | |
73 */ | |
74 static const float qcelp_hammsinc_table[4] = { -0.006822, 0.041249, -0.143459, 0.588863}; | |
75 | |
76 typedef struct { | |
77 uint8_t index; /*!< index into the QCELPContext structure */ | |
8092 | 78 uint8_t bitpos; /*!< position of the lowest bit in the value's byte */ |
8089 | 79 uint8_t bitlen; /*!< number of bits to read */ |
80 } QCELPBitmap; | |
81 | |
8230 | 82 #define QCELP_OF(variable, bit, len) {offsetof(QCELPFrame, variable), bit, len} |
8089 | 83 |
84 /** | |
85 * bitmap unpacking tables for RATE_FULL | |
86 * | |
87 * TIA/EIA/IS-733 Table 2.4.7.1-1 | |
88 */ | |
89 static const QCELPBitmap qcelp_rate_full_bitmap[] = { | |
90 // start on bit | |
91 QCELP_OF(lspv [ 2], 0, 3), // 265 | |
92 QCELP_OF(lspv [ 1], 0, 7), // 262 | |
93 QCELP_OF(lspv [ 0], 0, 6), // 255 | |
94 QCELP_OF(lspv [ 4], 0, 6), // 249 | |
95 QCELP_OF(lspv [ 3], 0, 6), // 243 | |
96 QCELP_OF(lspv [ 2], 3, 4), // 237 | |
97 QCELP_OF(cbsign[ 0], 0, 1), // 233 | |
98 QCELP_OF(cbgain[ 0], 0, 4), // 232 | |
99 QCELP_OF(pfrac [ 0], 0, 1), // 228 | |
100 QCELP_OF(plag [ 0], 0, 7), // 227 | |
101 QCELP_OF(pgain [ 0], 0, 3), // 220 | |
102 QCELP_OF(cindex[ 1], 0, 4), // 217 | |
103 QCELP_OF(cbsign[ 1], 0, 1), // 213 | |
104 QCELP_OF(cbgain[ 1], 0, 4), // 212 | |
105 QCELP_OF(cindex[ 0], 0, 7), // 208 | |
106 QCELP_OF(cbgain[ 3], 0, 1), // 201 | |
107 QCELP_OF(cindex[ 2], 0, 7), // 200 | |
108 QCELP_OF(cbsign[ 2], 0, 1), // 193 | |
109 QCELP_OF(cbgain[ 2], 0, 4), // 192 | |
110 QCELP_OF(cindex[ 1], 4, 3), // 188 | |
111 QCELP_OF(plag [ 1], 0, 3), // 185 | |
112 QCELP_OF(pgain [ 1], 0, 3), // 182 | |
113 QCELP_OF(cindex[ 3], 0, 7), // 179 | |
114 QCELP_OF(cbsign[ 3], 0, 1), // 172 | |
115 QCELP_OF(cbgain[ 3], 1, 2), // 171 | |
116 QCELP_OF(cindex[ 4], 0, 6), // 169 | |
117 QCELP_OF(cbsign[ 4], 0, 1), // 163 | |
118 QCELP_OF(cbgain[ 4], 0, 4), // 162 | |
119 QCELP_OF(pfrac [ 1], 0, 1), // 158 | |
120 QCELP_OF(plag [ 1], 3, 4), // 157 | |
121 QCELP_OF(cbgain[ 6], 0, 3), // 153 | |
122 QCELP_OF(cindex[ 5], 0, 7), // 150 | |
123 QCELP_OF(cbsign[ 5], 0, 1), // 143 | |
124 QCELP_OF(cbgain[ 5], 0, 4), // 142 | |
125 QCELP_OF(cindex[ 4], 6, 1), // 138 | |
126 QCELP_OF(cindex[ 7], 0, 3), // 137 | |
127 QCELP_OF(cbsign[ 7], 0, 1), // 134 | |
128 QCELP_OF(cbgain[ 7], 0, 3), // 133 | |
129 QCELP_OF(cindex[ 6], 0, 7), // 130 | |
130 QCELP_OF(cbsign[ 6], 0, 1), // 123 | |
131 QCELP_OF(cbgain[ 6], 3, 1), // 122 | |
132 QCELP_OF(cbgain[ 8], 0, 1), // 121 | |
133 QCELP_OF(pfrac [ 2], 0, 1), // 120 | |
134 QCELP_OF(plag [ 2], 0, 7), // 119 | |
135 QCELP_OF(pgain [ 2], 0, 3), // 112 | |
136 QCELP_OF(cindex[ 7], 3, 4), // 109 | |
137 QCELP_OF(cbsign[ 9], 0, 1), // 105 | |
138 QCELP_OF(cbgain[ 9], 0, 4), // 104 | |
139 QCELP_OF(cindex[ 8], 0, 7), // 100 | |
140 QCELP_OF(cbsign[ 8], 0, 1), // 93 | |
141 QCELP_OF(cbgain[ 8], 1, 3), // 92 | |
142 QCELP_OF(cindex[10], 0, 4), // 89 | |
143 QCELP_OF(cbsign[10], 0, 1), // 85 | |
144 QCELP_OF(cbgain[10], 0, 4), // 84 | |
145 QCELP_OF(cindex[ 9], 0, 7), // 80 | |
146 QCELP_OF(pgain [ 3], 0, 2), // 73 | |
147 QCELP_OF(cindex[11], 0, 7), // 71 | |
148 QCELP_OF(cbsign[11], 0, 1), // 64 | |
149 QCELP_OF(cbgain[11], 0, 3), // 63 | |
150 QCELP_OF(cindex[10], 4, 3), // 60 | |
151 QCELP_OF(cindex[12], 0, 2), // 57 | |
152 QCELP_OF(cbsign[12], 0, 1), // 55 | |
153 QCELP_OF(cbgain[12], 0, 4), // 54 | |
154 QCELP_OF(pfrac [ 3], 0, 1), // 50 | |
155 QCELP_OF(plag [ 3], 0, 7), // 49 | |
156 QCELP_OF(pgain [ 3], 2, 1), // 42 | |
157 QCELP_OF(cindex[13], 0, 6), // 41 | |
158 QCELP_OF(cbsign[13], 0, 1), // 35 | |
159 QCELP_OF(cbgain[13], 0, 4), // 34 | |
160 QCELP_OF(cindex[12], 2, 5), // 30 | |
161 QCELP_OF(cbgain[15], 0, 3), // 25 | |
162 QCELP_OF(cindex[14], 0, 7), // 22 | |
163 QCELP_OF(cbsign[14], 0, 1), // 15 | |
164 QCELP_OF(cbgain[14], 0, 4), // 14 | |
165 QCELP_OF(cindex[13], 6, 1), // 10 | |
166 QCELP_OF(reserved, 0, 2), // 9 | |
167 QCELP_OF(cindex[15], 0, 7), // 7 | |
168 QCELP_OF(cbsign[15], 0, 1) // 0 | |
169 }; | |
170 | |
171 /** | |
172 * bitmap unpacking tables for RATE_HALF | |
173 * | |
174 * TIA/EIA/IS-733 Table 2.4.7.2-1 | |
175 */ | |
176 static const QCELPBitmap qcelp_rate_half_bitmap[] = { | |
177 // start on bit | |
178 QCELP_OF(lspv [2], 0, 3), // 123 | |
179 QCELP_OF(lspv [1], 0, 7), // 120 | |
180 QCELP_OF(lspv [0], 0, 6), // 113 | |
181 QCELP_OF(lspv [4], 0, 6), // 107 | |
182 QCELP_OF(lspv [3], 0, 6), // 101 | |
183 QCELP_OF(lspv [2], 3, 4), // 95 | |
184 QCELP_OF(cbsign[0], 0, 1), // 91 | |
185 QCELP_OF(cbgain[0], 0, 4), // 90 | |
186 QCELP_OF(pfrac [0], 0, 1), // 86 | |
187 QCELP_OF(plag [0], 0, 7), // 85 | |
188 QCELP_OF(pgain [0], 0, 3), // 78 | |
189 QCELP_OF(plag [1], 0, 6), // 75 | |
190 QCELP_OF(pgain [1], 0, 3), // 69 | |
191 QCELP_OF(cindex[0], 0, 7), // 66 | |
192 QCELP_OF(pgain [2], 0, 2), // 59 | |
193 QCELP_OF(cindex[1], 0, 7), // 57 | |
194 QCELP_OF(cbsign[1], 0, 1), // 50 | |
195 QCELP_OF(cbgain[1], 0, 4), // 49 | |
196 QCELP_OF(pfrac [1], 0, 1), // 45 | |
197 QCELP_OF(plag [1], 6, 1), // 44 | |
198 QCELP_OF(cindex[2], 0, 2), // 43 | |
199 QCELP_OF(cbsign[2], 0, 1), // 41 | |
200 QCELP_OF(cbgain[2], 0, 4), // 40 | |
201 QCELP_OF(pfrac [2], 0, 1), // 36 | |
202 QCELP_OF(plag [2], 0, 7), // 35 | |
203 QCELP_OF(pgain [2], 2, 1), // 28 | |
204 QCELP_OF(pfrac [3], 0, 1), // 27 | |
205 QCELP_OF(plag [3], 0, 7), // 26 | |
206 QCELP_OF(pgain [3], 0, 3), // 19 | |
207 QCELP_OF(cindex[2], 2, 5), // 16 | |
208 QCELP_OF(cindex[3], 0, 7), // 11 | |
209 QCELP_OF(cbsign[3], 0, 1), // 4 | |
210 QCELP_OF(cbgain[3], 0, 4) // 3 | |
211 }; | |
212 | |
213 /** | |
214 * bitmap unpacking tables for RATE_QUARTER | |
215 * | |
216 * TIA/EIA/IS-733 Table 2.4.7.3-1 | |
217 */ | |
218 static const QCELPBitmap qcelp_rate_quarter_bitmap[] = { | |
219 // start on bit | |
220 QCELP_OF(lspv [2], 0, 3), // 53 | |
221 QCELP_OF(lspv [1], 0, 7), // 50 | |
222 QCELP_OF(lspv [0], 0, 6), // 43 | |
223 QCELP_OF(lspv [4], 0, 6), // 37 | |
224 QCELP_OF(lspv [3], 0, 6), // 31 | |
225 QCELP_OF(lspv [2], 3, 4), // 25 | |
226 QCELP_OF(cbgain[3], 0, 4), // 21 | |
227 QCELP_OF(cbgain[2], 0, 4), // 17 | |
228 QCELP_OF(cbgain[1], 0, 4), // 13 | |
229 QCELP_OF(cbgain[0], 0, 4), // 9 | |
230 QCELP_OF(reserved, 0, 2), // 5 | |
231 QCELP_OF(cbgain[4], 0, 4) // 3 | |
232 }; | |
233 | |
234 /** | |
235 * bitmap unpacking tables for RATE_OCTAVE | |
236 * | |
237 * trick: CBSEED is written into QCELPContext.cbsign[15], | |
238 * which is not used for RATE_OCTAVE. | |
239 * CBSEED is only used to ensure the occurrence of random bit | |
240 * patterns in the 16 first bits that are used as the seed. | |
241 * | |
242 * TIA/EIA/IS-733 Table 2.4.7.4-1 | |
243 */ | |
244 static const QCELPBitmap qcelp_rate_octave_bitmap[] = { | |
245 // start on bit | |
246 QCELP_OF(cbsign[15], 3, 1), // 19 | |
247 QCELP_OF(lspv [0], 0, 1), // 18 | |
248 QCELP_OF(lspv [1], 0, 1), // 17 | |
249 QCELP_OF(lspv [2], 0, 1), // 16 | |
250 QCELP_OF(cbsign[15], 2, 1), // 15 | |
251 QCELP_OF(lspv [3], 0, 1), // 14 | |
252 QCELP_OF(lspv [4], 0, 1), // 13 | |
253 QCELP_OF(lspv [5], 0, 1), // 12 | |
254 QCELP_OF(cbsign[15], 1, 1), // 11 | |
255 QCELP_OF(lspv [6], 0, 1), // 10 | |
256 QCELP_OF(lspv [7], 0, 1), // 9 | |
257 QCELP_OF(lspv [8], 0, 1), // 8 | |
258 QCELP_OF(cbsign[15], 0, 1), // 7 | |
259 QCELP_OF(lspv [9], 0, 1), // 6 | |
260 QCELP_OF(cbgain [0], 0, 2), // 7 | |
261 QCELP_OF(reserved, 0, 4) // 3 | |
262 }; | |
263 | |
264 /** | |
265 * position of the bitmapping data for each packet type in | |
266 * the QCELPContext | |
267 */ | |
268 static const QCELPBitmap * const qcelp_unpacking_bitmaps_per_rate[5] = { | |
269 NULL, ///!< for SILENCE rate | |
270 qcelp_rate_octave_bitmap, | |
271 qcelp_rate_quarter_bitmap, | |
272 qcelp_rate_half_bitmap, | |
273 qcelp_rate_full_bitmap, | |
274 }; | |
275 | |
8230 | 276 static const uint16_t qcelp_unpacking_bitmaps_lengths[5] = { |
8089 | 277 0, ///!< for SILENCE rate |
278 FF_ARRAY_ELEMS(qcelp_rate_octave_bitmap), | |
279 FF_ARRAY_ELEMS(qcelp_rate_quarter_bitmap), | |
280 FF_ARRAY_ELEMS(qcelp_rate_half_bitmap), | |
281 FF_ARRAY_ELEMS(qcelp_rate_full_bitmap), | |
282 }; | |
283 | |
8081 | 284 typedef uint16_t qcelp_vector[2]; |
285 | |
286 /** | |
287 * LSP vector quantization tables in x*10000 form | |
288 * | |
289 * TIA/EIA/IS-733 tables 2.4.3.2.6.3-1 through 2.4.3.2.6.3-5 | |
290 */ | |
291 | |
292 static const qcelp_vector qcelp_lspvq1[64]= { | |
293 { 327, 118},{ 919, 111},{ 427, 440},{1327, 185}, | |
294 { 469, 50},{1272, 91},{ 892, 59},{1771, 193}, | |
295 { 222, 158},{1100, 127},{ 827, 55},{ 978, 791}, | |
296 { 665, 47},{ 700,1401},{ 670, 859},{1913,1048}, | |
297 { 471, 215},{1046, 125},{ 645, 298},{1599, 160}, | |
298 { 593, 39},{1187, 462},{ 749, 341},{1520, 511}, | |
299 { 290, 792},{ 909, 362},{ 753, 81},{1111,1058}, | |
300 { 519, 253},{ 828, 839},{ 685, 541},{1421,1258}, | |
301 { 386, 130},{ 962, 119},{ 542, 387},{1431, 185}, | |
302 { 526, 51},{1175, 260},{ 831, 167},{1728, 510}, | |
303 { 273, 437},{1172, 113},{ 771, 144},{1122, 751}, | |
304 { 619, 119},{ 492,1276},{ 658, 695},{1882, 615}, | |
305 { 415, 200},{1018, 88},{ 681, 339},{1436, 325}, | |
306 { 555, 122},{1042, 485},{ 826, 345},{1374, 743}, | |
307 { 383,1018},{1005, 358},{ 704, 86},{1301, 586}, | |
308 { 597, 241},{ 832, 621},{ 555, 573},{1504, 839}}; | |
309 | |
310 static const qcelp_vector qcelp_lspvq2[128]= { | |
311 { 255, 293},{ 904, 219},{ 151,1211},{1447, 498}, | |
312 { 470, 253},{1559, 177},{1547, 994},{2394, 242}, | |
313 { 91, 813},{ 857, 590},{ 934,1326},{1889, 282}, | |
314 { 813, 472},{1057,1494},{ 450,3315},{2163,1895}, | |
315 { 538, 532},{1399, 218},{ 146,1552},{1755, 626}, | |
316 { 822, 202},{1299, 663},{ 706,1732},{2656, 401}, | |
317 { 418, 745},{ 762,1038},{ 583,1748},{1746,1285}, | |
318 { 527,1169},{1314, 830},{ 556,2116},{1073,2321}, | |
319 { 297, 570},{ 981, 403},{ 468,1103},{1740, 243}, | |
320 { 725, 179},{1255, 474},{1374,1362},{1922, 912}, | |
321 { 285, 947},{ 930, 700},{ 593,1372},{1909, 576}, | |
322 { 588, 916},{1110,1116},{ 224,2719},{1633,2220}, | |
323 { 402, 520},{1061, 448},{ 402,1352},{1499, 775}, | |
324 { 664, 589},{1081, 727},{ 801,2206},{2165,1157}, | |
325 { 566, 802},{ 911,1116},{ 306,1703},{1792, 836}, | |
326 { 655, 999},{1061,1038},{ 298,2089},{1110,1753}, | |
327 { 361, 311},{ 970, 239},{ 265,1231},{1495, 573}, | |
328 { 566, 262},{1569, 293},{1341,1144},{2271, 544}, | |
329 { 214, 877},{ 847, 719},{ 794,1384},{2067, 274}, | |
330 { 703, 688},{1099,1306},{ 391,2947},{2024,1670}, | |
331 { 471, 525},{1245, 290},{ 264,1557},{1568, 807}, | |
332 { 718, 399},{1193, 685},{ 883,1594},{2729, 764}, | |
333 { 500, 754},{ 809,1108},{ 541,1648},{1523,1385}, | |
334 { 614,1196},{1209, 847},{ 345,2242},{1442,1747}, | |
335 { 199, 560},{1092, 194},{ 349,1253},{1653, 507}, | |
336 { 625, 354},{1376, 431},{1187,1465},{2164, 872}, | |
337 { 360, 974},{1008, 698},{ 704,1346},{2114, 452}, | |
338 { 720, 816},{1240,1089},{ 439,2475},{1498,2040}, | |
339 { 336, 718},{1213, 187},{ 451,1450},{1368, 885}, | |
340 { 592, 578},{1131, 531},{ 861,1855},{1764,1500}, | |
341 { 444, 970},{ 935, 903},{ 424,1687},{1633,1102}, | |
342 { 793, 897},{1060, 897},{ 185,2011},{1205,1855}}; | |
343 | |
344 static const qcelp_vector qcelp_lspvq3[128]= { | |
345 { 225, 283},{1296, 355},{ 543, 343},{2073, 274}, | |
346 { 204,1099},{1562, 523},{1388, 161},{2784, 274}, | |
347 { 112, 849},{1870, 175},{1189, 160},{1490,1088}, | |
348 { 969,1115},{ 659,3322},{1158,1073},{3183,1363}, | |
349 { 517, 223},{1740, 223},{ 704, 387},{2637, 234}, | |
350 { 692,1005},{1287,1610},{ 952, 532},{2393, 646}, | |
351 { 490, 552},{1619, 657},{ 845, 670},{1784,2280}, | |
352 { 191,1775},{ 272,2868},{ 942, 952},{2628,1479}, | |
353 { 278, 579},{1565, 218},{ 814, 180},{2379, 187}, | |
354 { 276,1444},{1199,1223},{1200, 349},{3009, 307}, | |
355 { 312, 844},{1898, 306},{ 863, 470},{1685,1241}, | |
356 { 513,1727},{ 711,2233},{1085, 864},{3398, 527}, | |
357 { 414, 440},{1356, 612},{ 964, 147},{2173, 738}, | |
358 { 465,1292},{ 877,1749},{1104, 689},{2105,1311}, | |
359 { 580, 864},{1895, 752},{ 652, 609},{1485,1699}, | |
360 { 514,1400},{ 386,2131},{ 933, 798},{2473, 986}, | |
361 { 334, 360},{1375, 398},{ 621, 276},{2183, 280}, | |
362 { 311,1114},{1382, 807},{1284, 175},{2605, 636}, | |
363 { 230, 816},{1739, 408},{1074, 176},{1619,1120}, | |
364 { 784,1371},{ 448,3050},{1189, 880},{3039,1165}, | |
365 { 424, 241},{1672, 186},{ 815, 333},{2432, 324}, | |
366 { 584,1029},{1137,1546},{1015, 585},{2198, 995}, | |
367 { 574, 581},{1746, 647},{ 733, 740},{1938,1737}, | |
368 { 347,1710},{ 373,2429},{ 787,1061},{2439,1438}, | |
369 { 185, 536},{1489, 178},{ 703, 216},{2178, 487}, | |
370 { 154,1421},{1414, 994},{1103, 352},{3072, 473}, | |
371 { 408, 819},{2055, 168},{ 998, 354},{1917,1140}, | |
372 { 665,1799},{ 993,2213},{1234, 631},{3003, 762}, | |
373 { 373, 620},{1518, 425},{ 913, 300},{1966, 836}, | |
374 { 402,1185},{ 948,1385},{1121, 555},{1802,1509}, | |
375 { 474, 886},{1888, 610},{ 739, 585},{1231,2379}, | |
376 { 661,1335},{ 205,2211},{ 823, 822},{2480,1179}}; | |
377 | |
378 static const qcelp_vector qcelp_lspvq4[64]= { | |
379 { 348, 311},{ 812,1145},{ 552, 461},{1826, 263}, | |
380 { 601, 675},{1730, 172},{1523, 193},{2449, 277}, | |
381 { 334, 668},{ 805,1441},{1319, 207},{1684, 910}, | |
382 { 582,1318},{1403,1098},{ 979, 832},{2700,1359}, | |
383 { 624, 228},{1292, 979},{ 800, 195},{2226, 285}, | |
384 { 730, 862},{1537, 601},{1115, 509},{2720, 354}, | |
385 { 218,1167},{1212,1538},{1074, 247},{1674,1710}, | |
386 { 322,2142},{1263, 777},{ 981, 556},{2119,1710}, | |
387 { 193, 596},{1035, 957},{ 694, 397},{1997, 253}, | |
388 { 743, 603},{1584, 321},{1346, 346},{2221, 708}, | |
389 { 451, 732},{1040,1415},{1184, 230},{1853, 919}, | |
390 { 310,1661},{1625, 706},{ 856, 843},{2902, 702}, | |
391 { 467, 348},{1108,1048},{ 859, 306},{1964, 463}, | |
392 { 560,1013},{1425, 533},{1142, 634},{2391, 879}, | |
393 { 397,1084},{1345,1700},{ 976, 248},{1887,1189}, | |
394 { 644,2087},{1262, 603},{ 877, 550},{2203,1307}}; | |
395 | |
396 static const qcelp_vector qcelp_lspvq5[64]= { | |
397 { 360, 222},{ 820,1097},{ 601, 319},{1656, 198}, | |
398 { 604, 513},{1552, 141},{1391, 155},{2474, 261}, | |
399 { 269, 785},{1463, 646},{1123, 191},{2015, 223}, | |
400 { 785, 844},{1202,1011},{ 980, 807},{3014, 793}, | |
401 { 570, 180},{1135,1382},{ 778, 256},{1901, 179}, | |
402 { 807, 622},{1461, 458},{1231, 178},{2028, 821}, | |
403 { 387, 927},{1496,1004},{ 888, 392},{2246, 341}, | |
404 { 295,1462},{1156, 694},{1022, 473},{2226,1364}, | |
405 { 210, 478},{1029,1020},{ 722, 181},{1730, 251}, | |
406 { 730, 488},{1465, 293},{1303, 326},{2595, 387}, | |
407 { 458, 584},{1569, 742},{1029, 173},{1910, 495}, | |
408 { 605,1159},{1268, 719},{ 973, 646},{2872, 428}, | |
409 { 443, 334},{ 835,1465},{ 912, 138},{1716, 442}, | |
410 { 620, 778},{1316, 450},{1186, 335},{1446,1665}, | |
411 { 486,1050},{1675,1019},{ 880, 278},{2214, 202}, | |
412 { 539,1564},{1142, 533},{ 984, 391},{2130,1089}}; | |
413 | |
414 static const qcelp_vector * const qcelp_lspvq[5] = { | |
415 qcelp_lspvq1, | |
416 qcelp_lspvq2, | |
417 qcelp_lspvq3, | |
418 qcelp_lspvq4, | |
419 qcelp_lspvq5 | |
420 }; | |
421 | |
8123 | 422 /** |
8127 | 423 * the final gain scalefactor before clipping into a usable output float |
424 */ | |
425 #define QCELP_SCALE 8192. | |
426 | |
427 /** | |
8191 | 428 * the upper boundary of the clipping, depends on QCELP_SCALE |
429 */ | |
430 #define QCELP_CLIP_UPPER_BOUND (8191.75/8192.) | |
431 | |
432 /** | |
433 * the lower boundary of the clipping, depends on QCELP_SCALE | |
434 */ | |
435 #define QCELP_CLIP_LOWER_BOUND -1. | |
436 | |
437 /** | |
8127 | 438 * table for computing Ga (decoded linear codebook gain magnitude) |
439 * | |
440 * @note The table could fit in int16_t in x*8 form, but it seems | |
441 * to be slower on x86 | |
442 * | |
443 * TIA/EIA/IS-733 2.4.6.2.1-3 | |
444 */ | |
445 | |
446 static const float qcelp_g12ga[61] = { | |
447 1.000/QCELP_SCALE, 1.125/QCELP_SCALE, 1.250/QCELP_SCALE, 1.375/QCELP_SCALE, | |
448 1.625/QCELP_SCALE, 1.750/QCELP_SCALE, 2.000/QCELP_SCALE, 2.250/QCELP_SCALE, | |
449 2.500/QCELP_SCALE, 2.875/QCELP_SCALE, 3.125/QCELP_SCALE, 3.500/QCELP_SCALE, | |
450 4.000/QCELP_SCALE, 4.500/QCELP_SCALE, 5.000/QCELP_SCALE, 5.625/QCELP_SCALE, | |
451 6.250/QCELP_SCALE, 7.125/QCELP_SCALE, 8.000/QCELP_SCALE, 8.875/QCELP_SCALE, | |
452 10.000/QCELP_SCALE, 11.250/QCELP_SCALE, 12.625/QCELP_SCALE, 14.125/QCELP_SCALE, | |
453 15.875/QCELP_SCALE, 17.750/QCELP_SCALE, 20.000/QCELP_SCALE, 22.375/QCELP_SCALE, | |
454 25.125/QCELP_SCALE, 28.125/QCELP_SCALE, 31.625/QCELP_SCALE, 35.500/QCELP_SCALE, | |
455 39.750/QCELP_SCALE, 44.625/QCELP_SCALE, 50.125/QCELP_SCALE, 56.250/QCELP_SCALE, | |
456 63.125/QCELP_SCALE, 70.750/QCELP_SCALE, 79.375/QCELP_SCALE, 89.125/QCELP_SCALE, | |
457 100.000/QCELP_SCALE, 112.250/QCELP_SCALE, 125.875/QCELP_SCALE, 141.250/QCELP_SCALE, | |
458 158.500/QCELP_SCALE, 177.875/QCELP_SCALE, 199.500/QCELP_SCALE, 223.875/QCELP_SCALE, | |
459 251.250/QCELP_SCALE, 281.875/QCELP_SCALE, 316.250/QCELP_SCALE, 354.875/QCELP_SCALE, | |
8182 | 460 398.125/QCELP_SCALE, 446.625/QCELP_SCALE, 501.125/QCELP_SCALE, 562.375/QCELP_SCALE, |
8127 | 461 631.000/QCELP_SCALE, 708.000/QCELP_SCALE, 794.375/QCELP_SCALE, 891.250/QCELP_SCALE, |
462 1000.000/QCELP_SCALE}; | |
463 | |
464 /** | |
8123 | 465 * circular codebook for rate 1 frames in x*100 form |
466 * | |
467 * TIA/EIA/IS-733 2.4.6.1-2 | |
468 */ | |
469 static const int16_t qcelp_rate_full_codebook[128] = { | |
470 10, -65, -59, 12, 110, 34, -134, 157, | |
471 104, -84, -34, -115, 23, -101, 3, 45, | |
472 -101, -16, -59, 28, -45, 134, -67, 22, | |
473 61, -29, 226, -26, -55, -179, 157, -51, | |
474 -220, -93, -37, 60, 118, 74, -48, -95, | |
475 -181, 111, 36, -52, -215, 78, -112, 39, | |
476 -17, -47, -223, 19, 12, -98, -142, 130, | |
477 54, -127, 21, -12, 39, -48, 12, 128, | |
478 6, -167, 82, -102, -79, 55, -44, 48, | |
479 -20, -53, 8, -61, 11, -70, -157, -168, | |
480 20, -56, -74, 78, 33, -63, -173, -2, | |
481 -75, -53, -146, 77, 66, -29, 9, -75, | |
482 65, 119, -43, 76, 233, 98, 125, -156, | |
483 -27, 78, -9, 170, 176, 143, -148, -7, | |
484 27, -136, 5, 27, 18, 139, 204, 7, | |
485 -184, -197, 52, -3, 78, -189, 8, -65 | |
486 }; | |
487 #define QCELP_RATE_FULL_CODEBOOK_RATIO .01 | |
488 | |
489 /** | |
490 * circular codebook for rate 1/2 frames in x*2 form | |
491 * | |
492 * TIA/EIA/IS-733 2.4.6.1-1 | |
493 */ | |
494 static const int8_t qcelp_rate_half_codebook[128] = { | |
495 0, -4, 0, -3, 0, 0, 0, 0, | |
496 0, 0, 0, 0, 0, 0, 0, 0, | |
497 0, -3, -2, 0, 0, 0, 0, 0, | |
498 0, 0, 0, 0, 0, 0, 0, 5, | |
499 0, 0, 0, 0, 0, 0, 4, 0, | |
500 0, 3, 2, 0, 3, 4, 0, 0, | |
501 0, 0, 0, 0, 0, 0, 0, 0, | |
502 0, 0, 0, 0, 0, 3, 0, 0, | |
503 -3, 3, 0, 0, -2, 0, 3, 0, | |
504 0, 0, 0, 0, 0, 0, -5, 0, | |
505 0, 0, 0, 3, 0, 0, 0, 3, | |
506 0, 0, 0, 0, 0, 0, 0, 4, | |
507 0, 0, 0, 0, 0, 0, 0, 0, | |
508 0, 3, 6, -3, -4, 0, -3, -3, | |
509 3, -3, 0, 0, 0, 0, 0, 0, | |
510 0, 0, 0, 0, 0, 0, 0, 0 | |
511 }; | |
512 #define QCELP_RATE_HALF_CODEBOOK_RATIO 0.5 | |
513 | |
8127 | 514 /** |
8145 | 515 * sqrt(1.887) is the maximum of the pseudorandom |
516 * white sequence used to generate the scaled codebook | |
8191 | 517 * vector for bitrate 1/4. |
8145 | 518 * |
519 * TIA/EIA/IS-733 2.4.8.1.2 | |
520 */ | |
521 #define QCELP_SQRT1887 1.373681186 | |
522 | |
523 /** | |
8127 | 524 * table for impulse response of BPF used to filter |
8191 | 525 * the white excitation for bitrate 1/4 synthesis |
8127 | 526 * |
8190 | 527 * Only half the tables are needed because of symmetry. |
8127 | 528 * |
529 * TIA/EIA/IS-733 2.4.8.1.2-1.1 | |
530 */ | |
531 static const double qcelp_rnd_fir_coefs[11] = { | |
532 -1.344519e-1, 1.735384e-2, -6.905826e-2, 2.434368e-2, | |
533 -8.210701e-2, 3.041388e-2, -9.251384e-2, 3.501983e-2, | |
534 -9.918777e-2, 3.749518e-2, 8.985137e-1 | |
535 }; | |
536 | |
8191 | 537 /** |
538 * This spread factor is used, for bitrate 1/8 and I_F_Q, | |
539 * to force the LSP frequencies to be at least 80 Hz apart. | |
540 * | |
541 * TIA/EIA/IS-733 2.4.3.3.2 | |
542 */ | |
543 #define QCELP_LSP_SPREAD_FACTOR 0.02 | |
544 | |
545 /** | |
546 * predictor coefficient for the conversion of LSP codes | |
547 * to LSP frequencies for 1/8 and I_F_Q | |
548 * | |
549 * TIA/EIA/IS-733 2.4.3.2.7-2 | |
550 */ | |
551 #define QCELP_LSP_OCTAVE_PREDICTOR 29.0/32 | |
552 | |
9123
36a5caff8540
Part 2 of 2 of Kenan Gillet's 'make ff_qcelp_lspf2lpc
reynaldo
parents:
8718
diff
changeset
|
553 /** |
36a5caff8540
Part 2 of 2 of Kenan Gillet's 'make ff_qcelp_lspf2lpc
reynaldo
parents:
8718
diff
changeset
|
554 * initial coefficient to perform bandwidth expansion on LPC |
36a5caff8540
Part 2 of 2 of Kenan Gillet's 'make ff_qcelp_lspf2lpc
reynaldo
parents:
8718
diff
changeset
|
555 * |
36a5caff8540
Part 2 of 2 of Kenan Gillet's 'make ff_qcelp_lspf2lpc
reynaldo
parents:
8718
diff
changeset
|
556 * @note: 0.9883 looks like an approximation of 253/256. |
36a5caff8540
Part 2 of 2 of Kenan Gillet's 'make ff_qcelp_lspf2lpc
reynaldo
parents:
8718
diff
changeset
|
557 * |
36a5caff8540
Part 2 of 2 of Kenan Gillet's 'make ff_qcelp_lspf2lpc
reynaldo
parents:
8718
diff
changeset
|
558 * TIA/EIA/IS-733 2.4.3.3.6 6 |
36a5caff8540
Part 2 of 2 of Kenan Gillet's 'make ff_qcelp_lspf2lpc
reynaldo
parents:
8718
diff
changeset
|
559 */ |
9629 | 560 #define QCELP_BANDWIDTH_EXPANSION_COEFF 0.9883 |
9123
36a5caff8540
Part 2 of 2 of Kenan Gillet's 'make ff_qcelp_lspf2lpc
reynaldo
parents:
8718
diff
changeset
|
561 |
8081 | 562 #endif /* AVCODEC_QCELPDATA_H */ |