Subversion Repositories Kolibri OS

Rev

Go to most recent revision | Blame | Last modification | View Log | RSS feed

  1. /*
  2.  * QCELP decoder
  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.  
  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.  */
  32.  
  33. #include <stddef.h>
  34. #include <stdint.h>
  35. #include "libavutil/common.h"
  36.  
  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. /// @}
  47.  
  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. /// @}
  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.  
  62.     /**
  63.      * reserved bits only present in bitrate 1, 1/4 and 1/8 packets
  64.      */
  65.     uint8_t reserved;
  66. } QCELPFrame;
  67.  
  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};
  75.  
  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;
  81.  
  82. #define QCELP_OF(variable, bit, len) {offsetof(QCELPFrame, variable), bit, len}
  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), //  5
  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.  
  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 the 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 */
  553.