]> git.sesse.net Git - ffmpeg/blob - libavcodec/msmpeg4dec.c
avcodec/siren: Factorize category 5 and 6 noise computation
[ffmpeg] / libavcodec / msmpeg4dec.c
1 /*
2  * MSMPEG4 backend for encoder and decoder
3  * Copyright (c) 2001 Fabrice Bellard
4  * Copyright (c) 2002-2013 Michael Niedermayer <michaelni@gmx.at>
5  *
6  * msmpeg4v1 & v2 stuff by Michael Niedermayer <michaelni@gmx.at>
7  *
8  * This file is part of FFmpeg.
9  *
10  * FFmpeg is free software; you can redistribute it and/or
11  * modify it under the terms of the GNU Lesser General Public
12  * License as published by the Free Software Foundation; either
13  * version 2.1 of the License, or (at your option) any later version.
14  *
15  * FFmpeg is distributed in the hope that it will be useful,
16  * but WITHOUT ANY WARRANTY; without even the implied warranty of
17  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
18  * Lesser General Public License for more details.
19  *
20  * You should have received a copy of the GNU Lesser General Public
21  * License along with FFmpeg; if not, write to the Free Software
22  * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
23  */
24
25 #include "avcodec.h"
26 #include "internal.h"
27 #include "mpegutils.h"
28 #include "mpegvideo.h"
29 #include "msmpeg4.h"
30 #include "libavutil/imgutils.h"
31 #include "h263.h"
32 #include "mpeg4video.h"
33 #include "msmpeg4data.h"
34 #include "vc1data.h"
35 #include "wmv2.h"
36
37 #define DC_VLC_BITS 9
38 #define V2_INTRA_CBPC_VLC_BITS 3
39 #define V2_MB_TYPE_VLC_BITS 7
40 #define MV_VLC_BITS 9
41 #define TEX_VLC_BITS 9
42
43 #define DEFAULT_INTER_INDEX 3
44
45 static inline int msmpeg4v1_pred_dc(MpegEncContext * s, int n,
46                                     int32_t **dc_val_ptr)
47 {
48     int i;
49
50     if (n < 4) {
51         i= 0;
52     } else {
53         i= n-3;
54     }
55
56     *dc_val_ptr= &s->last_dc[i];
57     return s->last_dc[i];
58 }
59
60 /****************************************/
61 /* decoding stuff */
62
63 VLC ff_mb_non_intra_vlc[4];
64 static VLC v2_dc_lum_vlc;
65 static VLC v2_dc_chroma_vlc;
66 static VLC v2_intra_cbpc_vlc;
67 static VLC v2_mb_type_vlc;
68 VLC ff_inter_intra_vlc;
69
70 /* This is identical to H.263 except that its range is multiplied by 2. */
71 static int msmpeg4v2_decode_motion(MpegEncContext * s, int pred, int f_code)
72 {
73     int code, val, sign, shift;
74
75     code = get_vlc2(&s->gb, ff_h263_mv_vlc.table, H263_MV_VLC_BITS, 2);
76     ff_dlog(s, "MV code %d at %d %d pred: %d\n", code, s->mb_x,s->mb_y, pred);
77     if (code < 0)
78         return 0xffff;
79
80     if (code == 0)
81         return pred;
82     sign = get_bits1(&s->gb);
83     shift = f_code - 1;
84     val = code;
85     if (shift) {
86         val = (val - 1) << shift;
87         val |= get_bits(&s->gb, shift);
88         val++;
89     }
90     if (sign)
91         val = -val;
92
93     val += pred;
94     if (val <= -64)
95         val += 64;
96     else if (val >= 64)
97         val -= 64;
98
99     return val;
100 }
101
102 static int msmpeg4v12_decode_mb(MpegEncContext *s, int16_t block[6][64])
103 {
104     int cbp, code, i;
105     uint32_t * const mb_type_ptr = &s->current_picture.mb_type[s->mb_x + s->mb_y*s->mb_stride];
106
107     if (s->pict_type == AV_PICTURE_TYPE_P) {
108         if (s->use_skip_mb_code) {
109             if (get_bits1(&s->gb)) {
110                 /* skip mb */
111                 s->mb_intra = 0;
112                 for(i=0;i<6;i++)
113                     s->block_last_index[i] = -1;
114                 s->mv_dir = MV_DIR_FORWARD;
115                 s->mv_type = MV_TYPE_16X16;
116                 s->mv[0][0][0] = 0;
117                 s->mv[0][0][1] = 0;
118                 s->mb_skipped = 1;
119                 *mb_type_ptr = MB_TYPE_SKIP | MB_TYPE_L0 | MB_TYPE_16x16;
120                 return 0;
121             }
122         }
123
124         if(s->msmpeg4_version==2)
125             code = get_vlc2(&s->gb, v2_mb_type_vlc.table, V2_MB_TYPE_VLC_BITS, 1);
126         else
127             code = get_vlc2(&s->gb, ff_h263_inter_MCBPC_vlc.table, INTER_MCBPC_VLC_BITS, 2);
128         if(code<0 || code>7){
129             av_log(s->avctx, AV_LOG_ERROR, "cbpc %d invalid at %d %d\n", code, s->mb_x, s->mb_y);
130             return -1;
131         }
132
133         s->mb_intra = code >>2;
134
135         cbp = code & 0x3;
136     } else {
137         s->mb_intra = 1;
138         if(s->msmpeg4_version==2)
139             cbp= get_vlc2(&s->gb, v2_intra_cbpc_vlc.table, V2_INTRA_CBPC_VLC_BITS, 1);
140         else
141             cbp= get_vlc2(&s->gb, ff_h263_intra_MCBPC_vlc.table, INTRA_MCBPC_VLC_BITS, 2);
142         if(cbp<0 || cbp>3){
143             av_log(s->avctx, AV_LOG_ERROR, "cbpc %d invalid at %d %d\n", cbp, s->mb_x, s->mb_y);
144             return -1;
145         }
146     }
147
148     if (!s->mb_intra) {
149         int mx, my, cbpy;
150
151         cbpy= get_vlc2(&s->gb, ff_h263_cbpy_vlc.table, CBPY_VLC_BITS, 1);
152         if(cbpy<0){
153             av_log(s->avctx, AV_LOG_ERROR, "cbpy %d invalid at %d %d\n", cbp, s->mb_x, s->mb_y);
154             return -1;
155         }
156
157         cbp|= cbpy<<2;
158         if(s->msmpeg4_version==1 || (cbp&3) != 3) cbp^= 0x3C;
159
160         ff_h263_pred_motion(s, 0, 0, &mx, &my);
161         mx= msmpeg4v2_decode_motion(s, mx, 1);
162         my= msmpeg4v2_decode_motion(s, my, 1);
163
164         s->mv_dir = MV_DIR_FORWARD;
165         s->mv_type = MV_TYPE_16X16;
166         s->mv[0][0][0] = mx;
167         s->mv[0][0][1] = my;
168         *mb_type_ptr = MB_TYPE_L0 | MB_TYPE_16x16;
169     } else {
170         int v;
171         if(s->msmpeg4_version==2){
172             s->ac_pred = get_bits1(&s->gb);
173             v = get_vlc2(&s->gb, ff_h263_cbpy_vlc.table, CBPY_VLC_BITS, 1);
174             if (v < 0) {
175                 av_log(s->avctx, AV_LOG_ERROR, "cbpy vlc invalid\n");
176                 return -1;
177             }
178             cbp|= v<<2;
179         } else{
180             s->ac_pred = 0;
181             v = get_vlc2(&s->gb, ff_h263_cbpy_vlc.table, CBPY_VLC_BITS, 1);
182             if (v < 0) {
183                 av_log(s->avctx, AV_LOG_ERROR, "cbpy vlc invalid\n");
184                 return -1;
185             }
186             cbp|= v<<2;
187             if(s->pict_type==AV_PICTURE_TYPE_P) cbp^=0x3C;
188         }
189         *mb_type_ptr = MB_TYPE_INTRA;
190     }
191
192     s->bdsp.clear_blocks(s->block[0]);
193     for (i = 0; i < 6; i++) {
194         if (ff_msmpeg4_decode_block(s, block[i], i, (cbp >> (5 - i)) & 1, NULL) < 0)
195         {
196              av_log(s->avctx, AV_LOG_ERROR, "\nerror while decoding block: %d x %d (%d)\n", s->mb_x, s->mb_y, i);
197              return -1;
198         }
199     }
200     return 0;
201 }
202
203 static int msmpeg4v34_decode_mb(MpegEncContext *s, int16_t block[6][64])
204 {
205     int cbp, code, i;
206     uint8_t *coded_val;
207     uint32_t * const mb_type_ptr = &s->current_picture.mb_type[s->mb_x + s->mb_y*s->mb_stride];
208
209     if (get_bits_left(&s->gb) <= 0)
210         return AVERROR_INVALIDDATA;
211
212     if (s->pict_type == AV_PICTURE_TYPE_P) {
213         if (s->use_skip_mb_code) {
214             if (get_bits1(&s->gb)) {
215                 /* skip mb */
216                 s->mb_intra = 0;
217                 for(i=0;i<6;i++)
218                     s->block_last_index[i] = -1;
219                 s->mv_dir = MV_DIR_FORWARD;
220                 s->mv_type = MV_TYPE_16X16;
221                 s->mv[0][0][0] = 0;
222                 s->mv[0][0][1] = 0;
223                 s->mb_skipped = 1;
224                 *mb_type_ptr = MB_TYPE_SKIP | MB_TYPE_L0 | MB_TYPE_16x16;
225
226                 return 0;
227             }
228         }
229
230         code = get_vlc2(&s->gb, ff_mb_non_intra_vlc[DEFAULT_INTER_INDEX].table, MB_NON_INTRA_VLC_BITS, 3);
231         if (code < 0)
232             return -1;
233         //s->mb_intra = (code & 0x40) ? 0 : 1;
234         s->mb_intra = (~code & 0x40) >> 6;
235
236         cbp = code & 0x3f;
237     } else {
238         s->mb_intra = 1;
239         code = get_vlc2(&s->gb, ff_msmp4_mb_i_vlc.table, MB_INTRA_VLC_BITS, 2);
240         if (code < 0)
241             return -1;
242         /* predict coded block pattern */
243         cbp = 0;
244         for(i=0;i<6;i++) {
245             int val = ((code >> (5 - i)) & 1);
246             if (i < 4) {
247                 int pred = ff_msmpeg4_coded_block_pred(s, i, &coded_val);
248                 val = val ^ pred;
249                 *coded_val = val;
250             }
251             cbp |= val << (5 - i);
252         }
253     }
254
255     if (!s->mb_intra) {
256         int mx, my;
257         if(s->per_mb_rl_table && cbp){
258             s->rl_table_index = decode012(&s->gb);
259             s->rl_chroma_table_index = s->rl_table_index;
260         }
261         ff_h263_pred_motion(s, 0, 0, &mx, &my);
262         if (ff_msmpeg4_decode_motion(s, &mx, &my) < 0)
263             return -1;
264         s->mv_dir = MV_DIR_FORWARD;
265         s->mv_type = MV_TYPE_16X16;
266         s->mv[0][0][0] = mx;
267         s->mv[0][0][1] = my;
268         *mb_type_ptr = MB_TYPE_L0 | MB_TYPE_16x16;
269     } else {
270         ff_dlog(s, "I at %d %d %d %06X\n", s->mb_x, s->mb_y,
271                 ((cbp & 3) ? 1 : 0) +((cbp & 0x3C)? 2 : 0),
272                 show_bits(&s->gb, 24));
273         s->ac_pred = get_bits1(&s->gb);
274         *mb_type_ptr = MB_TYPE_INTRA;
275         if(s->inter_intra_pred){
276             s->h263_aic_dir= get_vlc2(&s->gb, ff_inter_intra_vlc.table, INTER_INTRA_VLC_BITS, 1);
277             ff_dlog(s, "%d%d %d %d/",
278                     s->ac_pred, s->h263_aic_dir, s->mb_x, s->mb_y);
279         }
280         if(s->per_mb_rl_table && cbp){
281             s->rl_table_index = decode012(&s->gb);
282             s->rl_chroma_table_index = s->rl_table_index;
283         }
284     }
285
286     s->bdsp.clear_blocks(s->block[0]);
287     for (i = 0; i < 6; i++) {
288         if (ff_msmpeg4_decode_block(s, block[i], i, (cbp >> (5 - i)) & 1, NULL) < 0)
289         {
290             av_log(s->avctx, AV_LOG_ERROR, "\nerror while decoding block: %d x %d (%d)\n", s->mb_x, s->mb_y, i);
291             return -1;
292         }
293     }
294
295     return 0;
296 }
297
298 /* init all vlc decoding tables */
299 av_cold int ff_msmpeg4_decode_init(AVCodecContext *avctx)
300 {
301     MpegEncContext *s = avctx->priv_data;
302     static volatile int done = 0;
303     int i, ret;
304     MVTable *mv;
305
306     if ((ret = av_image_check_size(avctx->width, avctx->height, 0, avctx)) < 0)
307         return ret;
308
309     if (ff_h263_decode_init(avctx) < 0)
310         return -1;
311
312     ff_msmpeg4_common_init(s);
313
314     if (!done) {
315         for(i=0;i<NB_RL_TABLES;i++) {
316             ff_rl_init(&ff_rl_table[i], ff_static_rl_table_store[i]);
317         }
318         INIT_FIRST_VLC_RL(ff_rl_table[0], 642);
319         INIT_FIRST_VLC_RL(ff_rl_table[1], 1104);
320         INIT_FIRST_VLC_RL(ff_rl_table[2], 554);
321         INIT_VLC_RL(ff_rl_table[3], 940);
322         INIT_VLC_RL(ff_rl_table[4], 962);
323         /* ff_rl_table[5] coincides with ff_h263_rl_inter which has just been
324          * initialized in ff_h263_decode_init() above. So just copy the VLCs. */
325         av_assert1(ff_h263_rl_inter.rl_vlc[0]);
326         memcpy(ff_rl_table[5].rl_vlc, ff_h263_rl_inter.rl_vlc, sizeof(ff_rl_table[5].rl_vlc));
327
328         mv = &ff_mv_tables[0];
329         INIT_VLC_STATIC(&mv->vlc, MV_VLC_BITS, mv->n + 1,
330                     mv->table_mv_bits, 1, 1,
331                     mv->table_mv_code, 2, 2, 3714);
332         mv = &ff_mv_tables[1];
333         INIT_VLC_STATIC(&mv->vlc, MV_VLC_BITS, mv->n + 1,
334                     mv->table_mv_bits, 1, 1,
335                     mv->table_mv_code, 2, 2, 2694);
336
337         INIT_VLC_STATIC(&ff_msmp4_dc_luma_vlc[0], DC_VLC_BITS, 120,
338                  &ff_table0_dc_lum[0][1], 8, 4,
339                  &ff_table0_dc_lum[0][0], 8, 4, 1158);
340         INIT_VLC_STATIC(&ff_msmp4_dc_chroma_vlc[0], DC_VLC_BITS, 120,
341                  &ff_table0_dc_chroma[0][1], 8, 4,
342                  &ff_table0_dc_chroma[0][0], 8, 4, 1118);
343         INIT_VLC_STATIC(&ff_msmp4_dc_luma_vlc[1], DC_VLC_BITS, 120,
344                  &ff_table1_dc_lum[0][1], 8, 4,
345                  &ff_table1_dc_lum[0][0], 8, 4, 1476);
346         INIT_VLC_STATIC(&ff_msmp4_dc_chroma_vlc[1], DC_VLC_BITS, 120,
347                  &ff_table1_dc_chroma[0][1], 8, 4,
348                  &ff_table1_dc_chroma[0][0], 8, 4, 1216);
349
350         INIT_VLC_STATIC(&v2_dc_lum_vlc, DC_VLC_BITS, 512,
351                  &ff_v2_dc_lum_table[0][1], 8, 4,
352                  &ff_v2_dc_lum_table[0][0], 8, 4, 1472);
353         INIT_VLC_STATIC(&v2_dc_chroma_vlc, DC_VLC_BITS, 512,
354                  &ff_v2_dc_chroma_table[0][1], 8, 4,
355                  &ff_v2_dc_chroma_table[0][0], 8, 4, 1506);
356
357         INIT_VLC_STATIC(&v2_intra_cbpc_vlc, V2_INTRA_CBPC_VLC_BITS, 4,
358                  &ff_v2_intra_cbpc[0][1], 2, 1,
359                  &ff_v2_intra_cbpc[0][0], 2, 1, 8);
360         INIT_VLC_STATIC(&v2_mb_type_vlc, V2_MB_TYPE_VLC_BITS, 8,
361                  &ff_v2_mb_type[0][1], 2, 1,
362                  &ff_v2_mb_type[0][0], 2, 1, 128);
363
364         INIT_VLC_STATIC(&ff_mb_non_intra_vlc[0], MB_NON_INTRA_VLC_BITS, 128,
365                      &ff_wmv2_inter_table[0][0][1], 8, 4,
366                      &ff_wmv2_inter_table[0][0][0], 8, 4, 1636);
367         INIT_VLC_STATIC(&ff_mb_non_intra_vlc[1], MB_NON_INTRA_VLC_BITS, 128,
368                      &ff_wmv2_inter_table[1][0][1], 8, 4,
369                      &ff_wmv2_inter_table[1][0][0], 8, 4, 2648);
370         INIT_VLC_STATIC(&ff_mb_non_intra_vlc[2], MB_NON_INTRA_VLC_BITS, 128,
371                      &ff_wmv2_inter_table[2][0][1], 8, 4,
372                      &ff_wmv2_inter_table[2][0][0], 8, 4, 1532);
373         INIT_VLC_STATIC(&ff_mb_non_intra_vlc[3], MB_NON_INTRA_VLC_BITS, 128,
374                      &ff_wmv2_inter_table[3][0][1], 8, 4,
375                      &ff_wmv2_inter_table[3][0][0], 8, 4, 2488);
376
377         INIT_VLC_STATIC(&ff_msmp4_mb_i_vlc, MB_INTRA_VLC_BITS, 64,
378                  &ff_msmp4_mb_i_table[0][1], 4, 2,
379                  &ff_msmp4_mb_i_table[0][0], 4, 2, 536);
380
381         INIT_VLC_STATIC(&ff_inter_intra_vlc, INTER_INTRA_VLC_BITS, 4,
382                  &ff_table_inter_intra[0][1], 2, 1,
383                  &ff_table_inter_intra[0][0], 2, 1, 8);
384         done = 1;
385     }
386
387     switch(s->msmpeg4_version){
388     case 1:
389     case 2:
390         s->decode_mb= msmpeg4v12_decode_mb;
391         break;
392     case 3:
393     case 4:
394         s->decode_mb= msmpeg4v34_decode_mb;
395         break;
396     case 5:
397         if (CONFIG_WMV2_DECODER)
398             s->decode_mb= ff_wmv2_decode_mb;
399     case 6:
400         //FIXME + TODO VC1 decode mb
401         break;
402     }
403
404     s->slice_height= s->mb_height; //to avoid 1/0 if the first frame is not a keyframe
405
406     return 0;
407 }
408
409 int ff_msmpeg4_decode_picture_header(MpegEncContext * s)
410 {
411     int code;
412
413     // at minimum one bit per macroblock is required at least in a valid frame,
414     // we discard frames much smaller than this. Frames smaller than 1/8 of the
415     // smallest "black/skip" frame generally contain not much recoverable content
416     // while at the same time they have the highest computational requirements
417     // per byte
418     if (get_bits_left(&s->gb) * 8LL < (s->width+15)/16 * ((s->height+15)/16))
419         return AVERROR_INVALIDDATA;
420
421     if(s->msmpeg4_version==1){
422         int start_code = get_bits_long(&s->gb, 32);
423         if(start_code!=0x00000100){
424             av_log(s->avctx, AV_LOG_ERROR, "invalid startcode\n");
425             return -1;
426         }
427
428         skip_bits(&s->gb, 5); // frame number */
429     }
430
431     s->pict_type = get_bits(&s->gb, 2) + 1;
432     if (s->pict_type != AV_PICTURE_TYPE_I &&
433         s->pict_type != AV_PICTURE_TYPE_P){
434         av_log(s->avctx, AV_LOG_ERROR, "invalid picture type\n");
435         return -1;
436     }
437     s->chroma_qscale= s->qscale = get_bits(&s->gb, 5);
438     if(s->qscale==0){
439         av_log(s->avctx, AV_LOG_ERROR, "invalid qscale\n");
440         return -1;
441     }
442
443     if (s->pict_type == AV_PICTURE_TYPE_I) {
444         code = get_bits(&s->gb, 5);
445         if(s->msmpeg4_version==1){
446             if(code==0 || code>s->mb_height){
447                 av_log(s->avctx, AV_LOG_ERROR, "invalid slice height %d\n", code);
448                 return -1;
449             }
450
451             s->slice_height = code;
452         }else{
453             /* 0x17: one slice, 0x18: two slices, ... */
454             if (code < 0x17){
455                 av_log(s->avctx, AV_LOG_ERROR, "error, slice code was %X\n", code);
456                 return -1;
457             }
458
459             s->slice_height = s->mb_height / (code - 0x16);
460         }
461
462         switch(s->msmpeg4_version){
463         case 1:
464         case 2:
465             s->rl_chroma_table_index = 2;
466             s->rl_table_index = 2;
467
468             s->dc_table_index = 0; //not used
469             break;
470         case 3:
471             s->rl_chroma_table_index = decode012(&s->gb);
472             s->rl_table_index = decode012(&s->gb);
473
474             s->dc_table_index = get_bits1(&s->gb);
475             break;
476         case 4:
477             ff_msmpeg4_decode_ext_header(s, (2+5+5+17+7)/8);
478
479             if(s->bit_rate > MBAC_BITRATE) s->per_mb_rl_table= get_bits1(&s->gb);
480             else                           s->per_mb_rl_table= 0;
481
482             if(!s->per_mb_rl_table){
483                 s->rl_chroma_table_index = decode012(&s->gb);
484                 s->rl_table_index = decode012(&s->gb);
485             }
486
487             s->dc_table_index = get_bits1(&s->gb);
488             s->inter_intra_pred= 0;
489             break;
490         }
491         s->no_rounding = 1;
492         if(s->avctx->debug&FF_DEBUG_PICT_INFO)
493             av_log(s->avctx, AV_LOG_DEBUG, "qscale:%d rlc:%d rl:%d dc:%d mbrl:%d slice:%d   \n",
494                 s->qscale,
495                 s->rl_chroma_table_index,
496                 s->rl_table_index,
497                 s->dc_table_index,
498                 s->per_mb_rl_table,
499                 s->slice_height);
500     } else {
501         switch(s->msmpeg4_version){
502         case 1:
503         case 2:
504             if(s->msmpeg4_version==1)
505                 s->use_skip_mb_code = 1;
506             else
507                 s->use_skip_mb_code = get_bits1(&s->gb);
508             s->rl_table_index = 2;
509             s->rl_chroma_table_index = s->rl_table_index;
510             s->dc_table_index = 0; //not used
511             s->mv_table_index = 0;
512             break;
513         case 3:
514             s->use_skip_mb_code = get_bits1(&s->gb);
515             s->rl_table_index = decode012(&s->gb);
516             s->rl_chroma_table_index = s->rl_table_index;
517
518             s->dc_table_index = get_bits1(&s->gb);
519
520             s->mv_table_index = get_bits1(&s->gb);
521             break;
522         case 4:
523             s->use_skip_mb_code = get_bits1(&s->gb);
524
525             if(s->bit_rate > MBAC_BITRATE) s->per_mb_rl_table= get_bits1(&s->gb);
526             else                           s->per_mb_rl_table= 0;
527
528             if(!s->per_mb_rl_table){
529                 s->rl_table_index = decode012(&s->gb);
530                 s->rl_chroma_table_index = s->rl_table_index;
531             }
532
533             s->dc_table_index = get_bits1(&s->gb);
534
535             s->mv_table_index = get_bits1(&s->gb);
536             s->inter_intra_pred= (s->width*s->height < 320*240 && s->bit_rate<=II_BITRATE);
537             break;
538         }
539
540         if(s->avctx->debug&FF_DEBUG_PICT_INFO)
541             av_log(s->avctx, AV_LOG_DEBUG, "skip:%d rl:%d rlc:%d dc:%d mv:%d mbrl:%d qp:%d   \n",
542                 s->use_skip_mb_code,
543                 s->rl_table_index,
544                 s->rl_chroma_table_index,
545                 s->dc_table_index,
546                 s->mv_table_index,
547                 s->per_mb_rl_table,
548                 s->qscale);
549
550         if(s->flipflop_rounding){
551             s->no_rounding ^= 1;
552         }else{
553             s->no_rounding = 0;
554         }
555     }
556     ff_dlog(s->avctx, "%d %"PRId64" %d %d %d\n", s->pict_type, s->bit_rate,
557             s->inter_intra_pred, s->width, s->height);
558
559     s->esc3_level_length= 0;
560     s->esc3_run_length= 0;
561
562     return 0;
563 }
564
565 int ff_msmpeg4_decode_ext_header(MpegEncContext * s, int buf_size)
566 {
567     int left= buf_size*8 - get_bits_count(&s->gb);
568     int length= s->msmpeg4_version>=3 ? 17 : 16;
569     /* the alt_bitstream reader could read over the end so we need to check it */
570     if(left>=length && left<length+8)
571     {
572         skip_bits(&s->gb, 5); /* fps */
573         s->bit_rate= get_bits(&s->gb, 11)*1024;
574         if(s->msmpeg4_version>=3)
575             s->flipflop_rounding= get_bits1(&s->gb);
576         else
577             s->flipflop_rounding= 0;
578     }
579     else if(left<length+8)
580     {
581         s->flipflop_rounding= 0;
582         if(s->msmpeg4_version != 2)
583             av_log(s->avctx, AV_LOG_ERROR, "ext header missing, %d left\n", left);
584     }
585     else
586     {
587         av_log(s->avctx, AV_LOG_ERROR, "I-frame too long, ignoring ext header\n");
588     }
589
590     return 0;
591 }
592
593 static int msmpeg4_decode_dc(MpegEncContext * s, int n, int *dir_ptr)
594 {
595     int level, pred;
596
597     if(s->msmpeg4_version<=2){
598         if (n < 4) {
599             level = get_vlc2(&s->gb, v2_dc_lum_vlc.table, DC_VLC_BITS, 3);
600         } else {
601             level = get_vlc2(&s->gb, v2_dc_chroma_vlc.table, DC_VLC_BITS, 3);
602         }
603         if (level < 0) {
604             av_log(s->avctx, AV_LOG_ERROR, "illegal dc vlc\n");
605             *dir_ptr = 0;
606             return -1;
607         }
608         level-=256;
609     }else{  //FIXME optimize use unified tables & index
610         if (n < 4) {
611             level = get_vlc2(&s->gb, ff_msmp4_dc_luma_vlc[s->dc_table_index].table, DC_VLC_BITS, 3);
612         } else {
613             level = get_vlc2(&s->gb, ff_msmp4_dc_chroma_vlc[s->dc_table_index].table, DC_VLC_BITS, 3);
614         }
615         if (level < 0){
616             av_log(s->avctx, AV_LOG_ERROR, "illegal dc vlc\n");
617             *dir_ptr = 0;
618             return -1;
619         }
620
621         if (level == DC_MAX) {
622             level = get_bits(&s->gb, 8);
623             if (get_bits1(&s->gb))
624                 level = -level;
625         } else if (level != 0) {
626             if (get_bits1(&s->gb))
627                 level = -level;
628         }
629     }
630
631     if(s->msmpeg4_version==1){
632         int32_t *dc_val;
633         pred = msmpeg4v1_pred_dc(s, n, &dc_val);
634         level += pred;
635
636         /* update predictor */
637         *dc_val= level;
638     }else{
639         int16_t *dc_val;
640         pred   = ff_msmpeg4_pred_dc(s, n, &dc_val, dir_ptr);
641         level += pred;
642
643         /* update predictor */
644         if (n < 4) {
645             *dc_val = level * s->y_dc_scale;
646         } else {
647             *dc_val = level * s->c_dc_scale;
648         }
649     }
650
651     return level;
652 }
653
654 int ff_msmpeg4_decode_block(MpegEncContext * s, int16_t * block,
655                               int n, int coded, const uint8_t *scan_table)
656 {
657     int level, i, last, run, run_diff;
658     int av_uninit(dc_pred_dir);
659     RLTable *rl;
660     RL_VLC_ELEM *rl_vlc;
661     int qmul, qadd;
662
663     if (s->mb_intra) {
664         qmul=1;
665         qadd=0;
666
667         /* DC coef */
668         level = msmpeg4_decode_dc(s, n, &dc_pred_dir);
669
670         if (level < 0){
671             av_log(s->avctx, AV_LOG_ERROR, "dc overflow- block: %d qscale: %d//\n", n, s->qscale);
672             if(s->inter_intra_pred) level=0;
673         }
674         if (n < 4) {
675             rl = &ff_rl_table[s->rl_table_index];
676             if(level > 256*s->y_dc_scale){
677                 av_log(s->avctx, AV_LOG_ERROR, "dc overflow+ L qscale: %d//\n", s->qscale);
678                 if(!s->inter_intra_pred) return -1;
679             }
680         } else {
681             rl = &ff_rl_table[3 + s->rl_chroma_table_index];
682             if(level > 256*s->c_dc_scale){
683                 av_log(s->avctx, AV_LOG_ERROR, "dc overflow+ C qscale: %d//\n", s->qscale);
684                 if(!s->inter_intra_pred) return -1;
685             }
686         }
687         block[0] = level;
688
689         run_diff = s->msmpeg4_version >= 4;
690         i = 0;
691         if (!coded) {
692             goto not_coded;
693         }
694         if (s->ac_pred) {
695             if (dc_pred_dir == 0)
696                 scan_table = s->intra_v_scantable.permutated; /* left */
697             else
698                 scan_table = s->intra_h_scantable.permutated; /* top */
699         } else {
700             scan_table = s->intra_scantable.permutated;
701         }
702         rl_vlc= rl->rl_vlc[0];
703     } else {
704         qmul = s->qscale << 1;
705         qadd = (s->qscale - 1) | 1;
706         i = -1;
707         rl = &ff_rl_table[3 + s->rl_table_index];
708
709         if(s->msmpeg4_version==2)
710             run_diff = 0;
711         else
712             run_diff = 1;
713
714         if (!coded) {
715             s->block_last_index[n] = i;
716             return 0;
717         }
718         if(!scan_table)
719             scan_table = s->inter_scantable.permutated;
720         rl_vlc= rl->rl_vlc[s->qscale];
721     }
722   {
723     OPEN_READER(re, &s->gb);
724     for(;;) {
725         UPDATE_CACHE(re, &s->gb);
726         GET_RL_VLC(level, run, re, &s->gb, rl_vlc, TEX_VLC_BITS, 2, 0);
727         if (level==0) {
728             int cache;
729             cache= GET_CACHE(re, &s->gb);
730             /* escape */
731             if (s->msmpeg4_version==1 || (cache&0x80000000)==0) {
732                 if (s->msmpeg4_version==1 || (cache&0x40000000)==0) {
733                     /* third escape */
734                     if(s->msmpeg4_version!=1) LAST_SKIP_BITS(re, &s->gb, 2);
735                     UPDATE_CACHE(re, &s->gb);
736                     if(s->msmpeg4_version<=3){
737                         last=  SHOW_UBITS(re, &s->gb, 1); SKIP_CACHE(re, &s->gb, 1);
738                         run=   SHOW_UBITS(re, &s->gb, 6); SKIP_CACHE(re, &s->gb, 6);
739                         level= SHOW_SBITS(re, &s->gb, 8);
740                         SKIP_COUNTER(re, &s->gb, 1+6+8);
741                     }else{
742                         int sign;
743                         last=  SHOW_UBITS(re, &s->gb, 1); SKIP_BITS(re, &s->gb, 1);
744                         if(!s->esc3_level_length){
745                             int ll;
746                             ff_dlog(s->avctx, "ESC-3 %X at %d %d\n",
747                                     show_bits(&s->gb, 24), s->mb_x, s->mb_y);
748                             if(s->qscale<8){
749                                 ll= SHOW_UBITS(re, &s->gb, 3); SKIP_BITS(re, &s->gb, 3);
750                                 if(ll==0){
751                                     ll= 8+SHOW_UBITS(re, &s->gb, 1); SKIP_BITS(re, &s->gb, 1);
752                                 }
753                             }else{
754                                 ll=2;
755                                 while(ll<8 && SHOW_UBITS(re, &s->gb, 1)==0){
756                                     ll++;
757                                     SKIP_BITS(re, &s->gb, 1);
758                                 }
759                                 if(ll<8) SKIP_BITS(re, &s->gb, 1);
760                             }
761
762                             s->esc3_level_length= ll;
763                             s->esc3_run_length= SHOW_UBITS(re, &s->gb, 2) + 3; SKIP_BITS(re, &s->gb, 2);
764                             UPDATE_CACHE(re, &s->gb);
765                         }
766                         run=   SHOW_UBITS(re, &s->gb, s->esc3_run_length);
767                         SKIP_BITS(re, &s->gb, s->esc3_run_length);
768
769                         sign=  SHOW_UBITS(re, &s->gb, 1);
770                         SKIP_BITS(re, &s->gb, 1);
771
772                         level= SHOW_UBITS(re, &s->gb, s->esc3_level_length);
773                         SKIP_BITS(re, &s->gb, s->esc3_level_length);
774                         if(sign) level= -level;
775                     }
776
777                     //level = level * qmul + (level>0) * qadd - (level<=0) * qadd ;
778                     if (level>0) level= level * qmul + qadd;
779                     else         level= level * qmul - qadd;
780                     i+= run + 1;
781                     if(last) i+=192;
782                 } else {
783                     /* second escape */
784                     SKIP_BITS(re, &s->gb, 2);
785                     GET_RL_VLC(level, run, re, &s->gb, rl_vlc, TEX_VLC_BITS, 2, 1);
786                     i+= run + rl->max_run[run>>7][level/qmul] + run_diff; //FIXME opt indexing
787                     level = (level ^ SHOW_SBITS(re, &s->gb, 1)) - SHOW_SBITS(re, &s->gb, 1);
788                     LAST_SKIP_BITS(re, &s->gb, 1);
789                 }
790             } else {
791                 /* first escape */
792                 SKIP_BITS(re, &s->gb, 1);
793                 GET_RL_VLC(level, run, re, &s->gb, rl_vlc, TEX_VLC_BITS, 2, 1);
794                 i+= run;
795                 level = level + rl->max_level[run>>7][(run-1)&63] * qmul;//FIXME opt indexing
796                 level = (level ^ SHOW_SBITS(re, &s->gb, 1)) - SHOW_SBITS(re, &s->gb, 1);
797                 LAST_SKIP_BITS(re, &s->gb, 1);
798             }
799         } else {
800             i+= run;
801             level = (level ^ SHOW_SBITS(re, &s->gb, 1)) - SHOW_SBITS(re, &s->gb, 1);
802             LAST_SKIP_BITS(re, &s->gb, 1);
803         }
804         if (i > 62){
805             i-= 192;
806             if(i&(~63)){
807                 const int left= get_bits_left(&s->gb);
808                 if (((i + 192 == 64 && level / qmul == -1) ||
809                      !(s->avctx->err_recognition & (AV_EF_BITSTREAM|AV_EF_COMPLIANT))) &&
810                     left >= 0) {
811                     av_log(s->avctx, AV_LOG_ERROR, "ignoring overflow at %d %d\n", s->mb_x, s->mb_y);
812                     i = 63;
813                     break;
814                 }else{
815                     av_log(s->avctx, AV_LOG_ERROR, "ac-tex damaged at %d %d\n", s->mb_x, s->mb_y);
816                     return -1;
817                 }
818             }
819
820             block[scan_table[i]] = level;
821             break;
822         }
823
824         block[scan_table[i]] = level;
825     }
826     CLOSE_READER(re, &s->gb);
827   }
828  not_coded:
829     if (s->mb_intra) {
830         ff_mpeg4_pred_ac(s, block, n, dc_pred_dir);
831         if (s->ac_pred) {
832             i = 63; /* XXX: not optimal */
833         }
834     }
835     if(s->msmpeg4_version>=4 && i>0) i=63; //FIXME/XXX optimize
836     s->block_last_index[n] = i;
837
838     return 0;
839 }
840
841 int ff_msmpeg4_decode_motion(MpegEncContext * s,
842                                  int *mx_ptr, int *my_ptr)
843 {
844     MVTable *mv;
845     int code, mx, my;
846
847     mv = &ff_mv_tables[s->mv_table_index];
848
849     code = get_vlc2(&s->gb, mv->vlc.table, MV_VLC_BITS, 2);
850     if (code < 0){
851         av_log(s->avctx, AV_LOG_ERROR, "illegal MV code at %d %d\n", s->mb_x, s->mb_y);
852         return -1;
853     }
854     if (code == mv->n) {
855         mx = get_bits(&s->gb, 6);
856         my = get_bits(&s->gb, 6);
857     } else {
858         mx = mv->table_mvx[code];
859         my = mv->table_mvy[code];
860     }
861
862     mx += *mx_ptr - 32;
863     my += *my_ptr - 32;
864     /* WARNING : they do not do exactly modulo encoding */
865     if (mx <= -64)
866         mx += 64;
867     else if (mx >= 64)
868         mx -= 64;
869
870     if (my <= -64)
871         my += 64;
872     else if (my >= 64)
873         my -= 64;
874     *mx_ptr = mx;
875     *my_ptr = my;
876     return 0;
877 }
878
879 AVCodec ff_msmpeg4v1_decoder = {
880     .name           = "msmpeg4v1",
881     .long_name      = NULL_IF_CONFIG_SMALL("MPEG-4 part 2 Microsoft variant version 1"),
882     .type           = AVMEDIA_TYPE_VIDEO,
883     .id             = AV_CODEC_ID_MSMPEG4V1,
884     .priv_data_size = sizeof(MpegEncContext),
885     .init           = ff_msmpeg4_decode_init,
886     .close          = ff_h263_decode_end,
887     .decode         = ff_h263_decode_frame,
888     .capabilities   = AV_CODEC_CAP_DRAW_HORIZ_BAND | AV_CODEC_CAP_DR1,
889     .caps_internal  = FF_CODEC_CAP_SKIP_FRAME_FILL_PARAM | FF_CODEC_CAP_INIT_CLEANUP,
890     .max_lowres     = 3,
891     .pix_fmts       = (const enum AVPixelFormat[]) {
892         AV_PIX_FMT_YUV420P,
893         AV_PIX_FMT_NONE
894     },
895 };
896
897 AVCodec ff_msmpeg4v2_decoder = {
898     .name           = "msmpeg4v2",
899     .long_name      = NULL_IF_CONFIG_SMALL("MPEG-4 part 2 Microsoft variant version 2"),
900     .type           = AVMEDIA_TYPE_VIDEO,
901     .id             = AV_CODEC_ID_MSMPEG4V2,
902     .priv_data_size = sizeof(MpegEncContext),
903     .init           = ff_msmpeg4_decode_init,
904     .close          = ff_h263_decode_end,
905     .decode         = ff_h263_decode_frame,
906     .capabilities   = AV_CODEC_CAP_DRAW_HORIZ_BAND | AV_CODEC_CAP_DR1,
907     .caps_internal  = FF_CODEC_CAP_SKIP_FRAME_FILL_PARAM | FF_CODEC_CAP_INIT_CLEANUP,
908     .max_lowres     = 3,
909     .pix_fmts       = (const enum AVPixelFormat[]) {
910         AV_PIX_FMT_YUV420P,
911         AV_PIX_FMT_NONE
912     },
913 };
914
915 AVCodec ff_msmpeg4v3_decoder = {
916     .name           = "msmpeg4",
917     .long_name      = NULL_IF_CONFIG_SMALL("MPEG-4 part 2 Microsoft variant version 3"),
918     .type           = AVMEDIA_TYPE_VIDEO,
919     .id             = AV_CODEC_ID_MSMPEG4V3,
920     .priv_data_size = sizeof(MpegEncContext),
921     .init           = ff_msmpeg4_decode_init,
922     .close          = ff_h263_decode_end,
923     .decode         = ff_h263_decode_frame,
924     .capabilities   = AV_CODEC_CAP_DRAW_HORIZ_BAND | AV_CODEC_CAP_DR1,
925     .caps_internal  = FF_CODEC_CAP_SKIP_FRAME_FILL_PARAM | FF_CODEC_CAP_INIT_CLEANUP,
926     .max_lowres     = 3,
927     .pix_fmts       = (const enum AVPixelFormat[]) {
928         AV_PIX_FMT_YUV420P,
929         AV_PIX_FMT_NONE
930     },
931 };
932
933 AVCodec ff_wmv1_decoder = {
934     .name           = "wmv1",
935     .long_name      = NULL_IF_CONFIG_SMALL("Windows Media Video 7"),
936     .type           = AVMEDIA_TYPE_VIDEO,
937     .id             = AV_CODEC_ID_WMV1,
938     .priv_data_size = sizeof(MpegEncContext),
939     .init           = ff_msmpeg4_decode_init,
940     .close          = ff_h263_decode_end,
941     .decode         = ff_h263_decode_frame,
942     .capabilities   = AV_CODEC_CAP_DRAW_HORIZ_BAND | AV_CODEC_CAP_DR1,
943     .caps_internal  = FF_CODEC_CAP_SKIP_FRAME_FILL_PARAM | FF_CODEC_CAP_INIT_CLEANUP,
944     .max_lowres     = 3,
945     .pix_fmts       = (const enum AVPixelFormat[]) {
946         AV_PIX_FMT_YUV420P,
947         AV_PIX_FMT_NONE
948     },
949 };