libavcodec/mpeg12.c

   1 /*
   2  * MPEG1 codec / MPEG2 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 mpeg12.c
  25  * MPEG1/2 codec
  26  */
  27
  28 //#define DEBUG
  29 #include "avcodec.h"
  30 #include "dsputil.h"
  31 #include "mpegvideo.h"
  32
  33 #include "mpeg12data.h"
  34 #include "bytestream.h"
  35
  36 //#undef NDEBUG
  37 //#include <assert.h>
  38
  39
  40 /* Start codes. */
  41 #define SEQ_END_CODE            0x000001b7
  42 #define SEQ_START_CODE          0x000001b3
  43 #define GOP_START_CODE          0x000001b8
  44 #define PICTURE_START_CODE      0x00000100
  45 #define SLICE_MIN_START_CODE    0x00000101
  46 #define SLICE_MAX_START_CODE    0x000001af
  47 #define EXT_START_CODE          0x000001b5
  48 #define USER_START_CODE         0x000001b2
  49
  50 #define DC_VLC_BITS 9
  51 #define MV_VLC_BITS 9
  52 #define MBINCR_VLC_BITS 9
  53 #define MB_PAT_VLC_BITS 9
  54 #define MB_PTYPE_VLC_BITS 6
  55 #define MB_BTYPE_VLC_BITS 6
  56 #define TEX_VLC_BITS 9
  57
  58 #ifdef CONFIG_ENCODERS
  59 static void mpeg1_encode_block(MpegEncContext *s,
  60                          DCTELEM *block,
  61                          int component);
  62 static void mpeg1_encode_motion(MpegEncContext *s, int val, int f_or_b_code);    // RAL: f_code parameter added
  63 #endif //CONFIG_ENCODERS
  64 static inline int mpeg1_decode_block_inter(MpegEncContext *s,
  65                               DCTELEM *block,
  66                               int n);
  67 static inline int mpeg1_decode_block_intra(MpegEncContext *s,
  68                               DCTELEM *block,
  69                               int n);
  70 static inline int mpeg1_fast_decode_block_inter(MpegEncContext *s, DCTELEM *block, int n);
  71 static inline int mpeg2_decode_block_non_intra(MpegEncContext *s,
  72                                         DCTELEM *block,
  73                                         int n);
  74 static inline int mpeg2_decode_block_intra(MpegEncContext *s,
  75                                     DCTELEM *block,
  76                                     int n);
  77 static inline int mpeg2_fast_decode_block_non_intra(MpegEncContext *s, DCTELEM *block, int n);
  78 static inline int mpeg2_fast_decode_block_intra(MpegEncContext *s, DCTELEM *block, int n);
  79 static int mpeg_decode_motion(MpegEncContext *s, int fcode, int pred);
  80 static void exchange_uv(MpegEncContext *s);
  81
  82 #ifdef HAVE_XVMC
  83 extern int XVMC_field_start(MpegEncContext *s, AVCodecContext *avctx);
  84 extern int XVMC_field_end(MpegEncContext *s);
  85 extern void XVMC_pack_pblocks(MpegEncContext *s,int cbp);
  86 extern void XVMC_init_block(MpegEncContext *s);//set s->block
  87 #endif
  88
  89 static const enum PixelFormat pixfmt_yuv_420[]= {PIX_FMT_YUV420P,-1};
  90 static const enum PixelFormat pixfmt_yuv_422[]= {PIX_FMT_YUV422P,-1};
  91 static const enum PixelFormat pixfmt_yuv_444[]= {PIX_FMT_YUV444P,-1};
  92 static const enum PixelFormat pixfmt_xvmc_mpg2_420[] = {
  93                                            PIX_FMT_XVMC_MPEG2_IDCT,
  94                                            PIX_FMT_XVMC_MPEG2_MC,
  95                                            -1};
  96 #ifdef CONFIG_ENCODERS
  97 static uint8_t mv_penalty[MAX_FCODE+1][MAX_MV*2+1];
  98 static uint8_t fcode_tab[MAX_MV*2+1];
  99
 100 static uint8_t  uni_mpeg1_ac_vlc_len [64*64*2];
 101 static uint8_t  uni_mpeg2_ac_vlc_len [64*64*2];
 102
 103 /* simple include everything table for dc, first byte is bits number next 3 are code*/
 104 static uint32_t mpeg1_lum_dc_uni[512];
 105 static uint32_t mpeg1_chr_dc_uni[512];
 106
 107 static uint8_t mpeg1_index_run[2][64];
 108 static int8_t mpeg1_max_level[2][64];
 109 #endif //CONFIG_ENCODERS
 110
 111 static uint8_t static_rl_table_store[2][2][2*MAX_RUN + MAX_LEVEL + 3];
 112
 113 static void init_2d_vlc_rl(RLTable *rl, int use_static)
 114 {
 115     int i;
 116
 117     init_vlc(&rl->vlc, TEX_VLC_BITS, rl->n + 2,
 118              &rl->table_vlc[0][1], 4, 2,
 119              &rl->table_vlc[0][0], 4, 2, use_static);
 120
 121     if(use_static)
 122         rl->rl_vlc[0]= av_mallocz_static(rl->vlc.table_size*sizeof(RL_VLC_ELEM));
 123     else
 124         rl->rl_vlc[0]= av_malloc(rl->vlc.table_size*sizeof(RL_VLC_ELEM));
 125
 126     for(i=0; i<rl->vlc.table_size; i++){
 127         int code= rl->vlc.table[i][0];
 128         int len = rl->vlc.table[i][1];
 129         int level, run;
 130
 131         if(len==0){ // illegal code
 132             run= 65;
 133             level= MAX_LEVEL;
 134         }else if(len<0){ //more bits needed
 135             run= 0;
 136             level= code;
 137         }else{
 138             if(code==rl->n){ //esc
 139                 run= 65;
 140                 level= 0;
 141             }else if(code==rl->n+1){ //eob
 142                 run= 0;
 143                 level= 127;
 144             }else{
 145                 run=   rl->table_run  [code] + 1;
 146                 level= rl->table_level[code];
 147             }
 148         }
 149         rl->rl_vlc[0][i].len= len;
 150         rl->rl_vlc[0][i].level= level;
 151         rl->rl_vlc[0][i].run= run;
 152     }
 153 }
 154
 155 #ifdef CONFIG_ENCODERS
 156 static void init_uni_ac_vlc(RLTable *rl, uint8_t *uni_ac_vlc_len){
 157     int i;
 158
 159     for(i=0; i<128; i++){
 160         int level= i-64;
 161         int run;
 162         for(run=0; run<64; run++){
 163             int len, bits, code;
 164
 165             int alevel= FFABS(level);
 166             int sign= (level>>31)&1;
 167
 168             if (alevel > rl->max_level[0][run])
 169                 code= 111; /*rl->n*/
 170             else
 171                 code= rl->index_run[0][run] + alevel - 1;
 172
 173             if (code < 111 /* rl->n */) {
 174                 /* store the vlc & sign at once */
 175                 len=   rl->table_vlc[code][1]+1;
 176                 bits= (rl->table_vlc[code][0]<<1) + sign;
 177             } else {
 178                 len=  rl->table_vlc[111/*rl->n*/][1]+6;
 179                 bits= rl->table_vlc[111/*rl->n*/][0]<<6;
 180
 181                 bits|= run;
 182                 if (alevel < 128) {
 183                     bits<<=8; len+=8;
 184                     bits|= level & 0xff;
 185                 } else {
 186                     bits<<=16; len+=16;
 187                     bits|= level & 0xff;
 188                     if (level < 0) {
 189                         bits|= 0x8001 + level + 255;
 190                     } else {
 191                         bits|= level & 0xffff;
 192                     }
 193                 }
 194             }
 195
 196             uni_ac_vlc_len [UNI_AC_ENC_INDEX(run, i)]= len;
 197         }
 198     }
 199 }
 200
 201
 202 static int find_frame_rate_index(MpegEncContext *s){
 203     int i;
 204     int64_t dmin= INT64_MAX;
 205     int64_t d;
 206
 207     for(i=1;i<14;i++) {
 208         int64_t n0= 1001LL/ff_frame_rate_tab[i].den*ff_frame_rate_tab[i].num*s->avctx->time_base.num;
 209         int64_t n1= 1001LL*s->avctx->time_base.den;
 210         if(s->avctx->strict_std_compliance > FF_COMPLIANCE_INOFFICIAL && i>=9) break;
 211
 212         d = FFABS(n0 - n1);
 213         if(d < dmin){
 214             dmin=d;
 215             s->frame_rate_index= i;
 216         }
 217     }
 218     if(dmin)
 219         return -1;
 220     else
 221         return 0;
 222 }
 223
 224 static int encode_init(AVCodecContext *avctx)
 225 {
 226     MpegEncContext *s = avctx->priv_data;
 227
 228     if(MPV_encode_init(avctx) < 0)
 229         return -1;
 230
 231     if(find_frame_rate_index(s) < 0){
 232         if(s->strict_std_compliance > FF_COMPLIANCE_EXPERIMENTAL){
 233             av_log(avctx, AV_LOG_ERROR, "MPEG1/2 does not support %d/%d fps\n", avctx->time_base.den, avctx->time_base.num);
 234             return -1;
 235         }else{
 236             av_log(avctx, AV_LOG_INFO, "MPEG1/2 does not support %d/%d fps, there may be AV sync issues\n", avctx->time_base.den, avctx->time_base.num);
 237         }
 238     }
 239
 240     if(avctx->profile == FF_PROFILE_UNKNOWN){
 241         if(avctx->level != FF_LEVEL_UNKNOWN){
 242             av_log(avctx, AV_LOG_ERROR, "Set profile and level\n");
 243             return -1;
 244         }
 245         avctx->profile = s->chroma_format == CHROMA_420 ? 4 : 0; /* Main or 4:2:2 */
 246     }
 247
 248     if(avctx->level == FF_LEVEL_UNKNOWN){
 249         if(avctx->profile == 0){ /* 4:2:2 */
 250             if(avctx->width <= 720 && avctx->height <= 608) avctx->level = 5; /* Main */
 251             else                                            avctx->level = 2; /* High */
 252         }else{
 253             if(avctx->profile != 1 && s->chroma_format != CHROMA_420){
 254                 av_log(avctx, AV_LOG_ERROR, "Only High(1) and 4:2:2(0) profiles support 4:2:2 color sampling\n");
 255                 return -1;
 256             }
 257             if(avctx->width <= 720 && avctx->height <= 576) avctx->level = 8; /* Main */
 258             else if(avctx->width <= 1440)                   avctx->level = 6; /* High 1440 */
 259             else                                            avctx->level = 4; /* High */
 260         }
 261     }
 262
 263     if((avctx->flags2 & CODEC_FLAG2_DROP_FRAME_TIMECODE) && s->frame_rate_index != 4){
 264         av_log(avctx, AV_LOG_ERROR, "Drop frame time code only allowed with 1001/30000 fps\n");
 265         return -1;
 266     }
 267
 268     return 0;
 269 }
 270
 271 static void put_header(MpegEncContext *s, int header)
 272 {
 273     align_put_bits(&s->pb);
 274     put_bits(&s->pb, 16, header>>16);
 275     put_bits(&s->pb, 16, header&0xFFFF);
 276 }
 277
 278 /* put sequence header if needed */
 279 static void mpeg1_encode_sequence_header(MpegEncContext *s)
 280 {
 281         unsigned int vbv_buffer_size;
 282         unsigned int fps, v;
 283         int i;
 284         uint64_t time_code;
 285         float best_aspect_error= 1E10;
 286         float aspect_ratio= av_q2d(s->avctx->sample_aspect_ratio);
 287         int constraint_parameter_flag;
 288
 289         if(aspect_ratio==0.0) aspect_ratio= 1.0; //pixel aspect 1:1 (VGA)
 290
 291         if (s->current_picture.key_frame) {
 292             AVRational framerate= ff_frame_rate_tab[s->frame_rate_index];
 293
 294             /* mpeg1 header repeated every gop */
 295             put_header(s, SEQ_START_CODE);
 296
 297             put_bits(&s->pb, 12, s->width);
 298             put_bits(&s->pb, 12, s->height);
 299
 300             for(i=1; i<15; i++){
 301                 float error= aspect_ratio;
 302                 if(s->codec_id == CODEC_ID_MPEG1VIDEO || i <=1)
 303                     error-= 1.0/mpeg1_aspect[i];
 304                 else
 305                     error-= av_q2d(mpeg2_aspect[i])*s->height/s->width;
 306
 307                 error= FFABS(error);
 308
 309                 if(error < best_aspect_error){
 310                     best_aspect_error= error;
 311                     s->aspect_ratio_info= i;
 312                 }
 313             }
 314
 315             put_bits(&s->pb, 4, s->aspect_ratio_info);
 316             put_bits(&s->pb, 4, s->frame_rate_index);
 317
 318             if(s->avctx->rc_max_rate){
 319                 v = (s->avctx->rc_max_rate + 399) / 400;
 320                 if (v > 0x3ffff && s->codec_id == CODEC_ID_MPEG1VIDEO)
 321                     v = 0x3ffff;
 322             }else{
 323                 v= 0x3FFFF;
 324             }
 325
 326             if(s->avctx->rc_buffer_size)
 327                 vbv_buffer_size = s->avctx->rc_buffer_size;
 328             else
 329                 /* VBV calculation: Scaled so that a VCD has the proper VBV size of 40 kilobytes */
 330                 vbv_buffer_size = (( 20 * s->bit_rate) / (1151929 / 2)) * 8 * 1024;
 331             vbv_buffer_size= (vbv_buffer_size + 16383) / 16384;
 332
 333             put_bits(&s->pb, 18, v & 0x3FFFF);
 334             put_bits(&s->pb, 1, 1); /* marker */
 335             put_bits(&s->pb, 10, vbv_buffer_size & 0x3FF);
 336
 337             constraint_parameter_flag=
 338                 s->width <= 768 && s->height <= 576 &&
 339                 s->mb_width * s->mb_height <= 396 &&
 340                 s->mb_width * s->mb_height * framerate.num <= framerate.den*396*25 &&
 341                 framerate.num <= framerate.den*30 &&
 342                 s->avctx->me_range && s->avctx->me_range < 128 &&
 343                 vbv_buffer_size <= 20 &&
 344                 v <= 1856000/400 &&
 345                 s->codec_id == CODEC_ID_MPEG1VIDEO;
 346
 347             put_bits(&s->pb, 1, constraint_parameter_flag);
 348
 349             ff_write_quant_matrix(&s->pb, s->avctx->intra_matrix);
 350             ff_write_quant_matrix(&s->pb, s->avctx->inter_matrix);
 351
 352             if(s->codec_id == CODEC_ID_MPEG2VIDEO){
 353                 put_header(s, EXT_START_CODE);
 354                 put_bits(&s->pb, 4, 1); //seq ext
 355
 356                 put_bits(&s->pb, 1, s->avctx->profile == 0); //escx 1 for 4:2:2 profile */
 357
 358                 put_bits(&s->pb, 3, s->avctx->profile); //profile
 359                 put_bits(&s->pb, 4, s->avctx->level); //level
 360
 361                 put_bits(&s->pb, 1, s->progressive_sequence);
 362                 put_bits(&s->pb, 2, s->chroma_format);
 363                 put_bits(&s->pb, 2, 0); //horizontal size ext
 364                 put_bits(&s->pb, 2, 0); //vertical size ext
 365                 put_bits(&s->pb, 12, v>>18); //bitrate ext
 366                 put_bits(&s->pb, 1, 1); //marker
 367                 put_bits(&s->pb, 8, vbv_buffer_size >>10); //vbv buffer ext
 368                 put_bits(&s->pb, 1, s->low_delay);
 369                 put_bits(&s->pb, 2, 0); // frame_rate_ext_n
 370                 put_bits(&s->pb, 5, 0); // frame_rate_ext_d
 371             }
 372
 373             put_header(s, GOP_START_CODE);
 374             put_bits(&s->pb, 1, !!(s->avctx->flags2 & CODEC_FLAG2_DROP_FRAME_TIMECODE)); /* drop frame flag */
 375             /* time code : we must convert from the real frame rate to a
 376                fake mpeg frame rate in case of low frame rate */
 377             fps = (framerate.num + framerate.den/2)/ framerate.den;
 378             time_code = s->current_picture_ptr->coded_picture_number + s->avctx->timecode_frame_start;
 379
 380             s->gop_picture_number = s->current_picture_ptr->coded_picture_number;
 381             if (s->avctx->flags2 & CODEC_FLAG2_DROP_FRAME_TIMECODE) {
 382                 /* only works for NTSC 29.97 */
 383                 int d = time_code / 17982;
 384                 int m = time_code % 17982;
 385                 //if (m < 2) m += 2; /* not needed since -2,-1 / 1798 in C returns 0 */
 386                 time_code += 18 * d + 2 * ((m - 2) / 1798);
 387             }
 388             put_bits(&s->pb, 5, (uint32_t)((time_code / (fps * 3600)) % 24));
 389             put_bits(&s->pb, 6, (uint32_t)((time_code / (fps * 60)) % 60));
 390             put_bits(&s->pb, 1, 1);
 391             put_bits(&s->pb, 6, (uint32_t)((time_code / fps) % 60));
 392             put_bits(&s->pb, 6, (uint32_t)((time_code % fps)));
 393             put_bits(&s->pb, 1, !!(s->flags & CODEC_FLAG_CLOSED_GOP));
 394             put_bits(&s->pb, 1, 0); /* broken link */
 395         }
 396 }
 397
 398 static inline void encode_mb_skip_run(MpegEncContext *s, int run){
 399     while (run >= 33) {
 400         put_bits(&s->pb, 11, 0x008);
 401         run -= 33;
 402     }
 403     put_bits(&s->pb, mbAddrIncrTable[run][1],
 404              mbAddrIncrTable[run][0]);
 405 }
 406 #endif //CONFIG_ENCODERS
 407
 408 static void common_init(MpegEncContext *s)
 409 {
 410
 411     s->y_dc_scale_table=
 412     s->c_dc_scale_table= mpeg2_dc_scale_table[s->intra_dc_precision];
 413
 414 }
 415
 416 void ff_mpeg1_clean_buffers(MpegEncContext *s){
 417     s->last_dc[0] = 1 << (7 + s->intra_dc_precision);
 418     s->last_dc[1] = s->last_dc[0];
 419     s->last_dc[2] = s->last_dc[0];
 420     memset(s->last_mv, 0, sizeof(s->last_mv));
 421 }
 422
 423 #ifdef CONFIG_ENCODERS
 424
 425 static av_always_inline void put_qscale(MpegEncContext *s)
 426 {
 427     if(s->q_scale_type){
 428         assert(s->qscale>=1 && s->qscale <=12);
 429         put_bits(&s->pb, 5, inv_non_linear_qscale[s->qscale]);
 430     }else{
 431         put_bits(&s->pb, 5, s->qscale);
 432     }
 433 }
 434
 435 void ff_mpeg1_encode_slice_header(MpegEncContext *s){
 436     put_header(s, SLICE_MIN_START_CODE + s->mb_y);
 437     put_qscale(s);
 438     put_bits(&s->pb, 1, 0); /* slice extra information */
 439 }
 440
 441 void mpeg1_encode_picture_header(MpegEncContext *s, int picture_number)
 442 {
 443     mpeg1_encode_sequence_header(s);
 444
 445     /* mpeg1 picture header */
 446     put_header(s, PICTURE_START_CODE);
 447     /* temporal reference */
 448
 449     // RAL: s->picture_number instead of s->fake_picture_number
 450     put_bits(&s->pb, 10, (s->picture_number -
 451                           s->gop_picture_number) & 0x3ff);
 452     put_bits(&s->pb, 3, s->pict_type);
 453
 454     s->vbv_delay_ptr= s->pb.buf + put_bits_count(&s->pb)/8;
 455     put_bits(&s->pb, 16, 0xFFFF); /* vbv_delay */
 456
 457     // RAL: Forward f_code also needed for B frames
 458     if (s->pict_type == P_TYPE || s->pict_type == B_TYPE) {
 459         put_bits(&s->pb, 1, 0); /* half pel coordinates */
 460         if(s->codec_id == CODEC_ID_MPEG1VIDEO)
 461             put_bits(&s->pb, 3, s->f_code); /* forward_f_code */
 462         else
 463             put_bits(&s->pb, 3, 7); /* forward_f_code */
 464     }
 465
 466     // RAL: Backward f_code necessary for B frames
 467     if (s->pict_type == B_TYPE) {
 468         put_bits(&s->pb, 1, 0); /* half pel coordinates */
 469         if(s->codec_id == CODEC_ID_MPEG1VIDEO)
 470             put_bits(&s->pb, 3, s->b_code); /* backward_f_code */
 471         else
 472             put_bits(&s->pb, 3, 7); /* backward_f_code */
 473     }
 474
 475     put_bits(&s->pb, 1, 0); /* extra bit picture */
 476
 477     s->frame_pred_frame_dct = 1;
 478     if(s->codec_id == CODEC_ID_MPEG2VIDEO){
 479         put_header(s, EXT_START_CODE);
 480         put_bits(&s->pb, 4, 8); //pic ext
 481         if (s->pict_type == P_TYPE || s->pict_type == B_TYPE) {
 482             put_bits(&s->pb, 4, s->f_code);
 483             put_bits(&s->pb, 4, s->f_code);
 484         }else{
 485             put_bits(&s->pb, 8, 255);
 486         }
 487         if (s->pict_type == B_TYPE) {
 488             put_bits(&s->pb, 4, s->b_code);
 489             put_bits(&s->pb, 4, s->b_code);
 490         }else{
 491             put_bits(&s->pb, 8, 255);
 492         }
 493         put_bits(&s->pb, 2, s->intra_dc_precision);
 494
 495         assert(s->picture_structure == PICT_FRAME);
 496         put_bits(&s->pb, 2, s->picture_structure);
 497         if (s->progressive_sequence) {
 498             put_bits(&s->pb, 1, 0); /* no repeat */
 499         } else {
 500             put_bits(&s->pb, 1, s->current_picture_ptr->top_field_first);
 501         }
 502         /* XXX: optimize the generation of this flag with entropy
 503            measures */
 504         s->frame_pred_frame_dct = s->progressive_sequence;
 505
 506         put_bits(&s->pb, 1, s->frame_pred_frame_dct);
 507         put_bits(&s->pb, 1, s->concealment_motion_vectors);
 508         put_bits(&s->pb, 1, s->q_scale_type);
 509         put_bits(&s->pb, 1, s->intra_vlc_format);
 510         put_bits(&s->pb, 1, s->alternate_scan);
 511         put_bits(&s->pb, 1, s->repeat_first_field);
 512         s->progressive_frame = s->progressive_sequence;
 513         put_bits(&s->pb, 1, s->chroma_format == CHROMA_420 ? s->progressive_frame : 0); /* chroma_420_type */
 514         put_bits(&s->pb, 1, s->progressive_frame);
 515         put_bits(&s->pb, 1, 0); //composite_display_flag
 516     }
 517     if(s->flags & CODEC_FLAG_SVCD_SCAN_OFFSET){
 518         int i;
 519
 520         put_header(s, USER_START_CODE);
 521         for(i=0; i<sizeof(svcd_scan_offset_placeholder); i++){
 522             put_bits(&s->pb, 8, svcd_scan_offset_placeholder[i]);
 523         }
 524     }
 525
 526     s->mb_y=0;
 527     ff_mpeg1_encode_slice_header(s);
 528 }
 529
 530 static inline void put_mb_modes(MpegEncContext *s, int n, int bits,
 531                                 int has_mv, int field_motion)
 532 {
 533     put_bits(&s->pb, n, bits);
 534     if (!s->frame_pred_frame_dct) {
 535         if (has_mv)
 536             put_bits(&s->pb, 2, 2 - field_motion); /* motion_type: frame/field */
 537         put_bits(&s->pb, 1, s->interlaced_dct);
 538     }
 539 }
 540
 541 static av_always_inline void mpeg1_encode_mb_internal(MpegEncContext *s,
 542                                                    DCTELEM block[6][64],
 543                                                    int motion_x, int motion_y,
 544                                                    int mb_block_count)
 545 {
 546     int i, cbp;
 547     const int mb_x = s->mb_x;
 548     const int mb_y = s->mb_y;
 549     const int first_mb= mb_x == s->resync_mb_x && mb_y == s->resync_mb_y;
 550
 551     /* compute cbp */
 552     cbp = 0;
 553     for(i=0;i<mb_block_count;i++) {
 554         if (s->block_last_index[i] >= 0)
 555             cbp |= 1 << (mb_block_count - 1 - i);
 556     }
 557
 558     if (cbp == 0 && !first_mb && s->mv_type == MV_TYPE_16X16 &&
 559         (mb_x != s->mb_width - 1 || (mb_y != s->mb_height - 1 && s->codec_id == CODEC_ID_MPEG1VIDEO)) &&
 560         ((s->pict_type == P_TYPE && (motion_x | motion_y) == 0) ||
 561         (s->pict_type == B_TYPE && s->mv_dir == s->last_mv_dir && (((s->mv_dir & MV_DIR_FORWARD) ? ((s->mv[0][0][0] - s->last_mv[0][0][0])|(s->mv[0][0][1] - s->last_mv[0][0][1])) : 0) |
 562         ((s->mv_dir & MV_DIR_BACKWARD) ? ((s->mv[1][0][0] - s->last_mv[1][0][0])|(s->mv[1][0][1] - s->last_mv[1][0][1])) : 0)) == 0))) {
 563         s->mb_skip_run++;
 564         s->qscale -= s->dquant;
 565         s->skip_count++;
 566         s->misc_bits++;
 567         s->last_bits++;
 568         if(s->pict_type == P_TYPE){
 569             s->last_mv[0][1][0]= s->last_mv[0][0][0]=
 570             s->last_mv[0][1][1]= s->last_mv[0][0][1]= 0;
 571         }
 572     } else {
 573         if(first_mb){
 574             assert(s->mb_skip_run == 0);
 575             encode_mb_skip_run(s, s->mb_x);
 576         }else{
 577             encode_mb_skip_run(s, s->mb_skip_run);
 578         }
 579
 580         if (s->pict_type == I_TYPE) {
 581             if(s->dquant && cbp){
 582                 put_mb_modes(s, 2, 1, 0, 0); /* macroblock_type : macroblock_quant = 1 */
 583                 put_qscale(s);
 584             }else{
 585                 put_mb_modes(s, 1, 1, 0, 0); /* macroblock_type : macroblock_quant = 0 */
 586                 s->qscale -= s->dquant;
 587             }
 588             s->misc_bits+= get_bits_diff(s);
 589             s->i_count++;
 590         } else if (s->mb_intra) {
 591             if(s->dquant && cbp){
 592                 put_mb_modes(s, 6, 0x01, 0, 0);
 593                 put_qscale(s);
 594             }else{
 595                 put_mb_modes(s, 5, 0x03, 0, 0);
 596                 s->qscale -= s->dquant;
 597             }
 598             s->misc_bits+= get_bits_diff(s);
 599             s->i_count++;
 600             memset(s->last_mv, 0, sizeof(s->last_mv));
 601         } else if (s->pict_type == P_TYPE) {
 602             if(s->mv_type == MV_TYPE_16X16){
 603                 if (cbp != 0) {
 604                     if ((motion_x|motion_y) == 0) {
 605                         if(s->dquant){
 606                             put_mb_modes(s, 5, 1, 0, 0); /* macroblock_pattern & quant */
 607                             put_qscale(s);
 608                         }else{
 609                             put_mb_modes(s, 2, 1, 0, 0); /* macroblock_pattern only */
 610                         }
 611                         s->misc_bits+= get_bits_diff(s);
 612                     } else {
 613                         if(s->dquant){
 614                             put_mb_modes(s, 5, 2, 1, 0); /* motion + cbp */
 615                             put_qscale(s);
 616                         }else{
 617                             put_mb_modes(s, 1, 1, 1, 0); /* motion + cbp */
 618                         }
 619                         s->misc_bits+= get_bits_diff(s);
 620                         mpeg1_encode_motion(s, motion_x - s->last_mv[0][0][0], s->f_code);    // RAL: f_code parameter added
 621                         mpeg1_encode_motion(s, motion_y - s->last_mv[0][0][1], s->f_code);    // RAL: f_code parameter added
 622                         s->mv_bits+= get_bits_diff(s);
 623                     }
 624                 } else {
 625                     put_bits(&s->pb, 3, 1); /* motion only */
 626                     if (!s->frame_pred_frame_dct)
 627                         put_bits(&s->pb, 2, 2); /* motion_type: frame */
 628                     s->misc_bits+= get_bits_diff(s);
 629                     mpeg1_encode_motion(s, motion_x - s->last_mv[0][0][0], s->f_code);    // RAL: f_code parameter added
 630                     mpeg1_encode_motion(s, motion_y - s->last_mv[0][0][1], s->f_code);    // RAL: f_code parameter added
 631                     s->qscale -= s->dquant;
 632                     s->mv_bits+= get_bits_diff(s);
 633                 }
 634                 s->last_mv[0][1][0]= s->last_mv[0][0][0]= motion_x;
 635                 s->last_mv[0][1][1]= s->last_mv[0][0][1]= motion_y;
 636             }else{
 637                 assert(!s->frame_pred_frame_dct && s->mv_type == MV_TYPE_FIELD);
 638
 639                 if (cbp) {
 640                     if(s->dquant){
 641                         put_mb_modes(s, 5, 2, 1, 1); /* motion + cbp */
 642                         put_qscale(s);
 643                     }else{
 644                         put_mb_modes(s, 1, 1, 1, 1); /* motion + cbp */
 645                     }
 646                 } else {
 647                     put_bits(&s->pb, 3, 1); /* motion only */
 648                     put_bits(&s->pb, 2, 1); /* motion_type: field */
 649                     s->qscale -= s->dquant;
 650                 }
 651                 s->misc_bits+= get_bits_diff(s);
 652                 for(i=0; i<2; i++){
 653                     put_bits(&s->pb, 1, s->field_select[0][i]);
 654                     mpeg1_encode_motion(s, s->mv[0][i][0] -  s->last_mv[0][i][0]    , s->f_code);
 655                     mpeg1_encode_motion(s, s->mv[0][i][1] - (s->last_mv[0][i][1]>>1), s->f_code);
 656                     s->last_mv[0][i][0]=   s->mv[0][i][0];
 657                     s->last_mv[0][i][1]= 2*s->mv[0][i][1];
 658                 }
 659                 s->mv_bits+= get_bits_diff(s);
 660             }
 661             if(cbp) {
 662                 if (s->chroma_y_shift) {
 663                     put_bits(&s->pb, mbPatTable[cbp][1], mbPatTable[cbp][0]);
 664                 } else {
 665                     put_bits(&s->pb, mbPatTable[cbp>>2][1], mbPatTable[cbp>>2][0]);
 666                     put_bits(&s->pb, 2, cbp & 3);
 667                 }
 668             }
 669             s->f_count++;
 670         } else{
 671             static const int mb_type_len[4]={0,3,4,2}; //bak,for,bi
 672
 673             if(s->mv_type == MV_TYPE_16X16){
 674                 if (cbp){    // With coded bloc pattern
 675                     if (s->dquant) {
 676                         if(s->mv_dir == MV_DIR_FORWARD)
 677                             put_mb_modes(s, 6, 3, 1, 0);
 678                         else
 679                             put_mb_modes(s, mb_type_len[s->mv_dir]+3, 2, 1, 0);
 680                         put_qscale(s);
 681                     } else {
 682                         put_mb_modes(s, mb_type_len[s->mv_dir], 3, 1, 0);
 683                     }
 684                 }else{    // No coded bloc pattern
 685                     put_bits(&s->pb, mb_type_len[s->mv_dir], 2);
 686                     if (!s->frame_pred_frame_dct)
 687                         put_bits(&s->pb, 2, 2); /* motion_type: frame */
 688                     s->qscale -= s->dquant;
 689                 }
 690                 s->misc_bits += get_bits_diff(s);
 691                 if (s->mv_dir&MV_DIR_FORWARD){
 692                     mpeg1_encode_motion(s, s->mv[0][0][0] - s->last_mv[0][0][0], s->f_code);
 693                     mpeg1_encode_motion(s, s->mv[0][0][1] - s->last_mv[0][0][1], s->f_code);
 694                     s->last_mv[0][0][0]=s->last_mv[0][1][0]= s->mv[0][0][0];
 695                     s->last_mv[0][0][1]=s->last_mv[0][1][1]= s->mv[0][0][1];
 696                     s->f_count++;
 697                 }
 698                 if (s->mv_dir&MV_DIR_BACKWARD){
 699                     mpeg1_encode_motion(s, s->mv[1][0][0] - s->last_mv[1][0][0], s->b_code);
 700                     mpeg1_encode_motion(s, s->mv[1][0][1] - s->last_mv[1][0][1], s->b_code);
 701                     s->last_mv[1][0][0]=s->last_mv[1][1][0]= s->mv[1][0][0];
 702                     s->last_mv[1][0][1]=s->last_mv[1][1][1]= s->mv[1][0][1];
 703                     s->b_count++;
 704                 }
 705             }else{
 706                 assert(s->mv_type == MV_TYPE_FIELD);
 707                 assert(!s->frame_pred_frame_dct);
 708                 if (cbp){    // With coded bloc pattern
 709                     if (s->dquant) {
 710                         if(s->mv_dir == MV_DIR_FORWARD)
 711                             put_mb_modes(s, 6, 3, 1, 1);
 712                         else
 713                             put_mb_modes(s, mb_type_len[s->mv_dir]+3, 2, 1, 1);
 714                         put_qscale(s);
 715                     } else {
 716                         put_mb_modes(s, mb_type_len[s->mv_dir], 3, 1, 1);
 717                     }
 718                 }else{    // No coded bloc pattern
 719                     put_bits(&s->pb, mb_type_len[s->mv_dir], 2);
 720                     put_bits(&s->pb, 2, 1); /* motion_type: field */
 721                     s->qscale -= s->dquant;
 722                 }
 723                 s->misc_bits += get_bits_diff(s);
 724                 if (s->mv_dir&MV_DIR_FORWARD){
 725                     for(i=0; i<2; i++){
 726                         put_bits(&s->pb, 1, s->field_select[0][i]);
 727                         mpeg1_encode_motion(s, s->mv[0][i][0] -  s->last_mv[0][i][0]    , s->f_code);
 728                         mpeg1_encode_motion(s, s->mv[0][i][1] - (s->last_mv[0][i][1]>>1), s->f_code);
 729                         s->last_mv[0][i][0]=   s->mv[0][i][0];
 730                         s->last_mv[0][i][1]= 2*s->mv[0][i][1];
 731                     }
 732                     s->f_count++;
 733                 }
 734                 if (s->mv_dir&MV_DIR_BACKWARD){
 735                     for(i=0; i<2; i++){
 736                         put_bits(&s->pb, 1, s->field_select[1][i]);
 737                         mpeg1_encode_motion(s, s->mv[1][i][0] -  s->last_mv[1][i][0]    , s->b_code);
 738                         mpeg1_encode_motion(s, s->mv[1][i][1] - (s->last_mv[1][i][1]>>1), s->b_code);
 739                         s->last_mv[1][i][0]=   s->mv[1][i][0];
 740                         s->last_mv[1][i][1]= 2*s->mv[1][i][1];
 741                     }
 742                     s->b_count++;
 743                 }
 744             }
 745             s->mv_bits += get_bits_diff(s);
 746             if(cbp) {
 747                 if (s->chroma_y_shift) {
 748                     put_bits(&s->pb, mbPatTable[cbp][1], mbPatTable[cbp][0]);
 749                 } else {
 750                     put_bits(&s->pb, mbPatTable[cbp>>2][1], mbPatTable[cbp>>2][0]);
 751                     put_bits(&s->pb, 2, cbp & 3);
 752                 }
 753             }
 754         }
 755         for(i=0;i<mb_block_count;i++) {
 756             if (cbp & (1 << (mb_block_count - 1 - i))) {
 757                 mpeg1_encode_block(s, block[i], i);
 758             }
 759         }
 760         s->mb_skip_run = 0;
 761         if(s->mb_intra)
 762             s->i_tex_bits+= get_bits_diff(s);
 763         else
 764             s->p_tex_bits+= get_bits_diff(s);
 765     }
 766 }
 767
 768 void mpeg1_encode_mb(MpegEncContext *s, DCTELEM block[6][64], int motion_x, int motion_y)
 769 {
 770     if (s->chroma_format == CHROMA_420) mpeg1_encode_mb_internal(s, block, motion_x, motion_y, 6);
 771     else                                mpeg1_encode_mb_internal(s, block, motion_x, motion_y, 8);
 772 }
 773
 774 // RAL: Parameter added: f_or_b_code
 775 static void mpeg1_encode_motion(MpegEncContext *s, int val, int f_or_b_code)
 776 {
 777     int code, bit_size, l, bits, range, sign;
 778
 779     if (val == 0) {
 780         /* zero vector */
 781         code = 0;
 782         put_bits(&s->pb,
 783                  mbMotionVectorTable[0][1],
 784                  mbMotionVectorTable[0][0]);
 785     } else {
 786         bit_size = f_or_b_code - 1;
 787         range = 1 << bit_size;
 788         /* modulo encoding */
 789         l= INT_BIT - 5 - bit_size;
 790         val= (val<<l)>>l;
 791
 792         if (val >= 0) {
 793             val--;
 794             code = (val >> bit_size) + 1;
 795             bits = val & (range - 1);
 796             sign = 0;
 797         } else {
 798             val = -val;
 799             val--;
 800             code = (val >> bit_size) + 1;
 801             bits = val & (range - 1);
 802             sign = 1;
 803         }
 804
 805         assert(code > 0 && code <= 16);
 806
 807         put_bits(&s->pb,
 808                  mbMotionVectorTable[code][1],
 809                  mbMotionVectorTable[code][0]);
 810
 811         put_bits(&s->pb, 1, sign);
 812         if (bit_size > 0) {
 813             put_bits(&s->pb, bit_size, bits);
 814         }
 815     }
 816 }
 817
 818 void ff_mpeg1_encode_init(MpegEncContext *s)
 819 {
 820     static int done=0;
 821
 822     common_init(s);
 823
 824     if(!done){
 825         int f_code;
 826         int mv;
 827         int i;
 828
 829         done=1;
 830         init_rl(&rl_mpeg1, static_rl_table_store[0]);
 831         init_rl(&rl_mpeg2, static_rl_table_store[1]);
 832
 833         for(i=0; i<64; i++)
 834         {
 835                 mpeg1_max_level[0][i]= rl_mpeg1.max_level[0][i];
 836                 mpeg1_index_run[0][i]= rl_mpeg1.index_run[0][i];
 837         }
 838
 839         init_uni_ac_vlc(&rl_mpeg1, uni_mpeg1_ac_vlc_len);
 840         if(s->intra_vlc_format)
 841             init_uni_ac_vlc(&rl_mpeg2, uni_mpeg2_ac_vlc_len);
 842
 843         /* build unified dc encoding tables */
 844         for(i=-255; i<256; i++)
 845         {
 846                 int adiff, index;
 847                 int bits, code;
 848                 int diff=i;
 849
 850                 adiff = FFABS(diff);
 851                 if(diff<0) diff--;
 852                 index = av_log2(2*adiff);
 853
 854                 bits= vlc_dc_lum_bits[index] + index;
 855                 code= (vlc_dc_lum_code[index]<<index) + (diff & ((1 << index) - 1));
 856                 mpeg1_lum_dc_uni[i+255]= bits + (code<<8);
 857
 858                 bits= vlc_dc_chroma_bits[index] + index;
 859                 code= (vlc_dc_chroma_code[index]<<index) + (diff & ((1 << index) - 1));
 860                 mpeg1_chr_dc_uni[i+255]= bits + (code<<8);
 861         }
 862
 863         for(f_code=1; f_code<=MAX_FCODE; f_code++){
 864             for(mv=-MAX_MV; mv<=MAX_MV; mv++){
 865                 int len;
 866
 867                 if(mv==0) len= mbMotionVectorTable[0][1];
 868                 else{
 869                     int val, bit_size, range, code;
 870
 871                     bit_size = f_code - 1;
 872                     range = 1 << bit_size;
 873
 874                     val=mv;
 875                     if (val < 0)
 876                         val = -val;
 877                     val--;
 878                     code = (val >> bit_size) + 1;
 879                     if(code<17){
 880                         len= mbMotionVectorTable[code][1] + 1 + bit_size;
 881                     }else{
 882                         len= mbMotionVectorTable[16][1] + 2 + bit_size;
 883                     }
 884                 }
 885
 886                 mv_penalty[f_code][mv+MAX_MV]= len;
 887             }
 888         }
 889
 890
 891         for(f_code=MAX_FCODE; f_code>0; f_code--){
 892             for(mv=-(8<<f_code); mv<(8<<f_code); mv++){
 893                 fcode_tab[mv+MAX_MV]= f_code;
 894             }
 895         }
 896     }
 897     s->me.mv_penalty= mv_penalty;
 898     s->fcode_tab= fcode_tab;
 899     if(s->codec_id == CODEC_ID_MPEG1VIDEO){
 900         s->min_qcoeff=-255;
 901         s->max_qcoeff= 255;
 902     }else{
 903         s->min_qcoeff=-2047;
 904         s->max_qcoeff= 2047;
 905     }
 906     if (s->intra_vlc_format) {
 907         s->intra_ac_vlc_length=
 908         s->intra_ac_vlc_last_length= uni_mpeg2_ac_vlc_len;
 909     } else {
 910         s->intra_ac_vlc_length=
 911         s->intra_ac_vlc_last_length= uni_mpeg1_ac_vlc_len;
 912     }
 913     s->inter_ac_vlc_length=
 914     s->inter_ac_vlc_last_length= uni_mpeg1_ac_vlc_len;
 915 }
 916
 917 static inline void encode_dc(MpegEncContext *s, int diff, int component)
 918 {
 919   if(((unsigned) (diff+255)) >= 511){
 920         int index;
 921
 922         if(diff<0){
 923             index= av_log2_16bit(-2*diff);
 924             diff--;
 925         }else{
 926             index= av_log2_16bit(2*diff);
 927         }
 928         if (component == 0) {
 929             put_bits(
 930                 &s->pb,
 931                 vlc_dc_lum_bits[index] + index,
 932                 (vlc_dc_lum_code[index]<<index) + (diff & ((1 << index) - 1)));
 933         }else{
 934             put_bits(
 935                 &s->pb,
 936                 vlc_dc_chroma_bits[index] + index,
 937                 (vlc_dc_chroma_code[index]<<index) + (diff & ((1 << index) - 1)));
 938         }
 939   }else{
 940     if (component == 0) {
 941         put_bits(
 942             &s->pb,
 943             mpeg1_lum_dc_uni[diff+255]&0xFF,
 944             mpeg1_lum_dc_uni[diff+255]>>8);
 945     } else {
 946         put_bits(
 947             &s->pb,
 948             mpeg1_chr_dc_uni[diff+255]&0xFF,
 949             mpeg1_chr_dc_uni[diff+255]>>8);
 950     }
 951   }
 952 }
 953
 954 static void mpeg1_encode_block(MpegEncContext *s,
 955                                DCTELEM *block,
 956                                int n)
 957 {
 958     int alevel, level, last_non_zero, dc, diff, i, j, run, last_index, sign;
 959     int code, component;
 960     const uint16_t (*table_vlc)[2] = rl_mpeg1.table_vlc;
 961
 962     last_index = s->block_last_index[n];
 963
 964     /* DC coef */
 965     if (s->mb_intra) {
 966         component = (n <= 3 ? 0 : (n&1) + 1);
 967         dc = block[0]; /* overflow is impossible */
 968         diff = dc - s->last_dc[component];
 969         encode_dc(s, diff, component);
 970         s->last_dc[component] = dc;
 971         i = 1;
 972         if (s->intra_vlc_format)
 973             table_vlc = rl_mpeg2.table_vlc;
 974     } else {
 975         /* encode the first coefficient : needs to be done here because
 976            it is handled slightly differently */
 977         level = block[0];
 978         if (abs(level) == 1) {
 979                 code = ((uint32_t)level >> 31); /* the sign bit */
 980                 put_bits(&s->pb, 2, code | 0x02);
 981                 i = 1;
 982         } else {
 983             i = 0;
 984             last_non_zero = -1;
 985             goto next_coef;
 986         }
 987     }
 988
 989     /* now quantify & encode AC coefs */
 990     last_non_zero = i - 1;
 991
 992     for(;i<=last_index;i++) {
 993         j = s->intra_scantable.permutated[i];
 994         level = block[j];
 995     next_coef:
 996 #if 0
 997         if (level != 0)
 998             dprintf(s->avctx, "level[%d]=%d\n", i, level);
 999 #endif
1000         /* encode using VLC */
1001         if (level != 0) {
1002             run = i - last_non_zero - 1;
1003
1004             alevel= level;
1005             MASK_ABS(sign, alevel)
1006             sign&=1;
1007
1008             if (alevel <= mpeg1_max_level[0][run]){
1009                 code= mpeg1_index_run[0][run] + alevel - 1;
1010                 /* store the vlc & sign at once */
1011                 put_bits(&s->pb, table_vlc[code][1]+1, (table_vlc[code][0]<<1) + sign);
1012             } else {
1013                 /* escape seems to be pretty rare <5% so i dont optimize it */
1014                 put_bits(&s->pb, table_vlc[111][1], table_vlc[111][0]);
1015                 /* escape: only clip in this case */
1016                 put_bits(&s->pb, 6, run);
1017                 if(s->codec_id == CODEC_ID_MPEG1VIDEO){
1018                     if (alevel < 128) {
1019                         put_bits(&s->pb, 8, level & 0xff);
1020                     } else {
1021                         if (level < 0) {
1022                             put_bits(&s->pb, 16, 0x8001 + level + 255);
1023                         } else {
1024                             put_bits(&s->pb, 16, level & 0xffff);
1025                         }
1026                     }
1027                 }else{
1028                     put_bits(&s->pb, 12, level & 0xfff);
1029                 }
1030             }
1031             last_non_zero = i;
1032         }
1033     }
1034     /* end of block */
1035     put_bits(&s->pb, table_vlc[112][1], table_vlc[112][0]);
1036 }
1037 #endif //CONFIG_ENCODERS
1038
1039 /******************************************/
1040 /* decoding */
1041
1042 static VLC dc_lum_vlc;
1043 static VLC dc_chroma_vlc;
1044 static VLC mv_vlc;
1045 static VLC mbincr_vlc;
1046 static VLC mb_ptype_vlc;
1047 static VLC mb_btype_vlc;
1048 static VLC mb_pat_vlc;
1049
1050 static void init_vlcs(void)
1051 {
1052     static int done = 0;
1053
1054     if (!done) {
1055         done = 1;
1056
1057         init_vlc(&dc_lum_vlc, DC_VLC_BITS, 12,
1058                  vlc_dc_lum_bits, 1, 1,
1059                  vlc_dc_lum_code, 2, 2, 1);
1060         init_vlc(&dc_chroma_vlc,  DC_VLC_BITS, 12,
1061                  vlc_dc_chroma_bits, 1, 1,
1062                  vlc_dc_chroma_code, 2, 2, 1);
1063         init_vlc(&mv_vlc, MV_VLC_BITS, 17,
1064                  &mbMotionVectorTable[0][1], 2, 1,
1065                  &mbMotionVectorTable[0][0], 2, 1, 1);
1066         init_vlc(&mbincr_vlc, MBINCR_VLC_BITS, 36,
1067                  &mbAddrIncrTable[0][1], 2, 1,
1068                  &mbAddrIncrTable[0][0], 2, 1, 1);
1069         init_vlc(&mb_pat_vlc, MB_PAT_VLC_BITS, 64,
1070                  &mbPatTable[0][1], 2, 1,
1071                  &mbPatTable[0][0], 2, 1, 1);
1072
1073         init_vlc(&mb_ptype_vlc, MB_PTYPE_VLC_BITS, 7,
1074                  &table_mb_ptype[0][1], 2, 1,
1075                  &table_mb_ptype[0][0], 2, 1, 1);
1076         init_vlc(&mb_btype_vlc, MB_BTYPE_VLC_BITS, 11,
1077                  &table_mb_btype[0][1], 2, 1,
1078                  &table_mb_btype[0][0], 2, 1, 1);
1079         init_rl(&rl_mpeg1, static_rl_table_store[0]);
1080         init_rl(&rl_mpeg2, static_rl_table_store[1]);
1081
1082         init_2d_vlc_rl(&rl_mpeg1, 1);
1083         init_2d_vlc_rl(&rl_mpeg2, 1);
1084     }
1085 }
1086
1087 static inline int get_dmv(MpegEncContext *s)
1088 {
1089     if(get_bits1(&s->gb))
1090         return 1 - (get_bits1(&s->gb) << 1);
1091     else
1092         return 0;
1093 }
1094
1095 static inline int get_qscale(MpegEncContext *s)
1096 {
1097     int qscale = get_bits(&s->gb, 5);
1098     if (s->q_scale_type) {
1099         return non_linear_qscale[qscale];
1100     } else {
1101         return qscale << 1;
1102     }
1103 }
1104
1105 /* motion type (for mpeg2) */
1106 #define MT_FIELD 1
1107 #define MT_FRAME 2
1108 #define MT_16X8  2
1109 #define MT_DMV   3
1110
1111 static int mpeg_decode_mb(MpegEncContext *s,
1112                           DCTELEM block[12][64])
1113 {
1114     int i, j, k, cbp, val, mb_type, motion_type;
1115     const int mb_block_count = 4 + (1<< s->chroma_format);
1116
1117     dprintf(s->avctx, "decode_mb: x=%d y=%d\n", s->mb_x, s->mb_y);
1118
1119     assert(s->mb_skipped==0);
1120
1121     if (s->mb_skip_run-- != 0) {
1122         if(s->pict_type == I_TYPE){
1123             av_log(s->avctx, AV_LOG_ERROR, "skipped MB in I frame at %d %d\n", s->mb_x, s->mb_y);
1124             return -1;
1125         }
1126
1127         /* skip mb */
1128         s->mb_intra = 0;
1129         for(i=0;i<12;i++)
1130             s->block_last_index[i] = -1;
1131         if(s->picture_structure == PICT_FRAME)
1132             s->mv_type = MV_TYPE_16X16;
1133         else
1134             s->mv_type = MV_TYPE_FIELD;
1135         if (s->pict_type == P_TYPE) {
1136             /* if P type, zero motion vector is implied */
1137             s->mv_dir = MV_DIR_FORWARD;
1138             s->mv[0][0][0] = s->mv[0][0][1] = 0;
1139             s->last_mv[0][0][0] = s->last_mv[0][0][1] = 0;
1140             s->last_mv[0][1][0] = s->last_mv[0][1][1] = 0;
1141             s->field_select[0][0]= s->picture_structure - 1;
1142             s->mb_skipped = 1;
1143             s->current_picture.mb_type[ s->mb_x + s->mb_y*s->mb_stride ]= MB_TYPE_SKIP | MB_TYPE_L0 | MB_TYPE_16x16;
1144         } else {
1145             int mb_type;
1146
1147             if(s->mb_x)
1148                 mb_type= s->current_picture.mb_type[ s->mb_x + s->mb_y*s->mb_stride - 1];
1149             else
1150                 mb_type= s->current_picture.mb_type[ s->mb_width + (s->mb_y-1)*s->mb_stride - 1]; // FIXME not sure if this is allowed in mpeg at all,
1151             if(IS_INTRA(mb_type))
1152                 return -1;
1153
1154             /* if B type, reuse previous vectors and directions */
1155             s->mv[0][0][0] = s->last_mv[0][0][0];
1156             s->mv[0][0][1] = s->last_mv[0][0][1];
1157             s->mv[1][0][0] = s->last_mv[1][0][0];
1158             s->mv[1][0][1] = s->last_mv[1][0][1];
1159
1160             s->current_picture.mb_type[ s->mb_x + s->mb_y*s->mb_stride ]=
1161                 mb_type | MB_TYPE_SKIP;
1162 //            assert(s->current_picture.mb_type[ s->mb_x + s->mb_y*s->mb_stride - 1]&(MB_TYPE_16x16|MB_TYPE_16x8));
1163
1164             if((s->mv[0][0][0]|s->mv[0][0][1]|s->mv[1][0][0]|s->mv[1][0][1])==0)
1165                 s->mb_skipped = 1;
1166         }
1167
1168         return 0;
1169     }
1170
1171     switch(s->pict_type) {
1172     default:
1173     case I_TYPE:
1174         if (get_bits1(&s->gb) == 0) {
1175             if (get_bits1(&s->gb) == 0){
1176                 av_log(s->avctx, AV_LOG_ERROR, "invalid mb type in I Frame at %d %d\n", s->mb_x, s->mb_y);
1177                 return -1;
1178             }
1179             mb_type = MB_TYPE_QUANT | MB_TYPE_INTRA;
1180         } else {
1181             mb_type = MB_TYPE_INTRA;
1182         }
1183         break;
1184     case P_TYPE:
1185         mb_type = get_vlc2(&s->gb, mb_ptype_vlc.table, MB_PTYPE_VLC_BITS, 1);
1186         if (mb_type < 0){
1187             av_log(s->avctx, AV_LOG_ERROR, "invalid mb type in P Frame at %d %d\n", s->mb_x, s->mb_y);
1188             return -1;
1189         }
1190         mb_type = ptype2mb_type[ mb_type ];
1191         break;
1192     case B_TYPE:
1193         mb_type = get_vlc2(&s->gb, mb_btype_vlc.table, MB_BTYPE_VLC_BITS, 1);
1194         if (mb_type < 0){
1195             av_log(s->avctx, AV_LOG_ERROR, "invalid mb type in B Frame at %d %d\n", s->mb_x, s->mb_y);
1196             return -1;
1197         }
1198         mb_type = btype2mb_type[ mb_type ];
1199         break;
1200     }
1201     dprintf(s->avctx, "mb_type=%x\n", mb_type);
1202 //    motion_type = 0; /* avoid warning */
1203     if (IS_INTRA(mb_type)) {
1204         s->dsp.clear_blocks(s->block[0]);
1205
1206         if(!s->chroma_y_shift){
1207             s->dsp.clear_blocks(s->block[6]);
1208         }
1209
1210         /* compute dct type */
1211         if (s->picture_structure == PICT_FRAME && //FIXME add a interlaced_dct coded var?
1212             !s->frame_pred_frame_dct) {
1213             s->interlaced_dct = get_bits1(&s->gb);
1214         }
1215
1216         if (IS_QUANT(mb_type))
1217             s->qscale = get_qscale(s);
1218
1219         if (s->concealment_motion_vectors) {
1220             /* just parse them */
1221             if (s->picture_structure != PICT_FRAME)
1222                 skip_bits1(&s->gb); /* field select */
1223
1224             s->mv[0][0][0]= s->last_mv[0][0][0]= s->last_mv[0][1][0] =
1225                 mpeg_decode_motion(s, s->mpeg_f_code[0][0], s->last_mv[0][0][0]);
1226             s->mv[0][0][1]= s->last_mv[0][0][1]= s->last_mv[0][1][1] =
1227                 mpeg_decode_motion(s, s->mpeg_f_code[0][1], s->last_mv[0][0][1]);
1228
1229             skip_bits1(&s->gb); /* marker */
1230         }else
1231             memset(s->last_mv, 0, sizeof(s->last_mv)); /* reset mv prediction */
1232         s->mb_intra = 1;
1233 #ifdef HAVE_XVMC
1234         //one 1 we memcpy blocks in xvmcvideo
1235         if(s->avctx->xvmc_acceleration > 1){
1236             XVMC_pack_pblocks(s,-1);//inter are always full blocks
1237             if(s->swap_uv){
1238                 exchange_uv(s);
1239             }
1240         }
1241 #endif
1242
1243         if (s->codec_id == CODEC_ID_MPEG2VIDEO) {
1244             if(s->flags2 & CODEC_FLAG2_FAST){
1245                 for(i=0;i<6;i++) {
1246                     mpeg2_fast_decode_block_intra(s, s->pblocks[i], i);
1247                 }
1248             }else{
1249                 for(i=0;i<mb_block_count;i++) {
1250                     if (mpeg2_decode_block_intra(s, s->pblocks[i], i) < 0)
1251                         return -1;
1252                 }
1253             }
1254         } else {
1255             for(i=0;i<6;i++) {
1256                 if (mpeg1_decode_block_intra(s, s->pblocks[i], i) < 0)
1257                     return -1;
1258             }
1259         }
1260     } else {
1261         if (mb_type & MB_TYPE_ZERO_MV){
1262             assert(mb_type & MB_TYPE_CBP);
1263
1264             /* compute dct type */
1265             if (s->picture_structure == PICT_FRAME && //FIXME add a interlaced_dct coded var?
1266                 !s->frame_pred_frame_dct) {
1267                 s->interlaced_dct = get_bits1(&s->gb);
1268             }
1269
1270             if (IS_QUANT(mb_type))
1271                 s->qscale = get_qscale(s);
1272
1273             s->mv_dir = MV_DIR_FORWARD;
1274             if(s->picture_structure == PICT_FRAME)
1275                 s->mv_type = MV_TYPE_16X16;
1276             else{
1277                 s->mv_type = MV_TYPE_FIELD;
1278                 mb_type |= MB_TYPE_INTERLACED;
1279                 s->field_select[0][0]= s->picture_structure - 1;
1280             }
1281             s->last_mv[0][0][0] = 0;
1282             s->last_mv[0][0][1] = 0;
1283             s->last_mv[0][1][0] = 0;
1284             s->last_mv[0][1][1] = 0;
1285             s->mv[0][0][0] = 0;
1286             s->mv[0][0][1] = 0;
1287         }else{
1288             assert(mb_type & MB_TYPE_L0L1);
1289 //FIXME decide if MBs in field pictures are MB_TYPE_INTERLACED
1290             /* get additional motion vector type */
1291             if (s->frame_pred_frame_dct)
1292                 motion_type = MT_FRAME;
1293             else{
1294                 motion_type = get_bits(&s->gb, 2);
1295             }
1296
1297             /* compute dct type */
1298             if (s->picture_structure == PICT_FRAME && //FIXME add a interlaced_dct coded var?
1299                 !s->frame_pred_frame_dct && HAS_CBP(mb_type)) {
1300                 s->interlaced_dct = get_bits1(&s->gb);
1301             }
1302
1303             if (IS_QUANT(mb_type))
1304                 s->qscale = get_qscale(s);
1305
1306             /* motion vectors */
1307             s->mv_dir = 0;
1308             for(i=0;i<2;i++) {
1309                 if (USES_LIST(mb_type, i)) {
1310                     s->mv_dir |= (MV_DIR_FORWARD >> i);
1311                     dprintf(s->avctx, "motion_type=%d\n", motion_type);
1312                     switch(motion_type) {
1313                     case MT_FRAME: /* or MT_16X8 */
1314                         if (s->picture_structure == PICT_FRAME) {
1315                             /* MT_FRAME */
1316                             mb_type |= MB_TYPE_16x16;
1317                             s->mv_type = MV_TYPE_16X16;
1318                             s->mv[i][0][0]= s->last_mv[i][0][0]= s->last_mv[i][1][0] =
1319                                 mpeg_decode_motion(s, s->mpeg_f_code[i][0], s->last_mv[i][0][0]);
1320                             s->mv[i][0][1]= s->last_mv[i][0][1]= s->last_mv[i][1][1] =
1321                                 mpeg_decode_motion(s, s->mpeg_f_code[i][1], s->last_mv[i][0][1]);
1322                             /* full_pel: only for mpeg1 */
1323                             if (s->full_pel[i]){
1324                                 s->mv[i][0][0] <<= 1;
1325                                 s->mv[i][0][1] <<= 1;
1326                             }
1327                         } else {
1328                             /* MT_16X8 */
1329                             mb_type |= MB_TYPE_16x8 | MB_TYPE_INTERLACED;
1330                             s->mv_type = MV_TYPE_16X8;
1331                             for(j=0;j<2;j++) {
1332                                 s->field_select[i][j] = get_bits1(&s->gb);
1333                                 for(k=0;k<2;k++) {
1334                                     val = mpeg_decode_motion(s, s->mpeg_f_code[i][k],
1335                                                              s->last_mv[i][j][k]);
1336                                     s->last_mv[i][j][k] = val;
1337                                     s->mv[i][j][k] = val;
1338                                 }
1339                             }
1340                         }
1341                         break;
1342                     case MT_FIELD:
1343                         s->mv_type = MV_TYPE_FIELD;
1344                         if (s->picture_structure == PICT_FRAME) {
1345                             mb_type |= MB_TYPE_16x8 | MB_TYPE_INTERLACED;
1346                             for(j=0;j<2;j++) {
1347                                 s->field_select[i][j] = get_bits1(&s->gb);
1348                                 val = mpeg_decode_motion(s, s->mpeg_f_code[i][0],
1349                                                          s->last_mv[i][j][0]);
1350                                 s->last_mv[i][j][0] = val;
1351                                 s->mv[i][j][0] = val;
1352                                 dprintf(s->avctx, "fmx=%d\n", val);
1353                                 val = mpeg_decode_motion(s, s->mpeg_f_code[i][1],
1354                                                          s->last_mv[i][j][1] >> 1);
1355                                 s->last_mv[i][j][1] = val << 1;
1356                                 s->mv[i][j][1] = val;
1357                                 dprintf(s->avctx, "fmy=%d\n", val);
1358                             }
1359                         } else {
1360                             mb_type |= MB_TYPE_16x16 | MB_TYPE_INTERLACED;
1361                             s->field_select[i][0] = get_bits1(&s->gb);
1362                             for(k=0;k<2;k++) {
1363                                 val = mpeg_decode_motion(s, s->mpeg_f_code[i][k],
1364                                                          s->last_mv[i][0][k]);
1365                                 s->last_mv[i][0][k] = val;
1366                                 s->last_mv[i][1][k] = val;
1367                                 s->mv[i][0][k] = val;
1368                             }
1369                         }
1370                         break;
1371                     case MT_DMV:
1372                         {
1373                             int dmx, dmy, mx, my, m;
1374
1375                             mx = mpeg_decode_motion(s, s->mpeg_f_code[i][0],
1376                                                     s->last_mv[i][0][0]);
1377                             s->last_mv[i][0][0] = mx;
1378                             s->last_mv[i][1][0] = mx;
1379                             dmx = get_dmv(s);
1380                             my = mpeg_decode_motion(s, s->mpeg_f_code[i][1],
1381                                                     s->last_mv[i][0][1] >> 1);
1382                             dmy = get_dmv(s);
1383                             s->mv_type = MV_TYPE_DMV;
1384
1385
1386                             s->last_mv[i][0][1] = my<<1;
1387                             s->last_mv[i][1][1] = my<<1;
1388
1389                             s->mv[i][0][0] = mx;
1390                             s->mv[i][0][1] = my;
1391                             s->mv[i][1][0] = mx;//not used
1392                             s->mv[i][1][1] = my;//not used
1393
1394                             if (s->picture_structure == PICT_FRAME) {
1395                                 mb_type |= MB_TYPE_16x16 | MB_TYPE_INTERLACED;
1396
1397                                 //m = 1 + 2 * s->top_field_first;
1398                                 m = s->top_field_first ? 1 : 3;
1399
1400                                 /* top -> top pred */
1401                                 s->mv[i][2][0] = ((mx * m + (mx > 0)) >> 1) + dmx;
1402                                 s->mv[i][2][1] = ((my * m + (my > 0)) >> 1) + dmy - 1;
1403                                 m = 4 - m;
1404                                 s->mv[i][3][0] = ((mx * m + (mx > 0)) >> 1) + dmx;
1405                                 s->mv[i][3][1] = ((my * m + (my > 0)) >> 1) + dmy + 1;
1406                             } else {
1407                                 mb_type |= MB_TYPE_16x16;
1408
1409                                 s->mv[i][2][0] = ((mx + (mx > 0)) >> 1) + dmx;
1410                                 s->mv[i][2][1] = ((my + (my > 0)) >> 1) + dmy;
1411                                 if(s->picture_structure == PICT_TOP_FIELD)
1412                                     s->mv[i][2][1]--;
1413                                 else
1414                                     s->mv[i][2][1]++;
1415                             }
1416                         }
1417                         break;
1418                     default:
1419                         av_log(s->avctx, AV_LOG_ERROR, "00 motion_type at %d %d\n", s->mb_x, s->mb_y);
1420                         return -1;
1421                     }
1422                 }
1423             }
1424         }
1425
1426         s->mb_intra = 0;
1427         if (HAS_CBP(mb_type)) {
1428             s->dsp.clear_blocks(s->block[0]);
1429
1430             if(!s->chroma_y_shift){
1431                 s->dsp.clear_blocks(s->block[6]);
1432             }
1433
1434             cbp = get_vlc2(&s->gb, mb_pat_vlc.table, MB_PAT_VLC_BITS, 1);
1435             if (cbp < 0 || ((cbp == 0) && (s->chroma_format < 2)) ){
1436                 av_log(s->avctx, AV_LOG_ERROR, "invalid cbp at %d %d\n", s->mb_x, s->mb_y);
1437                 return -1;
1438             }
1439             if(mb_block_count > 6){
1440                  cbp<<= mb_block_count-6;
1441                  cbp |= get_bits(&s->gb, mb_block_count-6);
1442             }
1443
1444 #ifdef HAVE_XVMC
1445             //on 1 we memcpy blocks in xvmcvideo
1446             if(s->avctx->xvmc_acceleration > 1){
1447                 XVMC_pack_pblocks(s,cbp);
1448                 if(s->swap_uv){
1449                     exchange_uv(s);
1450                 }
1451             }
1452 #endif
1453
1454             if (s->codec_id == CODEC_ID_MPEG2VIDEO) {
1455                 if(s->flags2 & CODEC_FLAG2_FAST){
1456                     for(i=0;i<6;i++) {
1457                         if(cbp & 32) {
1458                             mpeg2_fast_decode_block_non_intra(s, s->pblocks[i], i);
1459                         } else {
1460                             s->block_last_index[i] = -1;
1461                         }
1462                         cbp+=cbp;
1463                     }
1464                 }else{
1465                     cbp<<= 12-mb_block_count;
1466
1467                     for(i=0;i<mb_block_count;i++) {
1468                         if ( cbp & (1<<11) ) {
1469                             if (mpeg2_decode_block_non_intra(s, s->pblocks[i], i) < 0)
1470                                 return -1;
1471                         } else {
1472                             s->block_last_index[i] = -1;
1473                         }
1474                         cbp+=cbp;
1475                     }
1476                 }
1477             } else {
1478                 if(s->flags2 & CODEC_FLAG2_FAST){
1479                     for(i=0;i<6;i++) {
1480                         if (cbp & 32) {
1481                             mpeg1_fast_decode_block_inter(s, s->pblocks[i], i);
1482                         } else {
1483                             s->block_last_index[i] = -1;
1484                         }
1485                         cbp+=cbp;
1486                     }
1487                 }else{
1488                     for(i=0;i<6;i++) {
1489                         if (cbp & 32) {
1490                             if (mpeg1_decode_block_inter(s, s->pblocks[i], i) < 0)
1491                                 return -1;
1492                         } else {
1493                             s->block_last_index[i] = -1;
1494                         }
1495                         cbp+=cbp;
1496                     }
1497                 }
1498             }
1499         }else{
1500             for(i=0;i<12;i++)
1501                 s->block_last_index[i] = -1;
1502         }
1503     }
1504
1505     s->current_picture.mb_type[ s->mb_x + s->mb_y*s->mb_stride ]= mb_type;
1506
1507     return 0;
1508 }
1509
1510 /* as h263, but only 17 codes */
1511 static int mpeg_decode_motion(MpegEncContext *s, int fcode, int pred)
1512 {
1513     int code, sign, val, l, shift;
1514
1515     code = get_vlc2(&s->gb, mv_vlc.table, MV_VLC_BITS, 2);
1516     if (code == 0) {
1517         return pred;
1518     }
1519     if (code < 0) {
1520         return 0xffff;
1521     }
1522
1523     sign = get_bits1(&s->gb);
1524     shift = fcode - 1;
1525     val = code;
1526     if (shift) {
1527         val = (val - 1) << shift;
1528         val |= get_bits(&s->gb, shift);
1529         val++;
1530     }
1531     if (sign)
1532         val = -val;
1533     val += pred;
1534
1535     /* modulo decoding */
1536     l= INT_BIT - 5 - shift;
1537     val = (val<<l)>>l;
1538     return val;
1539 }
1540
1541 static inline int decode_dc(GetBitContext *gb, int component)
1542 {
1543     int code, diff;
1544
1545     if (component == 0) {
1546         code = get_vlc2(gb, dc_lum_vlc.table, DC_VLC_BITS, 2);
1547     } else {
1548         code = get_vlc2(gb, dc_chroma_vlc.table, DC_VLC_BITS, 2);
1549     }
1550     if (code < 0){
1551         av_log(NULL, AV_LOG_ERROR, "invalid dc code at\n");
1552         return 0xffff;
1553     }
1554     if (code == 0) {
1555         diff = 0;
1556     } else {
1557         diff = get_xbits(gb, code);
1558     }
1559     return diff;
1560 }
1561
1562 static inline int mpeg1_decode_block_intra(MpegEncContext *s,
1563                                DCTELEM *block,
1564                                int n)
1565 {
1566     int level, dc, diff, i, j, run;
1567     int component;
1568     RLTable *rl = &rl_mpeg1;
1569     uint8_t * const scantable= s->intra_scantable.permutated;
1570     const uint16_t *quant_matrix= s->intra_matrix;
1571     const int qscale= s->qscale;
1572
1573     /* DC coef */
1574     component = (n <= 3 ? 0 : n - 4 + 1);
1575     diff = decode_dc(&s->gb, component);
1576     if (diff >= 0xffff)
1577         return -1;
1578     dc = s->last_dc[component];
1579     dc += diff;
1580     s->last_dc[component] = dc;
1581     block[0] = dc<<3;
1582     dprintf(s->avctx, "dc=%d diff=%d\n", dc, diff);
1583     i = 0;
1584     {
1585         OPEN_READER(re, &s->gb);
1586         /* now quantify & encode AC coefs */
1587         for(;;) {
1588             UPDATE_CACHE(re, &s->gb);
1589             GET_RL_VLC(level, run, re, &s->gb, rl->rl_vlc[0], TEX_VLC_BITS, 2, 0);
1590
1591             if(level == 127){
1592                 break;
1593             } else if(level != 0) {
1594                 i += run;
1595                 j = scantable[i];
1596                 level= (level*qscale*quant_matrix[j])>>4;
1597                 level= (level-1)|1;
1598                 level = (level ^ SHOW_SBITS(re, &s->gb, 1)) - SHOW_SBITS(re, &s->gb, 1);
1599                 LAST_SKIP_BITS(re, &s->gb, 1);
1600             } else {
1601                 /* escape */
1602                 run = SHOW_UBITS(re, &s->gb, 6)+1; LAST_SKIP_BITS(re, &s->gb, 6);
1603                 UPDATE_CACHE(re, &s->gb);
1604                 level = SHOW_SBITS(re, &s->gb, 8); SKIP_BITS(re, &s->gb, 8);
1605                 if (level == -128) {
1606                     level = SHOW_UBITS(re, &s->gb, 8) - 256; LAST_SKIP_BITS(re, &s->gb, 8);
1607                 } else if (level == 0) {
1608                     level = SHOW_UBITS(re, &s->gb, 8)      ; LAST_SKIP_BITS(re, &s->gb, 8);
1609                 }
1610                 i += run;
1611                 j = scantable[i];
1612                 if(level<0){
1613                     level= -level;
1614                     level= (level*qscale*quant_matrix[j])>>4;
1615                     level= (level-1)|1;
1616                     level= -level;
1617                 }else{
1618                     level= (level*qscale*quant_matrix[j])>>4;
1619                     level= (level-1)|1;
1620                 }
1621             }
1622             if (i > 63){
1623                 av_log(s->avctx, AV_LOG_ERROR, "ac-tex damaged at %d %d\n", s->mb_x, s->mb_y);
1624                 return -1;
1625             }
1626
1627             block[j] = level;
1628         }
1629         CLOSE_READER(re, &s->gb);
1630     }
1631     s->block_last_index[n] = i;
1632    return 0;
1633 }
1634
1635 static inline int mpeg1_decode_block_inter(MpegEncContext *s,
1636                                DCTELEM *block,
1637                                int n)
1638 {
1639     int level, i, j, run;
1640     RLTable *rl = &rl_mpeg1;
1641     uint8_t * const scantable= s->intra_scantable.permutated;
1642     const uint16_t *quant_matrix= s->inter_matrix;
1643     const int qscale= s->qscale;
1644
1645     {
1646         OPEN_READER(re, &s->gb);
1647         i = -1;
1648         /* special case for the first coef. no need to add a second vlc table */
1649         UPDATE_CACHE(re, &s->gb);
1650         if (((int32_t)GET_CACHE(re, &s->gb)) < 0) {
1651             level= (3*qscale*quant_matrix[0])>>5;
1652             level= (level-1)|1;
1653             if(GET_CACHE(re, &s->gb)&0x40000000)
1654                 level= -level;
1655             block[0] = level;
1656             i++;
1657             SKIP_BITS(re, &s->gb, 2);
1658             if(((int32_t)GET_CACHE(re, &s->gb)) <= (int32_t)0xBFFFFFFF)
1659                 goto end;
1660         }
1661
1662         /* now quantify & encode AC coefs */
1663         for(;;) {
1664             GET_RL_VLC(level, run, re, &s->gb, rl->rl_vlc[0], TEX_VLC_BITS, 2, 0);
1665
1666             if(level != 0) {
1667                 i += run;
1668                 j = scantable[i];
1669                 level= ((level*2+1)*qscale*quant_matrix[j])>>5;
1670                 level= (level-1)|1;
1671                 level = (level ^ SHOW_SBITS(re, &s->gb, 1)) - SHOW_SBITS(re, &s->gb, 1);
1672                 SKIP_BITS(re, &s->gb, 1);
1673             } else {
1674                 /* escape */
1675                 run = SHOW_UBITS(re, &s->gb, 6)+1; LAST_SKIP_BITS(re, &s->gb, 6);
1676                 UPDATE_CACHE(re, &s->gb);
1677                 level = SHOW_SBITS(re, &s->gb, 8); SKIP_BITS(re, &s->gb, 8);
1678                 if (level == -128) {
1679                     level = SHOW_UBITS(re, &s->gb, 8) - 256; SKIP_BITS(re, &s->gb, 8);
1680                 } else if (level == 0) {
1681                     level = SHOW_UBITS(re, &s->gb, 8)      ; SKIP_BITS(re, &s->gb, 8);
1682                 }
1683                 i += run;
1684                 j = scantable[i];
1685                 if(level<0){
1686                     level= -level;
1687                     level= ((level*2+1)*qscale*quant_matrix[j])>>5;
1688                     level= (level-1)|1;
1689                     level= -level;
1690                 }else{
1691                     level= ((level*2+1)*qscale*quant_matrix[j])>>5;
1692                     level= (level-1)|1;
1693                 }
1694             }
1695             if (i > 63){
1696                 av_log(s->avctx, AV_LOG_ERROR, "ac-tex damaged at %d %d\n", s->mb_x, s->mb_y);
1697                 return -1;
1698             }
1699
1700             block[j] = level;
1701             if(((int32_t)GET_CACHE(re, &s->gb)) <= (int32_t)0xBFFFFFFF)
1702                 break;
1703             UPDATE_CACHE(re, &s->gb);
1704         }
1705 end:
1706         LAST_SKIP_BITS(re, &s->gb, 2);
1707         CLOSE_READER(re, &s->gb);
1708     }
1709     s->block_last_index[n] = i;
1710     return 0;
1711 }
1712
1713 static inline int mpeg1_fast_decode_block_inter(MpegEncContext *s, DCTELEM *block, int n)
1714 {
1715     int level, i, j, run;
1716     RLTable *rl = &rl_mpeg1;
1717     uint8_t * const scantable= s->intra_scantable.permutated;
1718     const int qscale= s->qscale;
1719
1720     {
1721         OPEN_READER(re, &s->gb);
1722         i = -1;
1723         /* special case for the first coef. no need to add a second vlc table */
1724         UPDATE_CACHE(re, &s->gb);
1725         if (((int32_t)GET_CACHE(re, &s->gb)) < 0) {
1726             level= (3*qscale)>>1;
1727             level= (level-1)|1;
1728             if(GET_CACHE(re, &s->gb)&0x40000000)
1729                 level= -level;
1730             block[0] = level;
1731             i++;
1732             SKIP_BITS(re, &s->gb, 2);
1733             if(((int32_t)GET_CACHE(re, &s->gb)) <= (int32_t)0xBFFFFFFF)
1734                 goto end;
1735         }
1736
1737         /* now quantify & encode AC coefs */
1738         for(;;) {
1739             GET_RL_VLC(level, run, re, &s->gb, rl->rl_vlc[0], TEX_VLC_BITS, 2, 0);
1740
1741             if(level != 0) {
1742                 i += run;
1743                 j = scantable[i];
1744                 level= ((level*2+1)*qscale)>>1;
1745                 level= (level-1)|1;
1746                 level = (level ^ SHOW_SBITS(re, &s->gb, 1)) - SHOW_SBITS(re, &s->gb, 1);
1747                 SKIP_BITS(re, &s->gb, 1);
1748             } else {
1749                 /* escape */
1750                 run = SHOW_UBITS(re, &s->gb, 6)+1; LAST_SKIP_BITS(re, &s->gb, 6);
1751                 UPDATE_CACHE(re, &s->gb);
1752                 level = SHOW_SBITS(re, &s->gb, 8); SKIP_BITS(re, &s->gb, 8);
1753                 if (level == -128) {
1754                     level = SHOW_UBITS(re, &s->gb, 8) - 256; SKIP_BITS(re, &s->gb, 8);
1755                 } else if (level == 0) {
1756                     level = SHOW_UBITS(re, &s->gb, 8)      ; SKIP_BITS(re, &s->gb, 8);
1757                 }
1758                 i += run;
1759                 j = scantable[i];
1760                 if(level<0){
1761                     level= -level;
1762                     level= ((level*2+1)*qscale)>>1;
1763                     level= (level-1)|1;
1764                     level= -level;
1765                 }else{
1766                     level= ((level*2+1)*qscale)>>1;
1767                     level= (level-1)|1;
1768                 }
1769             }
1770
1771             block[j] = level;
1772             if(((int32_t)GET_CACHE(re, &s->gb)) <= (int32_t)0xBFFFFFFF)
1773                 break;
1774             UPDATE_CACHE(re, &s->gb);
1775         }
1776 end:
1777         LAST_SKIP_BITS(re, &s->gb, 2);
1778         CLOSE_READER(re, &s->gb);
1779     }
1780     s->block_last_index[n] = i;
1781     return 0;
1782 }
1783
1784
1785 static inline int mpeg2_decode_block_non_intra(MpegEncContext *s,
1786                                DCTELEM *block,
1787                                int n)
1788 {
1789     int level, i, j, run;
1790     RLTable *rl = &rl_mpeg1;
1791     uint8_t * const scantable= s->intra_scantable.permutated;
1792     const uint16_t *quant_matrix;
1793     const int qscale= s->qscale;
1794     int mismatch;
1795
1796     mismatch = 1;
1797
1798     {
1799         OPEN_READER(re, &s->gb);
1800         i = -1;
1801         if (n < 4)
1802             quant_matrix = s->inter_matrix;
1803         else
1804             quant_matrix = s->chroma_inter_matrix;
1805
1806         /* special case for the first coef. no need to add a second vlc table */
1807         UPDATE_CACHE(re, &s->gb);
1808         if (((int32_t)GET_CACHE(re, &s->gb)) < 0) {
1809             level= (3*qscale*quant_matrix[0])>>5;
1810             if(GET_CACHE(re, &s->gb)&0x40000000)
1811                 level= -level;
1812             block[0] = level;
1813             mismatch ^= level;
1814             i++;
1815             SKIP_BITS(re, &s->gb, 2);
1816             if(((int32_t)GET_CACHE(re, &s->gb)) <= (int32_t)0xBFFFFFFF)
1817                 goto end;
1818         }
1819
1820         /* now quantify & encode AC coefs */
1821         for(;;) {
1822             GET_RL_VLC(level, run, re, &s->gb, rl->rl_vlc[0], TEX_VLC_BITS, 2, 0);
1823
1824             if(level != 0) {
1825                 i += run;
1826                 j = scantable[i];
1827                 level= ((level*2+1)*qscale*quant_matrix[j])>>5;
1828                 level = (level ^ SHOW_SBITS(re, &s->gb, 1)) - SHOW_SBITS(re, &s->gb, 1);
1829                 SKIP_BITS(re, &s->gb, 1);
1830             } else {
1831                 /* escape */
1832                 run = SHOW_UBITS(re, &s->gb, 6)+1; LAST_SKIP_BITS(re, &s->gb, 6);
1833                 UPDATE_CACHE(re, &s->gb);
1834                 level = SHOW_SBITS(re, &s->gb, 12); SKIP_BITS(re, &s->gb, 12);
1835
1836                 i += run;
1837                 j = scantable[i];
1838                 if(level<0){
1839                     level= ((-level*2+1)*qscale*quant_matrix[j])>>5;
1840                     level= -level;
1841                 }else{
1842                     level= ((level*2+1)*qscale*quant_matrix[j])>>5;
1843                 }
1844             }
1845             if (i > 63){
1846                 av_log(s->avctx, AV_LOG_ERROR, "ac-tex damaged at %d %d\n", s->mb_x, s->mb_y);
1847                 return -1;
1848             }
1849
1850             mismatch ^= level;
1851             block[j] = level;
1852             if(((int32_t)GET_CACHE(re, &s->gb)) <= (int32_t)0xBFFFFFFF)
1853                 break;
1854             UPDATE_CACHE(re, &s->gb);
1855         }
1856 end:
1857         LAST_SKIP_BITS(re, &s->gb, 2);
1858         CLOSE_READER(re, &s->gb);
1859     }
1860     block[63] ^= (mismatch & 1);
1861
1862     s->block_last_index[n] = i;
1863     return 0;
1864 }
1865
1866 static inline int mpeg2_fast_decode_block_non_intra(MpegEncContext *s,
1867                                DCTELEM *block,
1868                                int n)
1869 {
1870     int level, i, j, run;
1871     RLTable *rl = &rl_mpeg1;
1872     uint8_t * const scantable= s->intra_scantable.permutated;
1873     const int qscale= s->qscale;
1874     OPEN_READER(re, &s->gb);
1875     i = -1;
1876
1877     /* special case for the first coef. no need to add a second vlc table */
1878     UPDATE_CACHE(re, &s->gb);
1879     if (((int32_t)GET_CACHE(re, &s->gb)) < 0) {
1880         level= (3*qscale)>>1;
1881         if(GET_CACHE(re, &s->gb)&0x40000000)
1882             level= -level;
1883         block[0] = level;
1884         i++;
1885         SKIP_BITS(re, &s->gb, 2);
1886         if(((int32_t)GET_CACHE(re, &s->gb)) <= (int32_t)0xBFFFFFFF)
1887             goto end;
1888     }
1889
1890     /* now quantify & encode AC coefs */
1891     for(;;) {
1892         GET_RL_VLC(level, run, re, &s->gb, rl->rl_vlc[0], TEX_VLC_BITS, 2, 0);
1893
1894         if(level != 0) {
1895             i += run;
1896             j = scantable[i];
1897             level= ((level*2+1)*qscale)>>1;
1898             level = (level ^ SHOW_SBITS(re, &s->gb, 1)) - SHOW_SBITS(re, &s->gb, 1);
1899             SKIP_BITS(re, &s->gb, 1);
1900         } else {
1901             /* escape */
1902             run = SHOW_UBITS(re, &s->gb, 6)+1; LAST_SKIP_BITS(re, &s->gb, 6);
1903             UPDATE_CACHE(re, &s->gb);
1904             level = SHOW_SBITS(re, &s->gb, 12); SKIP_BITS(re, &s->gb, 12);
1905
1906             i += run;
1907             j = scantable[i];
1908             if(level<0){
1909                 level= ((-level*2+1)*qscale)>>1;
1910                 level= -level;
1911             }else{
1912                 level= ((level*2+1)*qscale)>>1;
1913             }
1914         }
1915
1916         block[j] = level;
1917         if(((int32_t)GET_CACHE(re, &s->gb)) <= (int32_t)0xBFFFFFFF)
1918             break;
1919         UPDATE_CACHE(re, &s->gb);
1920     }
1921 end:
1922     LAST_SKIP_BITS(re, &s->gb, 2);
1923     CLOSE_READER(re, &s->gb);
1924     s->block_last_index[n] = i;
1925     return 0;
1926 }
1927
1928
1929 static inline int mpeg2_decode_block_intra(MpegEncContext *s,
1930                                DCTELEM *block,
1931                                int n)
1932 {
1933     int level, dc, diff, i, j, run;
1934     int component;
1935     RLTable *rl;
1936     uint8_t * const scantable= s->intra_scantable.permutated;
1937     const uint16_t *quant_matrix;
1938     const int qscale= s->qscale;
1939     int mismatch;
1940
1941     /* DC coef */
1942     if (n < 4){
1943         quant_matrix = s->intra_matrix;
1944         component = 0;
1945     }else{
1946         quant_matrix = s->chroma_intra_matrix;
1947         component = (n&1) + 1;
1948     }
1949     diff = decode_dc(&s->gb, component);
1950     if (diff >= 0xffff)
1951         return -1;
1952     dc = s->last_dc[component];
1953     dc += diff;
1954     s->last_dc[component] = dc;
1955     block[0] = dc << (3 - s->intra_dc_precision);
1956     dprintf(s->avctx, "dc=%d\n", block[0]);
1957     mismatch = block[0] ^ 1;
1958     i = 0;
1959     if (s->intra_vlc_format)
1960         rl = &rl_mpeg2;
1961     else
1962         rl = &rl_mpeg1;
1963
1964     {
1965         OPEN_READER(re, &s->gb);
1966         /* now quantify & encode AC coefs */
1967         for(;;) {
1968             UPDATE_CACHE(re, &s->gb);
1969             GET_RL_VLC(level, run, re, &s->gb, rl->rl_vlc[0], TEX_VLC_BITS, 2, 0);
1970
1971             if(level == 127){
1972                 break;
1973             } else if(level != 0) {
1974                 i += run;
1975                 j = scantable[i];
1976                 level= (level*qscale*quant_matrix[j])>>4;
1977                 level = (level ^ SHOW_SBITS(re, &s->gb, 1)) - SHOW_SBITS(re, &s->gb, 1);
1978                 LAST_SKIP_BITS(re, &s->gb, 1);
1979             } else {
1980                 /* escape */
1981                 run = SHOW_UBITS(re, &s->gb, 6)+1; LAST_SKIP_BITS(re, &s->gb, 6);
1982                 UPDATE_CACHE(re, &s->gb);
1983                 level = SHOW_SBITS(re, &s->gb, 12); SKIP_BITS(re, &s->gb, 12);
1984                 i += run;
1985                 j = scantable[i];
1986                 if(level<0){
1987                     level= (-level*qscale*quant_matrix[j])>>4;
1988                     level= -level;
1989                 }else{
1990                     level= (level*qscale*quant_matrix[j])>>4;
1991                 }
1992             }
1993             if (i > 63){
1994                 av_log(s->avctx, AV_LOG_ERROR, "ac-tex damaged at %d %d\n", s->mb_x, s->mb_y);
1995                 return -1;
1996             }
1997
1998             mismatch^= level;
1999             block[j] = level;
2000         }
2001         CLOSE_READER(re, &s->gb);
2002     }
2003     block[63]^= mismatch&1;
2004
2005     s->block_last_index[n] = i;
2006     return 0;
2007 }
2008
2009 static inline int mpeg2_fast_decode_block_intra(MpegEncContext *s,
2010                                DCTELEM *block,
2011                                int n)
2012 {
2013     int level, dc, diff, j, run;
2014     int component;
2015     RLTable *rl;
2016     uint8_t * scantable= s->intra_scantable.permutated;
2017     const uint16_t *quant_matrix;
2018     const int qscale= s->qscale;
2019
2020     /* DC coef */
2021     if (n < 4){
2022         quant_matrix = s->intra_matrix;
2023         component = 0;
2024     }else{
2025         quant_matrix = s->chroma_intra_matrix;
2026         component = (n&1) + 1;
2027     }
2028     diff = decode_dc(&s->gb, component);
2029     if (diff >= 0xffff)
2030         return -1;
2031     dc = s->last_dc[component];
2032     dc += diff;
2033     s->last_dc[component] = dc;
2034     block[0] = dc << (3 - s->intra_dc_precision);
2035     if (s->intra_vlc_format)
2036         rl = &rl_mpeg2;
2037     else
2038         rl = &rl_mpeg1;
2039
2040     {
2041         OPEN_READER(re, &s->gb);
2042         /* now quantify & encode AC coefs */
2043         for(;;) {
2044             UPDATE_CACHE(re, &s->gb);
2045             GET_RL_VLC(level, run, re, &s->gb, rl->rl_vlc[0], TEX_VLC_BITS, 2, 0);
2046
2047             if(level == 127){
2048                 break;
2049             } else if(level != 0) {
2050                 scantable += run;
2051                 j = *scantable;
2052                 level= (level*qscale*quant_matrix[j])>>4;
2053                 level = (level ^ SHOW_SBITS(re, &s->gb, 1)) - SHOW_SBITS(re, &s->gb, 1);
2054                 LAST_SKIP_BITS(re, &s->gb, 1);
2055             } else {
2056                 /* escape */
2057                 run = SHOW_UBITS(re, &s->gb, 6)+1; LAST_SKIP_BITS(re, &s->gb, 6);
2058                 UPDATE_CACHE(re, &s->gb);
2059                 level = SHOW_SBITS(re, &s->gb, 12); SKIP_BITS(re, &s->gb, 12);
2060                 scantable += run;
2061                 j = *scantable;
2062                 if(level<0){
2063                     level= (-level*qscale*quant_matrix[j])>>4;
2064                     level= -level;
2065                 }else{
2066                     level= (level*qscale*quant_matrix[j])>>4;
2067                 }
2068             }
2069
2070             block[j] = level;
2071         }
2072         CLOSE_READER(re, &s->gb);
2073     }
2074
2075     s->block_last_index[n] = scantable - s->intra_scantable.permutated;
2076     return 0;
2077 }
2078
2079 typedef struct Mpeg1Context {
2080     MpegEncContext mpeg_enc_ctx;
2081     int mpeg_enc_ctx_allocated; /* true if decoding context allocated */
2082     int repeat_field; /* true if we must repeat the field */
2083     AVPanScan pan_scan; /** some temporary storage for the panscan */
2084     int slice_count;
2085     int swap_uv;//indicate VCR2
2086     int save_aspect_info;
2087     int save_width, save_height;
2088     AVRational frame_rate_ext;       ///< MPEG-2 specific framerate modificator
2089
2090 } Mpeg1Context;
2091
2092 static int mpeg_decode_init(AVCodecContext *avctx)
2093 {
2094     Mpeg1Context *s = avctx->priv_data;
2095     MpegEncContext *s2 = &s->mpeg_enc_ctx;
2096     int i;
2097
2098     //we need some parmutation to store
2099     //matrixes, until MPV_common_init()
2100     //set the real permutatuon
2101     for(i=0;i<64;i++)
2102        s2->dsp.idct_permutation[i]=i;
2103
2104     MPV_decode_defaults(s2);
2105
2106     s->mpeg_enc_ctx.avctx= avctx;
2107     s->mpeg_enc_ctx.flags= avctx->flags;
2108     s->mpeg_enc_ctx.flags2= avctx->flags2;
2109     common_init(&s->mpeg_enc_ctx);
2110     init_vlcs();
2111
2112     s->mpeg_enc_ctx_allocated = 0;
2113     s->mpeg_enc_ctx.picture_number = 0;
2114     s->repeat_field = 0;
2115     s->mpeg_enc_ctx.codec_id= avctx->codec->id;
2116     return 0;
2117 }
2118
2119 static void quant_matrix_rebuild(uint16_t *matrix, const uint8_t *old_perm,
2120                                      const uint8_t *new_perm){
2121     uint16_t temp_matrix[64];
2122     int i;
2123
2124     memcpy(temp_matrix,matrix,64*sizeof(uint16_t));
2125
2126     for(i=0;i<64;i++){
2127         matrix[new_perm[i]] = temp_matrix[old_perm[i]];
2128     }
2129 }
2130
2131 //Call this function when we know all parameters
2132 //it may be called in different places for mpeg1 and mpeg2
2133 static int mpeg_decode_postinit(AVCodecContext *avctx){
2134     Mpeg1Context *s1 = avctx->priv_data;
2135     MpegEncContext *s = &s1->mpeg_enc_ctx;
2136     uint8_t old_permutation[64];
2137
2138     if (
2139         (s1->mpeg_enc_ctx_allocated == 0)||
2140         avctx->coded_width  != s->width ||
2141         avctx->coded_height != s->height||
2142         s1->save_width != s->width ||
2143         s1->save_height != s->height ||
2144         s1->save_aspect_info != s->aspect_ratio_info||
2145         0)
2146     {
2147
2148         if (s1->mpeg_enc_ctx_allocated) {
2149             ParseContext pc= s->parse_context;
2150             s->parse_context.buffer=0;
2151             MPV_common_end(s);
2152             s->parse_context= pc;
2153         }
2154
2155         if( (s->width == 0 )||(s->height == 0))
2156             return -2;
2157
2158         avcodec_set_dimensions(avctx, s->width, s->height);
2159         avctx->bit_rate = s->bit_rate;
2160         s1->save_aspect_info = s->aspect_ratio_info;
2161         s1->save_width = s->width;
2162         s1->save_height = s->height;
2163
2164      //low_delay may be forced, in this case we will have B frames
2165      //that behave like P frames
2166         avctx->has_b_frames = !(s->low_delay);
2167
2168         if(avctx->sub_id==1){//s->codec_id==avctx->codec_id==CODEC_ID
2169             //mpeg1 fps
2170             avctx->time_base.den= ff_frame_rate_tab[s->frame_rate_index].num;
2171             avctx->time_base.num= ff_frame_rate_tab[s->frame_rate_index].den;
2172             //mpeg1 aspect
2173             avctx->sample_aspect_ratio= av_d2q(
2174                     1.0/mpeg1_aspect[s->aspect_ratio_info], 255);
2175
2176         }else{//mpeg2
2177         //mpeg2 fps
2178             av_reduce(
2179                 &s->avctx->time_base.den,
2180                 &s->avctx->time_base.num,
2181                 ff_frame_rate_tab[s->frame_rate_index].num * s1->frame_rate_ext.num,
2182                 ff_frame_rate_tab[s->frame_rate_index].den * s1->frame_rate_ext.den,
2183                 1<<30);
2184         //mpeg2 aspect
2185             if(s->aspect_ratio_info > 1){
2186                 if( (s1->pan_scan.width == 0 )||(s1->pan_scan.height == 0) ){
2187                     s->avctx->sample_aspect_ratio=
2188                         av_div_q(
2189                          mpeg2_aspect[s->aspect_ratio_info],
2190                          (AVRational){s->width, s->height}
2191                          );
2192                 }else{
2193                     s->avctx->sample_aspect_ratio=
2194                         av_div_q(
2195                          mpeg2_aspect[s->aspect_ratio_info],
2196                          (AVRational){s1->pan_scan.width, s1->pan_scan.height}
2197                         );
2198                 }
2199             }else{
2200                 s->avctx->sample_aspect_ratio=
2201                     mpeg2_aspect[s->aspect_ratio_info];
2202             }
2203         }//mpeg2
2204
2205         if(avctx->xvmc_acceleration){
2206             avctx->pix_fmt = avctx->get_format(avctx,pixfmt_xvmc_mpg2_420);
2207         }else{
2208             if(s->chroma_format <  2){
2209                 avctx->pix_fmt = avctx->get_format(avctx,pixfmt_yuv_420);
2210             }else
2211             if(s->chroma_format == 2){
2212                 avctx->pix_fmt = avctx->get_format(avctx,pixfmt_yuv_422);
2213             }else
2214             if(s->chroma_format >  2){
2215                 avctx->pix_fmt = avctx->get_format(avctx,pixfmt_yuv_444);
2216             }
2217         }
2218         //until then pix_fmt may be changed right after codec init
2219         if( avctx->pix_fmt == PIX_FMT_XVMC_MPEG2_IDCT )
2220             if( avctx->idct_algo == FF_IDCT_AUTO )
2221                 avctx->idct_algo = FF_IDCT_SIMPLE;
2222
2223         //quantization matrixes may need reordering
2224         //if dct permutation is changed
2225         memcpy(old_permutation,s->dsp.idct_permutation,64*sizeof(uint8_t));
2226
2227         if (MPV_common_init(s) < 0)
2228             return -2;
2229
2230         quant_matrix_rebuild(s->intra_matrix,       old_permutation,s->dsp.idct_permutation);
2231         quant_matrix_rebuild(s->inter_matrix,       old_permutation,s->dsp.idct_permutation);
2232         quant_matrix_rebuild(s->chroma_intra_matrix,old_permutation,s->dsp.idct_permutation);
2233         quant_matrix_rebuild(s->chroma_inter_matrix,old_permutation,s->dsp.idct_permutation);
2234
2235         s1->mpeg_enc_ctx_allocated = 1;
2236     }
2237     return 0;
2238 }
2239
2240 static int mpeg1_decode_picture(AVCodecContext *avctx,
2241                                 const uint8_t *buf, int buf_size)
2242 {
2243     Mpeg1Context *s1 = avctx->priv_data;
2244     MpegEncContext *s = &s1->mpeg_enc_ctx;
2245     int ref, f_code, vbv_delay;
2246
2247     if(mpeg_decode_postinit(s->avctx) < 0)
2248        return -2;
2249
2250     init_get_bits(&s->gb, buf, buf_size*8);
2251
2252     ref = get_bits(&s->gb, 10); /* temporal ref */
2253     s->pict_type = get_bits(&s->gb, 3);
2254     if(s->pict_type == 0 || s->pict_type > 3)
2255         return -1;
2256
2257     vbv_delay= get_bits(&s->gb, 16);
2258     if (s->pict_type == P_TYPE || s->pict_type == B_TYPE) {
2259         s->full_pel[0] = get_bits1(&s->gb);
2260         f_code = get_bits(&s->gb, 3);
2261         if (f_code == 0 && avctx->error_resilience >= FF_ER_COMPLIANT)
2262             return -1;
2263         s->mpeg_f_code[0][0] = f_code;
2264         s->mpeg_f_code[0][1] = f_code;
2265     }
2266     if (s->pict_type == B_TYPE) {
2267         s->full_pel[1] = get_bits1(&s->gb);
2268         f_code = get_bits(&s->gb, 3);
2269         if (f_code == 0 && avctx->error_resilience >= FF_ER_COMPLIANT)
2270             return -1;
2271         s->mpeg_f_code[1][0] = f_code;
2272         s->mpeg_f_code[1][1] = f_code;
2273     }
2274     s->current_picture.pict_type= s->pict_type;
2275     s->current_picture.key_frame= s->pict_type == I_TYPE;
2276
2277     if(avctx->debug & FF_DEBUG_PICT_INFO)
2278         av_log(avctx, AV_LOG_DEBUG, "vbv_delay %d, ref %d type:%d\n", vbv_delay, ref, s->pict_type);
2279
2280     s->y_dc_scale = 8;
2281     s->c_dc_scale = 8;
2282     s->first_slice = 1;
2283     return 0;
2284 }
2285
2286 static void mpeg_decode_sequence_extension(Mpeg1Context *s1)
2287 {
2288     MpegEncContext *s= &s1->mpeg_enc_ctx;
2289     int horiz_size_ext, vert_size_ext;
2290     int bit_rate_ext;
2291
2292     skip_bits(&s->gb, 1); /* profil and level esc*/
2293     s->avctx->profile= get_bits(&s->gb, 3);
2294     s->avctx->level= get_bits(&s->gb, 4);
2295     s->progressive_sequence = get_bits1(&s->gb); /* progressive_sequence */
2296     s->chroma_format = get_bits(&s->gb, 2); /* chroma_format 1=420, 2=422, 3=444 */
2297     horiz_size_ext = get_bits(&s->gb, 2);
2298     vert_size_ext = get_bits(&s->gb, 2);
2299     s->width |= (horiz_size_ext << 12);
2300     s->height |= (vert_size_ext << 12);
2301     bit_rate_ext = get_bits(&s->gb, 12);  /* XXX: handle it */
2302     s->bit_rate += (bit_rate_ext << 18) * 400;
2303     skip_bits1(&s->gb); /* marker */
2304     s->avctx->rc_buffer_size += get_bits(&s->gb, 8)*1024*16<<10;
2305
2306     s->low_delay = get_bits1(&s->gb);
2307     if(s->flags & CODEC_FLAG_LOW_DELAY) s->low_delay=1;
2308
2309     s1->frame_rate_ext.num = get_bits(&s->gb, 2)+1;
2310     s1->frame_rate_ext.den = get_bits(&s->gb, 5)+1;
2311
2312     dprintf(s->avctx, "sequence extension\n");
2313     s->codec_id= s->avctx->codec_id= CODEC_ID_MPEG2VIDEO;
2314     s->avctx->sub_id = 2; /* indicates mpeg2 found */
2315
2316     if(s->avctx->debug & FF_DEBUG_PICT_INFO)
2317         av_log(s->avctx, AV_LOG_DEBUG, "profile: %d, level: %d vbv buffer: %d, bitrate:%d\n",
2318                s->avctx->profile, s->avctx->level, s->avctx->rc_buffer_size, s->bit_rate);
2319
2320 }
2321
2322 static void mpeg_decode_sequence_display_extension(Mpeg1Context *s1)
2323 {
2324     MpegEncContext *s= &s1->mpeg_enc_ctx;
2325     int color_description, w, h;
2326
2327     skip_bits(&s->gb, 3); /* video format */
2328     color_description= get_bits1(&s->gb);
2329     if(color_description){
2330         skip_bits(&s->gb, 8); /* color primaries */
2331         skip_bits(&s->gb, 8); /* transfer_characteristics */
2332         skip_bits(&s->gb, 8); /* matrix_coefficients */
2333     }
2334     w= get_bits(&s->gb, 14);
2335     skip_bits(&s->gb, 1); //marker
2336     h= get_bits(&s->gb, 14);
2337     skip_bits(&s->gb, 1); //marker
2338
2339     s1->pan_scan.width= 16*w;
2340     s1->pan_scan.height=16*h;
2341
2342     if(s->avctx->debug & FF_DEBUG_PICT_INFO)
2343         av_log(s->avctx, AV_LOG_DEBUG, "sde w:%d, h:%d\n", w, h);
2344 }
2345
2346 static void mpeg_decode_picture_display_extension(Mpeg1Context *s1)
2347 {
2348     MpegEncContext *s= &s1->mpeg_enc_ctx;
2349     int i,nofco;
2350
2351     nofco = 1;
2352     if(s->progressive_sequence){
2353         if(s->repeat_first_field){
2354             nofco++;
2355             if(s->top_field_first)
2356                 nofco++;
2357         }
2358     }else{
2359         if(s->picture_structure == PICT_FRAME){
2360             nofco++;
2361             if(s->repeat_first_field)
2362                 nofco++;
2363         }
2364     }
2365     for(i=0; i<nofco; i++){
2366         s1->pan_scan.position[i][0]= get_sbits(&s->gb, 16);
2367         skip_bits(&s->gb, 1); //marker
2368         s1->pan_scan.position[i][1]= get_sbits(&s->gb, 16);
2369         skip_bits(&s->gb, 1); //marker
2370     }
2371
2372     if(s->avctx->debug & FF_DEBUG_PICT_INFO)
2373         av_log(s->avctx, AV_LOG_DEBUG, "pde (%d,%d) (%d,%d) (%d,%d)\n",
2374             s1->pan_scan.position[0][0], s1->pan_scan.position[0][1],
2375             s1->pan_scan.position[1][0], s1->pan_scan.position[1][1],
2376             s1->pan_scan.position[2][0], s1->pan_scan.position[2][1]
2377         );
2378 }
2379
2380 static void mpeg_decode_quant_matrix_extension(MpegEncContext *s)
2381 {
2382     int i, v, j;
2383
2384     dprintf(s->avctx, "matrix extension\n");
2385
2386     if (get_bits1(&s->gb)) {
2387         for(i=0;i<64;i++) {
2388             v = get_bits(&s->gb, 8);
2389             j= s->dsp.idct_permutation[ ff_zigzag_direct[i] ];
2390             s->intra_matrix[j] = v;
2391             s->chroma_intra_matrix[j] = v;
2392         }
2393     }
2394     if (get_bits1(&s->gb)) {
2395         for(i=0;i<64;i++) {
2396             v = get_bits(&s->gb, 8);
2397             j= s->dsp.idct_permutation[ ff_zigzag_direct[i] ];
2398             s->inter_matrix[j] = v;
2399             s->chroma_inter_matrix[j] = v;
2400         }
2401     }
2402     if (get_bits1(&s->gb)) {
2403         for(i=0;i<64;i++) {
2404             v = get_bits(&s->gb, 8);
2405             j= s->dsp.idct_permutation[ ff_zigzag_direct[i] ];
2406             s->chroma_intra_matrix[j] = v;
2407         }
2408     }
2409     if (get_bits1(&s->gb)) {
2410         for(i=0;i<64;i++) {
2411             v = get_bits(&s->gb, 8);
2412             j= s->dsp.idct_permutation[ ff_zigzag_direct[i] ];
2413             s->chroma_inter_matrix[j] = v;
2414         }
2415     }
2416 }
2417
2418 static void mpeg_decode_picture_coding_extension(MpegEncContext *s)
2419 {
2420     s->full_pel[0] = s->full_pel[1] = 0;
2421     s->mpeg_f_code[0][0] = get_bits(&s->gb, 4);
2422     s->mpeg_f_code[0][1] = get_bits(&s->gb, 4);
2423     s->mpeg_f_code[1][0] = get_bits(&s->gb, 4);
2424     s->mpeg_f_code[1][1] = get_bits(&s->gb, 4);
2425     s->intra_dc_precision = get_bits(&s->gb, 2);
2426     s->picture_structure = get_bits(&s->gb, 2);
2427     s->top_field_first = get_bits1(&s->gb);
2428     s->frame_pred_frame_dct = get_bits1(&s->gb);
2429     s->concealment_motion_vectors = get_bits1(&s->gb);
2430     s->q_scale_type = get_bits1(&s->gb);
2431     s->intra_vlc_format = get_bits1(&s->gb);
2432     s->alternate_scan = get_bits1(&s->gb);
2433     s->repeat_first_field = get_bits1(&s->gb);
2434     s->chroma_420_type = get_bits1(&s->gb);
2435     s->progressive_frame = get_bits1(&s->gb);
2436
2437     if(s->picture_structure == PICT_FRAME){
2438         s->first_field=0;
2439         s->v_edge_pos= 16*s->mb_height;
2440     }else{
2441         s->first_field ^= 1;
2442         s->v_edge_pos=  8*s->mb_height;
2443         memset(s->mbskip_table, 0, s->mb_stride*s->mb_height);
2444     }
2445
2446     if(s->alternate_scan){
2447         ff_init_scantable(s->dsp.idct_permutation, &s->inter_scantable  , ff_alternate_vertical_scan);
2448         ff_init_scantable(s->dsp.idct_permutation, &s->intra_scantable  , ff_alternate_vertical_scan);
2449     }else{
2450         ff_init_scantable(s->dsp.idct_permutation, &s->inter_scantable  , ff_zigzag_direct);
2451         ff_init_scantable(s->dsp.idct_permutation, &s->intra_scantable  , ff_zigzag_direct);
2452     }
2453
2454     /* composite display not parsed */
2455     dprintf(s->avctx, "intra_dc_precision=%d\n", s->intra_dc_precision);
2456     dprintf(s->avctx, "picture_structure=%d\n", s->picture_structure);
2457     dprintf(s->avctx, "top field first=%d\n", s->top_field_first);
2458     dprintf(s->avctx, "repeat first field=%d\n", s->repeat_first_field);
2459     dprintf(s->avctx, "conceal=%d\n", s->concealment_motion_vectors);
2460     dprintf(s->avctx, "intra_vlc_format=%d\n", s->intra_vlc_format);
2461     dprintf(s->avctx, "alternate_scan=%d\n", s->alternate_scan);
2462     dprintf(s->avctx, "frame_pred_frame_dct=%d\n", s->frame_pred_frame_dct);
2463     dprintf(s->avctx, "progressive_frame=%d\n", s->progressive_frame);
2464 }
2465
2466 static void mpeg_decode_extension(AVCodecContext *avctx,
2467                                   const uint8_t *buf, int buf_size)
2468 {
2469     Mpeg1Context *s1 = avctx->priv_data;
2470     MpegEncContext *s = &s1->mpeg_enc_ctx;
2471     int ext_type;
2472
2473     init_get_bits(&s->gb, buf, buf_size*8);
2474
2475     ext_type = get_bits(&s->gb, 4);
2476     switch(ext_type) {
2477     case 0x1:
2478         mpeg_decode_sequence_extension(s1);
2479         break;
2480     case 0x2:
2481         mpeg_decode_sequence_display_extension(s1);
2482         break;
2483     case 0x3:
2484         mpeg_decode_quant_matrix_extension(s);
2485         break;
2486     case 0x7:
2487         mpeg_decode_picture_display_extension(s1);
2488         break;
2489     case 0x8:
2490         mpeg_decode_picture_coding_extension(s);
2491         break;
2492     }
2493 }
2494
2495 static void exchange_uv(MpegEncContext *s){
2496     short * tmp = s->pblocks[4];
2497     s->pblocks[4] = s->pblocks[5];
2498     s->pblocks[5] = tmp;
2499 }
2500
2501 static int mpeg_field_start(MpegEncContext *s){
2502     AVCodecContext *avctx= s->avctx;
2503     Mpeg1Context *s1 = (Mpeg1Context*)s;
2504
2505     /* start frame decoding */
2506     if(s->first_field || s->picture_structure==PICT_FRAME){
2507         if(MPV_frame_start(s, avctx) < 0)
2508             return -1;
2509
2510         ff_er_frame_start(s);
2511
2512         /* first check if we must repeat the frame */
2513         s->current_picture_ptr->repeat_pict = 0;
2514         if (s->repeat_first_field) {
2515             if (s->progressive_sequence) {
2516                 if (s->top_field_first)
2517                     s->current_picture_ptr->repeat_pict = 4;
2518                 else
2519                     s->current_picture_ptr->repeat_pict = 2;
2520             } else if (s->progressive_frame) {
2521                 s->current_picture_ptr->repeat_pict = 1;
2522             }
2523         }
2524
2525         *s->current_picture_ptr->pan_scan= s1->pan_scan;
2526     }else{ //second field
2527             int i;
2528
2529             if(!s->current_picture_ptr){
2530                 av_log(s->avctx, AV_LOG_ERROR, "first field missing\n");
2531                 return -1;
2532             }
2533
2534             for(i=0; i<4; i++){
2535                 s->current_picture.data[i] = s->current_picture_ptr->data[i];
2536                 if(s->picture_structure == PICT_BOTTOM_FIELD){
2537                     s->current_picture.data[i] += s->current_picture_ptr->linesize[i];
2538                 }
2539             }
2540     }
2541 #ifdef HAVE_XVMC
2542 // MPV_frame_start will call this function too,
2543 // but we need to call it on every field
2544     if(s->avctx->xvmc_acceleration)
2545          XVMC_field_start(s,avctx);
2546 #endif
2547
2548     return 0;
2549 }
2550
2551 #define DECODE_SLICE_ERROR -1
2552 #define DECODE_SLICE_OK 0
2553
2554 /**
2555  * decodes a slice. MpegEncContext.mb_y must be set to the MB row from the startcode
2556  * @return DECODE_SLICE_ERROR if the slice is damaged<br>
2557  *         DECODE_SLICE_OK if this slice is ok<br>
2558  */
2559 static int mpeg_decode_slice(Mpeg1Context *s1, int mb_y,
2560                              const uint8_t **buf, int buf_size)
2561 {
2562     MpegEncContext *s = &s1->mpeg_enc_ctx;
2563     AVCodecContext *avctx= s->avctx;
2564     int ret;
2565     const int field_pic= s->picture_structure != PICT_FRAME;
2566     const int lowres= s->avctx->lowres;
2567
2568     s->resync_mb_x=
2569     s->resync_mb_y= -1;
2570
2571     if (mb_y<<field_pic >= s->mb_height){
2572         av_log(s->avctx, AV_LOG_ERROR, "slice below image (%d >= %d)\n", mb_y, s->mb_height);
2573         return -1;
2574     }
2575
2576     init_get_bits(&s->gb, *buf, buf_size*8);
2577
2578     ff_mpeg1_clean_buffers(s);
2579     s->interlaced_dct = 0;
2580
2581     s->qscale = get_qscale(s);
2582
2583     if(s->qscale == 0){
2584         av_log(s->avctx, AV_LOG_ERROR, "qscale == 0\n");
2585         return -1;
2586     }
2587
2588     /* extra slice info */
2589     while (get_bits1(&s->gb) != 0) {
2590         skip_bits(&s->gb, 8);
2591     }
2592
2593     s->mb_x=0;
2594
2595     for(;;) {
2596         int code = get_vlc2(&s->gb, mbincr_vlc.table, MBINCR_VLC_BITS, 2);
2597         if (code < 0){
2598             av_log(s->avctx, AV_LOG_ERROR, "first mb_incr damaged\n");
2599             return -1;
2600         }
2601         if (code >= 33) {
2602             if (code == 33) {
2603                 s->mb_x += 33;
2604             }
2605             /* otherwise, stuffing, nothing to do */
2606         } else {
2607             s->mb_x += code;
2608             break;
2609         }
2610     }
2611     if(s->mb_x >= (unsigned)s->mb_width){
2612         av_log(s->avctx, AV_LOG_ERROR, "initial skip overflow\n");
2613         return -1;
2614     }
2615
2616     s->resync_mb_x= s->mb_x;
2617     s->resync_mb_y= s->mb_y= mb_y;
2618     s->mb_skip_run= 0;
2619     ff_init_block_index(s);
2620
2621     if (s->mb_y==0 && s->mb_x==0 && (s->first_field || s->picture_structure==PICT_FRAME)) {
2622         if(s->avctx->debug&FF_DEBUG_PICT_INFO){
2623              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",
2624                  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],
2625                  s->pict_type == I_TYPE ? "I" : (s->pict_type == P_TYPE ? "P" : (s->pict_type == B_TYPE ? "B" : "S")),
2626                  s->progressive_sequence ? "ps" :"", s->progressive_frame ? "pf" : "", s->alternate_scan ? "alt" :"", s->top_field_first ? "top" :"",
2627                  s->intra_dc_precision, s->picture_structure, s->frame_pred_frame_dct, s->concealment_motion_vectors,
2628                  s->q_scale_type, s->intra_vlc_format, s->repeat_first_field, s->chroma_420_type ? "420" :"");
2629         }
2630     }
2631
2632     for(;;) {
2633 #ifdef HAVE_XVMC
2634         //one 1 we memcpy blocks in xvmcvideo
2635         if(s->avctx->xvmc_acceleration > 1)
2636             XVMC_init_block(s);//set s->block
2637 #endif
2638
2639         ret = mpeg_decode_mb(s, s->block);
2640         s->chroma_qscale= s->qscale;
2641
2642         dprintf(s->avctx, "ret=%d\n", ret);
2643         if (ret < 0)
2644             return -1;
2645
2646         if(s->current_picture.motion_val[0] && !s->encoding){ //note motion_val is normally NULL unless we want to extract the MVs
2647             const int wrap = field_pic ? 2*s->b8_stride : s->b8_stride;
2648             int xy = s->mb_x*2 + s->mb_y*2*wrap;
2649             int motion_x, motion_y, dir, i;
2650             if(field_pic && !s->first_field)
2651                 xy += wrap/2;
2652
2653             for(i=0; i<2; i++){
2654                 for(dir=0; dir<2; dir++){
2655                     if (s->mb_intra || (dir==1 && s->pict_type != B_TYPE)) {
2656                         motion_x = motion_y = 0;
2657                     }else if (s->mv_type == MV_TYPE_16X16 || (s->mv_type == MV_TYPE_FIELD && field_pic)){
2658                         motion_x = s->mv[dir][0][0];
2659                         motion_y = s->mv[dir][0][1];
2660                     } else /*if ((s->mv_type == MV_TYPE_FIELD) || (s->mv_type == MV_TYPE_16X8))*/ {
2661                         motion_x = s->mv[dir][i][0];
2662                         motion_y = s->mv[dir][i][1];
2663                     }
2664
2665                     s->current_picture.motion_val[dir][xy    ][0] = motion_x;
2666                     s->current_picture.motion_val[dir][xy    ][1] = motion_y;
2667                     s->current_picture.motion_val[dir][xy + 1][0] = motion_x;
2668                     s->current_picture.motion_val[dir][xy + 1][1] = motion_y;
2669                     s->current_picture.ref_index [dir][xy    ]=
2670                     s->current_picture.ref_index [dir][xy + 1]= s->field_select[dir][i];
2671                     assert(s->field_select[dir][i]==0 || s->field_select[dir][i]==1);
2672                 }
2673                 xy += wrap;
2674             }
2675         }
2676
2677         s->dest[0] += 16 >> lowres;
2678         s->dest[1] += 16 >> (s->chroma_x_shift + lowres);
2679         s->dest[2] += 16 >> (s->chroma_x_shift + lowres);
2680
2681         MPV_decode_mb(s, s->block);
2682
2683         if (++s->mb_x >= s->mb_width) {
2684             const int mb_size= 16>>s->avctx->lowres;
2685
2686             ff_draw_horiz_band(s, mb_size*s->mb_y, mb_size);
2687
2688             s->mb_x = 0;
2689             s->mb_y++;
2690
2691             if(s->mb_y<<field_pic >= s->mb_height){
2692                 int left= s->gb.size_in_bits - get_bits_count(&s->gb);
2693                 int is_d10= s->chroma_format==2 && s->pict_type==I_TYPE && avctx->profile==0 && avctx->level==5
2694                             && s->intra_dc_precision == 2 && s->q_scale_type == 1 && s->alternate_scan == 0
2695                             && s->progressive_frame == 0 /* vbv_delay == 0xBBB || 0xE10*/;
2696
2697                 if(left < 0 || (left && show_bits(&s->gb, FFMIN(left, 23)) && !is_d10)
2698                    || (avctx->error_resilience >= FF_ER_AGGRESSIVE && left>8)){
2699                     av_log(avctx, AV_LOG_ERROR, "end mismatch left=%d %0X\n", left, show_bits(&s->gb, FFMIN(left, 23)));
2700                     return -1;
2701                 }else
2702                     goto eos;
2703             }
2704
2705             ff_init_block_index(s);
2706         }
2707
2708         /* skip mb handling */
2709         if (s->mb_skip_run == -1) {
2710             /* read again increment */
2711             s->mb_skip_run = 0;
2712             for(;;) {
2713                 int code = get_vlc2(&s->gb, mbincr_vlc.table, MBINCR_VLC_BITS, 2);
2714                 if (code < 0){
2715                     av_log(s->avctx, AV_LOG_ERROR, "mb incr damaged\n");
2716                     return -1;
2717                 }
2718                 if (code >= 33) {
2719                     if (code == 33) {
2720                         s->mb_skip_run += 33;
2721                     }else if(code == 35){
2722                         if(s->mb_skip_run != 0 || show_bits(&s->gb, 15) != 0){
2723                             av_log(s->avctx, AV_LOG_ERROR, "slice mismatch\n");
2724                             return -1;
2725                         }
2726                         goto eos; /* end of slice */
2727                     }
2728                     /* otherwise, stuffing, nothing to do */
2729                 } else {
2730                     s->mb_skip_run += code;
2731                     break;
2732                 }
2733             }
2734         }
2735     }
2736 eos: // end of slice
2737     *buf += get_bits_count(&s->gb)/8 - 1;
2738 //printf("y %d %d %d %d\n", s->resync_mb_x, s->resync_mb_y, s->mb_x, s->mb_y);
2739     return 0;
2740 }
2741
2742 static int slice_decode_thread(AVCodecContext *c, void *arg){
2743     MpegEncContext *s= arg;
2744     const uint8_t *buf= s->gb.buffer;
2745     int mb_y= s->start_mb_y;
2746
2747     s->error_count= 3*(s->end_mb_y - s->start_mb_y)*s->mb_width;
2748
2749     for(;;){
2750         uint32_t start_code;
2751         int ret;
2752
2753         ret= mpeg_decode_slice((Mpeg1Context*)s, mb_y, &buf, s->gb.buffer_end - buf);
2754         emms_c();
2755 //av_log(c, AV_LOG_DEBUG, "ret:%d resync:%d/%d mb:%d/%d ts:%d/%d ec:%d\n",
2756 //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);
2757         if(ret < 0){
2758             if(s->resync_mb_x>=0 && s->resync_mb_y>=0)
2759                 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);
2760         }else{
2761             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);
2762         }
2763
2764         if(s->mb_y == s->end_mb_y)
2765             return 0;
2766
2767         start_code= -1;
2768         buf = ff_find_start_code(buf, s->gb.buffer_end, &start_code);
2769         mb_y= start_code - SLICE_MIN_START_CODE;
2770         if(mb_y < 0 || mb_y >= s->end_mb_y)
2771             return -1;
2772     }
2773
2774     return 0; //not reached
2775 }
2776
2777 /**
2778  * handles slice ends.
2779  * @return 1 if it seems to be the last slice of
2780  */
2781 static int slice_end(AVCodecContext *avctx, AVFrame *pict)
2782 {
2783     Mpeg1Context *s1 = avctx->priv_data;
2784     MpegEncContext *s = &s1->mpeg_enc_ctx;
2785
2786     if (!s1->mpeg_enc_ctx_allocated || !s->current_picture_ptr)
2787         return 0;
2788
2789 #ifdef HAVE_XVMC
2790     if(s->avctx->xvmc_acceleration)
2791         XVMC_field_end(s);
2792 #endif
2793     /* end of slice reached */
2794     if (/*s->mb_y<<field_pic == s->mb_height &&*/ !s->first_field) {
2795         /* end of image */
2796
2797         s->current_picture_ptr->qscale_type= FF_QSCALE_TYPE_MPEG2;
2798
2799         ff_er_frame_end(s);
2800
2801         MPV_frame_end(s);
2802
2803         if (s->pict_type == B_TYPE || s->low_delay) {
2804             *pict= *(AVFrame*)s->current_picture_ptr;
2805             ff_print_debug_info(s, pict);
2806         } else {
2807             s->picture_number++;
2808             /* latency of 1 frame for I and P frames */
2809             /* XXX: use another variable than picture_number */
2810             if (s->last_picture_ptr != NULL) {
2811                 *pict= *(AVFrame*)s->last_picture_ptr;
2812                  ff_print_debug_info(s, pict);
2813             }
2814         }
2815
2816         return 1;
2817     } else {
2818         return 0;
2819     }
2820 }
2821
2822 static int mpeg1_decode_sequence(AVCodecContext *avctx,
2823                                  const uint8_t *buf, int buf_size)
2824 {
2825     Mpeg1Context *s1 = avctx->priv_data;
2826     MpegEncContext *s = &s1->mpeg_enc_ctx;
2827     int width,height;
2828     int i, v, j;
2829
2830     init_get_bits(&s->gb, buf, buf_size*8);
2831
2832     width = get_bits(&s->gb, 12);
2833     height = get_bits(&s->gb, 12);
2834     if (width <= 0 || height <= 0 ||
2835         (width % 2) != 0 || (height % 2) != 0)
2836         return -1;
2837     s->aspect_ratio_info= get_bits(&s->gb, 4);
2838     if (s->aspect_ratio_info == 0)
2839         return -1;
2840     s->frame_rate_index = get_bits(&s->gb, 4);
2841     if (s->frame_rate_index == 0 || s->frame_rate_index > 13)
2842         return -1;
2843     s->bit_rate = get_bits(&s->gb, 18) * 400;
2844     if (get_bits1(&s->gb) == 0) /* marker */
2845         return -1;
2846     s->width = width;
2847     s->height = height;
2848
2849     s->avctx->rc_buffer_size= get_bits(&s->gb, 10) * 1024*16;
2850     skip_bits(&s->gb, 1);
2851
2852     /* get matrix */
2853     if (get_bits1(&s->gb)) {
2854         for(i=0;i<64;i++) {
2855             v = get_bits(&s->gb, 8);
2856             if(v==0){
2857                 av_log(s->avctx, AV_LOG_ERROR, "intra matrix damaged\n");
2858                 return -1;
2859             }
2860             j = s->dsp.idct_permutation[ ff_zigzag_direct[i] ];
2861             s->intra_matrix[j] = v;
2862             s->chroma_intra_matrix[j] = v;
2863         }
2864 #ifdef DEBUG
2865         dprintf(s->avctx, "intra matrix present\n");
2866         for(i=0;i<64;i++)
2867             dprintf(s->avctx, " %d", s->intra_matrix[s->dsp.idct_permutation[i]]);
2868         dprintf(s->avctx, "\n");
2869 #endif
2870     } else {
2871         for(i=0;i<64;i++) {
2872             j = s->dsp.idct_permutation[i];
2873             v = ff_mpeg1_default_intra_matrix[i];
2874             s->intra_matrix[j] = v;
2875             s->chroma_intra_matrix[j] = v;
2876         }
2877     }
2878     if (get_bits1(&s->gb)) {
2879         for(i=0;i<64;i++) {
2880             v = get_bits(&s->gb, 8);
2881             if(v==0){
2882                 av_log(s->avctx, AV_LOG_ERROR, "inter matrix damaged\n");
2883                 return -1;
2884             }
2885             j = s->dsp.idct_permutation[ ff_zigzag_direct[i] ];
2886             s->inter_matrix[j] = v;
2887             s->chroma_inter_matrix[j] = v;
2888         }
2889 #ifdef DEBUG
2890         dprintf(s->avctx, "non intra matrix present\n");
2891         for(i=0;i<64;i++)
2892             dprintf(s->avctx, " %d", s->inter_matrix[s->dsp.idct_permutation[i]]);
2893         dprintf(s->avctx, "\n");
2894 #endif
2895     } else {
2896         for(i=0;i<64;i++) {
2897             int j= s->dsp.idct_permutation[i];
2898             v = ff_mpeg1_default_non_intra_matrix[i];
2899             s->inter_matrix[j] = v;
2900             s->chroma_inter_matrix[j] = v;
2901         }
2902     }
2903
2904     if(show_bits(&s->gb, 23) != 0){
2905         av_log(s->avctx, AV_LOG_ERROR, "sequence header damaged\n");
2906         return -1;
2907     }
2908
2909     /* we set mpeg2 parameters so that it emulates mpeg1 */
2910     s->progressive_sequence = 1;
2911     s->progressive_frame = 1;
2912     s->picture_structure = PICT_FRAME;
2913     s->frame_pred_frame_dct = 1;
2914     s->chroma_format = 1;
2915     s->codec_id= s->avctx->codec_id= CODEC_ID_MPEG1VIDEO;
2916     avctx->sub_id = 1; /* indicates mpeg1 */
2917     s->out_format = FMT_MPEG1;
2918     s->swap_uv = 0;//AFAIK VCR2 don't have SEQ_HEADER
2919     if(s->flags & CODEC_FLAG_LOW_DELAY) s->low_delay=1;
2920
2921     if(s->avctx->debug & FF_DEBUG_PICT_INFO)
2922         av_log(s->avctx, AV_LOG_DEBUG, "vbv buffer: %d, bitrate:%d\n",
2923                s->avctx->rc_buffer_size, s->bit_rate);
2924
2925     return 0;
2926 }
2927
2928 static int vcr2_init_sequence(AVCodecContext *avctx)
2929 {
2930     Mpeg1Context *s1 = avctx->priv_data;
2931     MpegEncContext *s = &s1->mpeg_enc_ctx;
2932     int i, v;
2933
2934     /* start new mpeg1 context decoding */
2935     s->out_format = FMT_MPEG1;
2936     if (s1->mpeg_enc_ctx_allocated) {
2937         MPV_common_end(s);
2938     }
2939     s->width  = avctx->coded_width;
2940     s->height = avctx->coded_height;
2941     avctx->has_b_frames= 0; //true?
2942     s->low_delay= 1;
2943
2944     if(avctx->xvmc_acceleration){
2945         avctx->pix_fmt = avctx->get_format(avctx,pixfmt_xvmc_mpg2_420);
2946     }else{
2947         avctx->pix_fmt = avctx->get_format(avctx,pixfmt_yuv_420);
2948     }
2949
2950     if( avctx->pix_fmt == PIX_FMT_XVMC_MPEG2_IDCT )
2951         if( avctx->idct_algo == FF_IDCT_AUTO )
2952             avctx->idct_algo = FF_IDCT_SIMPLE;
2953
2954     if (MPV_common_init(s) < 0)
2955         return -1;
2956     exchange_uv(s);//common init reset pblocks, so we swap them here
2957     s->swap_uv = 1;// in case of xvmc we need to swap uv for each MB
2958     s1->mpeg_enc_ctx_allocated = 1;
2959
2960     for(i=0;i<64;i++) {
2961         int j= s->dsp.idct_permutation[i];
2962         v = ff_mpeg1_default_intra_matrix[i];
2963         s->intra_matrix[j] = v;
2964         s->chroma_intra_matrix[j] = v;
2965
2966         v = ff_mpeg1_default_non_intra_matrix[i];
2967         s->inter_matrix[j] = v;
2968         s->chroma_inter_matrix[j] = v;
2969     }
2970
2971     s->progressive_sequence = 1;
2972     s->progressive_frame = 1;
2973     s->picture_structure = PICT_FRAME;
2974     s->frame_pred_frame_dct = 1;
2975     s->chroma_format = 1;
2976     s->codec_id= s->avctx->codec_id= CODEC_ID_MPEG2VIDEO;
2977     avctx->sub_id = 2; /* indicates mpeg2 */
2978     return 0;
2979 }
2980
2981
2982 static void mpeg_decode_user_data(AVCodecContext *avctx,
2983                                   const uint8_t *buf, int buf_size)
2984 {
2985     const uint8_t *p;
2986     int len, flags;
2987     p = buf;
2988     len = buf_size;
2989
2990     /* we parse the DTG active format information */
2991     if (len >= 5 &&
2992         p[0] == 'D' && p[1] == 'T' && p[2] == 'G' && p[3] == '1') {
2993         flags = p[4];
2994         p += 5;
2995         len -= 5;
2996         if (flags & 0x80) {
2997             /* skip event id */
2998             if (len < 2)
2999                 return;
3000             p += 2;
3001             len -= 2;
3002         }
3003         if (flags & 0x40) {
3004             if (len < 1)
3005                 return;
3006             avctx->dtg_active_format = p[0] & 0x0f;
3007         }
3008     }
3009 }
3010
3011 static void mpeg_decode_gop(AVCodecContext *avctx,
3012                             const uint8_t *buf, int buf_size){
3013     Mpeg1Context *s1 = avctx->priv_data;
3014     MpegEncContext *s = &s1->mpeg_enc_ctx;
3015
3016     int drop_frame_flag;
3017     int time_code_hours, time_code_minutes;
3018     int time_code_seconds, time_code_pictures;
3019     int broken_link;
3020
3021     init_get_bits(&s->gb, buf, buf_size*8);
3022
3023     drop_frame_flag = get_bits1(&s->gb);
3024
3025     time_code_hours=get_bits(&s->gb,5);
3026     time_code_minutes = get_bits(&s->gb,6);
3027     skip_bits1(&s->gb);//marker bit
3028     time_code_seconds = get_bits(&s->gb,6);
3029     time_code_pictures = get_bits(&s->gb,6);
3030
3031     /*broken_link indicate that after editing the
3032       reference frames of the first B-Frames after GOP I-Frame
3033       are missing (open gop)*/
3034     broken_link = get_bits1(&s->gb);
3035
3036     if(s->avctx->debug & FF_DEBUG_PICT_INFO)
3037         av_log(s->avctx, AV_LOG_DEBUG, "GOP (%2d:%02d:%02d.[%02d]) broken_link=%d\n",
3038             time_code_hours, time_code_minutes, time_code_seconds,
3039             time_code_pictures, broken_link);
3040 }
3041 /**
3042  * finds the end of the current frame in the bitstream.
3043  * @return the position of the first byte of the next frame, or -1
3044  */
3045 static int mpeg1_find_frame_end(ParseContext *pc, const uint8_t *buf, int buf_size)
3046 {
3047     int i;
3048     uint32_t state= pc->state;
3049
3050     i=0;
3051     if(!pc->frame_start_found){
3052         for(i=0; i<buf_size; i++){
3053             i= ff_find_start_code(buf+i, buf+buf_size, &state) - buf - 1;
3054             if(state >= SLICE_MIN_START_CODE && state <= SLICE_MAX_START_CODE){
3055                 i++;
3056                 pc->frame_start_found=1;
3057                 break;
3058             }
3059         }
3060     }
3061
3062     if(pc->frame_start_found){
3063         /* EOF considered as end of frame */
3064         if (buf_size == 0)
3065             return 0;
3066         for(; i<buf_size; i++){
3067             i= ff_find_start_code(buf+i, buf+buf_size, &state) - buf - 1;
3068             if((state&0xFFFFFF00) == 0x100){
3069                 if(state < SLICE_MIN_START_CODE || state > SLICE_MAX_START_CODE){
3070                     pc->frame_start_found=0;
3071                     pc->state=-1;
3072                     return i-3;
3073                 }
3074             }
3075         }
3076     }
3077     pc->state= state;
3078     return END_NOT_FOUND;
3079 }
3080
3081 /* handle buffering and image synchronisation */
3082 static int mpeg_decode_frame(AVCodecContext *avctx,
3083                              void *data, int *data_size,
3084                              uint8_t *buf, int buf_size)
3085 {
3086     Mpeg1Context *s = avctx->priv_data;
3087     const uint8_t *buf_end;
3088     const uint8_t *buf_ptr;
3089     uint32_t start_code;
3090     int ret, input_size;
3091     AVFrame *picture = data;
3092     MpegEncContext *s2 = &s->mpeg_enc_ctx;
3093     dprintf(avctx, "fill_buffer\n");
3094
3095     if (buf_size == 0) {
3096         /* special case for last picture */
3097         if (s2->low_delay==0 && s2->next_picture_ptr) {
3098             *picture= *(AVFrame*)s2->next_picture_ptr;
3099             s2->next_picture_ptr= NULL;
3100
3101             *data_size = sizeof(AVFrame);
3102         }
3103         return 0;
3104     }
3105
3106     if(s2->flags&CODEC_FLAG_TRUNCATED){
3107         int next= mpeg1_find_frame_end(&s2->parse_context, buf, buf_size);
3108
3109         if( ff_combine_frame(&s2->parse_context, next, &buf, &buf_size) < 0 )
3110             return buf_size;
3111     }
3112
3113     buf_ptr = buf;
3114     buf_end = buf + buf_size;
3115
3116 #if 0
3117     if (s->repeat_field % 2 == 1) {
3118         s->repeat_field++;
3119         //fprintf(stderr,"\nRepeating last frame: %d -> %d! pict: %d %d", avctx->frame_number-1, avctx->frame_number,
3120         //        s2->picture_number, s->repeat_field);
3121         if (avctx->flags & CODEC_FLAG_REPEAT_FIELD) {
3122             *data_size = sizeof(AVPicture);
3123             goto the_end;
3124         }
3125     }
3126 #endif
3127
3128     if(s->mpeg_enc_ctx_allocated==0 && avctx->codec_tag == ff_get_fourcc("VCR2"))
3129         vcr2_init_sequence(avctx);
3130
3131     s->slice_count= 0;
3132
3133     for(;;) {
3134         /* find start next code */
3135         start_code = -1;
3136         buf_ptr = ff_find_start_code(buf_ptr,buf_end, &start_code);
3137         if (start_code > 0x1ff){
3138             if(s2->pict_type != B_TYPE || avctx->skip_frame <= AVDISCARD_DEFAULT){
3139                 if(avctx->thread_count > 1){
3140                     int i;
3141
3142                     avctx->execute(avctx, slice_decode_thread,  (void**)&(s2->thread_context[0]), NULL, s->slice_count);
3143                     for(i=0; i<s->slice_count; i++)
3144                         s2->error_count += s2->thread_context[i]->error_count;
3145                 }
3146                 if (slice_end(avctx, picture)) {
3147                     if(s2->last_picture_ptr || s2->low_delay) //FIXME merge with the stuff in mpeg_decode_slice
3148                         *data_size = sizeof(AVPicture);
3149                 }
3150             }
3151             return FFMAX(0, buf_ptr - buf - s2->parse_context.last_index);
3152         }
3153
3154         input_size = buf_end - buf_ptr;
3155
3156         if(avctx->debug & FF_DEBUG_STARTCODE){
3157             av_log(avctx, AV_LOG_DEBUG, "%3X at %zd left %d\n", start_code, buf_ptr-buf, input_size);
3158         }
3159
3160         /* prepare data for next start code */
3161         switch(start_code) {
3162         case SEQ_START_CODE:
3163             mpeg1_decode_sequence(avctx, buf_ptr,
3164                                     input_size);
3165             break;
3166
3167         case PICTURE_START_CODE:
3168             /* we have a complete image : we try to decompress it */
3169             mpeg1_decode_picture(avctx,
3170                                     buf_ptr, input_size);
3171             break;
3172         case EXT_START_CODE:
3173             mpeg_decode_extension(avctx,
3174                                     buf_ptr, input_size);
3175             break;
3176         case USER_START_CODE:
3177             mpeg_decode_user_data(avctx,
3178                                     buf_ptr, input_size);
3179             break;
3180         case GOP_START_CODE:
3181             s2->first_field=0;
3182             mpeg_decode_gop(avctx,
3183                                     buf_ptr, input_size);
3184             break;
3185         default:
3186             if (start_code >= SLICE_MIN_START_CODE &&
3187                 start_code <= SLICE_MAX_START_CODE) {
3188                 int mb_y= start_code - SLICE_MIN_START_CODE;
3189
3190                 if(s2->last_picture_ptr==NULL){
3191                 /* skip b frames if we dont have reference frames */
3192                     if(s2->pict_type==B_TYPE) break;
3193                 /* skip P frames if we dont have reference frame no valid header */
3194 //                    if(s2->pict_type==P_TYPE && s2->first_field && !s2->first_slice) break;
3195                 }
3196                 /* skip b frames if we are in a hurry */
3197                 if(avctx->hurry_up && s2->pict_type==B_TYPE) break;
3198                 if(  (avctx->skip_frame >= AVDISCARD_NONREF && s2->pict_type==B_TYPE)
3199                     ||(avctx->skip_frame >= AVDISCARD_NONKEY && s2->pict_type!=I_TYPE)
3200                     || avctx->skip_frame >= AVDISCARD_ALL)
3201                     break;
3202                 /* skip everything if we are in a hurry>=5 */
3203                 if(avctx->hurry_up>=5) break;
3204
3205                 if (!s->mpeg_enc_ctx_allocated) break;
3206
3207                 if(s2->codec_id == CODEC_ID_MPEG2VIDEO){
3208                     if(mb_y < avctx->skip_top || mb_y >= s2->mb_height - avctx->skip_bottom)
3209                         break;
3210                 }
3211
3212                 if(s2->first_slice){
3213                     s2->first_slice=0;
3214                             if(mpeg_field_start(s2) < 0)
3215                         return -1;
3216                     }
3217                 if(!s2->current_picture_ptr){
3218                     av_log(avctx, AV_LOG_ERROR, "current_picture not initalized\n");
3219                     return -1;
3220                 }
3221
3222                 if(avctx->thread_count > 1){
3223                     int threshold= (s2->mb_height*s->slice_count + avctx->thread_count/2) / avctx->thread_count;
3224                     if(threshold <= mb_y){
3225                         MpegEncContext *thread_context= s2->thread_context[s->slice_count];
3226
3227                         thread_context->start_mb_y= mb_y;
3228                         thread_context->end_mb_y  = s2->mb_height;
3229                         if(s->slice_count){
3230                             s2->thread_context[s->slice_count-1]->end_mb_y= mb_y;
3231                             ff_update_duplicate_context(thread_context, s2);
3232                         }
3233                         init_get_bits(&thread_context->gb, buf_ptr, input_size*8);
3234                         s->slice_count++;
3235                     }
3236                     buf_ptr += 2; //FIXME add minimum num of bytes per slice
3237                 }else{
3238                     ret = mpeg_decode_slice(s, mb_y, &buf_ptr, input_size);
3239                     emms_c();
3240
3241                     if(ret < 0){
3242                         if(s2->resync_mb_x>=0 && s2->resync_mb_y>=0)
3243                             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);
3244                     }else{
3245                         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);
3246                     }
3247                 }
3248             }
3249             break;
3250         }
3251     }
3252 }
3253
3254 static int mpeg_decode_end(AVCodecContext *avctx)
3255 {
3256     Mpeg1Context *s = avctx->priv_data;
3257
3258     if (s->mpeg_enc_ctx_allocated)
3259         MPV_common_end(&s->mpeg_enc_ctx);
3260     return 0;
3261 }
3262
3263 AVCodec mpeg1video_decoder = {
3264     "mpeg1video",
3265     CODEC_TYPE_VIDEO,
3266     CODEC_ID_MPEG1VIDEO,
3267     sizeof(Mpeg1Context),
3268     mpeg_decode_init,
3269     NULL,
3270     mpeg_decode_end,
3271     mpeg_decode_frame,
3272     CODEC_CAP_DRAW_HORIZ_BAND | CODEC_CAP_DR1 | CODEC_CAP_TRUNCATED | CODEC_CAP_DELAY,
3273     .flush= ff_mpeg_flush,
3274 };
3275
3276 AVCodec mpeg2video_decoder = {
3277     "mpeg2video",
3278     CODEC_TYPE_VIDEO,
3279     CODEC_ID_MPEG2VIDEO,
3280     sizeof(Mpeg1Context),
3281     mpeg_decode_init,
3282     NULL,
3283     mpeg_decode_end,
3284     mpeg_decode_frame,
3285     CODEC_CAP_DRAW_HORIZ_BAND | CODEC_CAP_DR1 | CODEC_CAP_TRUNCATED | CODEC_CAP_DELAY,
3286     .flush= ff_mpeg_flush,
3287 };
3288
3289 //legacy decoder
3290 AVCodec mpegvideo_decoder = {
3291     "mpegvideo",
3292     CODEC_TYPE_VIDEO,
3293     CODEC_ID_MPEG2VIDEO,
3294     sizeof(Mpeg1Context),
3295     mpeg_decode_init,
3296     NULL,
3297     mpeg_decode_end,
3298     mpeg_decode_frame,
3299     CODEC_CAP_DRAW_HORIZ_BAND | CODEC_CAP_DR1 | CODEC_CAP_TRUNCATED | CODEC_CAP_DELAY,
3300     .flush= ff_mpeg_flush,
3301 };
3302
3303 #ifdef CONFIG_ENCODERS
3304
3305 AVCodec mpeg1video_encoder = {
3306     "mpeg1video",
3307     CODEC_TYPE_VIDEO,
3308     CODEC_ID_MPEG1VIDEO,
3309     sizeof(MpegEncContext),
3310     encode_init,
3311     MPV_encode_picture,
3312     MPV_encode_end,
3313     .supported_framerates= ff_frame_rate_tab+1,
3314     .pix_fmts= (enum PixelFormat[]){PIX_FMT_YUV420P, -1},
3315     .capabilities= CODEC_CAP_DELAY,
3316 };
3317
3318 AVCodec mpeg2video_encoder = {
3319     "mpeg2video",
3320     CODEC_TYPE_VIDEO,
3321     CODEC_ID_MPEG2VIDEO,
3322     sizeof(MpegEncContext),
3323     encode_init,
3324     MPV_encode_picture,
3325     MPV_encode_end,
3326     .supported_framerates= ff_frame_rate_tab+1,
3327     .pix_fmts= (enum PixelFormat[]){PIX_FMT_YUV420P, PIX_FMT_YUV422P, -1},
3328     .capabilities= CODEC_CAP_DELAY,
3329 };
3330 #endif
3331
3332 #ifdef HAVE_XVMC
3333 static int mpeg_mc_decode_init(AVCodecContext *avctx){
3334     Mpeg1Context *s;
3335
3336     if( avctx->thread_count > 1)
3337         return -1;
3338     if( !(avctx->slice_flags & SLICE_FLAG_CODED_ORDER) )
3339         return -1;
3340     if( !(avctx->slice_flags & SLICE_FLAG_ALLOW_FIELD) ){
3341         dprintf(avctx, "mpeg12.c: XvMC decoder will work better if SLICE_FLAG_ALLOW_FIELD is set\n");
3342     }
3343     mpeg_decode_init(avctx);
3344     s = avctx->priv_data;
3345
3346     avctx->pix_fmt = PIX_FMT_XVMC_MPEG2_IDCT;
3347     avctx->xvmc_acceleration = 2;//2 - the blocks are packed!
3348
3349     return 0;
3350 }
3351
3352 AVCodec mpeg_xvmc_decoder = {
3353     "mpegvideo_xvmc",
3354     CODEC_TYPE_VIDEO,
3355     CODEC_ID_MPEG2VIDEO_XVMC,
3356     sizeof(Mpeg1Context),
3357     mpeg_mc_decode_init,
3358     NULL,
3359     mpeg_decode_end,
3360     mpeg_decode_frame,
3361     CODEC_CAP_DRAW_HORIZ_BAND | CODEC_CAP_DR1 | CODEC_CAP_TRUNCATED| CODEC_CAP_HWACCEL | CODEC_CAP_DELAY,
3362     .flush= ff_mpeg_flush,
3363 };
3364
3365 #endif
3366
3367 #ifdef CONFIG_MPEGVIDEO_PARSER
3368 static void mpegvideo_extract_headers(AVCodecParserContext *s,
3369                                       AVCodecContext *avctx,
3370                                       const uint8_t *buf, int buf_size)
3371 {
3372     ParseContext1 *pc = s->priv_data;
3373     const uint8_t *buf_end;
3374     uint32_t start_code;
3375     int frame_rate_index, ext_type, bytes_left;
3376     int frame_rate_ext_n, frame_rate_ext_d;
3377     int picture_structure, top_field_first, repeat_first_field, progressive_frame;
3378     int horiz_size_ext, vert_size_ext, bit_rate_ext;
3379 //FIXME replace the crap with get_bits()
3380     s->repeat_pict = 0;
3381     buf_end = buf + buf_size;
3382     while (buf < buf_end) {
3383         start_code= -1;
3384         buf= ff_find_start_code(buf, buf_end, &start_code);
3385         bytes_left = buf_end - buf;
3386         switch(start_code) {
3387         case PICTURE_START_CODE:
3388             if (bytes_left >= 2) {
3389                 s->pict_type = (buf[1] >> 3) & 7;
3390             }
3391             break;
3392         case SEQ_START_CODE:
3393             if (bytes_left >= 7) {
3394                 pc->width  = (buf[0] << 4) | (buf[1] >> 4);
3395                 pc->height = ((buf[1] & 0x0f) << 8) | buf[2];
3396                 avcodec_set_dimensions(avctx, pc->width, pc->height);
3397                 frame_rate_index = buf[3] & 0xf;
3398                 pc->frame_rate.den = avctx->time_base.den = ff_frame_rate_tab[frame_rate_index].num;
3399                 pc->frame_rate.num = avctx->time_base.num = ff_frame_rate_tab[frame_rate_index].den;
3400                 avctx->bit_rate = ((buf[4]<<10) | (buf[5]<<2) | (buf[6]>>6))*400;
3401                 avctx->codec_id = CODEC_ID_MPEG1VIDEO;
3402                 avctx->sub_id = 1;
3403             }
3404             break;
3405         case EXT_START_CODE:
3406             if (bytes_left >= 1) {
3407                 ext_type = (buf[0] >> 4);
3408                 switch(ext_type) {
3409                 case 0x1: /* sequence extension */
3410                     if (bytes_left >= 6) {
3411                         horiz_size_ext = ((buf[1] & 1) << 1) | (buf[2] >> 7);
3412                         vert_size_ext = (buf[2] >> 5) & 3;
3413                         bit_rate_ext = ((buf[2] & 0x1F)<<7) | (buf[3]>>1);
3414                         frame_rate_ext_n = (buf[5] >> 5) & 3;
3415                         frame_rate_ext_d = (buf[5] & 0x1f);
3416                         pc->progressive_sequence = buf[1] & (1 << 3);
3417                         avctx->has_b_frames= !(buf[5] >> 7);
3418
3419                         pc->width  |=(horiz_size_ext << 12);
3420                         pc->height |=( vert_size_ext << 12);
3421                         avctx->bit_rate += (bit_rate_ext << 18) * 400;
3422                         avcodec_set_dimensions(avctx, pc->width, pc->height);
3423                         avctx->time_base.den = pc->frame_rate.den * (frame_rate_ext_n + 1);
3424                         avctx->time_base.num = pc->frame_rate.num * (frame_rate_ext_d + 1);
3425                         avctx->codec_id = CODEC_ID_MPEG2VIDEO;
3426                         avctx->sub_id = 2; /* forces MPEG2 */
3427                     }
3428                     break;
3429                 case 0x8: /* picture coding extension */
3430                     if (bytes_left >= 5) {
3431                         picture_structure = buf[2]&3;
3432                         top_field_first = buf[3] & (1 << 7);
3433                         repeat_first_field = buf[3] & (1 << 1);
3434                         progressive_frame = buf[4] & (1 << 7);
3435
3436                         /* check if we must repeat the frame */
3437                         if (repeat_first_field) {
3438                             if (pc->progressive_sequence) {
3439                                 if (top_field_first)
3440                                     s->repeat_pict = 4;
3441                                 else
3442                                     s->repeat_pict = 2;
3443                             } else if (progressive_frame) {
3444                                 s->repeat_pict = 1;
3445                             }
3446                         }
3447
3448                         /* the packet only represents half a frame
3449                            XXX,FIXME maybe find a different solution */
3450                         if(picture_structure != 3)
3451                             s->repeat_pict = -1;
3452                     }
3453                     break;
3454                 }
3455             }
3456             break;
3457         case -1:
3458             goto the_end;
3459         default:
3460             /* we stop parsing when we encounter a slice. It ensures
3461                that this function takes a negligible amount of time */
3462             if (start_code >= SLICE_MIN_START_CODE &&
3463                 start_code <= SLICE_MAX_START_CODE)
3464                 goto the_end;
3465             break;
3466         }
3467     }
3468  the_end: ;
3469 }
3470
3471 static int mpegvideo_parse(AVCodecParserContext *s,
3472                            AVCodecContext *avctx,
3473                            uint8_t **poutbuf, int *poutbuf_size,
3474                            const uint8_t *buf, int buf_size)
3475 {
3476     ParseContext1 *pc1 = s->priv_data;
3477     ParseContext *pc= &pc1->pc;
3478     int next;
3479
3480     if(s->flags & PARSER_FLAG_COMPLETE_FRAMES){
3481         next= buf_size;
3482     }else{
3483         next= mpeg1_find_frame_end(pc, buf, buf_size);
3484
3485         if (ff_combine_frame(pc, next, (uint8_t **)&buf, &buf_size) < 0) {
3486             *poutbuf = NULL;
3487             *poutbuf_size = 0;
3488             return buf_size;
3489         }
3490
3491     }
3492     /* we have a full frame : we just parse the first few MPEG headers
3493        to have the full timing information. The time take by this
3494        function should be negligible for uncorrupted streams */
3495     mpegvideo_extract_headers(s, avctx, buf, buf_size);
3496 #if 0
3497     printf("pict_type=%d frame_rate=%0.3f repeat_pict=%d\n",
3498            s->pict_type, (double)avctx->time_base.den / avctx->time_base.num, s->repeat_pict);
3499 #endif
3500
3501     *poutbuf = (uint8_t *)buf;
3502     *poutbuf_size = buf_size;
3503     return next;
3504 }
3505
3506 static int mpegvideo_split(AVCodecContext *avctx,
3507                            const uint8_t *buf, int buf_size)
3508 {
3509     int i;
3510     uint32_t state= -1;
3511
3512     for(i=0; i<buf_size; i++){
3513         state= (state<<8) | buf[i];
3514         if(state != 0x1B3 && state != 0x1B5 && state < 0x200 && state >= 0x100)
3515             return i-3;
3516     }
3517     return 0;
3518 }
3519
3520 AVCodecParser mpegvideo_parser = {
3521     { CODEC_ID_MPEG1VIDEO, CODEC_ID_MPEG2VIDEO },
3522     sizeof(ParseContext1),
3523     NULL,
3524     mpegvideo_parse,
3525     ff_parse1_close,
3526     mpegvideo_split,
3527 };
3528 #endif /* !CONFIG_MPEGVIDEO_PARSER */
3529
3530 static int imx_dump_header(AVBitStreamFilterContext *bsfc, AVCodecContext *avctx, const char *args,
3531                            uint8_t **poutbuf, int *poutbuf_size,
3532                            const uint8_t *buf, int buf_size, int keyframe)
3533 {
3534     /* MXF essence element key */
3535     static const uint8_t imx_header[16] = { 0x06,0x0e,0x2b,0x34,0x01,0x02,0x01,0x01,0x0d,0x01,0x03,0x01,0x05,0x01,0x01,0x00 };
3536     uint8_t *poutbufp;
3537
3538     if (avctx->codec_id != CODEC_ID_MPEG2VIDEO) {
3539         av_log(avctx, AV_LOG_ERROR, "imx bitstream filter only applies to mpeg2video codec\n");
3540         return 0;
3541     }
3542
3543     *poutbuf = av_malloc(buf_size + 20 + FF_INPUT_BUFFER_PADDING_SIZE);
3544     poutbufp = *poutbuf;
3545     bytestream_put_buffer(&poutbufp, imx_header, 16);
3546     bytestream_put_byte(&poutbufp, 0x83); /* KLV BER long form */
3547     bytestream_put_be24(&poutbufp, buf_size);
3548     bytestream_put_buffer(&poutbufp, buf, buf_size);
3549     *poutbuf_size = poutbufp - *poutbuf;
3550     return 1;
3551 }
3552
3553 AVBitStreamFilter imx_dump_header_bsf = {
3554     "imxdump",
3555     0,
3556     imx_dump_header,
3557 };
3558
3559 /* this is ugly i know, but the alternative is too make
3560    hundreds of vars global and prefix them with ff_mpeg1_
3561    which is far uglier. */
3562 #include "mdec.c"