git.sesse.net Git - ffmpeg/blob - libavcodec/mpeg12.c

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