Imported Debian version 2.5.3~trusty1
[deb_ffmpeg.git] / ffmpeg / libavcodec / bgmc.c
1 /*
2 * Block Gilbert-Moore decoder
3 * Copyright (c) 2010 Thilo Borgmann <thilo.borgmann _at_ mail.de>
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_ mail.de>
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 */
460 av_cold int ff_bgmc_init(AVCodecContext *avctx,
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 */
498 void ff_bgmc_decode_end(GetBitContext *gb)
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 }