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