Mercurial > libavcodec.hg
annotate cabac.h @ 3446:d4992082c7f1 libavcodec
forgot to remove adx.o from objs always built
author | mru |
---|---|
date | Sat, 08 Jul 2006 22:38:35 +0000 |
parents | a224d9752912 |
children | 0efda682253c |
rev | line source |
---|---|
1287 | 1 /* |
2 * H.26L/H.264/AVC/JVT/14496-10/... encoder/decoder | |
3 * Copyright (c) 2003 Michael Niedermayer <michaelni@gmx.at> | |
4 * | |
5 * This library is free software; you can redistribute it and/or | |
6 * modify it under the terms of the GNU Lesser General Public | |
7 * License as published by the Free Software Foundation; either | |
8 * version 2 of the License, or (at your option) any later version. | |
9 * | |
10 * This library is distributed in the hope that it will be useful, | |
11 * but WITHOUT ANY WARRANTY; without even the implied warranty of | |
12 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU | |
13 * Lesser General Public License for more details. | |
14 * | |
15 * You should have received a copy of the GNU Lesser General Public | |
16 * License along with this library; if not, write to the Free Software | |
3036
0b546eab515d
Update licensing information: The FSF changed postal address.
diego
parents:
2967
diff
changeset
|
17 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA |
1287 | 18 * |
19 */ | |
2967 | 20 |
1287 | 21 /** |
22 * @file cabac.h | |
23 * Context Adaptive Binary Arithmetic Coder. | |
24 */ | |
25 | |
26 | |
3284
a224d9752912
don't force asserts in release builds. 2% faster h264.
lorenm
parents:
3036
diff
changeset
|
27 //#undef NDEBUG |
1287 | 28 #include <assert.h> |
29 | |
2323 | 30 #define CABAC_BITS 8 |
31 #define CABAC_MASK ((1<<CABAC_BITS)-1) | |
32 | |
1287 | 33 typedef struct CABACContext{ |
34 int low; | |
35 int range; | |
36 int outstanding_count; | |
37 #ifdef STRICT_LIMITS | |
38 int symCount; | |
39 #endif | |
2323 | 40 uint8_t lps_range[2*65][4]; ///< rangeTabLPS |
1287 | 41 uint8_t lps_state[2*64]; ///< transIdxLPS |
42 uint8_t mps_state[2*64]; ///< transIdxMPS | |
2024
f65d87bfdd5a
some of the warning fixes by (Michael Roitzsch <mroi at users dot sourceforge dot net>)
michael
parents:
1787
diff
changeset
|
43 const uint8_t *bytestream_start; |
f65d87bfdd5a
some of the warning fixes by (Michael Roitzsch <mroi at users dot sourceforge dot net>)
michael
parents:
1787
diff
changeset
|
44 const uint8_t *bytestream; |
2116 | 45 const uint8_t *bytestream_end; |
1287 | 46 PutBitContext pb; |
47 }CABACContext; | |
48 | |
1301 | 49 extern const uint8_t ff_h264_lps_range[64][4]; |
50 extern const uint8_t ff_h264_mps_state[64]; | |
51 extern const uint8_t ff_h264_lps_state[64]; | |
2323 | 52 extern const uint8_t ff_h264_norm_shift[256]; |
53 | |
1287 | 54 |
55 void ff_init_cabac_encoder(CABACContext *c, uint8_t *buf, int buf_size); | |
2024
f65d87bfdd5a
some of the warning fixes by (Michael Roitzsch <mroi at users dot sourceforge dot net>)
michael
parents:
1787
diff
changeset
|
56 void ff_init_cabac_decoder(CABACContext *c, const uint8_t *buf, int buf_size); |
2967 | 57 void ff_init_cabac_states(CABACContext *c, uint8_t const (*lps_range)[4], |
1287 | 58 uint8_t const *mps_state, uint8_t const *lps_state, int state_count); |
59 | |
60 | |
61 static inline void put_cabac_bit(CABACContext *c, int b){ | |
2967 | 62 put_bits(&c->pb, 1, b); |
63 for(;c->outstanding_count; c->outstanding_count--){ | |
1287 | 64 put_bits(&c->pb, 1, 1-b); |
65 } | |
66 } | |
67 | |
68 static inline void renorm_cabac_encoder(CABACContext *c){ | |
69 while(c->range < 0x100){ | |
70 //FIXME optimize | |
71 if(c->low<0x100){ | |
72 put_cabac_bit(c, 0); | |
73 }else if(c->low<0x200){ | |
74 c->outstanding_count++; | |
75 c->low -= 0x100; | |
76 }else{ | |
77 put_cabac_bit(c, 1); | |
78 c->low -= 0x200; | |
79 } | |
2967 | 80 |
1287 | 81 c->range+= c->range; |
82 c->low += c->low; | |
83 } | |
84 } | |
85 | |
86 static inline void put_cabac(CABACContext *c, uint8_t * const state, int bit){ | |
2323 | 87 int RangeLPS= c->lps_range[*state][c->range>>6]; |
2967 | 88 |
1287 | 89 if(bit == ((*state)&1)){ |
90 c->range -= RangeLPS; | |
91 *state= c->mps_state[*state]; | |
92 }else{ | |
93 c->low += c->range - RangeLPS; | |
94 c->range = RangeLPS; | |
95 *state= c->lps_state[*state]; | |
96 } | |
2967 | 97 |
1287 | 98 renorm_cabac_encoder(c); |
99 | |
100 #ifdef STRICT_LIMITS | |
101 c->symCount++; | |
102 #endif | |
103 } | |
104 | |
105 static inline void put_cabac_static(CABACContext *c, int RangeLPS, int bit){ | |
106 assert(c->range > RangeLPS); | |
107 | |
108 if(!bit){ | |
109 c->range -= RangeLPS; | |
110 }else{ | |
111 c->low += c->range - RangeLPS; | |
112 c->range = RangeLPS; | |
113 } | |
114 | |
115 renorm_cabac_encoder(c); | |
116 | |
117 #ifdef STRICT_LIMITS | |
118 c->symCount++; | |
119 #endif | |
120 } | |
121 | |
1290 | 122 /** |
123 * @param bit 0 -> write zero bit, !=0 write one bit | |
124 */ | |
1287 | 125 static inline void put_cabac_bypass(CABACContext *c, int bit){ |
126 c->low += c->low; | |
127 | |
128 if(bit){ | |
129 c->low += c->range; | |
130 } | |
131 //FIXME optimize | |
132 if(c->low<0x200){ | |
133 put_cabac_bit(c, 0); | |
134 }else if(c->low<0x400){ | |
135 c->outstanding_count++; | |
136 c->low -= 0x200; | |
137 }else{ | |
138 put_cabac_bit(c, 1); | |
139 c->low -= 0x400; | |
140 } | |
2967 | 141 |
1287 | 142 #ifdef STRICT_LIMITS |
143 c->symCount++; | |
144 #endif | |
145 } | |
146 | |
1300
e18667d1e94d
FFV1 codec (our very simple lossless intra only codec, compresses much better then huffyuv)
michaelni
parents:
1298
diff
changeset
|
147 /** |
e18667d1e94d
FFV1 codec (our very simple lossless intra only codec, compresses much better then huffyuv)
michaelni
parents:
1298
diff
changeset
|
148 * |
e18667d1e94d
FFV1 codec (our very simple lossless intra only codec, compresses much better then huffyuv)
michaelni
parents:
1298
diff
changeset
|
149 * @return the number of bytes written |
e18667d1e94d
FFV1 codec (our very simple lossless intra only codec, compresses much better then huffyuv)
michaelni
parents:
1298
diff
changeset
|
150 */ |
e18667d1e94d
FFV1 codec (our very simple lossless intra only codec, compresses much better then huffyuv)
michaelni
parents:
1298
diff
changeset
|
151 static inline int put_cabac_terminate(CABACContext *c, int bit){ |
1287 | 152 c->range -= 2; |
153 | |
154 if(!bit){ | |
155 renorm_cabac_encoder(c); | |
156 }else{ | |
157 c->low += c->range; | |
158 c->range= 2; | |
2967 | 159 |
1287 | 160 renorm_cabac_encoder(c); |
161 | |
162 assert(c->low <= 0x1FF); | |
163 put_cabac_bit(c, c->low>>9); | |
164 put_bits(&c->pb, 2, ((c->low>>7)&3)|1); | |
2967 | 165 |
1287 | 166 flush_put_bits(&c->pb); //FIXME FIXME FIXME XXX wrong |
167 } | |
2967 | 168 |
1287 | 169 #ifdef STRICT_LIMITS |
170 c->symCount++; | |
171 #endif | |
1300
e18667d1e94d
FFV1 codec (our very simple lossless intra only codec, compresses much better then huffyuv)
michaelni
parents:
1298
diff
changeset
|
172 |
1787 | 173 return (put_bits_count(&c->pb)+7)>>3; |
1287 | 174 } |
175 | |
1290 | 176 /** |
177 * put (truncated) unary binarization. | |
178 */ | |
179 static inline void put_cabac_u(CABACContext *c, uint8_t * state, int v, int max, int max_index, int truncated){ | |
180 int i; | |
2967 | 181 |
1290 | 182 assert(v <= max); |
2967 | 183 |
1290 | 184 #if 1 |
185 for(i=0; i<v; i++){ | |
186 put_cabac(c, state, 1); | |
187 if(i < max_index) state++; | |
188 } | |
189 if(truncated==0 || v<max) | |
190 put_cabac(c, state, 0); | |
191 #else | |
192 if(v <= max_index){ | |
193 for(i=0; i<v; i++){ | |
194 put_cabac(c, state+i, 1); | |
195 } | |
196 if(truncated==0 || v<max) | |
197 put_cabac(c, state+i, 0); | |
198 }else{ | |
199 for(i=0; i<=max_index; i++){ | |
200 put_cabac(c, state+i, 1); | |
201 } | |
202 for(; i<v; i++){ | |
203 put_cabac(c, state+max_index, 1); | |
204 } | |
205 if(truncated==0 || v<max) | |
206 put_cabac(c, state+max_index, 0); | |
207 } | |
208 #endif | |
209 } | |
210 | |
211 /** | |
212 * put unary exp golomb k-th order binarization. | |
213 */ | |
1298 | 214 static inline void put_cabac_ueg(CABACContext *c, uint8_t * state, int v, int max, int is_signed, int k, int max_index){ |
1290 | 215 int i; |
2967 | 216 |
1290 | 217 if(v==0) |
218 put_cabac(c, state, 0); | |
219 else{ | |
1298 | 220 const int sign= v < 0; |
2967 | 221 |
1298 | 222 if(is_signed) v= ABS(v); |
2967 | 223 |
1290 | 224 if(v<max){ |
225 for(i=0; i<v; i++){ | |
226 put_cabac(c, state, 1); | |
227 if(i < max_index) state++; | |
228 } | |
229 | |
230 put_cabac(c, state, 0); | |
231 }else{ | |
232 int m= 1<<k; | |
233 | |
234 for(i=0; i<max; i++){ | |
235 put_cabac(c, state, 1); | |
236 if(i < max_index) state++; | |
237 } | |
238 | |
239 v -= max; | |
240 while(v >= m){ //FIXME optimize | |
241 put_cabac_bypass(c, 1); | |
242 v-= m; | |
243 m+= m; | |
244 } | |
245 put_cabac_bypass(c, 0); | |
246 while(m>>=1){ | |
247 put_cabac_bypass(c, v&m); | |
248 } | |
249 } | |
250 | |
251 if(is_signed) | |
252 put_cabac_bypass(c, sign); | |
253 } | |
254 } | |
255 | |
2323 | 256 static void refill(CABACContext *c){ |
2736 | 257 if(c->bytestream <= c->bytestream_end) |
2323 | 258 #if CABAC_BITS == 16 |
259 c->low+= ((c->bytestream[0]<<9) + (c->bytestream[1])<<1); | |
260 #else | |
261 c->low+= c->bytestream[0]<<1; | |
262 #endif | |
263 c->low -= CABAC_MASK; | |
264 c->bytestream+= CABAC_BITS/8; | |
265 } | |
266 | |
2522
e25782262d7d
kill warnings patch by (M«©ns Rullg«©rd <mru inprovide com>)
michael
parents:
2323
diff
changeset
|
267 #if 0 /* all use commented */ |
2323 | 268 static void refill2(CABACContext *c){ |
269 int i, x; | |
270 | |
271 x= c->low ^ (c->low-1); | |
272 i= 8 - ff_h264_norm_shift[x>>(CABAC_BITS+1)]; | |
273 | |
274 x= -CABAC_MASK; | |
2967 | 275 |
2323 | 276 if(c->bytestream < c->bytestream_end) |
277 #if CABAC_BITS == 16 | |
278 x+= (c->bytestream[0]<<9) + (c->bytestream[1]<<1); | |
279 #else | |
280 x+= c->bytestream[0]<<1; | |
281 #endif | |
2967 | 282 |
2323 | 283 c->low += x<<i; |
284 c->bytestream+= CABAC_BITS/8; | |
285 } | |
2522
e25782262d7d
kill warnings patch by (M«©ns Rullg«©rd <mru inprovide com>)
michael
parents:
2323
diff
changeset
|
286 #endif |
2323 | 287 |
1287 | 288 static inline void renorm_cabac_decoder(CABACContext *c){ |
2323 | 289 while(c->range < (0x200 << CABAC_BITS)){ |
1287 | 290 c->range+= c->range; |
291 c->low+= c->low; | |
2323 | 292 if(!(c->low & CABAC_MASK)) |
293 refill(c); | |
1287 | 294 } |
295 } | |
296 | |
2323 | 297 static inline void renorm_cabac_decoder_once(CABACContext *c){ |
298 int mask= (c->range - (0x200 << CABAC_BITS))>>31; | |
299 c->range+= c->range&mask; | |
300 c->low += c->low &mask; | |
301 if(!(c->low & CABAC_MASK)) | |
302 refill(c); | |
303 } | |
304 | |
1287 | 305 static inline int get_cabac(CABACContext *c, uint8_t * const state){ |
2323 | 306 int RangeLPS= c->lps_range[*state][c->range>>(CABAC_BITS+7)]<<(CABAC_BITS+1); |
2522
e25782262d7d
kill warnings patch by (M«©ns Rullg«©rd <mru inprovide com>)
michael
parents:
2323
diff
changeset
|
307 int bit, lps_mask attribute_unused; |
2967 | 308 |
1287 | 309 c->range -= RangeLPS; |
2323 | 310 #if 1 |
1287 | 311 if(c->low < c->range){ |
312 bit= (*state)&1; | |
313 *state= c->mps_state[*state]; | |
2323 | 314 renorm_cabac_decoder_once(c); |
1287 | 315 }else{ |
2323 | 316 // int shift= ff_h264_norm_shift[RangeLPS>>17]; |
1287 | 317 bit= ((*state)&1)^1; |
318 c->low -= c->range; | |
2323 | 319 *state= c->lps_state[*state]; |
1287 | 320 c->range = RangeLPS; |
2323 | 321 renorm_cabac_decoder(c); |
322 /* c->range = RangeLPS<<shift; | |
323 c->low <<= shift; | |
324 if(!(c->low & 0xFFFF)){ | |
325 refill2(c); | |
326 }*/ | |
1287 | 327 } |
2323 | 328 #else |
329 lps_mask= (c->range - c->low)>>31; | |
2967 | 330 |
2323 | 331 c->low -= c->range & lps_mask; |
332 c->range += (RangeLPS - c->range) & lps_mask; | |
2967 | 333 |
2323 | 334 bit= ((*state)^lps_mask)&1; |
335 *state= c->mps_state[(*state) - (128&lps_mask)]; | |
2967 | 336 |
2323 | 337 lps_mask= ff_h264_norm_shift[c->range>>(CABAC_BITS+2)]; |
338 c->range<<= lps_mask; | |
339 c->low <<= lps_mask; | |
340 if(!(c->low & CABAC_MASK)) | |
341 refill2(c); | |
342 #endif | |
343 | |
2967 | 344 return bit; |
1287 | 345 } |
346 | |
347 static inline int get_cabac_bypass(CABACContext *c){ | |
348 c->low += c->low; | |
349 | |
2323 | 350 if(!(c->low & CABAC_MASK)) |
351 refill(c); | |
2967 | 352 |
1287 | 353 if(c->low < c->range){ |
354 return 0; | |
355 }else{ | |
356 c->low -= c->range; | |
357 return 1; | |
358 } | |
359 } | |
360 | |
1300
e18667d1e94d
FFV1 codec (our very simple lossless intra only codec, compresses much better then huffyuv)
michaelni
parents:
1298
diff
changeset
|
361 /** |
e18667d1e94d
FFV1 codec (our very simple lossless intra only codec, compresses much better then huffyuv)
michaelni
parents:
1298
diff
changeset
|
362 * |
e18667d1e94d
FFV1 codec (our very simple lossless intra only codec, compresses much better then huffyuv)
michaelni
parents:
1298
diff
changeset
|
363 * @return the number of bytes read or 0 if no end |
e18667d1e94d
FFV1 codec (our very simple lossless intra only codec, compresses much better then huffyuv)
michaelni
parents:
1298
diff
changeset
|
364 */ |
1287 | 365 static inline int get_cabac_terminate(CABACContext *c){ |
2323 | 366 c->range -= 4<<CABAC_BITS; |
1287 | 367 if(c->low < c->range){ |
2323 | 368 renorm_cabac_decoder_once(c); |
1287 | 369 return 0; |
370 }else{ | |
1300
e18667d1e94d
FFV1 codec (our very simple lossless intra only codec, compresses much better then huffyuv)
michaelni
parents:
1298
diff
changeset
|
371 return c->bytestream - c->bytestream_start; |
2967 | 372 } |
1287 | 373 } |
374 | |
1290 | 375 /** |
376 * get (truncated) unnary binarization. | |
377 */ | |
378 static inline int get_cabac_u(CABACContext *c, uint8_t * state, int max, int max_index, int truncated){ | |
379 int i; | |
2967 | 380 |
381 for(i=0; i<max; i++){ | |
1290 | 382 if(get_cabac(c, state)==0) |
383 return i; | |
2967 | 384 |
1290 | 385 if(i< max_index) state++; |
386 } | |
387 | |
388 return truncated ? max : -1; | |
389 } | |
390 | |
391 /** | |
392 * get unary exp golomb k-th order binarization. | |
393 */ | |
394 static inline int get_cabac_ueg(CABACContext *c, uint8_t * state, int max, int is_signed, int k, int max_index){ | |
395 int i, v; | |
396 int m= 1<<k; | |
2967 | 397 |
398 if(get_cabac(c, state)==0) | |
1290 | 399 return 0; |
2967 | 400 |
1290 | 401 if(0 < max_index) state++; |
2967 | 402 |
403 for(i=1; i<max; i++){ | |
1290 | 404 if(get_cabac(c, state)==0){ |
405 if(is_signed && get_cabac_bypass(c)){ | |
406 return -i; | |
407 }else | |
408 return i; | |
409 } | |
410 | |
411 if(i < max_index) state++; | |
412 } | |
2967 | 413 |
1290 | 414 while(get_cabac_bypass(c)){ |
415 i+= m; | |
416 m+= m; | |
417 } | |
2967 | 418 |
1290 | 419 v=0; |
420 while(m>>=1){ | |
421 v+= v + get_cabac_bypass(c); | |
422 } | |
423 i += v; | |
424 | |
425 if(is_signed && get_cabac_bypass(c)){ | |
426 return -i; | |
427 }else | |
428 return i; | |
429 } |