Statistics
| Branch: | Revision:

ffmpeg / libavcodec / golomb.h @ 8a01fc47

History | View | Annotate | Download (10.5 KB)

1 0da71265 Michael Niedermayer
/*
2
 * exp golomb vlc stuff
3
 * Copyright (c) 2003 Michael Niedermayer <michaelni@gmx.at>
4 4579adb0 Alex Beregszaszi
 * Copyright (c) 2004 Alex Beregszaszi
5 0da71265 Michael Niedermayer
 *
6
 * This library is free software; you can redistribute it and/or
7
 * modify it under the terms of the GNU Lesser General Public
8
 * License as published by the Free Software Foundation; either
9
 * version 2 of the License, or (at your option) any later version.
10
 *
11
 * This library is distributed in the hope that it will be useful,
12
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
13
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
14
 * Lesser General Public License for more details.
15
 *
16
 * You should have received a copy of the GNU Lesser General Public
17
 * License along with this library; if not, write to the Free Software
18
 * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
19
 *
20
 */
21
 
22
/**
23
 * @file golomb.h
24
 * @brief 
25
 *     exp golomb vlc stuff
26 4579adb0 Alex Beregszaszi
 * @author Michael Niedermayer <michaelni@gmx.at> and Alex Beregszaszi
27 0da71265 Michael Niedermayer
 */
28
29 8b82a956 Michael Niedermayer
#define INVALID_VLC           0x80000000
30
31 0da71265 Michael Niedermayer
extern const uint8_t ff_golomb_vlc_len[512];
32
extern const uint8_t ff_ue_golomb_vlc_code[512];
33
extern const  int8_t ff_se_golomb_vlc_code[512];
34
extern const uint8_t ff_ue_golomb_len[256];
35
36 4704097a Michael Niedermayer
extern const uint8_t ff_interleaved_golomb_vlc_len[256];
37
extern const uint8_t ff_interleaved_ue_golomb_vlc_code[256];
38
extern const  int8_t ff_interleaved_se_golomb_vlc_code[256];
39
40 0da71265 Michael Niedermayer
 
41
 /**
42
 * read unsigned exp golomb code.
43
 */
44
static inline int get_ue_golomb(GetBitContext *gb){
45
    unsigned int buf;
46
    int log;
47
    
48
    OPEN_READER(re, gb);
49
    UPDATE_CACHE(re, gb);
50
    buf=GET_CACHE(re, gb);
51
    
52
    if(buf >= (1<<27)){
53
        buf >>= 32 - 9;
54
        LAST_SKIP_BITS(re, gb, ff_golomb_vlc_len[buf]);
55
        CLOSE_READER(re, gb);
56
    
57
        return ff_ue_golomb_vlc_code[buf];
58
    }else{
59
        log= 2*av_log2(buf) - 31;
60
        buf>>= log;
61
        buf--;
62
        LAST_SKIP_BITS(re, gb, 32 - log);
63
        CLOSE_READER(re, gb);
64
    
65
        return buf;
66
    }
67
}
68
69 8b82a956 Michael Niedermayer
static inline int svq3_get_ue_golomb(GetBitContext *gb){
70 4704097a Michael Niedermayer
    uint32_t buf;
71 8b82a956 Michael Niedermayer
    int log;
72
73
    OPEN_READER(re, gb);
74
    UPDATE_CACHE(re, gb);
75 4704097a Michael Niedermayer
    buf=GET_CACHE(re, gb);
76
    
77
    if(buf&0xAA800000){
78
        buf >>= 32 - 8;
79
        LAST_SKIP_BITS(re, gb, ff_interleaved_golomb_vlc_len[buf]);
80
        CLOSE_READER(re, gb);
81
        
82
        return ff_interleaved_ue_golomb_vlc_code[buf];
83
    }else{
84 cceb990e Michael Niedermayer
        LAST_SKIP_BITS(re, gb, 8);
85
        UPDATE_CACHE(re, gb);
86
        buf |= 1 | (GET_CACHE(re, gb) >> 8);
87
88 4704097a Michael Niedermayer
        if((buf & 0xAAAAAAAA) == 0)
89
            return INVALID_VLC;
90 8b82a956 Michael Niedermayer
91 4704097a Michael Niedermayer
        for(log=31; (buf & 0x80000000) == 0; log--){
92
            buf = (buf << 2) - ((buf << log) >> (log - 1)) + (buf >> 30);
93
        }
94 8b82a956 Michael Niedermayer
95 cceb990e Michael Niedermayer
        LAST_SKIP_BITS(re, gb, 63 - 2*log - 8);
96 4704097a Michael Niedermayer
        CLOSE_READER(re, gb);
97 8b82a956 Michael Niedermayer
98 4704097a Michael Niedermayer
        return ((buf << log) >> log) - 1;
99
    }
100 8b82a956 Michael Niedermayer
}
101
102 0da71265 Michael Niedermayer
/**
103
 * read unsigned truncated exp golomb code.
104
 */
105
static inline int get_te0_golomb(GetBitContext *gb, int range){
106
    assert(range >= 1);
107
    
108
    if(range==1)      return 0;
109 1924f3ce Michael Niedermayer
    else if(range==2) return get_bits1(gb)^1;
110 0da71265 Michael Niedermayer
    else              return get_ue_golomb(gb);
111
}
112
113
/**
114
 * read unsigned truncated exp golomb code.
115
 */
116
static inline int get_te_golomb(GetBitContext *gb, int range){
117
    assert(range >= 1);
118
    
119 1924f3ce Michael Niedermayer
    if(range==2) return get_bits1(gb)^1;
120 0da71265 Michael Niedermayer
    else         return get_ue_golomb(gb);
121
}
122
123
124
/**
125
 * read signed exp golomb code.
126
 */
127
static inline int get_se_golomb(GetBitContext *gb){
128
    unsigned int buf;
129
    int log;
130
    
131
    OPEN_READER(re, gb);
132
    UPDATE_CACHE(re, gb);
133
    buf=GET_CACHE(re, gb);
134
    
135
    if(buf >= (1<<27)){
136
        buf >>= 32 - 9;
137
        LAST_SKIP_BITS(re, gb, ff_golomb_vlc_len[buf]);
138
        CLOSE_READER(re, gb);
139
    
140
        return ff_se_golomb_vlc_code[buf];
141
    }else{
142
        log= 2*av_log2(buf) - 31;
143
        buf>>= log;
144
        
145
        LAST_SKIP_BITS(re, gb, 32 - log);
146
        CLOSE_READER(re, gb);
147
    
148
        if(buf&1) buf= -(buf>>1);
149
        else      buf=  (buf>>1);
150
151
        return buf;
152
    }
153
}
154
155 8b82a956 Michael Niedermayer
static inline int svq3_get_se_golomb(GetBitContext *gb){
156
    unsigned int buf;
157
    int log;
158
159
    OPEN_READER(re, gb);
160
    UPDATE_CACHE(re, gb);
161 4704097a Michael Niedermayer
    buf=GET_CACHE(re, gb);
162 8b82a956 Michael Niedermayer
163 4704097a Michael Niedermayer
    if(buf&0xAA800000){
164
        buf >>= 32 - 8;
165
        LAST_SKIP_BITS(re, gb, ff_interleaved_golomb_vlc_len[buf]);
166
        CLOSE_READER(re, gb);
167
        
168
        return ff_interleaved_se_golomb_vlc_code[buf];
169
    }else{
170
        buf |=1;
171
        if((buf & 0xAAAAAAAA) == 0)
172
            return INVALID_VLC;
173 8b82a956 Michael Niedermayer
174 4704097a Michael Niedermayer
        for(log=31; (buf & 0x80000000) == 0; log--){
175
            buf = (buf << 2) - ((buf << log) >> (log - 1)) + (buf >> 30);
176
        }
177 8b82a956 Michael Niedermayer
178 4704097a Michael Niedermayer
        LAST_SKIP_BITS(re, gb, 63 - 2*log);
179
        CLOSE_READER(re, gb);
180 8b82a956 Michael Niedermayer
181 4704097a Michael Niedermayer
        return (signed) (((((buf << log) >> log) - 1) ^ -(buf & 0x1)) + 1) >> 1;
182
    }
183 8b82a956 Michael Niedermayer
}
184
185 11e659c2 Michael Niedermayer
/**
186 d9e6a6c6 Michael Niedermayer
 * read unsigned golomb rice code (ffv1).
187 11e659c2 Michael Niedermayer
 */
188
static inline int get_ur_golomb(GetBitContext *gb, int k, int limit, int esc_len){
189
    unsigned int buf;
190
    int log;
191
    
192
    OPEN_READER(re, gb);
193
    UPDATE_CACHE(re, gb);
194
    buf=GET_CACHE(re, gb);
195
196
    log= av_log2(buf);
197 d9e6a6c6 Michael Niedermayer
198 11e659c2 Michael Niedermayer
    if(log > 31-limit){
199
        buf >>= log - k;
200
        buf += (30-log)<<k;
201
        LAST_SKIP_BITS(re, gb, 32 + k - log);
202
        CLOSE_READER(re, gb);
203
    
204
        return buf;
205 d9e6a6c6 Michael Niedermayer
    }else{
206
        buf >>= 32 - limit - esc_len;
207
        LAST_SKIP_BITS(re, gb, esc_len + limit);
208 11e659c2 Michael Niedermayer
        CLOSE_READER(re, gb);
209
    
210 d9e6a6c6 Michael Niedermayer
        return buf + limit - 1;
211
    }
212
}
213
214
/**
215
 * read unsigned golomb rice code (jpegls).
216
 */
217
static inline int get_ur_golomb_jpegls(GetBitContext *gb, int k, int limit, int esc_len){
218
    unsigned int buf;
219
    int log;
220
    
221
    OPEN_READER(re, gb);
222
    UPDATE_CACHE(re, gb);
223
    buf=GET_CACHE(re, gb);
224
225
    log= av_log2(buf);
226
    
227 ce329952 Michael Niedermayer
    if(log > 31-11){
228 d9e6a6c6 Michael Niedermayer
        buf >>= log - k;
229
        buf += (30-log)<<k;
230
        LAST_SKIP_BITS(re, gb, 32 + k - log);
231
        CLOSE_READER(re, gb);
232
    
233
        return buf;
234
    }else{
235
        int i;
236
        for(i=0; SHOW_UBITS(re, gb, 1) == 0; i++){
237
            LAST_SKIP_BITS(re, gb, 1);
238
            UPDATE_CACHE(re, gb);
239
        }
240
        SKIP_BITS(re, gb, 1);
241
242
        if(i < limit - 1){
243
            if(k){
244
                buf = SHOW_UBITS(re, gb, k);
245
                LAST_SKIP_BITS(re, gb, k);
246
            }else{
247
                buf=0;
248
            }
249
250
            CLOSE_READER(re, gb);
251
            return buf + (i<<k);
252
        }else if(i == limit - 1){
253
            buf = SHOW_UBITS(re, gb, esc_len);
254
            LAST_SKIP_BITS(re, gb, esc_len);
255
            CLOSE_READER(re, gb);
256
    
257
            return buf + 1;
258
        }else
259
            return -1;
260
    }
261 11e659c2 Michael Niedermayer
}
262
263 ac2570a8 Michael Niedermayer
/**
264 4579adb0 Alex Beregszaszi
 * read signed golomb rice code (ffv1).
265
 */
266 a6c01275 Michael Niedermayer
static inline int get_sr_golomb(GetBitContext *gb, int k, int limit, int esc_len){
267 4579adb0 Alex Beregszaszi
    int v= get_ur_golomb(gb, k, limit, esc_len);
268
    
269
    v++;
270
    if (v&1) return v>>1;
271
    else return -(v>>1);
272
    
273
//    return (v>>1) ^ -(v&1);
274
}
275
276 a6c01275 Michael Niedermayer
/**
277 4579adb0 Alex Beregszaszi
 * read signed golomb rice code (flac).
278 ac2570a8 Michael Niedermayer
 */
279
static inline int get_sr_golomb_flac(GetBitContext *gb, int k, int limit, int esc_len){
280
    int v= get_ur_golomb_jpegls(gb, k, limit, esc_len);
281
    return (v>>1) ^ -(v&1);
282
}
283
284 0da71265 Michael Niedermayer
#ifdef TRACE
285
286
static inline int get_ue(GetBitContext *s, char *file, char *func, int line){
287
    int show= show_bits(s, 24);
288
    int pos= get_bits_count(s);
289
    int i= get_ue_golomb(s);
290
    int len= get_bits_count(s) - pos;
291
    int bits= show>>(24-len);
292
    
293
    print_bin(bits, len);
294
    
295 4579adb0 Alex Beregszaszi
    av_log(NULL, AV_LOG_DEBUG, "%5d %2d %3d ue  @%5d in %s %s:%d\n", bits, len, i, pos, file, func, line);
296 0da71265 Michael Niedermayer
    
297
    return i;
298
}
299
300
static inline int get_se(GetBitContext *s, char *file, char *func, int line){
301
    int show= show_bits(s, 24);
302
    int pos= get_bits_count(s);
303
    int i= get_se_golomb(s);
304
    int len= get_bits_count(s) - pos;
305
    int bits= show>>(24-len);
306
    
307
    print_bin(bits, len);
308
    
309 4579adb0 Alex Beregszaszi
    av_log(NULL, AV_LOG_DEBUG, "%5d %2d %3d se  @%5d in %s %s:%d\n", bits, len, i, pos, file, func, line);
310 0da71265 Michael Niedermayer
    
311
    return i;
312
}
313
314
static inline int get_te(GetBitContext *s, int r, char *file, char *func, int line){
315
    int show= show_bits(s, 24);
316
    int pos= get_bits_count(s);
317
    int i= get_te0_golomb(s, r);
318
    int len= get_bits_count(s) - pos;
319
    int bits= show>>(24-len);
320
    
321
    print_bin(bits, len);
322
    
323 4579adb0 Alex Beregszaszi
    av_log(NULL, AV_LOG_DEBUG, "%5d %2d %3d te  @%5d in %s %s:%d\n", bits, len, i, pos, file, func, line);
324 0da71265 Michael Niedermayer
    
325
    return i;
326
}
327
328
#define get_ue_golomb(a) get_ue(a, __FILE__, __PRETTY_FUNCTION__, __LINE__)
329
#define get_se_golomb(a) get_se(a, __FILE__, __PRETTY_FUNCTION__, __LINE__)
330
#define get_te_golomb(a, r) get_te(a, r, __FILE__, __PRETTY_FUNCTION__, __LINE__)
331
#define get_te0_golomb(a, r) get_te(a, r, __FILE__, __PRETTY_FUNCTION__, __LINE__)
332
333
#endif
334
335
/**
336
 * write unsigned exp golomb code.
337
 */
338
static inline void set_ue_golomb(PutBitContext *pb, int i){
339
    int e;
340
    
341
    assert(i>=0);
342
343
#if 0
344
    if(i=0){
345
        put_bits(pb, 1, 1);
346
        return;
347
    }
348
#endif
349
    if(i<256)
350
        put_bits(pb, ff_ue_golomb_len[i], i+1);
351
    else{
352
        e= av_log2(i+1);
353
    
354
        put_bits(pb, 2*e+1, i+1);
355
    }
356
}
357
358
/**
359
 * write truncated unsigned exp golomb code.
360
 */
361
static inline void set_te_golomb(PutBitContext *pb, int i, int range){
362
    assert(range >= 1);
363
    assert(i<=range);
364
365 1924f3ce Michael Niedermayer
    if(range==2) put_bits(pb, 1, i^1);
366 0da71265 Michael Niedermayer
    else         set_ue_golomb(pb, i);
367
}
368
369
/**
370
 * write signed exp golomb code.
371
 */
372
static inline void set_se_golomb(PutBitContext *pb, int i){
373
#if 0 
374
    if(i<=0) i= -2*i;
375
    else     i=  2*i-1;
376
#elif 1
377
    i= 2*i-1;
378
    if(i<0) i^= -1; //FIXME check if gcc does the right thing
379
#else
380
    i= 2*i-1;
381
    i^= (i>>31);
382
#endif
383
    set_ue_golomb(pb, i);
384
}
385 11e659c2 Michael Niedermayer
386
/**
387 d9e6a6c6 Michael Niedermayer
 * write unsigned golomb rice code (ffv1).
388 11e659c2 Michael Niedermayer
 */
389
static inline void set_ur_golomb(PutBitContext *pb, int i, int k, int limit, int esc_len){
390
    int e;
391
    
392
    assert(i>=0);
393
    
394
    e= i>>k;
395
    if(e<limit){
396
        put_bits(pb, e + k + 1, (1<<k) + (i&((1<<k)-1)));
397
    }else{
398 d9e6a6c6 Michael Niedermayer
        put_bits(pb, limit + esc_len, i - limit + 1);
399
    }
400
}
401
402
/**
403
 * write unsigned golomb rice code (jpegls).
404
 */
405
static inline void set_ur_golomb_jpegls(PutBitContext *pb, int i, int k, int limit, int esc_len){
406
    int e;
407
    
408
    assert(i>=0);
409
    
410
    e= (i>>k) + 1;
411
    if(e<limit){
412
        put_bits(pb, e, 1);
413 ce329952 Michael Niedermayer
        if(k)
414
            put_bits(pb, k, i&((1<<k)-1));
415 d9e6a6c6 Michael Niedermayer
    }else{
416
        put_bits(pb, limit  , 1);
417
        put_bits(pb, esc_len, i - 1);
418 11e659c2 Michael Niedermayer
    }
419
}
420 4579adb0 Alex Beregszaszi
421
/**
422
 * write signed golomb rice code (ffv1).
423
 */
424 a6c01275 Michael Niedermayer
static inline void set_sr_golomb(PutBitContext *pb, int i, int k, int limit, int esc_len){
425 4579adb0 Alex Beregszaszi
    int v;
426
427
    v = -2*i-1;
428
    v ^= (v>>31);
429
430
    set_ur_golomb(pb, v, k, limit, esc_len);
431
}
432
433
/**
434
 * write signed golomb rice code (flac).
435
 */
436
static inline void set_sr_golomb_flac(PutBitContext *pb, int i, int k, int limit, int esc_len){
437
    int v;
438
439
    v = -2*i-1;
440
    v ^= (v>>31);
441
442
    set_ur_golomb_jpegls(pb, v, k, limit, esc_len);
443
}