You can not select more than 25 topics Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.

382 lines
14KB

  1. /*
  2. * part of 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. #ifndef AVCODEC_QCELPDATA_H
  22. #define AVCODEC_QCELPDATA_H
  23. #include <stdint.h>
  24. /**
  25. * pre-calculated table for hammsinc function
  26. * Only half of the table is needed because of symmetry.
  27. *
  28. * TIA/EIA/IS-733 2.4.5.2-2/3
  29. */
  30. static const float qcelp_hammsinc_table[4] = { -0.006822, 0.041249, -0.143459, 0.588863};
  31. typedef struct {
  32. uint8_t index; /*!< index into the QCELPContext structure */
  33. uint8_t bitpos; /*!< position of the lowest bit in the value's byte */
  34. uint8_t bitlen; /*!< number of bits to read */
  35. } QCELPBitmap;
  36. #define QCELP_OF(variable, bit, len) {offsetof(QCELPContext, variable), bit, len}
  37. /**
  38. * bitmap unpacking tables for RATE_FULL
  39. *
  40. * TIA/EIA/IS-733 Table 2.4.7.1-1
  41. */
  42. static const QCELPBitmap qcelp_rate_full_bitmap[] = {
  43. // start on bit
  44. QCELP_OF(lspv [ 2], 0, 3), // 265
  45. QCELP_OF(lspv [ 1], 0, 7), // 262
  46. QCELP_OF(lspv [ 0], 0, 6), // 255
  47. QCELP_OF(lspv [ 4], 0, 6), // 249
  48. QCELP_OF(lspv [ 3], 0, 6), // 243
  49. QCELP_OF(lspv [ 2], 3, 4), // 237
  50. QCELP_OF(cbsign[ 0], 0, 1), // 233
  51. QCELP_OF(cbgain[ 0], 0, 4), // 232
  52. QCELP_OF(pfrac [ 0], 0, 1), // 228
  53. QCELP_OF(plag [ 0], 0, 7), // 227
  54. QCELP_OF(pgain [ 0], 0, 3), // 220
  55. QCELP_OF(cindex[ 1], 0, 4), // 217
  56. QCELP_OF(cbsign[ 1], 0, 1), // 213
  57. QCELP_OF(cbgain[ 1], 0, 4), // 212
  58. QCELP_OF(cindex[ 0], 0, 7), // 208
  59. QCELP_OF(cbgain[ 3], 0, 1), // 201
  60. QCELP_OF(cindex[ 2], 0, 7), // 200
  61. QCELP_OF(cbsign[ 2], 0, 1), // 193
  62. QCELP_OF(cbgain[ 2], 0, 4), // 192
  63. QCELP_OF(cindex[ 1], 4, 3), // 188
  64. QCELP_OF(plag [ 1], 0, 3), // 185
  65. QCELP_OF(pgain [ 1], 0, 3), // 182
  66. QCELP_OF(cindex[ 3], 0, 7), // 179
  67. QCELP_OF(cbsign[ 3], 0, 1), // 172
  68. QCELP_OF(cbgain[ 3], 1, 2), // 171
  69. QCELP_OF(cindex[ 4], 0, 6), // 169
  70. QCELP_OF(cbsign[ 4], 0, 1), // 163
  71. QCELP_OF(cbgain[ 4], 0, 4), // 162
  72. QCELP_OF(pfrac [ 1], 0, 1), // 158
  73. QCELP_OF(plag [ 1], 3, 4), // 157
  74. QCELP_OF(cbgain[ 6], 0, 3), // 153
  75. QCELP_OF(cindex[ 5], 0, 7), // 150
  76. QCELP_OF(cbsign[ 5], 0, 1), // 143
  77. QCELP_OF(cbgain[ 5], 0, 4), // 142
  78. QCELP_OF(cindex[ 4], 6, 1), // 138
  79. QCELP_OF(cindex[ 7], 0, 3), // 137
  80. QCELP_OF(cbsign[ 7], 0, 1), // 134
  81. QCELP_OF(cbgain[ 7], 0, 3), // 133
  82. QCELP_OF(cindex[ 6], 0, 7), // 130
  83. QCELP_OF(cbsign[ 6], 0, 1), // 123
  84. QCELP_OF(cbgain[ 6], 3, 1), // 122
  85. QCELP_OF(cbgain[ 8], 0, 1), // 121
  86. QCELP_OF(pfrac [ 2], 0, 1), // 120
  87. QCELP_OF(plag [ 2], 0, 7), // 119
  88. QCELP_OF(pgain [ 2], 0, 3), // 112
  89. QCELP_OF(cindex[ 7], 3, 4), // 109
  90. QCELP_OF(cbsign[ 9], 0, 1), // 105
  91. QCELP_OF(cbgain[ 9], 0, 4), // 104
  92. QCELP_OF(cindex[ 8], 0, 7), // 100
  93. QCELP_OF(cbsign[ 8], 0, 1), // 93
  94. QCELP_OF(cbgain[ 8], 1, 3), // 92
  95. QCELP_OF(cindex[10], 0, 4), // 89
  96. QCELP_OF(cbsign[10], 0, 1), // 85
  97. QCELP_OF(cbgain[10], 0, 4), // 84
  98. QCELP_OF(cindex[ 9], 0, 7), // 80
  99. QCELP_OF(pgain [ 3], 0, 2), // 73
  100. QCELP_OF(cindex[11], 0, 7), // 71
  101. QCELP_OF(cbsign[11], 0, 1), // 64
  102. QCELP_OF(cbgain[11], 0, 3), // 63
  103. QCELP_OF(cindex[10], 4, 3), // 60
  104. QCELP_OF(cindex[12], 0, 2), // 57
  105. QCELP_OF(cbsign[12], 0, 1), // 55
  106. QCELP_OF(cbgain[12], 0, 4), // 54
  107. QCELP_OF(pfrac [ 3], 0, 1), // 50
  108. QCELP_OF(plag [ 3], 0, 7), // 49
  109. QCELP_OF(pgain [ 3], 2, 1), // 42
  110. QCELP_OF(cindex[13], 0, 6), // 41
  111. QCELP_OF(cbsign[13], 0, 1), // 35
  112. QCELP_OF(cbgain[13], 0, 4), // 34
  113. QCELP_OF(cindex[12], 2, 5), // 30
  114. QCELP_OF(cbgain[15], 0, 3), // 25
  115. QCELP_OF(cindex[14], 0, 7), // 22
  116. QCELP_OF(cbsign[14], 0, 1), // 15
  117. QCELP_OF(cbgain[14], 0, 4), // 14
  118. QCELP_OF(cindex[13], 6, 1), // 10
  119. QCELP_OF(reserved, 0, 2), // 9
  120. QCELP_OF(cindex[15], 0, 7), // 7
  121. QCELP_OF(cbsign[15], 0, 1) // 0
  122. };
  123. /**
  124. * bitmap unpacking tables for RATE_HALF
  125. *
  126. * TIA/EIA/IS-733 Table 2.4.7.2-1
  127. */
  128. static const QCELPBitmap qcelp_rate_half_bitmap[] = {
  129. // start on bit
  130. QCELP_OF(lspv [2], 0, 3), // 123
  131. QCELP_OF(lspv [1], 0, 7), // 120
  132. QCELP_OF(lspv [0], 0, 6), // 113
  133. QCELP_OF(lspv [4], 0, 6), // 107
  134. QCELP_OF(lspv [3], 0, 6), // 101
  135. QCELP_OF(lspv [2], 3, 4), // 95
  136. QCELP_OF(cbsign[0], 0, 1), // 91
  137. QCELP_OF(cbgain[0], 0, 4), // 90
  138. QCELP_OF(pfrac [0], 0, 1), // 86
  139. QCELP_OF(plag [0], 0, 7), // 85
  140. QCELP_OF(pgain [0], 0, 3), // 78
  141. QCELP_OF(plag [1], 0, 6), // 75
  142. QCELP_OF(pgain [1], 0, 3), // 69
  143. QCELP_OF(cindex[0], 0, 7), // 66
  144. QCELP_OF(pgain [2], 0, 2), // 59
  145. QCELP_OF(cindex[1], 0, 7), // 57
  146. QCELP_OF(cbsign[1], 0, 1), // 50
  147. QCELP_OF(cbgain[1], 0, 4), // 49
  148. QCELP_OF(pfrac [1], 0, 1), // 45
  149. QCELP_OF(plag [1], 6, 1), // 44
  150. QCELP_OF(cindex[2], 0, 2), // 43
  151. QCELP_OF(cbsign[2], 0, 1), // 41
  152. QCELP_OF(cbgain[2], 0, 4), // 40
  153. QCELP_OF(pfrac [2], 0, 1), // 36
  154. QCELP_OF(plag [2], 0, 7), // 35
  155. QCELP_OF(pgain [2], 2, 1), // 28
  156. QCELP_OF(pfrac [3], 0, 1), // 27
  157. QCELP_OF(plag [3], 0, 7), // 26
  158. QCELP_OF(pgain [3], 0, 3), // 19
  159. QCELP_OF(cindex[2], 2, 5), // 16
  160. QCELP_OF(cindex[3], 0, 7), // 11
  161. QCELP_OF(cbsign[3], 0, 1), // 4
  162. QCELP_OF(cbgain[3], 0, 4) // 3
  163. };
  164. /**
  165. * bitmap unpacking tables for RATE_QUARTER
  166. *
  167. * TIA/EIA/IS-733 Table 2.4.7.3-1
  168. */
  169. static const QCELPBitmap qcelp_rate_quarter_bitmap[] = {
  170. // start on bit
  171. QCELP_OF(lspv [2], 0, 3), // 53
  172. QCELP_OF(lspv [1], 0, 7), // 50
  173. QCELP_OF(lspv [0], 0, 6), // 43
  174. QCELP_OF(lspv [4], 0, 6), // 37
  175. QCELP_OF(lspv [3], 0, 6), // 31
  176. QCELP_OF(lspv [2], 3, 4), // 25
  177. QCELP_OF(cbgain[3], 0, 4), // 21
  178. QCELP_OF(cbgain[2], 0, 4), // 17
  179. QCELP_OF(cbgain[1], 0, 4), // 13
  180. QCELP_OF(cbgain[0], 0, 4), // 9
  181. QCELP_OF(reserved, 0, 2), // 5
  182. QCELP_OF(cbgain[4], 0, 4) // 3
  183. };
  184. /**
  185. * bitmap unpacking tables for RATE_OCTAVE
  186. *
  187. * trick: CBSEED is written into QCELPContext.cbsign[15],
  188. * which is not used for RATE_OCTAVE.
  189. * CBSEED is only used to ensure the occurrence of random bit
  190. * patterns in the 16 first bits that are used as the seed.
  191. *
  192. * TIA/EIA/IS-733 Table 2.4.7.4-1
  193. */
  194. static const QCELPBitmap qcelp_rate_octave_bitmap[] = {
  195. // start on bit
  196. QCELP_OF(cbsign[15], 3, 1), // 19
  197. QCELP_OF(lspv [0], 0, 1), // 18
  198. QCELP_OF(lspv [1], 0, 1), // 17
  199. QCELP_OF(lspv [2], 0, 1), // 16
  200. QCELP_OF(cbsign[15], 2, 1), // 15
  201. QCELP_OF(lspv [3], 0, 1), // 14
  202. QCELP_OF(lspv [4], 0, 1), // 13
  203. QCELP_OF(lspv [5], 0, 1), // 12
  204. QCELP_OF(cbsign[15], 1, 1), // 11
  205. QCELP_OF(lspv [6], 0, 1), // 10
  206. QCELP_OF(lspv [7], 0, 1), // 9
  207. QCELP_OF(lspv [8], 0, 1), // 8
  208. QCELP_OF(cbsign[15], 0, 1), // 7
  209. QCELP_OF(lspv [9], 0, 1), // 6
  210. QCELP_OF(cbgain [0], 0, 2), // 7
  211. QCELP_OF(reserved, 0, 4) // 3
  212. };
  213. /**
  214. * position of the bitmapping data for each packet type in
  215. * the QCELPContext
  216. */
  217. static const QCELPBitmap * const qcelp_unpacking_bitmaps_per_rate[5] = {
  218. NULL, ///!< for SILENCE rate
  219. qcelp_rate_octave_bitmap,
  220. qcelp_rate_quarter_bitmap,
  221. qcelp_rate_half_bitmap,
  222. qcelp_rate_full_bitmap,
  223. };
  224. static const uint16_t qcelp_bits_per_rate[5] = {
  225. 0, ///!< for SILENCE rate
  226. FF_ARRAY_ELEMS(qcelp_rate_octave_bitmap),
  227. FF_ARRAY_ELEMS(qcelp_rate_quarter_bitmap),
  228. FF_ARRAY_ELEMS(qcelp_rate_half_bitmap),
  229. FF_ARRAY_ELEMS(qcelp_rate_full_bitmap),
  230. };
  231. typedef uint16_t qcelp_vector[2];
  232. /**
  233. * LSP vector quantization tables in x*10000 form
  234. *
  235. * TIA/EIA/IS-733 tables 2.4.3.2.6.3-1 through 2.4.3.2.6.3-5
  236. */
  237. static const qcelp_vector qcelp_lspvq1[64]= {
  238. { 327, 118},{ 919, 111},{ 427, 440},{1327, 185},
  239. { 469, 50},{1272, 91},{ 892, 59},{1771, 193},
  240. { 222, 158},{1100, 127},{ 827, 55},{ 978, 791},
  241. { 665, 47},{ 700,1401},{ 670, 859},{1913,1048},
  242. { 471, 215},{1046, 125},{ 645, 298},{1599, 160},
  243. { 593, 39},{1187, 462},{ 749, 341},{1520, 511},
  244. { 290, 792},{ 909, 362},{ 753, 81},{1111,1058},
  245. { 519, 253},{ 828, 839},{ 685, 541},{1421,1258},
  246. { 386, 130},{ 962, 119},{ 542, 387},{1431, 185},
  247. { 526, 51},{1175, 260},{ 831, 167},{1728, 510},
  248. { 273, 437},{1172, 113},{ 771, 144},{1122, 751},
  249. { 619, 119},{ 492,1276},{ 658, 695},{1882, 615},
  250. { 415, 200},{1018, 88},{ 681, 339},{1436, 325},
  251. { 555, 122},{1042, 485},{ 826, 345},{1374, 743},
  252. { 383,1018},{1005, 358},{ 704, 86},{1301, 586},
  253. { 597, 241},{ 832, 621},{ 555, 573},{1504, 839}};
  254. static const qcelp_vector qcelp_lspvq2[128]= {
  255. { 255, 293},{ 904, 219},{ 151,1211},{1447, 498},
  256. { 470, 253},{1559, 177},{1547, 994},{2394, 242},
  257. { 91, 813},{ 857, 590},{ 934,1326},{1889, 282},
  258. { 813, 472},{1057,1494},{ 450,3315},{2163,1895},
  259. { 538, 532},{1399, 218},{ 146,1552},{1755, 626},
  260. { 822, 202},{1299, 663},{ 706,1732},{2656, 401},
  261. { 418, 745},{ 762,1038},{ 583,1748},{1746,1285},
  262. { 527,1169},{1314, 830},{ 556,2116},{1073,2321},
  263. { 297, 570},{ 981, 403},{ 468,1103},{1740, 243},
  264. { 725, 179},{1255, 474},{1374,1362},{1922, 912},
  265. { 285, 947},{ 930, 700},{ 593,1372},{1909, 576},
  266. { 588, 916},{1110,1116},{ 224,2719},{1633,2220},
  267. { 402, 520},{1061, 448},{ 402,1352},{1499, 775},
  268. { 664, 589},{1081, 727},{ 801,2206},{2165,1157},
  269. { 566, 802},{ 911,1116},{ 306,1703},{1792, 836},
  270. { 655, 999},{1061,1038},{ 298,2089},{1110,1753},
  271. { 361, 311},{ 970, 239},{ 265,1231},{1495, 573},
  272. { 566, 262},{1569, 293},{1341,1144},{2271, 544},
  273. { 214, 877},{ 847, 719},{ 794,1384},{2067, 274},
  274. { 703, 688},{1099,1306},{ 391,2947},{2024,1670},
  275. { 471, 525},{1245, 290},{ 264,1557},{1568, 807},
  276. { 718, 399},{1193, 685},{ 883,1594},{2729, 764},
  277. { 500, 754},{ 809,1108},{ 541,1648},{1523,1385},
  278. { 614,1196},{1209, 847},{ 345,2242},{1442,1747},
  279. { 199, 560},{1092, 194},{ 349,1253},{1653, 507},
  280. { 625, 354},{1376, 431},{1187,1465},{2164, 872},
  281. { 360, 974},{1008, 698},{ 704,1346},{2114, 452},
  282. { 720, 816},{1240,1089},{ 439,2475},{1498,2040},
  283. { 336, 718},{1213, 187},{ 451,1450},{1368, 885},
  284. { 592, 578},{1131, 531},{ 861,1855},{1764,1500},
  285. { 444, 970},{ 935, 903},{ 424,1687},{1633,1102},
  286. { 793, 897},{1060, 897},{ 185,2011},{1205,1855}};
  287. static const qcelp_vector qcelp_lspvq3[128]= {
  288. { 225, 283},{1296, 355},{ 543, 343},{2073, 274},
  289. { 204,1099},{1562, 523},{1388, 161},{2784, 274},
  290. { 112, 849},{1870, 175},{1189, 160},{1490,1088},
  291. { 969,1115},{ 659,3322},{1158,1073},{3183,1363},
  292. { 517, 223},{1740, 223},{ 704, 387},{2637, 234},
  293. { 692,1005},{1287,1610},{ 952, 532},{2393, 646},
  294. { 490, 552},{1619, 657},{ 845, 670},{1784,2280},
  295. { 191,1775},{ 272,2868},{ 942, 952},{2628,1479},
  296. { 278, 579},{1565, 218},{ 814, 180},{2379, 187},
  297. { 276,1444},{1199,1223},{1200, 349},{3009, 307},
  298. { 312, 844},{1898, 306},{ 863, 470},{1685,1241},
  299. { 513,1727},{ 711,2233},{1085, 864},{3398, 527},
  300. { 414, 440},{1356, 612},{ 964, 147},{2173, 738},
  301. { 465,1292},{ 877,1749},{1104, 689},{2105,1311},
  302. { 580, 864},{1895, 752},{ 652, 609},{1485,1699},
  303. { 514,1400},{ 386,2131},{ 933, 798},{2473, 986},
  304. { 334, 360},{1375, 398},{ 621, 276},{2183, 280},
  305. { 311,1114},{1382, 807},{1284, 175},{2605, 636},
  306. { 230, 816},{1739, 408},{1074, 176},{1619,1120},
  307. { 784,1371},{ 448,3050},{1189, 880},{3039,1165},
  308. { 424, 241},{1672, 186},{ 815, 333},{2432, 324},
  309. { 584,1029},{1137,1546},{1015, 585},{2198, 995},
  310. { 574, 581},{1746, 647},{ 733, 740},{1938,1737},
  311. { 347,1710},{ 373,2429},{ 787,1061},{2439,1438},
  312. { 185, 536},{1489, 178},{ 703, 216},{2178, 487},
  313. { 154,1421},{1414, 994},{1103, 352},{3072, 473},
  314. { 408, 819},{2055, 168},{ 998, 354},{1917,1140},
  315. { 665,1799},{ 993,2213},{1234, 631},{3003, 762},
  316. { 373, 620},{1518, 425},{ 913, 300},{1966, 836},
  317. { 402,1185},{ 948,1385},{1121, 555},{1802,1509},
  318. { 474, 886},{1888, 610},{ 739, 585},{1231,2379},
  319. { 661,1335},{ 205,2211},{ 823, 822},{2480,1179}};
  320. static const qcelp_vector qcelp_lspvq4[64]= {
  321. { 348, 311},{ 812,1145},{ 552, 461},{1826, 263},
  322. { 601, 675},{1730, 172},{1523, 193},{2449, 277},
  323. { 334, 668},{ 805,1441},{1319, 207},{1684, 910},
  324. { 582,1318},{1403,1098},{ 979, 832},{2700,1359},
  325. { 624, 228},{1292, 979},{ 800, 195},{2226, 285},
  326. { 730, 862},{1537, 601},{1115, 509},{2720, 354},
  327. { 218,1167},{1212,1538},{1074, 247},{1674,1710},
  328. { 322,2142},{1263, 777},{ 981, 556},{2119,1710},
  329. { 193, 596},{1035, 957},{ 694, 397},{1997, 253},
  330. { 743, 603},{1584, 321},{1346, 346},{2221, 708},
  331. { 451, 732},{1040,1415},{1184, 230},{1853, 919},
  332. { 310,1661},{1625, 706},{ 856, 843},{2902, 702},
  333. { 467, 348},{1108,1048},{ 859, 306},{1964, 463},
  334. { 560,1013},{1425, 533},{1142, 634},{2391, 879},
  335. { 397,1084},{1345,1700},{ 976, 248},{1887,1189},
  336. { 644,2087},{1262, 603},{ 877, 550},{2203,1307}};
  337. static const qcelp_vector qcelp_lspvq5[64]= {
  338. { 360, 222},{ 820,1097},{ 601, 319},{1656, 198},
  339. { 604, 513},{1552, 141},{1391, 155},{2474, 261},
  340. { 269, 785},{1463, 646},{1123, 191},{2015, 223},
  341. { 785, 844},{1202,1011},{ 980, 807},{3014, 793},
  342. { 570, 180},{1135,1382},{ 778, 256},{1901, 179},
  343. { 807, 622},{1461, 458},{1231, 178},{2028, 821},
  344. { 387, 927},{1496,1004},{ 888, 392},{2246, 341},
  345. { 295,1462},{1156, 694},{1022, 473},{2226,1364},
  346. { 210, 478},{1029,1020},{ 722, 181},{1730, 251},
  347. { 730, 488},{1465, 293},{1303, 326},{2595, 387},
  348. { 458, 584},{1569, 742},{1029, 173},{1910, 495},
  349. { 605,1159},{1268, 719},{ 973, 646},{2872, 428},
  350. { 443, 334},{ 835,1465},{ 912, 138},{1716, 442},
  351. { 620, 778},{1316, 450},{1186, 335},{1446,1665},
  352. { 486,1050},{1675,1019},{ 880, 278},{2214, 202},
  353. { 539,1564},{1142, 533},{ 984, 391},{2130,1089}};
  354. static const qcelp_vector * const qcelp_lspvq[5] = {
  355. qcelp_lspvq1,
  356. qcelp_lspvq2,
  357. qcelp_lspvq3,
  358. qcelp_lspvq4,
  359. qcelp_lspvq5
  360. };
  361. #endif /* AVCODEC_QCELPDATA_H */