6c4882671a77f40bee70559f7a840039cfa2fbf7
[ffmpeg.git] / libavcodec / mpeg12.c
1 /*
2  * MPEG-1/2 decoder
3  * Copyright (c) 2000, 2001 Fabrice Bellard
4  * Copyright (c) 2002-2004 Michael Niedermayer <michaelni@gmx.at>
5  *
6  * This file is part of FFmpeg.
7  *
8  * FFmpeg is free software; you can redistribute it and/or
9  * modify it under the terms of the GNU Lesser General Public
10  * License as published by the Free Software Foundation; either
11  * version 2.1 of the License, or (at your option) any later version.
12  *
13  * FFmpeg is distributed in the hope that it will be useful,
14  * but WITHOUT ANY WARRANTY; without even the implied warranty of
15  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
16  * Lesser General Public License for more details.
17  *
18  * You should have received a copy of the GNU Lesser General Public
19  * License along with FFmpeg; if not, write to the Free Software
20  * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
21  */
22
23 /**
24  * @file
25  * MPEG-1/2 decoder
26  */
27
28 //#define DEBUG
29 #include "internal.h"
30 #include "avcodec.h"
31 #include "dsputil.h"
32 #include "mpegvideo.h"
33 #include "libavutil/avassert.h"
34
35 #include "mpeg12.h"
36 #include "mpeg12data.h"
37 #include "mpeg12decdata.h"
38 #include "bytestream.h"
39 #include "vdpau_internal.h"
40 #include "xvmc_internal.h"
41 #include "thread.h"
42
43 //#undef NDEBUG
44 //#include <assert.h>
45
46
47 #define MV_VLC_BITS 9
48 #define MBINCR_VLC_BITS 9
49 #define MB_PAT_VLC_BITS 9
50 #define MB_PTYPE_VLC_BITS 6
51 #define MB_BTYPE_VLC_BITS 6
52
53 static inline int mpeg1_decode_block_intra(MpegEncContext *s, DCTELEM *block, int n);
54 static inline int mpeg1_decode_block_inter(MpegEncContext *s, DCTELEM *block, int n);
55 static inline int mpeg1_fast_decode_block_inter(MpegEncContext *s, DCTELEM *block, int n);
56 static inline int mpeg2_decode_block_non_intra(MpegEncContext *s, DCTELEM *block, int n);
57 static inline int mpeg2_decode_block_intra(MpegEncContext *s, DCTELEM *block, int n);
58 static inline int mpeg2_fast_decode_block_non_intra(MpegEncContext *s, DCTELEM *block, int n);
59 static inline int mpeg2_fast_decode_block_intra(MpegEncContext *s, DCTELEM *block, int n);
60 static int mpeg_decode_motion(MpegEncContext *s, int fcode, int pred);
61 static void exchange_uv(MpegEncContext *s);
62
63 uint8_t ff_mpeg12_static_rl_table_store[2][2][2*MAX_RUN + MAX_LEVEL + 3];
64
65
66 #define INIT_2D_VLC_RL(rl, static_size)\
67 {\
68     static RL_VLC_ELEM rl_vlc_table[static_size];\
69     INIT_VLC_STATIC(&rl.vlc, TEX_VLC_BITS, rl.n + 2,\
70                     &rl.table_vlc[0][1], 4, 2,\
71                     &rl.table_vlc[0][0], 4, 2, static_size);\
72 \
73     rl.rl_vlc[0] = rl_vlc_table;\
74     init_2d_vlc_rl(&rl);\
75 }
76
77 static void init_2d_vlc_rl(RLTable *rl)
78 {
79     int i;
80
81     for (i = 0; i < rl->vlc.table_size; i++) {
82         int code = rl->vlc.table[i][0];
83         int len  = rl->vlc.table[i][1];
84         int level, run;
85
86         if (len == 0) { // illegal code
87             run   = 65;
88             level = MAX_LEVEL;
89         } else if (len<0) { //more bits needed
90             run   = 0;
91             level = code;
92         } else {
93             if (code == rl->n) { //esc
94                 run   = 65;
95                 level = 0;
96             } else if (code == rl->n+1) { //eob
97                 run   = 0;
98                 level = 127;
99             } else {
100                 run   = rl->table_run  [code] + 1;
101                 level = rl->table_level[code];
102             }
103         }
104         rl->rl_vlc[0][i].len   = len;
105         rl->rl_vlc[0][i].level = level;
106         rl->rl_vlc[0][i].run   = run;
107     }
108 }
109
110 void ff_mpeg12_common_init(MpegEncContext *s)
111 {
112
113     s->y_dc_scale_table =
114     s->c_dc_scale_table = ff_mpeg2_dc_scale_table[s->intra_dc_precision];
115
116 }
117
118 void ff_mpeg1_clean_buffers(MpegEncContext *s)
119 {
120     s->last_dc[0] = 1 << (7 + s->intra_dc_precision);
121     s->last_dc[1] = s->last_dc[0];
122     s->last_dc[2] = s->last_dc[0];
123     memset(s->last_mv, 0, sizeof(s->last_mv));
124 }
125
126
127 /******************************************/
128 /* decoding */
129
130 VLC ff_dc_lum_vlc;
131 VLC ff_dc_chroma_vlc;
132
133 static VLC mv_vlc;
134 static VLC mbincr_vlc;
135 static VLC mb_ptype_vlc;
136 static VLC mb_btype_vlc;
137 static VLC mb_pat_vlc;
138
139 av_cold void ff_mpeg12_init_vlcs(void)
140 {
141     static int done = 0;
142
143     if (!done) {
144         done = 1;
145
146         INIT_VLC_STATIC(&ff_dc_lum_vlc, DC_VLC_BITS, 12,
147                         ff_mpeg12_vlc_dc_lum_bits, 1, 1,
148                         ff_mpeg12_vlc_dc_lum_code, 2, 2, 512);
149         INIT_VLC_STATIC(&ff_dc_chroma_vlc,  DC_VLC_BITS, 12,
150                         ff_mpeg12_vlc_dc_chroma_bits, 1, 1,
151                         ff_mpeg12_vlc_dc_chroma_code, 2, 2, 514);
152         INIT_VLC_STATIC(&mv_vlc, MV_VLC_BITS, 17,
153                         &ff_mpeg12_mbMotionVectorTable[0][1], 2, 1,
154                         &ff_mpeg12_mbMotionVectorTable[0][0], 2, 1, 518);
155         INIT_VLC_STATIC(&mbincr_vlc, MBINCR_VLC_BITS, 36,
156                         &ff_mpeg12_mbAddrIncrTable[0][1], 2, 1,
157                         &ff_mpeg12_mbAddrIncrTable[0][0], 2, 1, 538);
158         INIT_VLC_STATIC(&mb_pat_vlc, MB_PAT_VLC_BITS, 64,
159                         &ff_mpeg12_mbPatTable[0][1], 2, 1,
160                         &ff_mpeg12_mbPatTable[0][0], 2, 1, 512);
161
162         INIT_VLC_STATIC(&mb_ptype_vlc, MB_PTYPE_VLC_BITS, 7,
163                         &table_mb_ptype[0][1], 2, 1,
164                         &table_mb_ptype[0][0], 2, 1, 64);
165         INIT_VLC_STATIC(&mb_btype_vlc, MB_BTYPE_VLC_BITS, 11,
166                         &table_mb_btype[0][1], 2, 1,
167                         &table_mb_btype[0][0], 2, 1, 64);
168         init_rl(&ff_rl_mpeg1, ff_mpeg12_static_rl_table_store[0]);
169         init_rl(&ff_rl_mpeg2, ff_mpeg12_static_rl_table_store[1]);
170
171         INIT_2D_VLC_RL(ff_rl_mpeg1, 680);
172         INIT_2D_VLC_RL(ff_rl_mpeg2, 674);
173     }
174 }
175
176 static inline int get_dmv(MpegEncContext *s)
177 {
178     if (get_bits1(&s->gb))
179         return 1 - (get_bits1(&s->gb) << 1);
180     else
181         return 0;
182 }
183
184 static inline int get_qscale(MpegEncContext *s)
185 {
186     int qscale = get_bits(&s->gb, 5);
187     if (s->q_scale_type) {
188         return non_linear_qscale[qscale];
189     } else {
190         return qscale << 1;
191     }
192 }
193
194 /* motion type (for MPEG-2) */
195 #define MT_FIELD 1
196 #define MT_FRAME 2
197 #define MT_16X8  2
198 #define MT_DMV   3
199
200 static int mpeg_decode_mb(MpegEncContext *s, DCTELEM block[12][64])
201 {
202     int i, j, k, cbp, val, mb_type, motion_type;
203     const int mb_block_count = 4 + (1 << s->chroma_format);
204
205     av_dlog(s->avctx, "decode_mb: x=%d y=%d\n", s->mb_x, s->mb_y);
206
207     assert(s->mb_skipped == 0);
208
209     if (s->mb_skip_run-- != 0) {
210         if (s->pict_type == AV_PICTURE_TYPE_P) {
211             s->mb_skipped = 1;
212             s->current_picture.f.mb_type[s->mb_x + s->mb_y * s->mb_stride] = MB_TYPE_SKIP | MB_TYPE_L0 | MB_TYPE_16x16;
213         } else {
214             int mb_type;
215
216             if (s->mb_x)
217                 mb_type = s->current_picture.f.mb_type[s->mb_x + s->mb_y * s->mb_stride - 1];
218             else
219                 mb_type = s->current_picture.f.mb_type[s->mb_width + (s->mb_y - 1) * s->mb_stride - 1]; // FIXME not sure if this is allowed in MPEG at all
220             if (IS_INTRA(mb_type))
221                 return -1;
222
223             s->current_picture.f.mb_type[s->mb_x + s->mb_y*s->mb_stride] =
224                 mb_type | MB_TYPE_SKIP;
225 //            assert(s->current_picture.f.mb_type[s->mb_x + s->mb_y * s->mb_stride - 1] & (MB_TYPE_16x16 | MB_TYPE_16x8));
226
227             if ((s->mv[0][0][0] | s->mv[0][0][1] | s->mv[1][0][0] | s->mv[1][0][1]) == 0)
228                 s->mb_skipped = 1;
229         }
230
231         return 0;
232     }
233
234     switch (s->pict_type) {
235     default:
236     case AV_PICTURE_TYPE_I:
237         if (get_bits1(&s->gb) == 0) {
238             if (get_bits1(&s->gb) == 0) {
239                 av_log(s->avctx, AV_LOG_ERROR, "invalid mb type in I Frame at %d %d\n", s->mb_x, s->mb_y);
240                 return -1;
241             }
242             mb_type = MB_TYPE_QUANT | MB_TYPE_INTRA;
243         } else {
244             mb_type = MB_TYPE_INTRA;
245         }
246         break;
247     case AV_PICTURE_TYPE_P:
248         mb_type = get_vlc2(&s->gb, mb_ptype_vlc.table, MB_PTYPE_VLC_BITS, 1);
249         if (mb_type < 0) {
250             av_log(s->avctx, AV_LOG_ERROR, "invalid mb type in P Frame at %d %d\n", s->mb_x, s->mb_y);
251             return -1;
252         }
253         mb_type = ptype2mb_type[mb_type];
254         break;
255     case AV_PICTURE_TYPE_B:
256         mb_type = get_vlc2(&s->gb, mb_btype_vlc.table, MB_BTYPE_VLC_BITS, 1);
257         if (mb_type < 0) {
258             av_log(s->avctx, AV_LOG_ERROR, "invalid mb type in B Frame at %d %d\n", s->mb_x, s->mb_y);
259             return -1;
260         }
261         mb_type = btype2mb_type[mb_type];
262         break;
263     }
264     av_dlog(s->avctx, "mb_type=%x\n", mb_type);
265 //    motion_type = 0; /* avoid warning */
266     if (IS_INTRA(mb_type)) {
267         s->dsp.clear_blocks(s->block[0]);
268
269         if (!s->chroma_y_shift) {
270             s->dsp.clear_blocks(s->block[6]);
271         }
272
273         /* compute DCT type */
274         if (s->picture_structure == PICT_FRAME && // FIXME add an interlaced_dct coded var?
275             !s->frame_pred_frame_dct) {
276             s->interlaced_dct = get_bits1(&s->gb);
277         }
278
279         if (IS_QUANT(mb_type))
280             s->qscale = get_qscale(s);
281
282         if (s->concealment_motion_vectors) {
283             /* just parse them */
284             if (s->picture_structure != PICT_FRAME)
285                 skip_bits1(&s->gb); /* field select */
286
287             s->mv[0][0][0]= s->last_mv[0][0][0]= s->last_mv[0][1][0] =
288                 mpeg_decode_motion(s, s->mpeg_f_code[0][0], s->last_mv[0][0][0]);
289             s->mv[0][0][1]= s->last_mv[0][0][1]= s->last_mv[0][1][1] =
290                 mpeg_decode_motion(s, s->mpeg_f_code[0][1], s->last_mv[0][0][1]);
291
292             skip_bits1(&s->gb); /* marker */
293         } else
294             memset(s->last_mv, 0, sizeof(s->last_mv)); /* reset mv prediction */
295         s->mb_intra = 1;
296         // if 1, we memcpy blocks in xvmcvideo
297         if (CONFIG_MPEG_XVMC_DECODER && s->avctx->xvmc_acceleration > 1) {
298             ff_xvmc_pack_pblocks(s, -1); // inter are always full blocks
299             if (s->swap_uv) {
300                 exchange_uv(s);
301             }
302         }
303
304         if (s->codec_id == CODEC_ID_MPEG2VIDEO) {
305             if (s->flags2 & CODEC_FLAG2_FAST) {
306                 for (i = 0; i < 6; i++) {
307                     mpeg2_fast_decode_block_intra(s, *s->pblocks[i], i);
308                 }
309             } else {
310                 for (i = 0; i < mb_block_count; i++) {
311                     if (mpeg2_decode_block_intra(s, *s->pblocks[i], i) < 0)
312                         return -1;
313                 }
314             }
315         } else {
316             for (i = 0; i < 6; i++) {
317                 if (mpeg1_decode_block_intra(s, *s->pblocks[i], i) < 0)
318                     return -1;
319             }
320         }
321     } else {
322         if (mb_type & MB_TYPE_ZERO_MV) {
323             assert(mb_type & MB_TYPE_CBP);
324
325             s->mv_dir = MV_DIR_FORWARD;
326             if (s->picture_structure == PICT_FRAME) {
327                 if (!s->frame_pred_frame_dct)
328                     s->interlaced_dct = get_bits1(&s->gb);
329                 s->mv_type = MV_TYPE_16X16;
330             } else {
331                 s->mv_type = MV_TYPE_FIELD;
332                 mb_type |= MB_TYPE_INTERLACED;
333                 s->field_select[0][0] = s->picture_structure - 1;
334             }
335
336             if (IS_QUANT(mb_type))
337                 s->qscale = get_qscale(s);
338
339             s->last_mv[0][0][0] = 0;
340             s->last_mv[0][0][1] = 0;
341             s->last_mv[0][1][0] = 0;
342             s->last_mv[0][1][1] = 0;
343             s->mv[0][0][0] = 0;
344             s->mv[0][0][1] = 0;
345         } else {
346             assert(mb_type & MB_TYPE_L0L1);
347             // FIXME decide if MBs in field pictures are MB_TYPE_INTERLACED
348             /* get additional motion vector type */
349             if (s->frame_pred_frame_dct)
350                 motion_type = MT_FRAME;
351             else {
352                 motion_type = get_bits(&s->gb, 2);
353                 if (s->picture_structure == PICT_FRAME && HAS_CBP(mb_type))
354                     s->interlaced_dct = get_bits1(&s->gb);
355             }
356
357             if (IS_QUANT(mb_type))
358                 s->qscale = get_qscale(s);
359
360             /* motion vectors */
361             s->mv_dir = (mb_type >> 13) & 3;
362             av_dlog(s->avctx, "motion_type=%d\n", motion_type);
363             switch (motion_type) {
364             case MT_FRAME: /* or MT_16X8 */
365                 if (s->picture_structure == PICT_FRAME) {
366                     mb_type |= MB_TYPE_16x16;
367                     s->mv_type = MV_TYPE_16X16;
368                     for (i = 0; i < 2; i++) {
369                         if (USES_LIST(mb_type, i)) {
370                             /* MT_FRAME */
371                             s->mv[i][0][0]= s->last_mv[i][0][0]= s->last_mv[i][1][0] =
372                                 mpeg_decode_motion(s, s->mpeg_f_code[i][0], s->last_mv[i][0][0]);
373                             s->mv[i][0][1]= s->last_mv[i][0][1]= s->last_mv[i][1][1] =
374                                 mpeg_decode_motion(s, s->mpeg_f_code[i][1], s->last_mv[i][0][1]);
375                             /* full_pel: only for MPEG-1 */
376                             if (s->full_pel[i]) {
377                                 s->mv[i][0][0] <<= 1;
378                                 s->mv[i][0][1] <<= 1;
379                             }
380                         }
381                     }
382                 } else {
383                     mb_type |= MB_TYPE_16x8 | MB_TYPE_INTERLACED;
384                     s->mv_type = MV_TYPE_16X8;
385                     for (i = 0; i < 2; i++) {
386                         if (USES_LIST(mb_type, i)) {
387                             /* MT_16X8 */
388                             for (j = 0; j < 2; j++) {
389                                 s->field_select[i][j] = get_bits1(&s->gb);
390                                 for (k = 0; k < 2; k++) {
391                                     val = mpeg_decode_motion(s, s->mpeg_f_code[i][k],
392                                                              s->last_mv[i][j][k]);
393                                     s->last_mv[i][j][k] = val;
394                                     s->mv[i][j][k]      = val;
395                                 }
396                             }
397                         }
398                     }
399                 }
400                 break;
401             case MT_FIELD:
402                 if(s->progressive_sequence){
403                     av_log(s->avctx, AV_LOG_ERROR, "MT_FIELD in progressive_sequence\n");
404                     return -1;
405                 }
406                 s->mv_type = MV_TYPE_FIELD;
407                 if (s->picture_structure == PICT_FRAME) {
408                     mb_type |= MB_TYPE_16x8 | MB_TYPE_INTERLACED;
409                     for (i = 0; i < 2; i++) {
410                         if (USES_LIST(mb_type, i)) {
411                             for (j = 0; j < 2; j++) {
412                                 s->field_select[i][j] = get_bits1(&s->gb);
413                                 val = mpeg_decode_motion(s, s->mpeg_f_code[i][0],
414                                                          s->last_mv[i][j][0]);
415                                 s->last_mv[i][j][0] = val;
416                                 s->mv[i][j][0]      = val;
417                                 av_dlog(s->avctx, "fmx=%d\n", val);
418                                 val = mpeg_decode_motion(s, s->mpeg_f_code[i][1],
419                                                          s->last_mv[i][j][1] >> 1);
420                                 s->last_mv[i][j][1] = val << 1;
421                                 s->mv[i][j][1]      = val;
422                                 av_dlog(s->avctx, "fmy=%d\n", val);
423                             }
424                         }
425                     }
426                 } else {
427                     mb_type |= MB_TYPE_16x16 | MB_TYPE_INTERLACED;
428                     for (i = 0; i < 2; i++) {
429                         if (USES_LIST(mb_type, i)) {
430                             s->field_select[i][0] = get_bits1(&s->gb);
431                             for (k = 0; k < 2; k++) {
432                                 val = mpeg_decode_motion(s, s->mpeg_f_code[i][k],
433                                                          s->last_mv[i][0][k]);
434                                 s->last_mv[i][0][k] = val;
435                                 s->last_mv[i][1][k] = val;
436                                 s->mv[i][0][k]      = val;
437                             }
438                         }
439                     }
440                 }
441                 break;
442             case MT_DMV:
443                 if(s->progressive_sequence){
444                     av_log(s->avctx, AV_LOG_ERROR, "MT_DMV in progressive_sequence\n");
445                     return -1;
446                 }
447                 s->mv_type = MV_TYPE_DMV;
448                 for (i = 0; i < 2; i++) {
449                     if (USES_LIST(mb_type, i)) {
450                         int dmx, dmy, mx, my, m;
451                         const int my_shift = s->picture_structure == PICT_FRAME;
452
453                         mx = mpeg_decode_motion(s, s->mpeg_f_code[i][0],
454                                                 s->last_mv[i][0][0]);
455                         s->last_mv[i][0][0] = mx;
456                         s->last_mv[i][1][0] = mx;
457                         dmx = get_dmv(s);
458                         my  = mpeg_decode_motion(s, s->mpeg_f_code[i][1],
459                                                  s->last_mv[i][0][1] >> my_shift);
460                         dmy = get_dmv(s);
461
462
463                         s->last_mv[i][0][1] = my << my_shift;
464                         s->last_mv[i][1][1] = my << my_shift;
465
466                         s->mv[i][0][0] = mx;
467                         s->mv[i][0][1] = my;
468                         s->mv[i][1][0] = mx; // not used
469                         s->mv[i][1][1] = my; // not used
470
471                         if (s->picture_structure == PICT_FRAME) {
472                             mb_type |= MB_TYPE_16x16 | MB_TYPE_INTERLACED;
473
474                             // m = 1 + 2 * s->top_field_first;
475                             m = s->top_field_first ? 1 : 3;
476
477                             /* top -> top pred */
478                             s->mv[i][2][0] = ((mx * m + (mx > 0)) >> 1) + dmx;
479                             s->mv[i][2][1] = ((my * m + (my > 0)) >> 1) + dmy - 1;
480                             m = 4 - m;
481                             s->mv[i][3][0] = ((mx * m + (mx > 0)) >> 1) + dmx;
482                             s->mv[i][3][1] = ((my * m + (my > 0)) >> 1) + dmy + 1;
483                         } else {
484                             mb_type |= MB_TYPE_16x16;
485
486                             s->mv[i][2][0] = ((mx + (mx > 0)) >> 1) + dmx;
487                             s->mv[i][2][1] = ((my + (my > 0)) >> 1) + dmy;
488                             if (s->picture_structure == PICT_TOP_FIELD)
489                                 s->mv[i][2][1]--;
490                             else
491                                 s->mv[i][2][1]++;
492                         }
493                     }
494                 }
495                 break;
496             default:
497                 av_log(s->avctx, AV_LOG_ERROR, "00 motion_type at %d %d\n", s->mb_x, s->mb_y);
498                 return -1;
499             }
500         }
501
502         s->mb_intra = 0;
503         if (HAS_CBP(mb_type)) {
504             s->dsp.clear_blocks(s->block[0]);
505
506             cbp = get_vlc2(&s->gb, mb_pat_vlc.table, MB_PAT_VLC_BITS, 1);
507             if (mb_block_count > 6) {
508                  cbp <<= mb_block_count - 6;
509                  cbp  |= get_bits(&s->gb, mb_block_count - 6);
510                  s->dsp.clear_blocks(s->block[6]);
511             }
512             if (cbp <= 0) {
513                 av_log(s->avctx, AV_LOG_ERROR, "invalid cbp at %d %d\n", s->mb_x, s->mb_y);
514                 return -1;
515             }
516
517             //if 1, we memcpy blocks in xvmcvideo
518             if (CONFIG_MPEG_XVMC_DECODER && s->avctx->xvmc_acceleration > 1) {
519                 ff_xvmc_pack_pblocks(s, cbp);
520                 if (s->swap_uv) {
521                     exchange_uv(s);
522                 }
523             }
524
525             if (s->codec_id == CODEC_ID_MPEG2VIDEO) {
526                 if (s->flags2 & CODEC_FLAG2_FAST) {
527                     for (i = 0; i < 6; i++) {
528                         if (cbp & 32) {
529                             mpeg2_fast_decode_block_non_intra(s, *s->pblocks[i], i);
530                         } else {
531                             s->block_last_index[i] = -1;
532                         }
533                         cbp += cbp;
534                     }
535                 } else {
536                     cbp <<= 12-mb_block_count;
537
538                     for (i = 0; i < mb_block_count; i++) {
539                         if (cbp & (1 << 11)) {
540                             if (mpeg2_decode_block_non_intra(s, *s->pblocks[i], i) < 0)
541                                 return -1;
542                         } else {
543                             s->block_last_index[i] = -1;
544                         }
545                         cbp += cbp;
546                     }
547                 }
548             } else {
549                 if (s->flags2 & CODEC_FLAG2_FAST) {
550                     for (i = 0; i < 6; i++) {
551                         if (cbp & 32) {
552                             mpeg1_fast_decode_block_inter(s, *s->pblocks[i], i);
553                         } else {
554                             s->block_last_index[i] = -1;
555                         }
556                         cbp += cbp;
557                     }
558                 } else {
559                     for (i = 0; i < 6; i++) {
560                         if (cbp & 32) {
561                             if (mpeg1_decode_block_inter(s, *s->pblocks[i], i) < 0)
562                                 return -1;
563                         } else {
564                             s->block_last_index[i] = -1;
565                         }
566                         cbp += cbp;
567                     }
568                 }
569             }
570         } else {
571             for (i = 0; i < 12; i++)
572                 s->block_last_index[i] = -1;
573         }
574     }
575
576     s->current_picture.f.mb_type[s->mb_x + s->mb_y * s->mb_stride] = mb_type;
577
578     return 0;
579 }
580
581 /* as H.263, but only 17 codes */
582 static int mpeg_decode_motion(MpegEncContext *s, int fcode, int pred)
583 {
584     int code, sign, val, l, shift;
585
586     code = get_vlc2(&s->gb, mv_vlc.table, MV_VLC_BITS, 2);
587     if (code == 0) {
588         return pred;
589     }
590     if (code < 0) {
591         return 0xffff;
592     }
593
594     sign  = get_bits1(&s->gb);
595     shift = fcode - 1;
596     val   = code;
597     if (shift) {
598         val  = (val - 1) << shift;
599         val |= get_bits(&s->gb, shift);
600         val++;
601     }
602     if (sign)
603         val = -val;
604     val += pred;
605
606     /* modulo decoding */
607     l   = INT_BIT - 5 - shift;
608     val = (val << l) >> l;
609     return val;
610 }
611
612 static inline int mpeg1_decode_block_intra(MpegEncContext *s, DCTELEM *block, int n)
613 {
614     int level, dc, diff, i, j, run;
615     int component;
616     RLTable *rl = &ff_rl_mpeg1;
617     uint8_t * const scantable    = s->intra_scantable.permutated;
618     const uint16_t *quant_matrix = s->intra_matrix;
619     const int qscale             = s->qscale;
620
621     /* DC coefficient */
622     component = (n <= 3 ? 0 : n - 4 + 1);
623     diff = decode_dc(&s->gb, component);
624     if (diff >= 0xffff)
625         return -1;
626     dc  = s->last_dc[component];
627     dc += diff;
628     s->last_dc[component] = dc;
629     block[0] = dc * quant_matrix[0];
630     av_dlog(s->avctx, "dc=%d diff=%d\n", dc, diff);
631     i = 0;
632     {
633         OPEN_READER(re, &s->gb);
634         /* now quantify & encode AC coefficients */
635         for (;;) {
636             UPDATE_CACHE(re, &s->gb);
637             GET_RL_VLC(level, run, re, &s->gb, rl->rl_vlc[0], TEX_VLC_BITS, 2, 0);
638
639             if (level == 127) {
640                 break;
641             } else if (level != 0) {
642                 i += run;
643                 j = scantable[i];
644                 level = (level * qscale * quant_matrix[j]) >> 4;
645                 level = (level - 1) | 1;
646                 level = (level ^ SHOW_SBITS(re, &s->gb, 1)) - SHOW_SBITS(re, &s->gb, 1);
647                 LAST_SKIP_BITS(re, &s->gb, 1);
648             } else {
649                 /* escape */
650                 run = SHOW_UBITS(re, &s->gb, 6) + 1; LAST_SKIP_BITS(re, &s->gb, 6);
651                 UPDATE_CACHE(re, &s->gb);
652                 level = SHOW_SBITS(re, &s->gb, 8); SKIP_BITS(re, &s->gb, 8);
653                 if (level == -128) {
654                     level = SHOW_UBITS(re, &s->gb, 8) - 256; LAST_SKIP_BITS(re, &s->gb, 8);
655                 } else if (level == 0) {
656                     level = SHOW_UBITS(re, &s->gb, 8)      ; LAST_SKIP_BITS(re, &s->gb, 8);
657                 }
658                 i += run;
659                 j = scantable[i];
660                 if (level < 0) {
661                     level = -level;
662                     level = (level * qscale * quant_matrix[j]) >> 4;
663                     level = (level - 1) | 1;
664                     level = -level;
665                 } else {
666                     level = (level * qscale * quant_matrix[j]) >> 4;
667                     level = (level - 1) | 1;
668                 }
669             }
670             if (i > 63) {
671                 av_log(s->avctx, AV_LOG_ERROR, "ac-tex damaged at %d %d\n", s->mb_x, s->mb_y);
672                 return -1;
673             }
674
675             block[j] = level;
676         }
677         CLOSE_READER(re, &s->gb);
678     }
679     s->block_last_index[n] = i;
680    return 0;
681 }
682
683 int ff_mpeg1_decode_block_intra(MpegEncContext *s, DCTELEM *block, int n)
684 {
685     return mpeg1_decode_block_intra(s, block, n);
686 }
687
688 static inline int mpeg1_decode_block_inter(MpegEncContext *s, DCTELEM *block, int n)
689 {
690     int level, i, j, run;
691     RLTable *rl = &ff_rl_mpeg1;
692     uint8_t * const scantable    = s->intra_scantable.permutated;
693     const uint16_t *quant_matrix = s->inter_matrix;
694     const int qscale             = s->qscale;
695
696     {
697         OPEN_READER(re, &s->gb);
698         i = -1;
699         // special case for first coefficient, no need to add second VLC table
700         UPDATE_CACHE(re, &s->gb);
701         if (((int32_t)GET_CACHE(re, &s->gb)) < 0) {
702             level = (3 * qscale * quant_matrix[0]) >> 5;
703             level = (level - 1) | 1;
704             if (GET_CACHE(re, &s->gb) & 0x40000000)
705                 level = -level;
706             block[0] = level;
707             i++;
708             SKIP_BITS(re, &s->gb, 2);
709             if (((int32_t)GET_CACHE(re, &s->gb)) <= (int32_t)0xBFFFFFFF)
710                 goto end;
711         }
712         /* now quantify & encode AC coefficients */
713         for (;;) {
714             GET_RL_VLC(level, run, re, &s->gb, rl->rl_vlc[0], TEX_VLC_BITS, 2, 0);
715
716             if (level != 0) {
717                 i += run;
718                 j = scantable[i];
719                 level = ((level * 2 + 1) * qscale * quant_matrix[j]) >> 5;
720                 level = (level - 1) | 1;
721                 level = (level ^ SHOW_SBITS(re, &s->gb, 1)) - SHOW_SBITS(re, &s->gb, 1);
722                 SKIP_BITS(re, &s->gb, 1);
723             } else {
724                 /* escape */
725                 run = SHOW_UBITS(re, &s->gb, 6) + 1; LAST_SKIP_BITS(re, &s->gb, 6);
726                 UPDATE_CACHE(re, &s->gb);
727                 level = SHOW_SBITS(re, &s->gb, 8); SKIP_BITS(re, &s->gb, 8);
728                 if (level == -128) {
729                     level = SHOW_UBITS(re, &s->gb, 8) - 256; SKIP_BITS(re, &s->gb, 8);
730                 } else if (level == 0) {
731                     level = SHOW_UBITS(re, &s->gb, 8)      ; SKIP_BITS(re, &s->gb, 8);
732                 }
733                 i += run;
734                 j = scantable[i];
735                 if (level < 0) {
736                     level = -level;
737                     level = ((level * 2 + 1) * qscale * quant_matrix[j]) >> 5;
738                     level = (level - 1) | 1;
739                     level = -level;
740                 } else {
741                     level = ((level * 2 + 1) * qscale * quant_matrix[j]) >> 5;
742                     level = (level - 1) | 1;
743                 }
744             }
745             if (i > 63) {
746                 av_log(s->avctx, AV_LOG_ERROR, "ac-tex damaged at %d %d\n", s->mb_x, s->mb_y);
747                 return -1;
748             }
749
750             block[j] = level;
751             if (((int32_t)GET_CACHE(re, &s->gb)) <= (int32_t)0xBFFFFFFF)
752                 break;
753             UPDATE_CACHE(re, &s->gb);
754         }
755 end:
756         LAST_SKIP_BITS(re, &s->gb, 2);
757         CLOSE_READER(re, &s->gb);
758     }
759     s->block_last_index[n] = i;
760     return 0;
761 }
762
763 static inline int mpeg1_fast_decode_block_inter(MpegEncContext *s, DCTELEM *block, int n)
764 {
765     int level, i, j, run;
766     RLTable *rl = &ff_rl_mpeg1;
767     uint8_t * const scantable = s->intra_scantable.permutated;
768     const int qscale          = s->qscale;
769
770     {
771         OPEN_READER(re, &s->gb);
772         i = -1;
773         // special case for first coefficient, no need to add second VLC table
774         UPDATE_CACHE(re, &s->gb);
775         if (((int32_t)GET_CACHE(re, &s->gb)) < 0) {
776             level = (3 * qscale) >> 1;
777             level = (level - 1) | 1;
778             if (GET_CACHE(re, &s->gb) & 0x40000000)
779                 level = -level;
780             block[0] = level;
781             i++;
782             SKIP_BITS(re, &s->gb, 2);
783             if (((int32_t)GET_CACHE(re, &s->gb)) <= (int32_t)0xBFFFFFFF)
784                 goto end;
785         }
786
787         /* now quantify & encode AC coefficients */
788         for (;;) {
789             GET_RL_VLC(level, run, re, &s->gb, rl->rl_vlc[0], TEX_VLC_BITS, 2, 0);
790
791             if (level != 0) {
792                 i += run;
793                 j = scantable[i];
794                 level = ((level * 2 + 1) * qscale) >> 1;
795                 level = (level - 1) | 1;
796                 level = (level ^ SHOW_SBITS(re, &s->gb, 1)) - SHOW_SBITS(re, &s->gb, 1);
797                 SKIP_BITS(re, &s->gb, 1);
798             } else {
799                 /* escape */
800                 run = SHOW_UBITS(re, &s->gb, 6)+1; LAST_SKIP_BITS(re, &s->gb, 6);
801                 UPDATE_CACHE(re, &s->gb);
802                 level = SHOW_SBITS(re, &s->gb, 8); SKIP_BITS(re, &s->gb, 8);
803                 if (level == -128) {
804                     level = SHOW_UBITS(re, &s->gb, 8) - 256; SKIP_BITS(re, &s->gb, 8);
805                 } else if (level == 0) {
806                     level = SHOW_UBITS(re, &s->gb, 8)      ; SKIP_BITS(re, &s->gb, 8);
807                 }
808                 i += run;
809                 j = scantable[i];
810                 if (level < 0) {
811                     level = -level;
812                     level = ((level * 2 + 1) * qscale) >> 1;
813                     level = (level - 1) | 1;
814                     level = -level;
815                 } else {
816                     level = ((level * 2 + 1) * qscale) >> 1;
817                     level = (level - 1) | 1;
818                 }
819             }
820
821             block[j] = level;
822             if (((int32_t)GET_CACHE(re, &s->gb)) <= (int32_t)0xBFFFFFFF)
823                 break;
824             UPDATE_CACHE(re, &s->gb);
825         }
826 end:
827         LAST_SKIP_BITS(re, &s->gb, 2);
828         CLOSE_READER(re, &s->gb);
829     }
830     s->block_last_index[n] = i;
831     return 0;
832 }
833
834
835 static inline int mpeg2_decode_block_non_intra(MpegEncContext *s, DCTELEM *block, int n)
836 {
837     int level, i, j, run;
838     RLTable *rl = &ff_rl_mpeg1;
839     uint8_t * const scantable = s->intra_scantable.permutated;
840     const uint16_t *quant_matrix;
841     const int qscale = s->qscale;
842     int mismatch;
843
844     mismatch = 1;
845
846     {
847         OPEN_READER(re, &s->gb);
848         i = -1;
849         if (n < 4)
850             quant_matrix = s->inter_matrix;
851         else
852             quant_matrix = s->chroma_inter_matrix;
853
854         // special case for first coefficient, no need to add second VLC table
855         UPDATE_CACHE(re, &s->gb);
856         if (((int32_t)GET_CACHE(re, &s->gb)) < 0) {
857             level= (3 * qscale * quant_matrix[0]) >> 5;
858             if (GET_CACHE(re, &s->gb) & 0x40000000)
859                 level = -level;
860             block[0]  = level;
861             mismatch ^= level;
862             i++;
863             SKIP_BITS(re, &s->gb, 2);
864             if (((int32_t)GET_CACHE(re, &s->gb)) <= (int32_t)0xBFFFFFFF)
865                 goto end;
866         }
867
868         /* now quantify & encode AC coefficients */
869         for (;;) {
870             GET_RL_VLC(level, run, re, &s->gb, rl->rl_vlc[0], TEX_VLC_BITS, 2, 0);
871
872             if (level != 0) {
873                 i += run;
874                 j = scantable[i];
875                 level = ((level * 2 + 1) * qscale * quant_matrix[j]) >> 5;
876                 level = (level ^ SHOW_SBITS(re, &s->gb, 1)) - SHOW_SBITS(re, &s->gb, 1);
877                 SKIP_BITS(re, &s->gb, 1);
878             } else {
879                 /* escape */
880                 run = SHOW_UBITS(re, &s->gb, 6) + 1; LAST_SKIP_BITS(re, &s->gb, 6);
881                 UPDATE_CACHE(re, &s->gb);
882                 level = SHOW_SBITS(re, &s->gb, 12); SKIP_BITS(re, &s->gb, 12);
883
884                 i += run;
885                 j = scantable[i];
886                 if (level < 0) {
887                     level = ((-level * 2 + 1) * qscale * quant_matrix[j]) >> 5;
888                     level = -level;
889                 } else {
890                     level = ((level * 2 + 1) * qscale * quant_matrix[j]) >> 5;
891                 }
892             }
893             if (i > 63) {
894                 av_log(s->avctx, AV_LOG_ERROR, "ac-tex damaged at %d %d\n", s->mb_x, s->mb_y);
895                 return -1;
896             }
897
898             mismatch ^= level;
899             block[j]  = level;
900             if (((int32_t)GET_CACHE(re, &s->gb)) <= (int32_t)0xBFFFFFFF)
901                 break;
902             UPDATE_CACHE(re, &s->gb);
903         }
904 end:
905         LAST_SKIP_BITS(re, &s->gb, 2);
906         CLOSE_READER(re, &s->gb);
907     }
908     block[63] ^= (mismatch & 1);
909
910     s->block_last_index[n] = i;
911     return 0;
912 }
913
914 static inline int mpeg2_fast_decode_block_non_intra(MpegEncContext *s,
915                                                     DCTELEM *block, int n)
916 {
917     int level, i, j, run;
918     RLTable *rl = &ff_rl_mpeg1;
919     uint8_t * const scantable = s->intra_scantable.permutated;
920     const int qscale          = s->qscale;
921     OPEN_READER(re, &s->gb);
922     i = -1;
923
924     // special case for first coefficient, no need to add second VLC table
925     UPDATE_CACHE(re, &s->gb);
926     if (((int32_t)GET_CACHE(re, &s->gb)) < 0) {
927         level = (3 * qscale) >> 1;
928         if (GET_CACHE(re, &s->gb) & 0x40000000)
929             level = -level;
930         block[0] = level;
931         i++;
932         SKIP_BITS(re, &s->gb, 2);
933         if (((int32_t)GET_CACHE(re, &s->gb)) <= (int32_t)0xBFFFFFFF)
934             goto end;
935     }
936
937     /* now quantify & encode AC coefficients */
938     for (;;) {
939         GET_RL_VLC(level, run, re, &s->gb, rl->rl_vlc[0], TEX_VLC_BITS, 2, 0);
940
941         if (level != 0) {
942             i += run;
943             j  = scantable[i];
944             level = ((level * 2 + 1) * qscale) >> 1;
945             level = (level ^ SHOW_SBITS(re, &s->gb, 1)) - SHOW_SBITS(re, &s->gb, 1);
946             SKIP_BITS(re, &s->gb, 1);
947         } else {
948             /* escape */
949             run = SHOW_UBITS(re, &s->gb, 6) + 1; LAST_SKIP_BITS(re, &s->gb, 6);
950             UPDATE_CACHE(re, &s->gb);
951             level = SHOW_SBITS(re, &s->gb, 12); SKIP_BITS(re, &s->gb, 12);
952
953             i += run;
954             j  = scantable[i];
955             if (level < 0) {
956                 level = ((-level * 2 + 1) * qscale) >> 1;
957                 level = -level;
958             } else {
959                 level = ((level * 2 + 1) * qscale) >> 1;
960             }
961         }
962
963         block[j] = level;
964         if (((int32_t)GET_CACHE(re, &s->gb)) <= (int32_t)0xBFFFFFFF)
965             break;
966         UPDATE_CACHE(re, &s->gb);
967     }
968 end:
969     LAST_SKIP_BITS(re, &s->gb, 2);
970     CLOSE_READER(re, &s->gb);
971     s->block_last_index[n] = i;
972     return 0;
973 }
974
975
976 static inline int mpeg2_decode_block_intra(MpegEncContext *s, DCTELEM *block, int n)
977 {
978     int level, dc, diff, i, j, run;
979     int component;
980     RLTable *rl;
981     uint8_t * const scantable = s->intra_scantable.permutated;
982     const uint16_t *quant_matrix;
983     const int qscale = s->qscale;
984     int mismatch;
985
986     /* DC coefficient */
987     if (n < 4) {
988         quant_matrix = s->intra_matrix;
989         component = 0;
990     } else {
991         quant_matrix = s->chroma_intra_matrix;
992         component = (n & 1) + 1;
993     }
994     diff = decode_dc(&s->gb, component);
995     if (diff >= 0xffff)
996         return -1;
997     dc  = s->last_dc[component];
998     dc += diff;
999     s->last_dc[component] = dc;
1000     block[0] = dc << (3 - s->intra_dc_precision);
1001     av_dlog(s->avctx, "dc=%d\n", block[0]);
1002     mismatch = block[0] ^ 1;
1003     i = 0;
1004     if (s->intra_vlc_format)
1005         rl = &ff_rl_mpeg2;
1006     else
1007         rl = &ff_rl_mpeg1;
1008
1009     {
1010         OPEN_READER(re, &s->gb);
1011         /* now quantify & encode AC coefficients */
1012         for (;;) {
1013             UPDATE_CACHE(re, &s->gb);
1014             GET_RL_VLC(level, run, re, &s->gb, rl->rl_vlc[0], TEX_VLC_BITS, 2, 0);
1015
1016             if (level == 127) {
1017                 break;
1018             } else if (level != 0) {
1019                 i += run;
1020                 j  = scantable[i];
1021                 level = (level * qscale * quant_matrix[j]) >> 4;
1022                 level = (level ^ SHOW_SBITS(re, &s->gb, 1)) - SHOW_SBITS(re, &s->gb, 1);
1023                 LAST_SKIP_BITS(re, &s->gb, 1);
1024             } else {
1025                 /* escape */
1026                 run = SHOW_UBITS(re, &s->gb, 6) + 1; LAST_SKIP_BITS(re, &s->gb, 6);
1027                 UPDATE_CACHE(re, &s->gb);
1028                 level = SHOW_SBITS(re, &s->gb, 12); SKIP_BITS(re, &s->gb, 12);
1029                 i += run;
1030                 j  = scantable[i];
1031                 if (level < 0) {
1032                     level = (-level * qscale * quant_matrix[j]) >> 4;
1033                     level = -level;
1034                 } else {
1035                     level = (level * qscale * quant_matrix[j]) >> 4;
1036                 }
1037             }
1038             if (i > 63) {
1039                 av_log(s->avctx, AV_LOG_ERROR, "ac-tex damaged at %d %d\n", s->mb_x, s->mb_y);
1040                 return -1;
1041             }
1042
1043             mismatch ^= level;
1044             block[j]  = level;
1045         }
1046         CLOSE_READER(re, &s->gb);
1047     }
1048     block[63] ^= mismatch & 1;
1049
1050     s->block_last_index[n] = i;
1051     return 0;
1052 }
1053
1054 static inline int mpeg2_fast_decode_block_intra(MpegEncContext *s, DCTELEM *block, int n)
1055 {
1056     int level, dc, diff, j, run;
1057     int component;
1058     RLTable *rl;
1059     uint8_t * scantable = s->intra_scantable.permutated;
1060     const uint16_t *quant_matrix;
1061     const int qscale = s->qscale;
1062
1063     /* DC coefficient */
1064     if (n < 4) {
1065         quant_matrix = s->intra_matrix;
1066         component = 0;
1067     } else {
1068         quant_matrix = s->chroma_intra_matrix;
1069         component = (n & 1) + 1;
1070     }
1071     diff = decode_dc(&s->gb, component);
1072     if (diff >= 0xffff)
1073         return -1;
1074     dc = s->last_dc[component];
1075     dc += diff;
1076     s->last_dc[component] = dc;
1077     block[0] = dc << (3 - s->intra_dc_precision);
1078     if (s->intra_vlc_format)
1079         rl = &ff_rl_mpeg2;
1080     else
1081         rl = &ff_rl_mpeg1;
1082
1083     {
1084         OPEN_READER(re, &s->gb);
1085         /* now quantify & encode AC coefficients */
1086         for (;;) {
1087             UPDATE_CACHE(re, &s->gb);
1088             GET_RL_VLC(level, run, re, &s->gb, rl->rl_vlc[0], TEX_VLC_BITS, 2, 0);
1089
1090             if (level == 127) {
1091                 break;
1092             } else if (level != 0) {
1093                 scantable += run;
1094                 j = *scantable;
1095                 level = (level * qscale * quant_matrix[j]) >> 4;
1096                 level = (level ^ SHOW_SBITS(re, &s->gb, 1)) - SHOW_SBITS(re, &s->gb, 1);
1097                 LAST_SKIP_BITS(re, &s->gb, 1);
1098             } else {
1099                 /* escape */
1100                 run = SHOW_UBITS(re, &s->gb, 6) + 1; LAST_SKIP_BITS(re, &s->gb, 6);
1101                 UPDATE_CACHE(re, &s->gb);
1102                 level = SHOW_SBITS(re, &s->gb, 12); SKIP_BITS(re, &s->gb, 12);
1103                 scantable += run;
1104                 j = *scantable;
1105                 if (level < 0) {
1106                     level = (-level * qscale * quant_matrix[j]) >> 4;
1107                     level = -level;
1108                 } else {
1109                     level = (level * qscale * quant_matrix[j]) >> 4;
1110                 }
1111             }
1112
1113             block[j] = level;
1114         }
1115         CLOSE_READER(re, &s->gb);
1116     }
1117
1118     s->block_last_index[n] = scantable - s->intra_scantable.permutated;
1119     return 0;
1120 }
1121
1122 typedef struct Mpeg1Context {
1123     MpegEncContext mpeg_enc_ctx;
1124     int mpeg_enc_ctx_allocated; /* true if decoding context allocated */
1125     int repeat_field; /* true if we must repeat the field */
1126     AVPanScan pan_scan;              /**< some temporary storage for the panscan */
1127     int slice_count;
1128     int swap_uv;//indicate VCR2
1129     int save_aspect_info;
1130     int save_width, save_height, save_progressive_seq;
1131     AVRational frame_rate_ext;       ///< MPEG-2 specific framerate modificator
1132     int sync;                        ///< Did we reach a sync point like a GOP/SEQ/KEYFrame?
1133 } Mpeg1Context;
1134
1135 static av_cold int mpeg_decode_init(AVCodecContext *avctx)
1136 {
1137     Mpeg1Context *s = avctx->priv_data;
1138     MpegEncContext *s2 = &s->mpeg_enc_ctx;
1139     int i;
1140
1141     /* we need some permutation to store matrices,
1142      * until MPV_common_init() sets the real permutation. */
1143     for (i = 0; i < 64; i++)
1144        s2->dsp.idct_permutation[i]=i;
1145
1146     MPV_decode_defaults(s2);
1147
1148     s->mpeg_enc_ctx.avctx  = avctx;
1149     s->mpeg_enc_ctx.flags  = avctx->flags;
1150     s->mpeg_enc_ctx.flags2 = avctx->flags2;
1151     ff_mpeg12_common_init(&s->mpeg_enc_ctx);
1152     ff_mpeg12_init_vlcs();
1153
1154     s->mpeg_enc_ctx_allocated      = 0;
1155     s->mpeg_enc_ctx.picture_number = 0;
1156     s->repeat_field                = 0;
1157     s->mpeg_enc_ctx.codec_id       = avctx->codec->id;
1158     avctx->color_range = AVCOL_RANGE_MPEG;
1159     if (avctx->codec->id == CODEC_ID_MPEG1VIDEO)
1160         avctx->chroma_sample_location = AVCHROMA_LOC_CENTER;
1161     else
1162         avctx->chroma_sample_location = AVCHROMA_LOC_LEFT;
1163     return 0;
1164 }
1165
1166 static int mpeg_decode_update_thread_context(AVCodecContext *avctx, const AVCodecContext *avctx_from)
1167 {
1168     Mpeg1Context *ctx = avctx->priv_data, *ctx_from = avctx_from->priv_data;
1169     MpegEncContext *s = &ctx->mpeg_enc_ctx, *s1 = &ctx_from->mpeg_enc_ctx;
1170     int err;
1171
1172     if (avctx == avctx_from || !ctx_from->mpeg_enc_ctx_allocated || !s1->context_initialized)
1173         return 0;
1174
1175     err = ff_mpeg_update_thread_context(avctx, avctx_from);
1176     if (err) return err;
1177
1178     if (!ctx->mpeg_enc_ctx_allocated)
1179         memcpy(s + 1, s1 + 1, sizeof(Mpeg1Context) - sizeof(MpegEncContext));
1180
1181     if (!(s->pict_type == FF_B_TYPE || s->low_delay))
1182         s->picture_number++;
1183
1184     return 0;
1185 }
1186
1187 static void quant_matrix_rebuild(uint16_t *matrix, const uint8_t *old_perm,
1188                                  const uint8_t *new_perm)
1189 {
1190     uint16_t temp_matrix[64];
1191     int i;
1192
1193     memcpy(temp_matrix, matrix, 64 * sizeof(uint16_t));
1194
1195     for (i = 0; i < 64; i++) {
1196         matrix[new_perm[i]] = temp_matrix[old_perm[i]];
1197     }
1198 }
1199
1200 static const enum PixelFormat mpeg1_hwaccel_pixfmt_list_420[] = {
1201 #if CONFIG_MPEG_XVMC_DECODER
1202     PIX_FMT_XVMC_MPEG2_IDCT,
1203     PIX_FMT_XVMC_MPEG2_MC,
1204 #endif
1205 #if CONFIG_MPEG1_VDPAU_HWACCEL
1206     PIX_FMT_VDPAU_MPEG1,
1207 #endif
1208     PIX_FMT_YUV420P,
1209     PIX_FMT_NONE
1210 };
1211
1212 static const enum PixelFormat mpeg2_hwaccel_pixfmt_list_420[] = {
1213 #if CONFIG_MPEG_XVMC_DECODER
1214     PIX_FMT_XVMC_MPEG2_IDCT,
1215     PIX_FMT_XVMC_MPEG2_MC,
1216 #endif
1217 #if CONFIG_MPEG2_VDPAU_HWACCEL
1218     PIX_FMT_VDPAU_MPEG2,
1219 #endif
1220 #if CONFIG_MPEG2_DXVA2_HWACCEL
1221     PIX_FMT_DXVA2_VLD,
1222 #endif
1223 #if CONFIG_MPEG2_VAAPI_HWACCEL
1224     PIX_FMT_VAAPI_VLD,
1225 #endif
1226     PIX_FMT_YUV420P,
1227     PIX_FMT_NONE
1228 };
1229
1230 static inline int uses_vdpau(AVCodecContext *avctx) {
1231     return avctx->pix_fmt == PIX_FMT_VDPAU_MPEG1 || avctx->pix_fmt == PIX_FMT_VDPAU_MPEG2;
1232 }
1233
1234 static enum PixelFormat mpeg_get_pixelformat(AVCodecContext *avctx)
1235 {
1236     Mpeg1Context *s1 = avctx->priv_data;
1237     MpegEncContext *s = &s1->mpeg_enc_ctx;
1238
1239     if(s->chroma_format < 2) {
1240         enum PixelFormat res;
1241         res = avctx->get_format(avctx,
1242                                 avctx->codec_id == CODEC_ID_MPEG1VIDEO ?
1243                                 mpeg1_hwaccel_pixfmt_list_420 :
1244                                 mpeg2_hwaccel_pixfmt_list_420);
1245         if (res != PIX_FMT_XVMC_MPEG2_IDCT && res != PIX_FMT_XVMC_MPEG2_MC) {
1246             avctx->xvmc_acceleration = 0;
1247         } else if (!avctx->xvmc_acceleration) {
1248             avctx->xvmc_acceleration = 2;
1249         }
1250         return res;
1251     } else if(s->chroma_format == 2)
1252         return PIX_FMT_YUV422P;
1253     else
1254         return PIX_FMT_YUV444P;
1255 }
1256
1257 /* Call this function when we know all parameters.
1258  * It may be called in different places for MPEG-1 and MPEG-2. */
1259 static int mpeg_decode_postinit(AVCodecContext *avctx)
1260 {
1261     Mpeg1Context *s1 = avctx->priv_data;
1262     MpegEncContext *s = &s1->mpeg_enc_ctx;
1263     uint8_t old_permutation[64];
1264
1265     if ((s1->mpeg_enc_ctx_allocated == 0) ||
1266         avctx->coded_width  != s->width   ||
1267         avctx->coded_height != s->height  ||
1268         s1->save_width           != s->width                ||
1269         s1->save_height          != s->height               ||
1270         s1->save_aspect_info     != s->aspect_ratio_info    ||
1271         s1->save_progressive_seq != s->progressive_sequence ||
1272         0)
1273     {
1274
1275         if (s1->mpeg_enc_ctx_allocated) {
1276             ParseContext pc = s->parse_context;
1277             s->parse_context.buffer = 0;
1278             MPV_common_end(s);
1279             s->parse_context = pc;
1280         }
1281
1282         if ((s->width == 0) || (s->height == 0))
1283             return -2;
1284
1285         avcodec_set_dimensions(avctx, s->width, s->height);
1286         avctx->bit_rate          = s->bit_rate;
1287         s1->save_aspect_info     = s->aspect_ratio_info;
1288         s1->save_width           = s->width;
1289         s1->save_height          = s->height;
1290         s1->save_progressive_seq = s->progressive_sequence;
1291
1292         /* low_delay may be forced, in this case we will have B-frames
1293          * that behave like P-frames. */
1294         avctx->has_b_frames = !(s->low_delay);
1295
1296         assert((avctx->sub_id == 1) == (avctx->codec_id == CODEC_ID_MPEG1VIDEO));
1297         if (avctx->codec_id == CODEC_ID_MPEG1VIDEO) {
1298             //MPEG-1 fps
1299             avctx->time_base.den = ff_frame_rate_tab[s->frame_rate_index].num;
1300             avctx->time_base.num = ff_frame_rate_tab[s->frame_rate_index].den;
1301             //MPEG-1 aspect
1302             avctx->sample_aspect_ratio = av_d2q(1.0/ff_mpeg1_aspect[s->aspect_ratio_info], 255);
1303             avctx->ticks_per_frame=1;
1304         } else {//MPEG-2
1305         //MPEG-2 fps
1306             av_reduce(&s->avctx->time_base.den,
1307                       &s->avctx->time_base.num,
1308                       ff_frame_rate_tab[s->frame_rate_index].num * s1->frame_rate_ext.num*2,
1309                       ff_frame_rate_tab[s->frame_rate_index].den * s1->frame_rate_ext.den,
1310                       1 << 30);
1311             avctx->ticks_per_frame = 2;
1312             //MPEG-2 aspect
1313             if (s->aspect_ratio_info > 1) {
1314                 AVRational dar =
1315                     av_mul_q(av_div_q(ff_mpeg2_aspect[s->aspect_ratio_info],
1316                                       (AVRational) {s1->pan_scan.width, s1->pan_scan.height}),
1317                              (AVRational) {s->width, s->height});
1318
1319                 // we ignore the spec here and guess a bit as reality does not match the spec, see for example
1320                 // res_change_ffmpeg_aspect.ts and sequence-display-aspect.mpg
1321                 // issue1613, 621, 562
1322                 if ((s1->pan_scan.width == 0) || (s1->pan_scan.height == 0) ||
1323                    (av_cmp_q(dar, (AVRational) {4, 3}) && av_cmp_q(dar, (AVRational) {16, 9}))) {
1324                     s->avctx->sample_aspect_ratio =
1325                         av_div_q(ff_mpeg2_aspect[s->aspect_ratio_info],
1326                                  (AVRational) {s->width, s->height});
1327                 } else {
1328                     s->avctx->sample_aspect_ratio =
1329                         av_div_q(ff_mpeg2_aspect[s->aspect_ratio_info],
1330                                  (AVRational) {s1->pan_scan.width, s1->pan_scan.height});
1331 //issue1613 4/3 16/9 -> 16/9
1332 //res_change_ffmpeg_aspect.ts 4/3 225/44 ->4/3
1333 //widescreen-issue562.mpg 4/3 16/9 -> 16/9
1334 //                    s->avctx->sample_aspect_ratio = av_mul_q(s->avctx->sample_aspect_ratio, (AVRational) {s->width, s->height});
1335 //av_log(NULL, AV_LOG_ERROR, "A %d/%d\n", ff_mpeg2_aspect[s->aspect_ratio_info].num, ff_mpeg2_aspect[s->aspect_ratio_info].den);
1336 //av_log(NULL, AV_LOG_ERROR, "B %d/%d\n", s->avctx->sample_aspect_ratio.num, s->avctx->sample_aspect_ratio.den);
1337                 }
1338             } else {
1339                 s->avctx->sample_aspect_ratio =
1340                     ff_mpeg2_aspect[s->aspect_ratio_info];
1341             }
1342         } // MPEG-2
1343
1344         avctx->pix_fmt = mpeg_get_pixelformat(avctx);
1345         avctx->hwaccel = ff_find_hwaccel(avctx->codec->id, avctx->pix_fmt);
1346         // until then pix_fmt may be changed right after codec init
1347         if (avctx->pix_fmt == PIX_FMT_XVMC_MPEG2_IDCT ||
1348             avctx->hwaccel )
1349             if (avctx->idct_algo == FF_IDCT_AUTO)
1350                 avctx->idct_algo = FF_IDCT_SIMPLE;
1351
1352         /* Quantization matrices may need reordering
1353          * if DCT permutation is changed. */
1354         memcpy(old_permutation, s->dsp.idct_permutation, 64 * sizeof(uint8_t));
1355
1356         if (MPV_common_init(s) < 0)
1357             return -2;
1358
1359         quant_matrix_rebuild(s->intra_matrix,        old_permutation, s->dsp.idct_permutation);
1360         quant_matrix_rebuild(s->inter_matrix,        old_permutation, s->dsp.idct_permutation);
1361         quant_matrix_rebuild(s->chroma_intra_matrix, old_permutation, s->dsp.idct_permutation);
1362         quant_matrix_rebuild(s->chroma_inter_matrix, old_permutation, s->dsp.idct_permutation);
1363
1364         s1->mpeg_enc_ctx_allocated = 1;
1365     }
1366     return 0;
1367 }
1368
1369 static int mpeg1_decode_picture(AVCodecContext *avctx,
1370                                 const uint8_t *buf, int buf_size)
1371 {
1372     Mpeg1Context *s1 = avctx->priv_data;
1373     MpegEncContext *s = &s1->mpeg_enc_ctx;
1374     int ref, f_code, vbv_delay;
1375
1376     init_get_bits(&s->gb, buf, buf_size*8);
1377
1378     ref = get_bits(&s->gb, 10); /* temporal ref */
1379     s->pict_type = get_bits(&s->gb, 3);
1380     if (s->pict_type == 0 || s->pict_type > 3)
1381         return -1;
1382
1383     vbv_delay = get_bits(&s->gb, 16);
1384     if (s->pict_type == AV_PICTURE_TYPE_P || s->pict_type == AV_PICTURE_TYPE_B) {
1385         s->full_pel[0] = get_bits1(&s->gb);
1386         f_code = get_bits(&s->gb, 3);
1387         if (f_code == 0 && avctx->error_recognition >= FF_ER_COMPLIANT)
1388             return -1;
1389         s->mpeg_f_code[0][0] = f_code;
1390         s->mpeg_f_code[0][1] = f_code;
1391     }
1392     if (s->pict_type == AV_PICTURE_TYPE_B) {
1393         s->full_pel[1] = get_bits1(&s->gb);
1394         f_code = get_bits(&s->gb, 3);
1395         if (f_code == 0 && avctx->error_recognition >= FF_ER_COMPLIANT)
1396             return -1;
1397         s->mpeg_f_code[1][0] = f_code;
1398         s->mpeg_f_code[1][1] = f_code;
1399     }
1400     s->current_picture.f.pict_type = s->pict_type;
1401     s->current_picture.f.key_frame = s->pict_type == AV_PICTURE_TYPE_I;
1402
1403     if (avctx->debug & FF_DEBUG_PICT_INFO)
1404         av_log(avctx, AV_LOG_DEBUG, "vbv_delay %d, ref %d type:%d\n", vbv_delay, ref, s->pict_type);
1405
1406     s->y_dc_scale = 8;
1407     s->c_dc_scale = 8;
1408     return 0;
1409 }
1410
1411 static void mpeg_decode_sequence_extension(Mpeg1Context *s1)
1412 {
1413     MpegEncContext *s= &s1->mpeg_enc_ctx;
1414     int horiz_size_ext, vert_size_ext;
1415     int bit_rate_ext;
1416
1417     skip_bits(&s->gb, 1); /* profile and level esc*/
1418     s->avctx->profile       = get_bits(&s->gb, 3);
1419     s->avctx->level         = get_bits(&s->gb, 4);
1420     s->progressive_sequence = get_bits1(&s->gb); /* progressive_sequence */
1421     s->chroma_format        = get_bits(&s->gb, 2); /* chroma_format 1=420, 2=422, 3=444 */
1422     horiz_size_ext          = get_bits(&s->gb, 2);
1423     vert_size_ext           = get_bits(&s->gb, 2);
1424     s->width  |= (horiz_size_ext << 12);
1425     s->height |= (vert_size_ext  << 12);
1426     bit_rate_ext = get_bits(&s->gb, 12);  /* XXX: handle it */
1427     s->bit_rate += (bit_rate_ext << 18) * 400;
1428     skip_bits1(&s->gb); /* marker */
1429     s->avctx->rc_buffer_size += get_bits(&s->gb, 8) * 1024 * 16 << 10;
1430
1431     s->low_delay = get_bits1(&s->gb);
1432     if (s->flags & CODEC_FLAG_LOW_DELAY)
1433         s->low_delay = 1;
1434
1435     s1->frame_rate_ext.num = get_bits(&s->gb, 2) + 1;
1436     s1->frame_rate_ext.den = get_bits(&s->gb, 5) + 1;
1437
1438     av_dlog(s->avctx, "sequence extension\n");
1439     s->codec_id      = s->avctx->codec_id = CODEC_ID_MPEG2VIDEO;
1440     s->avctx->sub_id = 2; /* indicates MPEG-2 found */
1441
1442     if (s->avctx->debug & FF_DEBUG_PICT_INFO)
1443         av_log(s->avctx, AV_LOG_DEBUG, "profile: %d, level: %d vbv buffer: %d, bitrate:%d\n",
1444                s->avctx->profile, s->avctx->level, s->avctx->rc_buffer_size, s->bit_rate);
1445
1446 }
1447
1448 static void mpeg_decode_sequence_display_extension(Mpeg1Context *s1)
1449 {
1450     MpegEncContext *s = &s1->mpeg_enc_ctx;
1451     int color_description, w, h;
1452
1453     skip_bits(&s->gb, 3); /* video format */
1454     color_description = get_bits1(&s->gb);
1455     if (color_description) {
1456         s->avctx->color_primaries = get_bits(&s->gb, 8);
1457         s->avctx->color_trc       = get_bits(&s->gb, 8);
1458         s->avctx->colorspace      = get_bits(&s->gb, 8);
1459     }
1460     w = get_bits(&s->gb, 14);
1461     skip_bits(&s->gb, 1); //marker
1462     h = get_bits(&s->gb, 14);
1463     // remaining 3 bits are zero padding
1464
1465     s1->pan_scan.width  = 16 * w;
1466     s1->pan_scan.height = 16 * h;
1467
1468     if (s->avctx->debug & FF_DEBUG_PICT_INFO)
1469         av_log(s->avctx, AV_LOG_DEBUG, "sde w:%d, h:%d\n", w, h);
1470 }
1471
1472 static void mpeg_decode_picture_display_extension(Mpeg1Context *s1)
1473 {
1474     MpegEncContext *s = &s1->mpeg_enc_ctx;
1475     int i, nofco;
1476
1477     nofco = 1;
1478     if (s->progressive_sequence) {
1479         if (s->repeat_first_field) {
1480             nofco++;
1481             if (s->top_field_first)
1482                 nofco++;
1483         }
1484     } else {
1485         if (s->picture_structure == PICT_FRAME) {
1486             nofco++;
1487             if (s->repeat_first_field)
1488                 nofco++;
1489         }
1490     }
1491     for (i = 0; i < nofco; i++) {
1492         s1->pan_scan.position[i][0] = get_sbits(&s->gb, 16);
1493         skip_bits(&s->gb, 1); // marker
1494         s1->pan_scan.position[i][1] = get_sbits(&s->gb, 16);
1495         skip_bits(&s->gb, 1); // marker
1496     }
1497
1498     if (s->avctx->debug & FF_DEBUG_PICT_INFO)
1499         av_log(s->avctx, AV_LOG_DEBUG, "pde (%d,%d) (%d,%d) (%d,%d)\n",
1500                s1->pan_scan.position[0][0], s1->pan_scan.position[0][1],
1501                s1->pan_scan.position[1][0], s1->pan_scan.position[1][1],
1502                s1->pan_scan.position[2][0], s1->pan_scan.position[2][1]);
1503 }
1504
1505 static int load_matrix(MpegEncContext *s, uint16_t matrix0[64], uint16_t matrix1[64], int intra)
1506 {
1507     int i;
1508
1509     for (i = 0; i < 64; i++) {
1510         int j = s->dsp.idct_permutation[ff_zigzag_direct[i]];
1511         int v = get_bits(&s->gb, 8);
1512         if (v == 0) {
1513             av_log(s->avctx, AV_LOG_ERROR, "matrix damaged\n");
1514             return -1;
1515         }
1516         if (intra && i == 0 && v != 8) {
1517             av_log(s->avctx, AV_LOG_ERROR, "intra matrix invalid, ignoring\n");
1518             v = 8; // needed by pink.mpg / issue1046
1519         }
1520         matrix0[j] = v;
1521         if (matrix1)
1522             matrix1[j] = v;
1523     }
1524     return 0;
1525 }
1526
1527 static void mpeg_decode_quant_matrix_extension(MpegEncContext *s)
1528 {
1529     av_dlog(s->avctx, "matrix extension\n");
1530
1531     if (get_bits1(&s->gb)) load_matrix(s, s->chroma_intra_matrix, s->intra_matrix, 1);
1532     if (get_bits1(&s->gb)) load_matrix(s, s->chroma_inter_matrix, s->inter_matrix, 0);
1533     if (get_bits1(&s->gb)) load_matrix(s, s->chroma_intra_matrix, NULL           , 1);
1534     if (get_bits1(&s->gb)) load_matrix(s, s->chroma_inter_matrix, NULL           , 0);
1535 }
1536
1537 static void mpeg_decode_picture_coding_extension(Mpeg1Context *s1)
1538 {
1539     MpegEncContext *s = &s1->mpeg_enc_ctx;
1540
1541     s->full_pel[0] = s->full_pel[1] = 0;
1542     s->mpeg_f_code[0][0] = get_bits(&s->gb, 4);
1543     s->mpeg_f_code[0][1] = get_bits(&s->gb, 4);
1544     s->mpeg_f_code[1][0] = get_bits(&s->gb, 4);
1545     s->mpeg_f_code[1][1] = get_bits(&s->gb, 4);
1546     if (!s->pict_type && s1->mpeg_enc_ctx_allocated) {
1547         av_log(s->avctx, AV_LOG_ERROR, "Missing picture start code, guessing missing values\n");
1548         if (s->mpeg_f_code[1][0] == 15 && s->mpeg_f_code[1][1] == 15) {
1549             if (s->mpeg_f_code[0][0] == 15 && s->mpeg_f_code[0][1] == 15)
1550                 s->pict_type = AV_PICTURE_TYPE_I;
1551             else
1552                 s->pict_type = AV_PICTURE_TYPE_P;
1553         } else
1554             s->pict_type = AV_PICTURE_TYPE_B;
1555         s->current_picture.f.pict_type = s->pict_type;
1556         s->current_picture.f.key_frame = s->pict_type == AV_PICTURE_TYPE_I;
1557     }
1558     s->intra_dc_precision         = get_bits(&s->gb, 2);
1559     s->picture_structure          = get_bits(&s->gb, 2);
1560     s->top_field_first            = get_bits1(&s->gb);
1561     s->frame_pred_frame_dct       = get_bits1(&s->gb);
1562     s->concealment_motion_vectors = get_bits1(&s->gb);
1563     s->q_scale_type               = get_bits1(&s->gb);
1564     s->intra_vlc_format           = get_bits1(&s->gb);
1565     s->alternate_scan             = get_bits1(&s->gb);
1566     s->repeat_first_field         = get_bits1(&s->gb);
1567     s->chroma_420_type            = get_bits1(&s->gb);
1568     s->progressive_frame          = get_bits1(&s->gb);
1569
1570     if (s->progressive_sequence && !s->progressive_frame) {
1571         s->progressive_frame = 1;
1572         av_log(s->avctx, AV_LOG_ERROR, "interlaced frame in progressive sequence, ignoring\n");
1573     }
1574
1575     if (s->picture_structure == 0 || (s->progressive_frame && s->picture_structure != PICT_FRAME)) {
1576         av_log(s->avctx, AV_LOG_ERROR, "picture_structure %d invalid, ignoring\n", s->picture_structure);
1577         s->picture_structure = PICT_FRAME;
1578     }
1579
1580     if (s->progressive_sequence && !s->frame_pred_frame_dct) {
1581         av_log(s->avctx, AV_LOG_ERROR, "invalid frame_pred_frame_dct\n");
1582     }
1583
1584     if (s->picture_structure == PICT_FRAME) {
1585         s->first_field = 0;
1586         s->v_edge_pos  = 16 * s->mb_height;
1587     } else {
1588         s->first_field ^= 1;
1589         s->v_edge_pos   = 8 * s->mb_height;
1590         memset(s->mbskip_table, 0, s->mb_stride * s->mb_height);
1591     }
1592
1593     if (s->alternate_scan) {
1594         ff_init_scantable(s->dsp.idct_permutation, &s->inter_scantable, ff_alternate_vertical_scan);
1595         ff_init_scantable(s->dsp.idct_permutation, &s->intra_scantable, ff_alternate_vertical_scan);
1596     } else {
1597         ff_init_scantable(s->dsp.idct_permutation, &s->inter_scantable, ff_zigzag_direct);
1598         ff_init_scantable(s->dsp.idct_permutation, &s->intra_scantable, ff_zigzag_direct);
1599     }
1600
1601     /* composite display not parsed */
1602     av_dlog(s->avctx, "intra_dc_precision=%d\n", s->intra_dc_precision);
1603     av_dlog(s->avctx, "picture_structure=%d\n", s->picture_structure);
1604     av_dlog(s->avctx, "top field first=%d\n", s->top_field_first);
1605     av_dlog(s->avctx, "repeat first field=%d\n", s->repeat_first_field);
1606     av_dlog(s->avctx, "conceal=%d\n", s->concealment_motion_vectors);
1607     av_dlog(s->avctx, "intra_vlc_format=%d\n", s->intra_vlc_format);
1608     av_dlog(s->avctx, "alternate_scan=%d\n", s->alternate_scan);
1609     av_dlog(s->avctx, "frame_pred_frame_dct=%d\n", s->frame_pred_frame_dct);
1610     av_dlog(s->avctx, "progressive_frame=%d\n", s->progressive_frame);
1611 }
1612
1613 static void exchange_uv(MpegEncContext *s)
1614 {
1615     DCTELEM (*tmp)[64];
1616
1617     tmp           = s->pblocks[4];
1618     s->pblocks[4] = s->pblocks[5];
1619     s->pblocks[5] = tmp;
1620 }
1621
1622 static int mpeg_field_start(MpegEncContext *s, const uint8_t *buf, int buf_size)
1623 {
1624     AVCodecContext *avctx = s->avctx;
1625     Mpeg1Context *s1 = (Mpeg1Context*)s;
1626
1627     /* start frame decoding */
1628     if (s->first_field || s->picture_structure == PICT_FRAME) {
1629         if (MPV_frame_start(s, avctx) < 0)
1630             return -1;
1631
1632         ff_er_frame_start(s);
1633
1634         /* first check if we must repeat the frame */
1635         s->current_picture_ptr->f.repeat_pict = 0;
1636         if (s->repeat_first_field) {
1637             if (s->progressive_sequence) {
1638                 if (s->top_field_first)
1639                     s->current_picture_ptr->f.repeat_pict = 4;
1640                 else
1641                     s->current_picture_ptr->f.repeat_pict = 2;
1642             } else if (s->progressive_frame) {
1643                 s->current_picture_ptr->f.repeat_pict = 1;
1644             }
1645         }
1646
1647         *s->current_picture_ptr->f.pan_scan = s1->pan_scan;
1648
1649         if (HAVE_PTHREADS && (avctx->active_thread_type & FF_THREAD_FRAME))
1650             ff_thread_finish_setup(avctx);
1651     } else { // second field
1652         int i;
1653
1654         if (!s->current_picture_ptr) {
1655             av_log(s->avctx, AV_LOG_ERROR, "first field missing\n");
1656             return -1;
1657         }
1658
1659         for (i = 0; i < 4; i++) {
1660             s->current_picture.f.data[i] = s->current_picture_ptr->f.data[i];
1661             if (s->picture_structure == PICT_BOTTOM_FIELD) {
1662                 s->current_picture.f.data[i] += s->current_picture_ptr->f.linesize[i];
1663             }
1664         }
1665     }
1666
1667     if (avctx->hwaccel) {
1668         if (avctx->hwaccel->start_frame(avctx, buf, buf_size) < 0)
1669             return -1;
1670     }
1671
1672 // MPV_frame_start will call this function too,
1673 // but we need to call it on every field
1674     if (CONFIG_MPEG_XVMC_DECODER && s->avctx->xvmc_acceleration)
1675         if (ff_xvmc_field_start(s, avctx) < 0)
1676             return -1;
1677
1678     return 0;
1679 }
1680
1681 #define DECODE_SLICE_ERROR -1
1682 #define DECODE_SLICE_OK     0
1683
1684 /**
1685  * decodes a slice. MpegEncContext.mb_y must be set to the MB row from the startcode
1686  * @return DECODE_SLICE_ERROR if the slice is damaged<br>
1687  *         DECODE_SLICE_OK if this slice is ok<br>
1688  */
1689 static int mpeg_decode_slice(Mpeg1Context *s1, int mb_y,
1690                              const uint8_t **buf, int buf_size)
1691 {
1692     MpegEncContext *s     = &s1->mpeg_enc_ctx;
1693     AVCodecContext *avctx = s->avctx;
1694     const int lowres      = s->avctx->lowres;
1695     const int field_pic   = s->picture_structure != PICT_FRAME;
1696
1697     s->resync_mb_x =
1698     s->resync_mb_y = -1;
1699
1700     assert(mb_y < s->mb_height);
1701
1702     init_get_bits(&s->gb, *buf, buf_size * 8);
1703
1704     ff_mpeg1_clean_buffers(s);
1705     s->interlaced_dct = 0;
1706
1707     s->qscale = get_qscale(s);
1708
1709     if (s->qscale == 0) {
1710         av_log(s->avctx, AV_LOG_ERROR, "qscale == 0\n");
1711         return -1;
1712     }
1713
1714     /* extra slice info */
1715     while (get_bits1(&s->gb) != 0) {
1716         skip_bits(&s->gb, 8);
1717     }
1718
1719     s->mb_x = 0;
1720
1721     if (mb_y == 0 && s->codec_tag == AV_RL32("SLIF")) {
1722         skip_bits1(&s->gb);
1723     } else {
1724         for (;;) {
1725             int code = get_vlc2(&s->gb, mbincr_vlc.table, MBINCR_VLC_BITS, 2);
1726             if (code < 0) {
1727                 av_log(s->avctx, AV_LOG_ERROR, "first mb_incr damaged\n");
1728                 return -1;
1729             }
1730             if (code >= 33) {
1731                 if (code == 33) {
1732                     s->mb_x += 33;
1733                 }
1734                 /* otherwise, stuffing, nothing to do */
1735             } else {
1736                 s->mb_x += code;
1737                 break;
1738             }
1739         }
1740     }
1741
1742     if (s->mb_x >= (unsigned)s->mb_width) {
1743         av_log(s->avctx, AV_LOG_ERROR, "initial skip overflow\n");
1744         return -1;
1745     }
1746
1747     if (avctx->hwaccel) {
1748         const uint8_t *buf_end, *buf_start = *buf - 4; /* include start_code */
1749         int start_code = -1;
1750         buf_end = ff_find_start_code(buf_start + 2, *buf + buf_size, &start_code);
1751         if (buf_end < *buf + buf_size)
1752             buf_end -= 4;
1753         s->mb_y = mb_y;
1754         if (avctx->hwaccel->decode_slice(avctx, buf_start, buf_end - buf_start) < 0)
1755             return DECODE_SLICE_ERROR;
1756         *buf = buf_end;
1757         return DECODE_SLICE_OK;
1758     }
1759
1760     s->resync_mb_x = s->mb_x;
1761     s->resync_mb_y = s->mb_y = mb_y;
1762     s->mb_skip_run = 0;
1763     ff_init_block_index(s);
1764
1765     if (s->mb_y == 0 && s->mb_x == 0 && (s->first_field || s->picture_structure == PICT_FRAME)) {
1766         if (s->avctx->debug & FF_DEBUG_PICT_INFO) {
1767              av_log(s->avctx, AV_LOG_DEBUG, "qp:%d fc:%2d%2d%2d%2d %s %s %s %s %s dc:%d pstruct:%d fdct:%d cmv:%d qtype:%d ivlc:%d rff:%d %s\n",
1768                     s->qscale, s->mpeg_f_code[0][0], s->mpeg_f_code[0][1], s->mpeg_f_code[1][0], s->mpeg_f_code[1][1],
1769                     s->pict_type == AV_PICTURE_TYPE_I ? "I" : (s->pict_type == AV_PICTURE_TYPE_P ? "P" : (s->pict_type == AV_PICTURE_TYPE_B ? "B" : "S")),
1770                     s->progressive_sequence ? "ps" :"", s->progressive_frame ? "pf" : "", s->alternate_scan ? "alt" :"", s->top_field_first ? "top" :"",
1771                     s->intra_dc_precision, s->picture_structure, s->frame_pred_frame_dct, s->concealment_motion_vectors,
1772                     s->q_scale_type, s->intra_vlc_format, s->repeat_first_field, s->chroma_420_type ? "420" :"");
1773         }
1774     }
1775
1776     for (;;) {
1777         // If 1, we memcpy blocks in xvmcvideo.
1778         if (CONFIG_MPEG_XVMC_DECODER && s->avctx->xvmc_acceleration > 1)
1779             ff_xvmc_init_block(s); // set s->block
1780
1781         if (mpeg_decode_mb(s, s->block) < 0)
1782             return -1;
1783
1784         if (s->current_picture.f.motion_val[0] && !s->encoding) { // note motion_val is normally NULL unless we want to extract the MVs
1785             const int wrap = s->b8_stride;
1786             int xy         = s->mb_x * 2 + s->mb_y * 2 * wrap;
1787             int b8_xy      = 4 * (s->mb_x + s->mb_y * s->mb_stride);
1788             int motion_x, motion_y, dir, i;
1789
1790             for (i = 0; i < 2; i++) {
1791                 for (dir = 0; dir < 2; dir++) {
1792                     if (s->mb_intra || (dir == 1 && s->pict_type != AV_PICTURE_TYPE_B)) {
1793                         motion_x = motion_y = 0;
1794                     } else if (s->mv_type == MV_TYPE_16X16 || (s->mv_type == MV_TYPE_FIELD && field_pic)) {
1795                         motion_x = s->mv[dir][0][0];
1796                         motion_y = s->mv[dir][0][1];
1797                     } else /*if ((s->mv_type == MV_TYPE_FIELD) || (s->mv_type == MV_TYPE_16X8))*/ {
1798                         motion_x = s->mv[dir][i][0];
1799                         motion_y = s->mv[dir][i][1];
1800                     }
1801
1802                     s->current_picture.f.motion_val[dir][xy    ][0] = motion_x;
1803                     s->current_picture.f.motion_val[dir][xy    ][1] = motion_y;
1804                     s->current_picture.f.motion_val[dir][xy + 1][0] = motion_x;
1805                     s->current_picture.f.motion_val[dir][xy + 1][1] = motion_y;
1806                     s->current_picture.f.ref_index [dir][b8_xy    ] =
1807                     s->current_picture.f.ref_index [dir][b8_xy + 1] = s->field_select[dir][i];
1808                     assert(s->field_select[dir][i] == 0 || s->field_select[dir][i] == 1);
1809                 }
1810                 xy += wrap;
1811                 b8_xy +=2;
1812             }
1813         }
1814
1815         s->dest[0] += 16 >> lowres;
1816         s->dest[1] +=(16 >> lowres) >> s->chroma_x_shift;
1817         s->dest[2] +=(16 >> lowres) >> s->chroma_x_shift;
1818
1819         MPV_decode_mb(s, s->block);
1820
1821         if (++s->mb_x >= s->mb_width) {
1822             const int mb_size = 16 >> s->avctx->lowres;
1823
1824             ff_draw_horiz_band(s, mb_size*(s->mb_y >> field_pic), mb_size);
1825             MPV_report_decode_progress(s);
1826
1827             s->mb_x = 0;
1828             s->mb_y += 1 << field_pic;
1829
1830             if (s->mb_y >= s->mb_height) {
1831                 int left   = get_bits_left(&s->gb);
1832                 int is_d10 = s->chroma_format == 2 && s->pict_type == AV_PICTURE_TYPE_I && avctx->profile == 0 && avctx->level == 5
1833                              && s->intra_dc_precision == 2 && s->q_scale_type == 1 && s->alternate_scan == 0
1834                              && s->progressive_frame == 0 /* vbv_delay == 0xBBB || 0xE10*/;
1835
1836                 if (left < 0 || (left && show_bits(&s->gb, FFMIN(left, 23)) && !is_d10)
1837                     || (avctx->error_recognition >= FF_ER_AGGRESSIVE && left > 8)) {
1838                     av_log(avctx, AV_LOG_ERROR, "end mismatch left=%d %0X\n", left, show_bits(&s->gb, FFMIN(left, 23)));
1839                     return -1;
1840                 } else
1841                     goto eos;
1842             }
1843
1844             ff_init_block_index(s);
1845         }
1846
1847         /* skip mb handling */
1848         if (s->mb_skip_run == -1) {
1849             /* read increment again */
1850             s->mb_skip_run = 0;
1851             for (;;) {
1852                 int code = get_vlc2(&s->gb, mbincr_vlc.table, MBINCR_VLC_BITS, 2);
1853                 if (code < 0) {
1854                     av_log(s->avctx, AV_LOG_ERROR, "mb incr damaged\n");
1855                     return -1;
1856                 }
1857                 if (code >= 33) {
1858                     if (code == 33) {
1859                         s->mb_skip_run += 33;
1860                     } else if (code == 35) {
1861                         if (s->mb_skip_run != 0 || show_bits(&s->gb, 15) != 0) {
1862                             av_log(s->avctx, AV_LOG_ERROR, "slice mismatch\n");
1863                             return -1;
1864                         }
1865                         goto eos; /* end of slice */
1866                     }
1867                     /* otherwise, stuffing, nothing to do */
1868                 } else {
1869                     s->mb_skip_run += code;
1870                     break;
1871                 }
1872             }
1873             if (s->mb_skip_run) {
1874                 int i;
1875                 if (s->pict_type == AV_PICTURE_TYPE_I) {
1876                     av_log(s->avctx, AV_LOG_ERROR, "skipped MB in I frame at %d %d\n", s->mb_x, s->mb_y);
1877                     return -1;
1878                 }
1879
1880                 /* skip mb */
1881                 s->mb_intra = 0;
1882                 for (i = 0; i < 12; i++)
1883                     s->block_last_index[i] = -1;
1884                 if (s->picture_structure == PICT_FRAME)
1885                     s->mv_type = MV_TYPE_16X16;
1886                 else
1887                     s->mv_type = MV_TYPE_FIELD;
1888                 if (s->pict_type == AV_PICTURE_TYPE_P) {
1889                     /* if P type, zero motion vector is implied */
1890                     s->mv_dir             = MV_DIR_FORWARD;
1891                     s->mv[0][0][0]        = s->mv[0][0][1]      = 0;
1892                     s->last_mv[0][0][0]   = s->last_mv[0][0][1] = 0;
1893                     s->last_mv[0][1][0]   = s->last_mv[0][1][1] = 0;
1894                     s->field_select[0][0] = (s->picture_structure - 1) & 1;
1895                 } else {
1896                     /* if B type, reuse previous vectors and directions */
1897                     s->mv[0][0][0] = s->last_mv[0][0][0];
1898                     s->mv[0][0][1] = s->last_mv[0][0][1];
1899                     s->mv[1][0][0] = s->last_mv[1][0][0];
1900                     s->mv[1][0][1] = s->last_mv[1][0][1];
1901                 }
1902             }
1903         }
1904     }
1905 eos: // end of slice
1906     *buf += (get_bits_count(&s->gb)-1)/8;
1907 //printf("y %d %d %d %d\n", s->resync_mb_x, s->resync_mb_y, s->mb_x, s->mb_y);
1908     return 0;
1909 }
1910
1911 static int slice_decode_thread(AVCodecContext *c, void *arg)
1912 {
1913     MpegEncContext *s   = *(void**)arg;
1914     const uint8_t *buf  = s->gb.buffer;
1915     int mb_y            = s->start_mb_y;
1916     const int field_pic = s->picture_structure != PICT_FRAME;
1917
1918     s->error_count = (3 * (s->end_mb_y - s->start_mb_y) * s->mb_width) >> field_pic;
1919
1920     for (;;) {
1921         uint32_t start_code;
1922         int ret;
1923
1924         ret = mpeg_decode_slice((Mpeg1Context*)s, mb_y, &buf, s->gb.buffer_end - buf);
1925         emms_c();
1926 //av_log(c, AV_LOG_DEBUG, "ret:%d resync:%d/%d mb:%d/%d ts:%d/%d ec:%d\n",
1927 //ret, s->resync_mb_x, s->resync_mb_y, s->mb_x, s->mb_y, s->start_mb_y, s->end_mb_y, s->error_count);
1928         if (ret < 0) {
1929             if (c->error_recognition >= FF_ER_EXPLODE)
1930                 return ret;
1931             if (s->resync_mb_x >= 0 && s->resync_mb_y >= 0)
1932                 ff_er_add_slice(s, s->resync_mb_x, s->resync_mb_y, s->mb_x, s->mb_y, AC_ERROR | DC_ERROR | MV_ERROR);
1933         } else {
1934             ff_er_add_slice(s, s->resync_mb_x, s->resync_mb_y, s->mb_x-1, s->mb_y, AC_END | DC_END | MV_END);
1935         }
1936
1937         if (s->mb_y == s->end_mb_y)
1938             return 0;
1939
1940         start_code = -1;
1941         buf = ff_find_start_code(buf, s->gb.buffer_end, &start_code);
1942         mb_y= (start_code - SLICE_MIN_START_CODE) << field_pic;
1943         if (s->picture_structure == PICT_BOTTOM_FIELD)
1944             mb_y++;
1945         if (mb_y < 0 || mb_y >= s->end_mb_y)
1946             return -1;
1947     }
1948 }
1949
1950 /**
1951  * Handle slice ends.
1952  * @return 1 if it seems to be the last slice
1953  */
1954 static int slice_end(AVCodecContext *avctx, AVFrame *pict)
1955 {
1956     Mpeg1Context *s1 = avctx->priv_data;
1957     MpegEncContext *s = &s1->mpeg_enc_ctx;
1958
1959     if (!s1->mpeg_enc_ctx_allocated || !s->current_picture_ptr)
1960         return 0;
1961
1962     if (s->avctx->hwaccel) {
1963         if (s->avctx->hwaccel->end_frame(s->avctx) < 0)
1964             av_log(avctx, AV_LOG_ERROR, "hardware accelerator failed to decode picture\n");
1965     }
1966
1967     if (CONFIG_MPEG_XVMC_DECODER && s->avctx->xvmc_acceleration)
1968         ff_xvmc_field_end(s);
1969
1970     /* end of slice reached */
1971     if (/*s->mb_y << field_pic == s->mb_height &&*/ !s->first_field && !s->first_slice) {
1972         /* end of image */
1973
1974         s->current_picture_ptr->f.qscale_type = FF_QSCALE_TYPE_MPEG2;
1975
1976         ff_er_frame_end(s);
1977
1978         MPV_frame_end(s);
1979
1980         if (s->pict_type == AV_PICTURE_TYPE_B || s->low_delay) {
1981             *pict = *(AVFrame*)s->current_picture_ptr;
1982             ff_print_debug_info(s, pict);
1983         } else {
1984             if (avctx->active_thread_type & FF_THREAD_FRAME)
1985                 s->picture_number++;
1986             /* latency of 1 frame for I- and P-frames */
1987             /* XXX: use another variable than picture_number */
1988             if (s->last_picture_ptr != NULL) {
1989                 *pict = *(AVFrame*)s->last_picture_ptr;
1990                  ff_print_debug_info(s, pict);
1991             }
1992         }
1993
1994         return 1;
1995     } else {
1996         return 0;
1997     }
1998 }
1999
2000 static int mpeg1_decode_sequence(AVCodecContext *avctx,
2001                                  const uint8_t *buf, int buf_size)
2002 {
2003     Mpeg1Context *s1 = avctx->priv_data;
2004     MpegEncContext *s = &s1->mpeg_enc_ctx;
2005     int width, height;
2006     int i, v, j;
2007
2008     init_get_bits(&s->gb, buf, buf_size*8);
2009
2010     width  = get_bits(&s->gb, 12);
2011     height = get_bits(&s->gb, 12);
2012     if (width <= 0 || height <= 0)
2013         return -1;
2014     s->aspect_ratio_info = get_bits(&s->gb, 4);
2015     if (s->aspect_ratio_info == 0) {
2016         av_log(avctx, AV_LOG_ERROR, "aspect ratio has forbidden 0 value\n");
2017         if (avctx->error_recognition >= FF_ER_COMPLIANT)
2018             return -1;
2019     }
2020     s->frame_rate_index = get_bits(&s->gb, 4);
2021     if (s->frame_rate_index == 0 || s->frame_rate_index > 13)
2022         return -1;
2023     s->bit_rate = get_bits(&s->gb, 18) * 400;
2024     if (get_bits1(&s->gb) == 0) /* marker */
2025         return -1;
2026     s->width  = width;
2027     s->height = height;
2028
2029     s->avctx->rc_buffer_size = get_bits(&s->gb, 10) * 1024 * 16;
2030     skip_bits(&s->gb, 1);
2031
2032     /* get matrix */
2033     if (get_bits1(&s->gb)) {
2034         load_matrix(s, s->chroma_intra_matrix, s->intra_matrix, 1);
2035     } else {
2036         for (i = 0; i < 64; i++) {
2037             j = s->dsp.idct_permutation[i];
2038             v = ff_mpeg1_default_intra_matrix[i];
2039             s->intra_matrix[j]        = v;
2040             s->chroma_intra_matrix[j] = v;
2041         }
2042     }
2043     if (get_bits1(&s->gb)) {
2044         load_matrix(s, s->chroma_inter_matrix, s->inter_matrix, 0);
2045     } else {
2046         for (i = 0; i < 64; i++) {
2047             int j = s->dsp.idct_permutation[i];
2048             v = ff_mpeg1_default_non_intra_matrix[i];
2049             s->inter_matrix[j]        = v;
2050             s->chroma_inter_matrix[j] = v;
2051         }
2052     }
2053
2054     if (show_bits(&s->gb, 23) != 0) {
2055         av_log(s->avctx, AV_LOG_ERROR, "sequence header damaged\n");
2056         return -1;
2057     }
2058
2059     /* we set MPEG-2 parameters so that it emulates MPEG-1 */
2060     s->progressive_sequence = 1;
2061     s->progressive_frame    = 1;
2062     s->picture_structure    = PICT_FRAME;
2063     s->frame_pred_frame_dct = 1;
2064     s->chroma_format        = 1;
2065     s->codec_id             = s->avctx->codec_id = CODEC_ID_MPEG1VIDEO;
2066     avctx->sub_id           = 1; /* indicates MPEG-1 */
2067     s->out_format           = FMT_MPEG1;
2068     s->swap_uv              = 0; // AFAIK VCR2 does not have SEQ_HEADER
2069     if (s->flags & CODEC_FLAG_LOW_DELAY)
2070         s->low_delay = 1;
2071
2072     if (s->avctx->debug & FF_DEBUG_PICT_INFO)
2073         av_log(s->avctx, AV_LOG_DEBUG, "vbv buffer: %d, bitrate:%d\n",
2074                s->avctx->rc_buffer_size, s->bit_rate);
2075
2076     return 0;
2077 }
2078
2079 static int vcr2_init_sequence(AVCodecContext *avctx)
2080 {
2081     Mpeg1Context *s1 = avctx->priv_data;
2082     MpegEncContext *s = &s1->mpeg_enc_ctx;
2083     int i, v;
2084
2085     /* start new MPEG-1 context decoding */
2086     s->out_format = FMT_MPEG1;
2087     if (s1->mpeg_enc_ctx_allocated) {
2088         MPV_common_end(s);
2089     }
2090     s->width  = avctx->coded_width;
2091     s->height = avctx->coded_height;
2092     avctx->has_b_frames = 0; // true?
2093     s->low_delay = 1;
2094
2095     avctx->pix_fmt = mpeg_get_pixelformat(avctx);
2096     avctx->hwaccel = ff_find_hwaccel(avctx->codec->id, avctx->pix_fmt);
2097
2098     if( avctx->pix_fmt == PIX_FMT_XVMC_MPEG2_IDCT || avctx->hwaccel )
2099         if (avctx->idct_algo == FF_IDCT_AUTO)
2100             avctx->idct_algo = FF_IDCT_SIMPLE;
2101
2102     if (MPV_common_init(s) < 0)
2103         return -1;
2104     exchange_uv(s); // common init reset pblocks, so we swap them here
2105     s->swap_uv = 1; // in case of xvmc we need to swap uv for each MB
2106     s1->mpeg_enc_ctx_allocated = 1;
2107
2108     for (i = 0; i < 64; i++) {
2109         int j = s->dsp.idct_permutation[i];
2110         v = ff_mpeg1_default_intra_matrix[i];
2111         s->intra_matrix[j]        = v;
2112         s->chroma_intra_matrix[j] = v;
2113
2114         v = ff_mpeg1_default_non_intra_matrix[i];
2115         s->inter_matrix[j]        = v;
2116         s->chroma_inter_matrix[j] = v;
2117     }
2118
2119     s->progressive_sequence  = 1;
2120     s->progressive_frame     = 1;
2121     s->picture_structure     = PICT_FRAME;
2122     s->frame_pred_frame_dct  = 1;
2123     s->chroma_format         = 1;
2124     s->codec_id              = s->avctx->codec_id = CODEC_ID_MPEG2VIDEO;
2125     avctx->sub_id            = 2; /* indicates MPEG-2 */
2126     s1->save_width           = s->width;
2127     s1->save_height          = s->height;
2128     s1->save_progressive_seq = s->progressive_sequence;
2129     return 0;
2130 }
2131
2132
2133 static void mpeg_decode_user_data(AVCodecContext *avctx,
2134                                   const uint8_t *p, int buf_size)
2135 {
2136     const uint8_t *buf_end = p + buf_size;
2137
2138     /* we parse the DTG active format information */
2139     if (buf_end - p >= 5 &&
2140         p[0] == 'D' && p[1] == 'T' && p[2] == 'G' && p[3] == '1') {
2141         int flags = p[4];
2142         p += 5;
2143         if (flags & 0x80) {
2144             /* skip event id */
2145             p += 2;
2146         }
2147         if (flags & 0x40) {
2148             if (buf_end - p < 1)
2149                 return;
2150             avctx->dtg_active_format = p[0] & 0x0f;
2151         }
2152     }
2153 }
2154
2155 static void mpeg_decode_gop(AVCodecContext *avctx,
2156                             const uint8_t *buf, int buf_size)
2157 {
2158     Mpeg1Context *s1  = avctx->priv_data;
2159     MpegEncContext *s = &s1->mpeg_enc_ctx;
2160
2161     int time_code_hours, time_code_minutes;
2162     int time_code_seconds, time_code_pictures;
2163     int broken_link;
2164
2165     init_get_bits(&s->gb, buf, buf_size*8);
2166
2167     skip_bits1(&s->gb); /* drop_frame_flag */
2168
2169     time_code_hours   = get_bits(&s->gb, 5);
2170     time_code_minutes = get_bits(&s->gb, 6);
2171     skip_bits1(&s->gb); // marker bit
2172     time_code_seconds  = get_bits(&s->gb, 6);
2173     time_code_pictures = get_bits(&s->gb, 6);
2174
2175     s->closed_gop = get_bits1(&s->gb);
2176     /*broken_link indicate that after editing the
2177       reference frames of the first B-Frames after GOP I-Frame
2178       are missing (open gop)*/
2179     broken_link = get_bits1(&s->gb);
2180
2181     if (s->avctx->debug & FF_DEBUG_PICT_INFO)
2182         av_log(s->avctx, AV_LOG_DEBUG, "GOP (%2d:%02d:%02d.[%02d]) closed_gop=%d broken_link=%d\n",
2183                time_code_hours, time_code_minutes, time_code_seconds,
2184                time_code_pictures, s->closed_gop, broken_link);
2185 }
2186 /**
2187  * Find the end of the current frame in the bitstream.
2188  * @return the position of the first byte of the next frame, or -1
2189  */
2190 int ff_mpeg1_find_frame_end(ParseContext *pc, const uint8_t *buf, int buf_size, AVCodecParserContext *s)
2191 {
2192     int i;
2193     uint32_t state = pc->state;
2194
2195     /* EOF considered as end of frame */
2196     if (buf_size == 0)
2197         return 0;
2198
2199 /*
2200  0  frame start         -> 1/4
2201  1  first_SEQEXT        -> 0/2
2202  2  first field start   -> 3/0
2203  3  second_SEQEXT       -> 2/0
2204  4  searching end
2205 */
2206
2207     for (i = 0; i < buf_size; i++) {
2208         assert(pc->frame_start_found >= 0 && pc->frame_start_found <= 4);
2209         if (pc->frame_start_found & 1) {
2210             if (state == EXT_START_CODE && (buf[i] & 0xF0) != 0x80)
2211                 pc->frame_start_found--;
2212             else if (state == EXT_START_CODE + 2) {
2213                 if ((buf[i] & 3) == 3)
2214                     pc->frame_start_found = 0;
2215                 else
2216                     pc->frame_start_found = (pc->frame_start_found + 1) & 3;
2217             }
2218             state++;
2219         } else {
2220             i = ff_find_start_code(buf + i, buf + buf_size, &state) - buf - 1;
2221             if (pc->frame_start_found == 0 && state >= SLICE_MIN_START_CODE && state <= SLICE_MAX_START_CODE) {
2222                 i++;
2223                 pc->frame_start_found = 4;
2224             }
2225             if (state == SEQ_END_CODE) {
2226                 pc->state=-1;
2227                 return i+1;
2228             }
2229             if (pc->frame_start_found == 2 && state == SEQ_START_CODE)
2230                 pc->frame_start_found = 0;
2231             if (pc->frame_start_found  < 4 && state == EXT_START_CODE)
2232                 pc->frame_start_found++;
2233             if (pc->frame_start_found == 4 && (state & 0xFFFFFF00) == 0x100) {
2234                 if (state < SLICE_MIN_START_CODE || state > SLICE_MAX_START_CODE) {
2235                     pc->frame_start_found = 0;
2236                     pc->state             = -1;
2237                     return i - 3;
2238                 }
2239             }
2240             if (pc->frame_start_found == 0 && s && state == PICTURE_START_CODE) {
2241                 ff_fetch_timestamp(s, i - 3, 1);
2242             }
2243         }
2244     }
2245     pc->state = state;
2246     return END_NOT_FOUND;
2247 }
2248
2249 static int decode_chunks(AVCodecContext *avctx,
2250                          AVFrame *picture, int *data_size,
2251                          const uint8_t *buf, int buf_size);
2252
2253 /* handle buffering and image synchronisation */
2254 static int mpeg_decode_frame(AVCodecContext *avctx,
2255                              void *data, int *data_size,
2256                              AVPacket *avpkt)
2257 {
2258     const uint8_t *buf = avpkt->data;
2259     int buf_size = avpkt->size;
2260     Mpeg1Context *s = avctx->priv_data;
2261     AVFrame *picture = data;
2262     MpegEncContext *s2 = &s->mpeg_enc_ctx;
2263     av_dlog(avctx, "fill_buffer\n");
2264
2265     if (buf_size == 0 || (buf_size == 4 && AV_RB32(buf) == SEQ_END_CODE)) {
2266         /* special case for last picture */
2267         if (s2->low_delay == 0 && s2->next_picture_ptr) {
2268             *picture = *(AVFrame*)s2->next_picture_ptr;
2269             s2->next_picture_ptr = NULL;
2270
2271             *data_size = sizeof(AVFrame);
2272         }
2273         return buf_size;
2274     }
2275
2276     if (s2->flags & CODEC_FLAG_TRUNCATED) {
2277         int next = ff_mpeg1_find_frame_end(&s2->parse_context, buf, buf_size, NULL);
2278
2279         if (ff_combine_frame(&s2->parse_context, next, (const uint8_t **)&buf, &buf_size) < 0)
2280             return buf_size;
2281     }
2282
2283     if (s->mpeg_enc_ctx_allocated == 0 && avctx->codec_tag == AV_RL32("VCR2"))
2284         vcr2_init_sequence(avctx);
2285
2286     s->slice_count = 0;
2287
2288     if (avctx->extradata && !avctx->frame_number) {
2289         int ret = decode_chunks(avctx, picture, data_size, avctx->extradata, avctx->extradata_size);
2290         if (ret < 0 && avctx->error_recognition >= FF_ER_EXPLODE)
2291             return ret;
2292     }
2293
2294     return decode_chunks(avctx, picture, data_size, buf, buf_size);
2295 }
2296
2297 static int decode_chunks(AVCodecContext *avctx,
2298                          AVFrame *picture, int *data_size,
2299                          const uint8_t *buf, int buf_size)
2300 {
2301     Mpeg1Context *s = avctx->priv_data;
2302     MpegEncContext *s2 = &s->mpeg_enc_ctx;
2303     const uint8_t *buf_ptr = buf;
2304     const uint8_t *buf_end = buf + buf_size;
2305     int ret, input_size;
2306     int last_code = 0;
2307
2308     for (;;) {
2309         /* find next start code */
2310         uint32_t start_code = -1;
2311         buf_ptr = ff_find_start_code(buf_ptr, buf_end, &start_code);
2312         if (start_code > 0x1ff) {
2313             if (s2->pict_type != AV_PICTURE_TYPE_B || avctx->skip_frame <= AVDISCARD_DEFAULT) {
2314                 if (HAVE_THREADS && (avctx->active_thread_type & FF_THREAD_SLICE)) {
2315                     int i;
2316                     av_assert0(avctx->thread_count > 1);
2317
2318                     avctx->execute(avctx, slice_decode_thread,  &s2->thread_context[0], NULL, s->slice_count, sizeof(void*));
2319                     for (i = 0; i < s->slice_count; i++)
2320                         s2->error_count += s2->thread_context[i]->error_count;
2321                 }
2322
2323                 if (CONFIG_VDPAU && uses_vdpau(avctx))
2324                     ff_vdpau_mpeg_picture_complete(s2, buf, buf_size, s->slice_count);
2325
2326                 if (slice_end(avctx, picture)) {
2327                     if (s2->last_picture_ptr || s2->low_delay) //FIXME merge with the stuff in mpeg_decode_slice
2328                         *data_size = sizeof(AVPicture);
2329                 }
2330             }
2331             s2->pict_type = 0;
2332             return FFMAX(0, buf_ptr - buf - s2->parse_context.last_index);
2333         }
2334
2335         input_size = buf_end - buf_ptr;
2336
2337         if (avctx->debug & FF_DEBUG_STARTCODE) {
2338             av_log(avctx, AV_LOG_DEBUG, "%3X at %td left %d\n", start_code, buf_ptr-buf, input_size);
2339         }
2340
2341         /* prepare data for next start code */
2342         switch (start_code) {
2343         case SEQ_START_CODE:
2344             if (last_code == 0) {
2345                 mpeg1_decode_sequence(avctx, buf_ptr, input_size);
2346                 s->sync=1;
2347             } else {
2348                 av_log(avctx, AV_LOG_ERROR, "ignoring SEQ_START_CODE after %X\n", last_code);
2349                 if (avctx->error_recognition >= FF_ER_EXPLODE)
2350                     return AVERROR_INVALIDDATA;
2351             }
2352             break;
2353
2354         case PICTURE_START_CODE:
2355             if (HAVE_THREADS && (avctx->active_thread_type & FF_THREAD_SLICE) && s->slice_count) {
2356                 int i;
2357
2358                 avctx->execute(avctx, slice_decode_thread,
2359                                s2->thread_context, NULL,
2360                                s->slice_count, sizeof(void*));
2361                 for (i = 0; i < s->slice_count; i++)
2362                     s2->error_count += s2->thread_context[i]->error_count;
2363                 s->slice_count = 0;
2364             }
2365             if (last_code == 0 || last_code == SLICE_MIN_START_CODE) {
2366                 ret = mpeg_decode_postinit(avctx);
2367                 if (ret < 0) {
2368                     av_log(avctx, AV_LOG_ERROR, "mpeg_decode_postinit() failure\n");
2369                     return ret;
2370                 }
2371
2372                 /* we have a complete image: we try to decompress it */
2373                 if (mpeg1_decode_picture(avctx, buf_ptr, input_size) < 0)
2374                     s2->pict_type = 0;
2375                 s2->first_slice = 1;
2376                 last_code = PICTURE_START_CODE;
2377             } else {
2378                 av_log(avctx, AV_LOG_ERROR, "ignoring pic after %X\n", last_code);
2379                 if (avctx->error_recognition >= FF_ER_EXPLODE)
2380                     return AVERROR_INVALIDDATA;
2381             }
2382             break;
2383         case EXT_START_CODE:
2384             init_get_bits(&s2->gb, buf_ptr, input_size*8);
2385
2386             switch (get_bits(&s2->gb, 4)) {
2387             case 0x1:
2388                 if (last_code == 0) {
2389                 mpeg_decode_sequence_extension(s);
2390                 } else {
2391                     av_log(avctx, AV_LOG_ERROR, "ignoring seq ext after %X\n", last_code);
2392                     if (avctx->error_recognition >= FF_ER_EXPLODE)
2393                         return AVERROR_INVALIDDATA;
2394                 }
2395                 break;
2396             case 0x2:
2397                 mpeg_decode_sequence_display_extension(s);
2398                 break;
2399             case 0x3:
2400                 mpeg_decode_quant_matrix_extension(s2);
2401                 break;
2402             case 0x7:
2403                 mpeg_decode_picture_display_extension(s);
2404                 break;
2405             case 0x8:
2406                 if (last_code == PICTURE_START_CODE) {
2407                     mpeg_decode_picture_coding_extension(s);
2408                 } else {
2409                     av_log(avctx, AV_LOG_ERROR, "ignoring pic cod ext after %X\n", last_code);
2410                     if (avctx->error_recognition >= FF_ER_EXPLODE)
2411                         return AVERROR_INVALIDDATA;
2412                 }
2413                 break;
2414             }
2415             break;
2416         case USER_START_CODE:
2417             mpeg_decode_user_data(avctx, buf_ptr, input_size);
2418             break;
2419         case GOP_START_CODE:
2420             if (last_code == 0) {
2421                 s2->first_field=0;
2422                 mpeg_decode_gop(avctx, buf_ptr, input_size);
2423                 s->sync=1;
2424             } else {
2425                 av_log(avctx, AV_LOG_ERROR, "ignoring GOP_START_CODE after %X\n", last_code);
2426                 if (avctx->error_recognition >= FF_ER_EXPLODE)
2427                     return AVERROR_INVALIDDATA;
2428             }
2429             break;
2430         default:
2431             if (start_code >= SLICE_MIN_START_CODE &&
2432                 start_code <= SLICE_MAX_START_CODE && last_code != 0) {
2433                 const int field_pic = s2->picture_structure != PICT_FRAME;
2434                 int mb_y = (start_code - SLICE_MIN_START_CODE) << field_pic;
2435                 last_code = SLICE_MIN_START_CODE;
2436
2437                 if (s2->picture_structure == PICT_BOTTOM_FIELD)
2438                     mb_y++;
2439
2440                 if (mb_y >= s2->mb_height) {
2441                     av_log(s2->avctx, AV_LOG_ERROR, "slice below image (%d >= %d)\n", mb_y, s2->mb_height);
2442                     return -1;
2443                 }
2444
2445                 if (s2->last_picture_ptr == NULL) {
2446                 /* Skip B-frames if we do not have reference frames and gop is not closed */
2447                     if (s2->pict_type == AV_PICTURE_TYPE_B) {
2448                         if (!s2->closed_gop)
2449                             break;
2450                     }
2451                 }
2452                 if (s2->pict_type == AV_PICTURE_TYPE_I)
2453                     s->sync=1;
2454                 if (s2->next_picture_ptr == NULL) {
2455                 /* Skip P-frames if we do not have a reference frame or we have an invalid header. */
2456                     if (s2->pict_type == AV_PICTURE_TYPE_P && !s->sync) break;
2457                 }
2458                 if ((avctx->skip_frame >= AVDISCARD_NONREF && s2->pict_type == AV_PICTURE_TYPE_B) ||
2459                     (avctx->skip_frame >= AVDISCARD_NONKEY && s2->pict_type != AV_PICTURE_TYPE_I) ||
2460                      avctx->skip_frame >= AVDISCARD_ALL)
2461                     break;
2462
2463                 if (!s->mpeg_enc_ctx_allocated)
2464                     break;
2465
2466                 if (s2->codec_id == CODEC_ID_MPEG2VIDEO) {
2467                     if (mb_y < avctx->skip_top || mb_y >= s2->mb_height - avctx->skip_bottom)
2468                         break;
2469                 }
2470
2471                 if (!s2->pict_type) {
2472                     av_log(avctx, AV_LOG_ERROR, "Missing picture start code\n");
2473                     if (avctx->error_recognition >= FF_ER_EXPLODE)
2474                         return AVERROR_INVALIDDATA;
2475                     break;
2476                 }
2477
2478                 if (s2->first_slice) {
2479                     s2->first_slice = 0;
2480                     if (mpeg_field_start(s2, buf, buf_size) < 0)
2481                         return -1;
2482                 }
2483                 if (!s2->current_picture_ptr) {
2484                     av_log(avctx, AV_LOG_ERROR, "current_picture not initialized\n");
2485                     return AVERROR_INVALIDDATA;
2486                 }
2487
2488                 if (uses_vdpau(avctx)) {
2489                     s->slice_count++;
2490                     break;
2491                 }
2492
2493                 if (HAVE_THREADS && (avctx->active_thread_type & FF_THREAD_SLICE)) {
2494                     int threshold= (s2->mb_height*s->slice_count + avctx->thread_count/2) / avctx->thread_count;
2495                     av_assert0(avctx->thread_count > 1);
2496                     if (threshold <= mb_y) {
2497                         MpegEncContext *thread_context = s2->thread_context[s->slice_count];
2498
2499                         thread_context->start_mb_y = mb_y;
2500                         thread_context->end_mb_y   = s2->mb_height;
2501                         if (s->slice_count) {
2502                             s2->thread_context[s->slice_count-1]->end_mb_y = mb_y;
2503                             ff_update_duplicate_context(thread_context, s2);
2504                         }
2505                         init_get_bits(&thread_context->gb, buf_ptr, input_size*8);
2506                         s->slice_count++;
2507                     }
2508                     buf_ptr += 2; // FIXME add minimum number of bytes per slice
2509                 } else {
2510                     ret = mpeg_decode_slice(s, mb_y, &buf_ptr, input_size);
2511                     emms_c();
2512
2513                     if (ret < 0) {
2514                         if (avctx->error_recognition >= FF_ER_EXPLODE)
2515                             return ret;
2516                         if (s2->resync_mb_x >= 0 && s2->resync_mb_y >= 0)
2517                             ff_er_add_slice(s2, s2->resync_mb_x, s2->resync_mb_y, s2->mb_x, s2->mb_y, AC_ERROR | DC_ERROR | MV_ERROR);
2518                     } else {
2519                         ff_er_add_slice(s2, s2->resync_mb_x, s2->resync_mb_y, s2->mb_x-1, s2->mb_y, AC_END | DC_END | MV_END);
2520                     }
2521                 }
2522             }
2523             break;
2524         }
2525     }
2526 }
2527
2528 static void flush(AVCodecContext *avctx)
2529 {
2530     Mpeg1Context *s = avctx->priv_data;
2531
2532     s->sync=0;
2533
2534     ff_mpeg_flush(avctx);
2535 }
2536
2537 static int mpeg_decode_end(AVCodecContext *avctx)
2538 {
2539     Mpeg1Context *s = avctx->priv_data;
2540
2541     if (s->mpeg_enc_ctx_allocated)
2542         MPV_common_end(&s->mpeg_enc_ctx);
2543     return 0;
2544 }
2545
2546 static const AVProfile mpeg2_video_profiles[] = {
2547     { FF_PROFILE_MPEG2_422,          "4:2:2"              },
2548     { FF_PROFILE_MPEG2_HIGH,         "High"               },
2549     { FF_PROFILE_MPEG2_SS,           "Spatially Scalable" },
2550     { FF_PROFILE_MPEG2_SNR_SCALABLE, "SNR Scalable"       },
2551     { FF_PROFILE_MPEG2_MAIN,         "Main"               },
2552     { FF_PROFILE_MPEG2_SIMPLE,       "Simple"             },
2553     { FF_PROFILE_RESERVED,           "Reserved"           },
2554     { FF_PROFILE_RESERVED,           "Reserved"           },
2555     { FF_PROFILE_UNKNOWN },
2556 };
2557
2558
2559 AVCodec ff_mpeg1video_decoder = {
2560     .name           = "mpeg1video",
2561     .type           = AVMEDIA_TYPE_VIDEO,
2562     .id             = CODEC_ID_MPEG1VIDEO,
2563     .priv_data_size = sizeof(Mpeg1Context),
2564     .init           = mpeg_decode_init,
2565     .close          = mpeg_decode_end,
2566     .decode         = mpeg_decode_frame,
2567     .capabilities   = CODEC_CAP_DRAW_HORIZ_BAND | CODEC_CAP_DR1 | CODEC_CAP_TRUNCATED | CODEC_CAP_DELAY | CODEC_CAP_SLICE_THREADS,
2568     .flush          = flush,
2569     .max_lowres     = 3,
2570     .long_name      = NULL_IF_CONFIG_SMALL("MPEG-1 video"),
2571     .update_thread_context = ONLY_IF_THREADS_ENABLED(mpeg_decode_update_thread_context)
2572 };
2573
2574 AVCodec ff_mpeg2video_decoder = {
2575     .name           = "mpeg2video",
2576     .type           = AVMEDIA_TYPE_VIDEO,
2577     .id             = CODEC_ID_MPEG2VIDEO,
2578     .priv_data_size = sizeof(Mpeg1Context),
2579     .init           = mpeg_decode_init,
2580     .close          = mpeg_decode_end,
2581     .decode         = mpeg_decode_frame,
2582     .capabilities   = CODEC_CAP_DRAW_HORIZ_BAND | CODEC_CAP_DR1 | CODEC_CAP_TRUNCATED | CODEC_CAP_DELAY | CODEC_CAP_SLICE_THREADS,
2583     .flush          = flush,
2584     .max_lowres     = 3,
2585     .long_name      = NULL_IF_CONFIG_SMALL("MPEG-2 video"),
2586     .profiles       = NULL_IF_CONFIG_SMALL(mpeg2_video_profiles),
2587 };
2588
2589 //legacy decoder
2590 AVCodec ff_mpegvideo_decoder = {
2591     .name           = "mpegvideo",
2592     .type           = AVMEDIA_TYPE_VIDEO,
2593     .id             = CODEC_ID_MPEG2VIDEO,
2594     .priv_data_size = sizeof(Mpeg1Context),
2595     .init           = mpeg_decode_init,
2596     .close          = mpeg_decode_end,
2597     .decode         = mpeg_decode_frame,
2598     .capabilities   = CODEC_CAP_DRAW_HORIZ_BAND | CODEC_CAP_DR1 | CODEC_CAP_TRUNCATED | CODEC_CAP_DELAY | CODEC_CAP_SLICE_THREADS,
2599     .flush          = flush,
2600     .max_lowres     = 3,
2601     .long_name      = NULL_IF_CONFIG_SMALL("MPEG-1 video"),
2602 };
2603
2604 #if CONFIG_MPEG_XVMC_DECODER
2605 static av_cold int mpeg_mc_decode_init(AVCodecContext *avctx)
2606 {
2607     if (avctx->active_thread_type & FF_THREAD_SLICE)
2608         return -1;
2609     if (!(avctx->slice_flags & SLICE_FLAG_CODED_ORDER))
2610         return -1;
2611     if (!(avctx->slice_flags & SLICE_FLAG_ALLOW_FIELD)) {
2612         av_dlog(avctx, "mpeg12.c: XvMC decoder will work better if SLICE_FLAG_ALLOW_FIELD is set\n");
2613     }
2614     mpeg_decode_init(avctx);
2615
2616     avctx->pix_fmt           = PIX_FMT_XVMC_MPEG2_IDCT;
2617     avctx->xvmc_acceleration = 2; // 2 - the blocks are packed!
2618
2619     return 0;
2620 }
2621
2622 AVCodec ff_mpeg_xvmc_decoder = {
2623     .name           = "mpegvideo_xvmc",
2624     .type           = AVMEDIA_TYPE_VIDEO,
2625     .id             = CODEC_ID_MPEG2VIDEO_XVMC,
2626     .priv_data_size = sizeof(Mpeg1Context),
2627     .init           = mpeg_mc_decode_init,
2628     .close          = mpeg_decode_end,
2629     .decode         = mpeg_decode_frame,
2630     .capabilities   = CODEC_CAP_DRAW_HORIZ_BAND | CODEC_CAP_DR1 | CODEC_CAP_TRUNCATED| CODEC_CAP_HWACCEL | CODEC_CAP_DELAY,
2631     .flush          = flush,
2632     .long_name      = NULL_IF_CONFIG_SMALL("MPEG-1/2 video XvMC (X-Video Motion Compensation)"),
2633 };
2634
2635 #endif
2636
2637 #if CONFIG_MPEG_VDPAU_DECODER
2638 AVCodec ff_mpeg_vdpau_decoder = {
2639     .name           = "mpegvideo_vdpau",
2640     .type           = AVMEDIA_TYPE_VIDEO,
2641     .id             = CODEC_ID_MPEG2VIDEO,
2642     .priv_data_size = sizeof(Mpeg1Context),
2643     .init           = mpeg_decode_init,
2644     .close          = mpeg_decode_end,
2645     .decode         = mpeg_decode_frame,
2646     .capabilities   = CODEC_CAP_DR1 | CODEC_CAP_TRUNCATED | CODEC_CAP_HWACCEL_VDPAU | CODEC_CAP_DELAY,
2647     .flush          = flush,
2648     .long_name      = NULL_IF_CONFIG_SMALL("MPEG-1/2 video (VDPAU acceleration)"),
2649 };
2650 #endif
2651
2652 #if CONFIG_MPEG1_VDPAU_DECODER
2653 AVCodec ff_mpeg1_vdpau_decoder = {
2654     .name           = "mpeg1video_vdpau",
2655     .type           = AVMEDIA_TYPE_VIDEO,
2656     .id             = CODEC_ID_MPEG1VIDEO,
2657     .priv_data_size = sizeof(Mpeg1Context),
2658     .init           = mpeg_decode_init,
2659     .close          = mpeg_decode_end,
2660     .decode         = mpeg_decode_frame,
2661     .capabilities   = CODEC_CAP_DR1 | CODEC_CAP_TRUNCATED | CODEC_CAP_HWACCEL_VDPAU | CODEC_CAP_DELAY,
2662     .flush          = flush,
2663     .long_name      = NULL_IF_CONFIG_SMALL("MPEG-1 video (VDPAU acceleration)"),
2664 };
2665 #endif
2666