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