FFmpeg
 All Data Structures Namespaces Files Functions Variables Typedefs Enumerations Enumerator Macros Groups Pages
bgmc.c
Go to the documentation of this file.
1 /*
2  * Block Gilbert-Moore decoder
3  * Copyright (c) 2010 Thilo Borgmann <thilo.borgmann _at_ googlemail.com>
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 /**
23  * @file
24  * Block Gilbert-Moore decoder as used by MPEG-4 ALS
25  * @author Thilo Borgmann <thilo.borgmann _at_ googlemail.com>
26  */
27 
28 #include "libavutil/attributes.h"
29 #include "bgmc.h"
30 
31 #define FREQ_BITS 14 // bits used by frequency counters
32 #define VALUE_BITS 18 // bits used to represent the values
33 #define TOP_VALUE ((1 << VALUE_BITS) - 1) // maximum value
34 #define FIRST_QTR (TOP_VALUE / 4 + 1) // first quarter of values maximum value
35 #define HALF (2 * FIRST_QTR) // first half of values maximum value
36 #define THIRD_QTR (3 * FIRST_QTR) // third quarter of values maximum value
37 
38 #define LUT_BITS (FREQ_BITS - 8) // number of bits used to index lookup tables
39 #define LUT_SIZE (1 << LUT_BITS) // size of the lookup tables
40 #define LUT_BUFF 4 // number of buffered lookup tables
41 
42 
43 /** Cumulative frequency tables for block Gilbert-Moore coding. */
44 static const uint16_t cf_tables_1[3][129] = {
45  {
46  16384, 16066, 15748, 15431, 15114, 14799, 14485, 14173, 13861, 13552,
47  13243, 12939, 12635, 12336, 12038, 11745, 11452, 11161, 10870, 10586,
48  10303, 10027, 9751, 9483, 9215, 8953, 8692, 8440, 8189, 7946,
49  7704, 7472, 7240, 7008, 6776, 6554, 6333, 6122, 5912, 5711,
50  5512, 5320, 5128, 4947, 4766, 4595, 4425, 4264, 4104, 3946,
51  3788, 3640, 3493, 3355, 3218, 3090, 2963, 2842, 2721, 2609,
52  2498, 2395, 2292, 2196, 2100, 2004, 1908, 1820, 1732, 1651,
53  1570, 1497, 1424, 1355, 1287, 1223, 1161, 1100, 1044, 988,
54  938, 888, 839, 790, 746, 702, 662, 623, 588, 553,
55  520, 488, 459, 431, 405, 380, 357, 334, 311, 288,
56  268, 248, 230, 213, 197, 182, 168, 154, 142, 130,
57  119, 108, 99, 90, 81, 72, 64, 56, 49, 42,
58  36, 30, 25, 20, 15, 11, 7, 3, 0
59  },
60  {
61  16384, 16080, 15776, 15473, 15170, 14868, 14567, 14268, 13970, 13674,
62  13378, 13086, 12794, 12505, 12218, 11936, 11654, 11373, 11092, 10818,
63  10544, 10276, 10008, 9749, 9490, 9236, 8982, 8737, 8492, 8256,
64  8020, 7792, 7564, 7336, 7108, 6888, 6669, 6459, 6249, 6050,
65  5852, 5660, 5468, 5286, 5104, 4931, 4760, 4598, 4436, 4275,
66  4115, 3965, 3816, 3674, 3534, 3403, 3272, 3147, 3023, 2907,
67  2792, 2684, 2577, 2476, 2375, 2274, 2173, 2079, 1986, 1897,
68  1810, 1724, 1645, 1567, 1493, 1419, 1351, 1284, 1222, 1161,
69  1105, 1050, 995, 941, 891, 842, 797, 753, 713, 673,
70  636, 599, 566, 533, 503, 473, 446, 419, 392, 365,
71  340, 316, 294, 272, 253, 234, 216, 199, 184, 169,
72  155, 142, 130, 118, 106, 95, 85, 75, 66, 57,
73  49, 41, 34, 27, 21, 15, 10, 5, 0
74  },
75  {
76  16384, 16092, 15801, 15510, 15219, 14930, 14641, 14355, 14069, 13785,
77  13501, 13219, 12938, 12661, 12384, 12112, 11841, 11571, 11301, 11037,
78  10773, 10514, 10256, 10005, 9754, 9508, 9263, 9025, 8787, 8557,
79  8327, 8103, 7879, 7655, 7431, 7215, 7000, 6792, 6585, 6387,
80  6190, 5998, 5807, 5625, 5445, 5272, 5100, 4937, 4774, 4613,
81  4452, 4301, 4150, 4007, 3865, 3731, 3597, 3469, 3341, 3218,
82  3099, 2981, 2869, 2758, 2652, 2546, 2440, 2334, 2234, 2134,
83  2041, 1949, 1864, 1779, 1699, 1620, 1547, 1474, 1407, 1340,
84  1278, 1217, 1157, 1097, 1043, 989, 940, 891, 846, 801,
85  759, 718, 680, 643, 609, 575, 543, 511, 479, 447,
86  418, 389, 363, 337, 314, 291, 270, 249, 230, 212,
87  195, 179, 164, 149, 135, 121, 108, 96, 85, 74,
88  64, 54, 45, 36, 28, 20, 13, 6, 0
89  }
90 };
91 
92 
93 static const uint16_t cf_tables_2[8][193] = {
94  {
95  16384, 16104, 15825, 15546, 15268, 14991, 14714, 14439, 14164, 13891,
96  13620, 13350, 13081, 12815, 12549, 12287, 12025, 11765, 11505, 11250,
97  10996, 10746, 10497, 10254, 10011, 9772, 9534, 9303, 9072, 8848,
98  8624, 8406, 8188, 7970, 7752, 7539, 7327, 7123, 6919, 6724,
99  6529, 6339, 6150, 5970, 5790, 5618, 5446, 5282, 5119, 4957,
100  4795, 4642, 4490, 4345, 4201, 4065, 3929, 3798, 3669, 3547,
101  3425, 3310, 3196, 3086, 2976, 2866, 2756, 2650, 2545, 2447,
102  2350, 2260, 2170, 2085, 2000, 1921, 1843, 1770, 1698, 1632,
103  1566, 1501, 1436, 1376, 1316, 1261, 1207, 1157, 1108, 1061,
104  1015, 973, 931, 893, 855, 819, 783, 747, 711, 677,
105  644, 614, 584, 557, 530, 505, 480, 458, 436, 416,
106  396, 378, 360, 343, 326, 310, 295, 281, 267, 255,
107  243, 232, 221, 211, 201, 192, 183, 174, 166, 158,
108  150, 142, 134, 126, 119, 112, 106, 100, 95, 90,
109  85, 80, 76, 72, 69, 66, 63, 60, 57, 54,
110  51, 48, 46, 44, 42, 40, 38, 36, 34, 33,
111  32, 31, 30, 29, 28, 27, 26, 25, 24, 23,
112  22, 21, 20, 19, 18, 17, 16, 15, 14, 13,
113  12, 11, 10, 9, 8, 7, 6, 5, 4, 3,
114  2, 1, 0
115  },
116  {
117  16384, 16116, 15849, 15582, 15316, 15050, 14785, 14521, 14257, 13995,
118  13734, 13476, 13218, 12963, 12708, 12457, 12206, 11956, 11706, 11460,
119  11215, 10975, 10735, 10500, 10265, 10034, 9803, 9579, 9355, 9136,
120  8917, 8703, 8489, 8275, 8061, 7853, 7645, 7444, 7244, 7051,
121  6858, 6671, 6484, 6305, 6127, 5956, 5785, 5622, 5459, 5298,
122  5137, 4983, 4830, 4684, 4539, 4401, 4263, 4131, 3999, 3874,
123  3750, 3632, 3515, 3401, 3287, 3173, 3059, 2949, 2840, 2737,
124  2635, 2539, 2444, 2354, 2264, 2181, 2098, 2020, 1943, 1872,
125  1801, 1731, 1661, 1596, 1532, 1472, 1412, 1357, 1303, 1251,
126  1200, 1153, 1106, 1063, 1020, 979, 938, 897, 856, 818,
127  780, 746, 712, 681, 650, 621, 592, 566, 540, 517,
128  494, 473, 452, 431, 410, 391, 373, 356, 340, 325,
129  310, 296, 282, 270, 258, 247, 236, 225, 214, 203,
130  192, 182, 172, 162, 153, 144, 136, 128, 121, 114,
131  108, 102, 97, 92, 87, 82, 77, 73, 69, 65,
132  62, 59, 56, 53, 50, 47, 45, 43, 41, 39,
133  37, 35, 33, 31, 29, 27, 26, 25, 24, 23,
134  22, 21, 20, 19, 18, 17, 16, 15, 14, 13,
135  12, 11, 10, 9, 8, 7, 6, 5, 4, 3,
136  2, 1, 0
137  },
138  {
139  16384, 16128, 15872, 15617, 15362, 15107, 14853, 14600, 14347, 14096,
140  13846, 13597, 13350, 13105, 12860, 12618, 12376, 12135, 11894, 11657,
141  11421, 11189, 10957, 10730, 10503, 10279, 10056, 9838, 9620, 9407,
142  9195, 8987, 8779, 8571, 8363, 8159, 7955, 7758, 7561, 7371,
143  7182, 6997, 6812, 6635, 6459, 6289, 6120, 5957, 5795, 5634,
144  5473, 5319, 5165, 5018, 4871, 4732, 4593, 4458, 4324, 4197,
145  4071, 3951, 3831, 3714, 3597, 3480, 3363, 3250, 3138, 3032,
146  2927, 2828, 2729, 2635, 2541, 2453, 2366, 2284, 2202, 2126,
147  2050, 1975, 1900, 1830, 1761, 1697, 1633, 1574, 1515, 1459,
148  1403, 1351, 1300, 1252, 1205, 1160, 1115, 1070, 1025, 982,
149  939, 899, 860, 824, 789, 756, 723, 693, 663, 636,
150  609, 584, 559, 535, 511, 489, 467, 447, 427, 409,
151  391, 374, 358, 343, 328, 313, 300, 287, 274, 261,
152  248, 235, 223, 211, 200, 189, 179, 169, 160, 151,
153  143, 135, 128, 121, 115, 109, 103, 97, 92, 87,
154  82, 77, 73, 69, 65, 61, 58, 55, 52, 49,
155  46, 43, 40, 37, 35, 33, 31, 29, 27, 25,
156  23, 21, 20, 19, 18, 17, 16, 15, 14, 13,
157  12, 11, 10, 9, 8, 7, 6, 5, 4, 3,
158  2, 1, 0
159  },
160  {
161  16384, 16139, 15894, 15649, 15405, 15162, 14919, 14677, 14435, 14195,
162  13955, 13717, 13479, 13243, 13008, 12775, 12542, 12310, 12079, 11851,
163  11623, 11399, 11176, 10956, 10737, 10521, 10305, 10094, 9883, 9677,
164  9471, 9268, 9065, 8862, 8659, 8459, 8260, 8067, 7874, 7688,
165  7502, 7321, 7140, 6965, 6790, 6621, 6452, 6290, 6128, 5968,
166  5808, 5655, 5503, 5356, 5209, 5069, 4929, 4794, 4660, 4532,
167  4404, 4282, 4160, 4041, 3922, 3803, 3684, 3568, 3452, 3343,
168  3234, 3131, 3029, 2931, 2833, 2741, 2649, 2563, 2477, 2396,
169  2316, 2236, 2157, 2083, 2009, 1940, 1871, 1807, 1743, 1683,
170  1623, 1567, 1511, 1459, 1407, 1357, 1307, 1257, 1207, 1159,
171  1111, 1067, 1023, 983, 943, 905, 868, 834, 800, 769,
172  738, 709, 681, 653, 625, 600, 575, 552, 529, 508,
173  487, 466, 447, 428, 410, 392, 376, 360, 344, 328,
174  313, 298, 283, 268, 255, 242, 230, 218, 207, 196,
175  186, 176, 167, 158, 150, 142, 135, 128, 121, 114,
176  108, 102, 97, 92, 87, 82, 78, 74, 70, 66,
177  62, 58, 54, 50, 47, 44, 41, 38, 35, 32,
178  30, 28, 26, 24, 22, 20, 18, 16, 14, 13,
179  12, 11, 10, 9, 8, 7, 6, 5, 4, 3,
180  2, 1, 0
181  },
182  {
183  16384, 16149, 15915, 15681, 15447, 15214, 14981, 14749, 14517, 14286,
184  14055, 13827, 13599, 13373, 13147, 12923, 12699, 12476, 12253, 12034,
185  11815, 11599, 11383, 11171, 10959, 10750, 10541, 10337, 10133, 9933,
186  9733, 9536, 9339, 9142, 8945, 8751, 8557, 8369, 8181, 7998,
187  7816, 7638, 7460, 7288, 7116, 6950, 6785, 6625, 6465, 6306,
188  6147, 5995, 5843, 5697, 5551, 5411, 5271, 5135, 5000, 4871,
189  4742, 4618, 4495, 4374, 4253, 4132, 4011, 3893, 3775, 3663,
190  3552, 3446, 3340, 3239, 3138, 3043, 2948, 2858, 2768, 2684,
191  2600, 2516, 2433, 2355, 2278, 2205, 2133, 2065, 1997, 1932,
192  1867, 1807, 1747, 1690, 1634, 1580, 1526, 1472, 1418, 1366,
193  1314, 1266, 1218, 1174, 1130, 1088, 1047, 1009, 971, 936,
194  901, 868, 836, 804, 772, 743, 714, 685, 658, 631,
195  606, 582, 559, 536, 515, 494, 475, 456, 437, 418,
196  399, 380, 362, 344, 328, 312, 297, 283, 270, 257,
197  245, 233, 222, 211, 201, 191, 181, 172, 163, 155,
198  147, 139, 132, 125, 119, 113, 107, 101, 96, 91,
199  86, 81, 76, 71, 66, 62, 58, 54, 50, 46,
200  43, 40, 37, 34, 31, 28, 26, 24, 22, 20,
201  18, 16, 14, 12, 10, 8, 6, 5, 4, 3,
202  2, 1, 0
203  },
204  {
205  16384, 16159, 15934, 15709, 15485, 15261, 15038, 14816, 14594, 14373,
206  14152, 13933, 13714, 13497, 13280, 13065, 12850, 12636, 12422, 12211,
207  12000, 11791, 11583, 11378, 11173, 10971, 10769, 10571, 10373, 10179,
208  9985, 9793, 9601, 9409, 9217, 9029, 8842, 8658, 8475, 8297,
209  8120, 7946, 7773, 7604, 7435, 7271, 7108, 6950, 6792, 6634,
210  6477, 6326, 6175, 6029, 5883, 5742, 5602, 5466, 5330, 5199,
211  5068, 4943, 4818, 4696, 4574, 4452, 4330, 4211, 4093, 3979,
212  3866, 3759, 3652, 3549, 3446, 3348, 3250, 3157, 3065, 2977,
213  2889, 2802, 2716, 2634, 2553, 2476, 2399, 2326, 2254, 2185,
214  2117, 2052, 1987, 1926, 1866, 1808, 1750, 1692, 1634, 1578,
215  1522, 1470, 1418, 1369, 1321, 1275, 1229, 1187, 1145, 1105,
216  1066, 1027, 991, 955, 919, 883, 850, 817, 786, 756,
217  728, 700, 674, 648, 624, 600, 578, 556, 534, 512,
218  490, 468, 447, 426, 407, 388, 371, 354, 338, 322,
219  307, 293, 280, 267, 255, 243, 231, 219, 209, 199,
220  189, 179, 170, 161, 153, 145, 138, 131, 124, 117,
221  111, 105, 99, 93, 87, 81, 76, 71, 66, 61,
222  57, 53, 49, 45, 42, 39, 36, 33, 30, 27,
223  24, 21, 19, 17, 15, 13, 11, 9, 7, 5,
224  3, 1, 0
225  },
226  {
227  16384, 16169, 15954, 15739, 15524, 15310, 15096, 14883, 14670, 14458,
228  14246, 14035, 13824, 13614, 13405, 13198, 12991, 12785, 12579, 12376,
229  12173, 11972, 11772, 11574, 11377, 11182, 10987, 10795, 10603, 10414,
230  10226, 10040, 9854, 9668, 9482, 9299, 9116, 8937, 8759, 8585,
231  8411, 8241, 8071, 7906, 7741, 7580, 7419, 7263, 7107, 6952,
232  6797, 6647, 6497, 6353, 6209, 6070, 5931, 5796, 5661, 5531,
233  5401, 5275, 5150, 5027, 4904, 4781, 4658, 4538, 4419, 4304,
234  4190, 4081, 3972, 3867, 3762, 3662, 3562, 3467, 3372, 3281,
235  3191, 3101, 3012, 2928, 2844, 2764, 2684, 2608, 2533, 2460,
236  2387, 2318, 2250, 2185, 2121, 2059, 1997, 1935, 1873, 1813,
237  1754, 1698, 1642, 1588, 1535, 1483, 1433, 1384, 1338, 1292,
238  1249, 1206, 1165, 1125, 1085, 1045, 1008, 971, 937, 903,
239  871, 840, 810, 780, 752, 724, 698, 672, 647, 622,
240  597, 572, 548, 524, 502, 480, 460, 440, 421, 403,
241  386, 369, 353, 337, 323, 309, 295, 281, 268, 255,
242  243, 231, 220, 209, 199, 189, 180, 171, 163, 155,
243  147, 139, 131, 123, 116, 109, 102, 95, 89, 83,
244  77, 72, 67, 62, 57, 52, 48, 44, 40, 36,
245  32, 28, 25, 22, 19, 16, 13, 10, 8, 6,
246  4, 2, 0
247  },
248  {
249  16384, 16177, 15970, 15764, 15558, 15353, 15148, 14944, 14740, 14537,
250  14334, 14132, 13930, 13729, 13529, 13330, 13131, 12933, 12735, 12539,
251  12343, 12150, 11957, 11766, 11576, 11388, 11200, 11015, 10830, 10647,
252  10465, 10285, 10105, 9925, 9745, 9568, 9391, 9218, 9045, 8876,
253  8707, 8541, 8375, 8213, 8051, 7894, 7737, 7583, 7429, 7277,
254  7125, 6977, 6830, 6687, 6544, 6406, 6268, 6133, 5998, 5868,
255  5738, 5612, 5487, 5364, 5241, 5118, 4995, 4875, 4755, 4640,
256  4525, 4414, 4304, 4198, 4092, 3990, 3888, 3790, 3693, 3600,
257  3507, 3415, 3323, 3235, 3147, 3064, 2981, 2902, 2823, 2746,
258  2670, 2594, 2522, 2450, 2382, 2314, 2248, 2182, 2116, 2050,
259  1987, 1924, 1864, 1804, 1748, 1692, 1638, 1585, 1534, 1484,
260  1437, 1390, 1346, 1302, 1258, 1215, 1174, 1133, 1095, 1057,
261  1021, 986, 952, 918, 887, 856, 827, 798, 770, 742,
262  714, 686, 659, 632, 607, 582, 559, 536, 514, 492,
263  472, 452, 433, 415, 398, 381, 364, 348, 333, 318,
264  304, 290, 277, 264, 252, 240, 229, 218, 208, 198,
265  188, 178, 168, 158, 149, 140, 132, 124, 116, 108,
266  101, 94, 87, 81, 75, 69, 64, 59, 54, 49,
267  44, 39, 35, 31, 27, 23, 19, 15, 12, 9,
268  6, 3, 0
269  }
270 };
271 
272 
273 static const uint16_t cf_tables_3[5][257] = {
274  {
275  16384, 16187, 15990, 15793, 15597, 15401, 15205, 15009, 14813, 14618,
276  14423, 14230, 14037, 13845, 13653, 13463, 13273, 13083, 12894, 12706,
277  12518, 12332, 12146, 11962, 11778, 11597, 11416, 11237, 11059, 10882,
278  10706, 10532, 10358, 10184, 10010, 9838, 9666, 9497, 9328, 9163,
279  8999, 8837, 8675, 8517, 8359, 8205, 8051, 7901, 7751, 7602,
280  7453, 7308, 7163, 7022, 6882, 6745, 6609, 6476, 6343, 6214,
281  6085, 5960, 5835, 5712, 5589, 5466, 5343, 5223, 5103, 4987,
282  4872, 4761, 4650, 4542, 4435, 4332, 4229, 4130, 4031, 3936,
283  3841, 3747, 3653, 3563, 3473, 3387, 3302, 3220, 3138, 3059,
284  2980, 2905, 2830, 2759, 2688, 2619, 2550, 2481, 2412, 2345,
285  2278, 2215, 2152, 2092, 2032, 1974, 1917, 1863, 1809, 1758,
286  1707, 1659, 1611, 1564, 1517, 1473, 1429, 1387, 1346, 1307,
287  1268, 1230, 1193, 1158, 1123, 1090, 1058, 1026, 994, 962,
288  930, 899, 869, 841, 813, 786, 760, 735, 710, 687,
289  664, 643, 622, 602, 582, 562, 543, 525, 507, 490,
290  473, 457, 442, 427, 412, 398, 385, 373, 361, 349,
291  337, 325, 313, 301, 290, 279, 269, 259, 249, 240,
292  231, 222, 214, 206, 199, 192, 185, 178, 171, 165,
293  159, 153, 148, 143, 138, 133, 128, 123, 119, 115,
294  111, 107, 103, 99, 95, 91, 87, 83, 80, 77,
295  74, 71, 68, 65, 63, 61, 59, 57, 55, 53,
296  51, 49, 47, 45, 43, 41, 40, 39, 38, 37,
297  36, 35, 34, 33, 32, 31, 30, 29, 28, 27,
298  26, 25, 24, 23, 22, 21, 20, 19, 18, 17,
299  16, 15, 14, 13, 12, 11, 10, 9, 8, 7,
300  6, 5, 4, 3, 2, 1, 0
301  },
302  {
303  16384, 16195, 16006, 15817, 15629, 15441, 15253, 15065, 14878, 14692,
304  14506, 14321, 14136, 13952, 13768, 13585, 13402, 13219, 13037, 12857,
305  12677, 12499, 12321, 12144, 11967, 11792, 11617, 11444, 11271, 11100,
306  10930, 10762, 10594, 10426, 10258, 10091, 9925, 9761, 9598, 9438,
307  9278, 9120, 8963, 8809, 8655, 8504, 8354, 8207, 8060, 7914,
308  7769, 7627, 7485, 7347, 7209, 7074, 6939, 6807, 6676, 6548,
309  6420, 6296, 6172, 6050, 5928, 5806, 5684, 5564, 5444, 5328,
310  5212, 5100, 4988, 4879, 4771, 4667, 4563, 4462, 4362, 4265,
311  4169, 4073, 3978, 3886, 3795, 3707, 3619, 3535, 3451, 3369,
312  3288, 3210, 3133, 3059, 2985, 2913, 2841, 2769, 2697, 2627,
313  2557, 2490, 2424, 2360, 2297, 2237, 2177, 2119, 2062, 2007,
314  1953, 1901, 1849, 1798, 1748, 1700, 1652, 1607, 1562, 1519,
315  1476, 1435, 1394, 1355, 1317, 1281, 1245, 1210, 1175, 1140,
316  1105, 1071, 1037, 1005, 973, 943, 913, 885, 857, 830,
317  804, 779, 754, 731, 708, 685, 663, 642, 621, 601,
318  581, 563, 545, 528, 511, 495, 479, 463, 448, 433,
319  419, 405, 391, 377, 364, 351, 338, 326, 314, 302,
320  291, 280, 270, 260, 251, 242, 234, 226, 218, 210,
321  202, 195, 188, 181, 174, 168, 162, 156, 150, 144,
322  139, 134, 129, 124, 119, 114, 109, 104, 100, 96,
323  92, 88, 84, 80, 77, 74, 71, 68, 65, 62,
324  59, 56, 54, 52, 50, 48, 46, 44, 42, 40,
325  38, 36, 34, 33, 32, 31, 30, 29, 28, 27,
326  26, 25, 24, 23, 22, 21, 20, 19, 18, 17,
327  16, 15, 14, 13, 12, 11, 10, 9, 8, 7,
328  6, 5, 4, 3, 2, 1, 0
329  },
330  {
331  16384, 16203, 16022, 15842, 15662, 15482, 15302, 15122, 14942, 14763,
332  14584, 14406, 14228, 14051, 13874, 13698, 13522, 13347, 13172, 12998,
333  12824, 12652, 12480, 12310, 12140, 11971, 11803, 11637, 11471, 11307,
334  11143, 10980, 10817, 10654, 10491, 10330, 10169, 10011, 9853, 9697,
335  9542, 9389, 9236, 9086, 8936, 8789, 8642, 8498, 8355, 8212,
336  8070, 7931, 7792, 7656, 7520, 7388, 7256, 7126, 6996, 6870,
337  6744, 6621, 6498, 6377, 6256, 6135, 6014, 5895, 5776, 5660,
338  5545, 5433, 5321, 5212, 5104, 4999, 4895, 4793, 4692, 4594,
339  4496, 4400, 4304, 4211, 4118, 4028, 3939, 3853, 3767, 3684,
340  3601, 3521, 3441, 3364, 3287, 3212, 3137, 3062, 2987, 2915,
341  2843, 2773, 2704, 2638, 2572, 2508, 2445, 2384, 2324, 2266,
342  2208, 2153, 2098, 2044, 1990, 1939, 1888, 1839, 1791, 1745,
343  1699, 1655, 1611, 1569, 1527, 1487, 1448, 1409, 1370, 1331,
344  1292, 1255, 1218, 1183, 1148, 1115, 1082, 1051, 1020, 990,
345  960, 932, 904, 878, 852, 826, 801, 777, 753, 731,
346  709, 687, 666, 645, 625, 605, 586, 567, 550, 533,
347  516, 499, 482, 465, 449, 433, 418, 403, 389, 375,
348  362, 349, 337, 325, 314, 303, 293, 283, 273, 263,
349  254, 245, 236, 227, 219, 211, 204, 197, 190, 183,
350  177, 171, 165, 159, 153, 147, 141, 135, 130, 125,
351  120, 115, 110, 105, 101, 97, 93, 89, 85, 81,
352  77, 74, 71, 68, 65, 62, 59, 56, 53, 51,
353  49, 47, 45, 43, 41, 39, 37, 35, 33, 31,
354  29, 27, 25, 23, 22, 21, 20, 19, 18, 17,
355  16, 15, 14, 13, 12, 11, 10, 9, 8, 7,
356  6, 5, 4, 3, 2, 1, 0
357  },
358  {
359  16384, 16210, 16036, 15863, 15690, 15517, 15344, 15172, 15000, 14828,
360  14656, 14485, 14314, 14145, 13976, 13808, 13640, 13472, 13304, 13137,
361  12970, 12804, 12639, 12475, 12312, 12149, 11987, 11827, 11667, 11508,
362  11349, 11192, 11035, 10878, 10721, 10565, 10410, 10257, 10104, 9953,
363  9802, 9654, 9506, 9359, 9213, 9070, 8927, 8787, 8647, 8508,
364  8369, 8233, 8097, 7964, 7831, 7700, 7570, 7442, 7315, 7190,
365  7065, 6943, 6821, 6701, 6581, 6461, 6341, 6223, 6105, 5990,
366  5876, 5764, 5653, 5545, 5437, 5331, 5226, 5124, 5022, 4924,
367  4826, 4729, 4632, 4538, 4444, 4353, 4262, 4174, 4087, 4002,
368  3917, 3835, 3753, 3674, 3595, 3518, 3441, 3364, 3287, 3212,
369  3138, 3066, 2995, 2926, 2858, 2792, 2726, 2662, 2599, 2538,
370  2478, 2420, 2362, 2305, 2249, 2195, 2141, 2089, 2037, 1988,
371  1939, 1891, 1844, 1799, 1754, 1711, 1668, 1626, 1584, 1542,
372  1500, 1459, 1418, 1380, 1342, 1305, 1269, 1234, 1199, 1166,
373  1133, 1102, 1071, 1041, 1012, 983, 954, 926, 899, 872,
374  847, 822, 798, 774, 751, 728, 707, 686, 666, 646,
375  627, 608, 589, 570, 552, 534, 517, 500, 484, 468,
376  453, 438, 424, 410, 397, 384, 372, 360, 348, 336,
377  325, 314, 303, 293, 283, 273, 264, 255, 246, 237,
378  229, 221, 213, 205, 197, 189, 181, 174, 167, 160,
379  154, 148, 142, 136, 131, 126, 121, 116, 111, 106,
380  101, 97, 93, 89, 85, 81, 77, 73, 70, 67,
381  64, 61, 58, 55, 52, 49, 46, 43, 40, 37,
382  35, 33, 31, 29, 27, 25, 23, 21, 19, 17,
383  16, 15, 14, 13, 12, 11, 10, 9, 8, 7,
384  6, 5, 4, 3, 2, 1, 0
385  },
386  {
387  16384, 16218, 16052, 15886, 15720, 15554, 15389, 15224, 15059, 14895,
388  14731, 14567, 14403, 14240, 14077, 13915, 13753, 13591, 13429, 13269,
389  13109, 12950, 12791, 12633, 12476, 12320, 12164, 12009, 11854, 11701,
390  11548, 11396, 11244, 11092, 10940, 10790, 10640, 10492, 10344, 10198,
391  10052, 9908, 9764, 9622, 9481, 9342, 9203, 9066, 8929, 8793,
392  8657, 8524, 8391, 8261, 8131, 8003, 7875, 7749, 7624, 7502,
393  7380, 7260, 7140, 7022, 6904, 6786, 6668, 6551, 6435, 6322,
394  6209, 6099, 5989, 5881, 5773, 5668, 5563, 5461, 5359, 5260,
395  5161, 5063, 4965, 4871, 4777, 4686, 4595, 4506, 4417, 4331,
396  4245, 4162, 4079, 3999, 3919, 3841, 3763, 3685, 3607, 3530,
397  3454, 3380, 3307, 3236, 3166, 3097, 3029, 2963, 2897, 2834,
398  2771, 2710, 2650, 2591, 2532, 2475, 2418, 2363, 2309, 2257,
399  2205, 2155, 2105, 2057, 2009, 1963, 1918, 1873, 1828, 1783,
400  1738, 1694, 1650, 1607, 1565, 1524, 1484, 1445, 1407, 1369,
401  1333, 1297, 1263, 1229, 1197, 1165, 1134, 1103, 1073, 1043,
402  1015, 987, 960, 933, 907, 882, 858, 834, 811, 788,
403  766, 744, 722, 700, 679, 658, 638, 618, 599, 581,
404  563, 545, 528, 511, 495, 480, 465, 451, 437, 423,
405  410, 397, 384, 372, 360, 348, 337, 326, 315, 305,
406  295, 285, 275, 265, 255, 245, 236, 227, 219, 211,
407  203, 195, 188, 181, 174, 167, 161, 155, 149, 143,
408  137, 131, 126, 121, 116, 111, 106, 101, 97, 93,
409  89, 85, 81, 77, 73, 69, 65, 61, 58, 55,
410  52, 49, 46, 43, 40, 37, 34, 32, 30, 28,
411  26, 24, 22, 20, 18, 16, 14, 12, 10, 8,
412  6, 5, 4, 3, 2, 1, 0
413  }
414 };
415 
416 
417 static const uint16_t *const cf_table[16] = {
418  cf_tables_1[0], cf_tables_1[1], cf_tables_1[2], cf_tables_2[0],
419  cf_tables_2[1], cf_tables_2[2], cf_tables_2[3], cf_tables_2[4],
420  cf_tables_2[5], cf_tables_2[6], cf_tables_2[7], cf_tables_3[0],
421  cf_tables_3[1], cf_tables_3[2], cf_tables_3[3], cf_tables_3[4]
422 };
423 
424 
425 /** Initialize a given lookup table using a given delta */
426 static void bgmc_lut_fillp(uint8_t *lut, int *lut_status, int delta)
427 {
428  unsigned int sx, i;
429 
430  for (sx = 0; sx < 16; sx++)
431  for (i = 0; i < LUT_SIZE; i++) {
432  unsigned int target = (i + 1) << (FREQ_BITS - LUT_BITS);
433  unsigned int symbol = 1 << delta;
434 
435  while (cf_table[sx][symbol] > target)
436  symbol += 1 << delta;
437 
438  *lut++ = symbol >> delta;
439  }
440 
441  *lut_status = delta;
442 }
443 
444 
445 /** Retune the index of a suitable lookup table for a given delta */
446 static uint8_t *bgmc_lut_getp(uint8_t *lut, int *lut_status, int delta)
447 {
448  unsigned int i = av_clip(delta, 0, LUT_BUFF - 1);
449 
450  lut += (i * LUT_SIZE) << 4;
451 
452  if (lut_status[i] != delta)
453  bgmc_lut_fillp(lut, &lut_status[i], delta);
454 
455  return lut;
456 }
457 
458 
459 /** Initialize the lookup table arrays */
461  uint8_t **cf_lut, int **cf_lut_status)
462 {
463  *cf_lut = av_malloc(sizeof(**cf_lut) * LUT_BUFF * 16 * LUT_SIZE);
464  *cf_lut_status = av_malloc(sizeof(**cf_lut_status) * LUT_BUFF);
465 
466  if (!*cf_lut || !*cf_lut_status) {
467  ff_bgmc_end(cf_lut, cf_lut_status);
468  av_log(avctx, AV_LOG_ERROR, "Allocating buffer memory failed.\n");
469  return AVERROR(ENOMEM);
470  } else {
471  // initialize lut_status buffer to a value never used to compare against
472  memset(*cf_lut_status, -1, sizeof(**cf_lut_status) * LUT_BUFF);
473  }
474 
475  return 0;
476 }
477 
478 
479 /** Release the lookup table arrays */
480 av_cold void ff_bgmc_end(uint8_t **cf_lut, int **cf_lut_status)
481 {
482  av_freep(cf_lut);
483  av_freep(cf_lut_status);
484 }
485 
486 
487 /** Initialize decoding and reads the first value */
488 void ff_bgmc_decode_init(GetBitContext *gb, unsigned int *h,
489  unsigned int *l, unsigned int *v)
490 {
491  *h = TOP_VALUE;
492  *l = 0;
493  *v = get_bits_long(gb, VALUE_BITS);
494 }
495 
496 
497 /** Finish decoding */
499 {
500  skip_bits_long(gb, -(VALUE_BITS - 2));
501 }
502 
503 
504 /** Read and decode a block Gilbert-Moore coded symbol */
505 void ff_bgmc_decode(GetBitContext *gb, unsigned int num, int32_t *dst,
506  int delta, unsigned int sx,
507  unsigned int *h, unsigned int *l, unsigned int *v,
508  uint8_t *cf_lut, int *cf_lut_status)
509 {
510  unsigned int i;
511  uint8_t *lut = bgmc_lut_getp(cf_lut, cf_lut_status, delta);
512 
513  // read current state
514  unsigned int high = *h;
515  unsigned int low = *l;
516  unsigned int value = *v;
517 
518  lut += sx * LUT_SIZE;
519 
520  // decode num samples
521  for (i = 0; i < num; i++) {
522  unsigned int range = high - low + 1;
523  unsigned int target = (((value - low + 1) << FREQ_BITS) - 1) / range;
524  unsigned int symbol = lut[target >> (FREQ_BITS - LUT_BITS)] << delta;
525 
526  while (cf_table[sx][symbol] > target)
527  symbol += 1 << delta;
528 
529  symbol = (symbol >> delta) - 1;
530 
531  high = low + ((range * cf_table[sx][(symbol) << delta] - (1 << FREQ_BITS)) >> FREQ_BITS);
532  low = low + ((range * cf_table[sx][(symbol + 1) << delta]) >> FREQ_BITS);
533 
534  while (1) {
535  if (high >= HALF) {
536  if (low >= HALF) {
537  value -= HALF;
538  low -= HALF;
539  high -= HALF;
540  } else if (low >= FIRST_QTR && high < THIRD_QTR) {
541  value -= FIRST_QTR;
542  low -= FIRST_QTR;
543  high -= FIRST_QTR;
544  } else
545  break;
546  }
547 
548  low *= 2;
549  high = 2 * high + 1;
550  value = 2 * value + get_bits1(gb);
551  }
552 
553  *dst++ = symbol;
554  }
555 
556  // save current state
557  *h = high;
558  *l = low;
559  *v = value;
560 }