FFmpeg: libavcodec/qcelpdata.h Source File

FFmpeg

[フレーム]

libavcodec

qcelpdata.h

Go to the documentation of this file.

1 /*

2 * QCELP decoder

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 */

22 #ifndef AVCODEC_QCELPDATA_H

23 #define AVCODEC_QCELPDATA_H

25 /**

26 * @file

27 * Data tables for the QCELP decoder

28 * @author Reynaldo H. Verdejo Pinochet

29 * @remark FFmpeg merging spearheaded by Kenan Gillet

30 * @remark Development mentored by Benjamin Larson

31 */

33 #include <stddef.h>

34 #include <stdint.h>

35 #include "libavutil/common.h"

37 /**

38 * QCELP unpacked data frame

39 */

40 typedef struct QCELPFrame {

41 /// @name 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 /// @}

48 /// @name 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 /// @}

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];

62 /**

63 * reserved bits only present in bitrate 1, 1/4 and 1/8 packets

64 */

65 uint8_t reserved;

66 } QCELPFrame;

68 /**

69 * Pre-calculated table for hammsinc function.

70 * Only half of the table is needed because of symmetry.

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};

76 typedef struct QCELPBitmap {

77 uint8_t index; /**< index into the QCELPContext structure */

78 uint8_t bitpos; /**< position of the lowest bit in the value's byte */

79 uint8_t bitlen; /**< number of bits to read */

80 } QCELPBitmap;

82 #define QCELP_OF(variable, bit, len) {offsetof(QCELPFrame, variable), bit, len}

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), // 5

261 QCELP_OF(reserved, 0, 4) // 3

262 };

263

264 /**

265 * Bitmapping data position 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

276 static const uint16_t qcelp_unpacking_bitmaps_lengths[5] = {

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

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

422 /**

423 * The final gain scalefactor before clipping into a usable output float

424 */

425 #define QCELP_SCALE 8192.

426

427 /**

428 * Table for computing Ga (decoded linear codebook gain magnitude)

429 *

430 * @note The table could fit in int16_t in x*8 form, but it seems

431 * to be slower on x86

432 *

433 * TIA/EIA/IS-733 2.4.6.2.1-3

434 */

435

436 static const float qcelp_g12ga[61] = {

437 1.000/QCELP_SCALE, 1.125/QCELP_SCALE, 1.250/QCELP_SCALE, 1.375/QCELP_SCALE,

438 1.625/QCELP_SCALE, 1.750/QCELP_SCALE, 2.000/QCELP_SCALE, 2.250/QCELP_SCALE,

439 2.500/QCELP_SCALE, 2.875/QCELP_SCALE, 3.125/QCELP_SCALE, 3.500/QCELP_SCALE,

440 4.000/QCELP_SCALE, 4.500/QCELP_SCALE, 5.000/QCELP_SCALE, 5.625/QCELP_SCALE,

441 6.250/QCELP_SCALE, 7.125/QCELP_SCALE, 8.000/QCELP_SCALE, 8.875/QCELP_SCALE,

442 10.000/QCELP_SCALE, 11.250/QCELP_SCALE, 12.625/QCELP_SCALE, 14.125/QCELP_SCALE,

443 15.875/QCELP_SCALE, 17.750/QCELP_SCALE, 20.000/QCELP_SCALE, 22.375/QCELP_SCALE,

444 25.125/QCELP_SCALE, 28.125/QCELP_SCALE, 31.625/QCELP_SCALE, 35.500/QCELP_SCALE,

445 39.750/QCELP_SCALE, 44.625/QCELP_SCALE, 50.125/QCELP_SCALE, 56.250/QCELP_SCALE,

446 63.125/QCELP_SCALE, 70.750/QCELP_SCALE, 79.375/QCELP_SCALE, 89.125/QCELP_SCALE,

447 100.000/QCELP_SCALE, 112.250/QCELP_SCALE, 125.875/QCELP_SCALE, 141.250/QCELP_SCALE,

448 158.500/QCELP_SCALE, 177.875/QCELP_SCALE, 199.500/QCELP_SCALE, 223.875/QCELP_SCALE,

449 251.250/QCELP_SCALE, 281.875/QCELP_SCALE, 316.250/QCELP_SCALE, 354.875/QCELP_SCALE,

450 398.125/QCELP_SCALE, 446.625/QCELP_SCALE, 501.125/QCELP_SCALE, 562.375/QCELP_SCALE,

451 631.000/QCELP_SCALE, 708.000/QCELP_SCALE, 794.375/QCELP_SCALE, 891.250/QCELP_SCALE,

452 1000.000/QCELP_SCALE};

453

454 /**

455 * Circular codebook for rate 1 frames in x*100 form

456 *

457 * TIA/EIA/IS-733 2.4.6.1-2

458 */

459 static const int16_t qcelp_rate_full_codebook[128] = {

460 10, -65, -59, 12, 110, 34, -134, 157,

461 104, -84, -34, -115, 23, -101, 3, 45,

462 -101, -16, -59, 28, -45, 134, -67, 22,

463 61, -29, 226, -26, -55, -179, 157, -51,

464 -220, -93, -37, 60, 118, 74, -48, -95,

465 -181, 111, 36, -52, -215, 78, -112, 39,

466 -17, -47, -223, 19, 12, -98, -142, 130,

467 54, -127, 21, -12, 39, -48, 12, 128,

468 6, -167, 82, -102, -79, 55, -44, 48,

469 -20, -53, 8, -61, 11, -70, -157, -168,

470 20, -56, -74, 78, 33, -63, -173, -2,

471 -75, -53, -146, 77, 66, -29, 9, -75,

472 65, 119, -43, 76, 233, 98, 125, -156,

473 -27, 78, -9, 170, 176, 143, -148, -7,

474 27, -136, 5, 27, 18, 139, 204, 7,

475 -184, -197, 52, -3, 78, -189, 8, -65

476 };

477 #define QCELP_RATE_FULL_CODEBOOK_RATIO .01

478

479 /**

480 * Circular codebook for rate 1/2 frames in x*2 form

481 *

482 * TIA/EIA/IS-733 2.4.6.1-1

483 */

484 static const int8_t qcelp_rate_half_codebook[128] = {

485 0, -4, 0, -3, 0, 0, 0, 0,

486 0, 0, 0, 0, 0, 0, 0, 0,

487 0, -3, -2, 0, 0, 0, 0, 0,

488 0, 0, 0, 0, 0, 0, 0, 5,

489 0, 0, 0, 0, 0, 0, 4, 0,

490 0, 3, 2, 0, 3, 4, 0, 0,

491 0, 0, 0, 0, 0, 0, 0, 0,

492 0, 0, 0, 0, 0, 3, 0, 0,

493 -3, 3, 0, 0, -2, 0, 3, 0,

494 0, 0, 0, 0, 0, 0, -5, 0,

495 0, 0, 0, 3, 0, 0, 0, 3,

496 0, 0, 0, 0, 0, 0, 0, 4,

497 0, 0, 0, 0, 0, 0, 0, 0,

498 0, 3, 6, -3, -4, 0, -3, -3,

499 3, -3, 0, 0, 0, 0, 0, 0,

500 0, 0, 0, 0, 0, 0, 0, 0

501 };

502 #define QCELP_RATE_HALF_CODEBOOK_RATIO 0.5

503

504 /**

505 * sqrt(1.887) is the maximum of the pseudorandom

506 * white sequence used to generate the scaled codebook

507 * vector for bitrate 1/4.

508 *

509 * TIA/EIA/IS-733 2.4.8.1.2

510 */

511 #define QCELP_SQRT1887 1.373681186

512

513 /**

514 * Table for impulse response of BPF used to filter

515 * the white excitation for bitrate 1/4 synthesis

516 *

517 * Only half the tables are needed because of symmetry.

518 *

519 * TIA/EIA/IS-733 2.4.8.1.2-1.1

520 */

521 static const double qcelp_rnd_fir_coefs[11] = {

522 -1.344519e-1, 1.735384e-2, -6.905826e-2, 2.434368e-2,

523 -8.210701e-2, 3.041388e-2, -9.251384e-2, 3.501983e-2,

524 -9.918777e-2, 3.749518e-2, 8.985137e-1

525 };

526

527 /**

528 * This spread factor is used, for bitrate 1/8 and I_F_Q,

529 * to force LSP frequencies to be at least 80 Hz apart.

530 *

531 * TIA/EIA/IS-733 2.4.3.3.2

532 */

533 #define QCELP_LSP_SPREAD_FACTOR 0.02

534

535 /**

536 * Predictor coefficient for the conversion of LSP codes

537 * to LSP frequencies for 1/8 and I_F_Q

538 *

539 * TIA/EIA/IS-733 2.4.3.2.7-2

540 */

541 #define QCELP_LSP_OCTAVE_PREDICTOR 29.0/32

542

543 /**

544 * Initial coefficient to perform bandwidth expansion on LPC

545 *

546 * @note: 0.9883 looks like an approximation of 253/256.

547 *

548 * TIA/EIA/IS-733 2.4.3.3.6 6

549 */

550 #define QCELP_BANDWIDTH_EXPANSION_COEFF 0.9883

551

552 #endif /* AVCODEC_QCELPDATA_H */

Generated on Sun Sep 14 2014 18:56:02 for FFmpeg by doxygen 1.8.2