mpeg12.c 94.1 KB
Newer Older
F
Fabrice Bellard 已提交
1
/*
M
Michael Niedermayer 已提交
2
 * MPEG1 codec / MPEG2 decoder
F
Fabrice Bellard 已提交
3
 * Copyright (c) 2000,2001 Fabrice Bellard.
F
Fabrice Bellard 已提交
4
 *
F
Fabrice Bellard 已提交
5 6 7 8
 * This library is free software; you can redistribute it and/or
 * modify it under the terms of the GNU Lesser General Public
 * License as published by the Free Software Foundation; either
 * version 2 of the License, or (at your option) any later version.
F
Fabrice Bellard 已提交
9
 *
F
Fabrice Bellard 已提交
10
 * This library is distributed in the hope that it will be useful,
F
Fabrice Bellard 已提交
11
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
F
Fabrice Bellard 已提交
12 13
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
 * Lesser General Public License for more details.
F
Fabrice Bellard 已提交
14
 *
F
Fabrice Bellard 已提交
15 16 17
 * You should have received a copy of the GNU Lesser General Public
 * License along with this library; if not, write to the Free Software
 * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
F
Fabrice Bellard 已提交
18
 */
M
Michael Niedermayer 已提交
19 20 21
 
/**
 * @file mpeg12.c
M
Michael Niedermayer 已提交
22
 * MPEG1/2 codec
M
Michael Niedermayer 已提交
23 24
 */
 
F
Fabrice Bellard 已提交
25
//#define DEBUG
F
Fabrice Bellard 已提交
26 27 28 29 30 31
#include "avcodec.h"
#include "dsputil.h"
#include "mpegvideo.h"

#include "mpeg12data.h"

32

F
Fabrice Bellard 已提交
33 34 35 36 37 38 39 40 41 42
/* Start codes. */
#define SEQ_END_CODE		0x000001b7
#define SEQ_START_CODE		0x000001b3
#define GOP_START_CODE		0x000001b8
#define PICTURE_START_CODE	0x00000100
#define SLICE_MIN_START_CODE	0x00000101
#define SLICE_MAX_START_CODE	0x000001af
#define EXT_START_CODE		0x000001b5
#define USER_START_CODE		0x000001b2

43 44 45 46 47 48 49 50
#define DC_VLC_BITS 9
#define MV_VLC_BITS 9
#define MBINCR_VLC_BITS 9
#define MB_PAT_VLC_BITS 9
#define MB_PTYPE_VLC_BITS 6
#define MB_BTYPE_VLC_BITS 6
#define TEX_VLC_BITS 9

51
#ifdef CONFIG_ENCODERS
F
Fabrice Bellard 已提交
52 53 54
static void mpeg1_encode_block(MpegEncContext *s, 
                         DCTELEM *block, 
                         int component);
55
static void mpeg1_encode_motion(MpegEncContext *s, int val, int f_or_b_code);    // RAL: f_code parameter added
56
#endif //CONFIG_ENCODERS
F
Fabrice Bellard 已提交
57
static void mpeg1_skip_picture(MpegEncContext *s, int pict_num);
58 59 60 61
static inline int mpeg1_decode_block_inter(MpegEncContext *s, 
                              DCTELEM *block, 
                              int n);
static inline int mpeg1_decode_block_intra(MpegEncContext *s, 
F
Fabrice Bellard 已提交
62 63
                              DCTELEM *block, 
                              int n);
64
static inline int mpeg2_decode_block_non_intra(MpegEncContext *s, 
F
Fabrice Bellard 已提交
65 66
                                        DCTELEM *block, 
                                        int n);
67
static inline int mpeg2_decode_block_intra(MpegEncContext *s, 
F
Fabrice Bellard 已提交
68 69 70
                                    DCTELEM *block, 
                                    int n);
static int mpeg_decode_motion(MpegEncContext *s, int fcode, int pred);
I
Ivan Kalvachev 已提交
71
static void exchange_uv(MpegEncContext *s);
F
Fabrice Bellard 已提交
72

I
Ivan Kalvachev 已提交
73 74 75
#ifdef HAVE_XVMC
extern int XVMC_field_start(MpegEncContext *s, AVCodecContext *avctx);
extern int XVMC_field_end(MpegEncContext *s);
76
extern void XVMC_pack_pblocks(MpegEncContext *s,int cbp);
I
Ivan Kalvachev 已提交
77
extern void XVMC_init_block(MpegEncContext *s);//set s->block
I
Ivan Kalvachev 已提交
78 79
#endif

80
#ifdef CONFIG_ENCODERS
81
static uint8_t (*mv_penalty)[MAX_MV*2+1]= NULL;
82
static uint8_t fcode_tab[MAX_MV*2+1];
83

84 85
static uint32_t uni_mpeg1_ac_vlc_bits[64*64*2];
static uint8_t  uni_mpeg1_ac_vlc_len [64*64*2];
M
Michael Niedermayer 已提交
86 87 88 89 90 91 92

/* simple include everything table for dc, first byte is bits number next 3 are code*/
static uint32_t mpeg1_lum_dc_uni[512];
static uint32_t mpeg1_chr_dc_uni[512];

static uint8_t mpeg1_index_run[2][64];
static int8_t mpeg1_max_level[2][64];
93
#endif //CONFIG_ENCODERS
94

95 96
static void init_2d_vlc_rl(RLTable *rl)
{
97
    int i;
98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120
    
    init_vlc(&rl->vlc, TEX_VLC_BITS, rl->n + 2, 
             &rl->table_vlc[0][1], 4, 2,
             &rl->table_vlc[0][0], 4, 2);

    
    rl->rl_vlc[0]= av_malloc(rl->vlc.table_size*sizeof(RL_VLC_ELEM));
    for(i=0; i<rl->vlc.table_size; i++){
        int code= rl->vlc.table[i][0];
        int len = rl->vlc.table[i][1];
        int level, run;
    
        if(len==0){ // illegal code
            run= 65;
            level= MAX_LEVEL;
        }else if(len<0){ //more bits needed
            run= 0;
            level= code;
        }else{
            if(code==rl->n){ //esc
                run= 65;
                level= 0;
            }else if(code==rl->n+1){ //eob
121 122
                run= 0;
                level= 127;
123 124 125 126 127 128 129 130 131 132 133
            }else{
                run=   rl->table_run  [code] + 1;
                level= rl->table_level[code];
            }
        }
        rl->rl_vlc[0][i].len= len;
        rl->rl_vlc[0][i].level= level;
        rl->rl_vlc[0][i].run= run;
    }
}

M
Michael Niedermayer 已提交
134
#ifdef CONFIG_ENCODERS
135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179
static void init_uni_ac_vlc(RLTable *rl, uint32_t *uni_ac_vlc_bits, uint8_t *uni_ac_vlc_len){
    int i;

    for(i=0; i<128; i++){
        int level= i-64;
        int run;
        for(run=0; run<64; run++){
            int len, bits, code;
            
            int alevel= ABS(level);
            int sign= (level>>31)&1;

            if (alevel > rl->max_level[0][run])
                code= 111; /*rl->n*/
            else
                code= rl->index_run[0][run] + alevel - 1;

            if (code < 111 /* rl->n */) {
	    	/* store the vlc & sign at once */
                len=   mpeg1_vlc[code][1]+1;
                bits= (mpeg1_vlc[code][0]<<1) + sign;
            } else {
                len=  mpeg1_vlc[111/*rl->n*/][1]+6;
                bits= mpeg1_vlc[111/*rl->n*/][0]<<6;

                bits|= run;
                if (alevel < 128) {
                    bits<<=8; len+=8;
                    bits|= level & 0xff;
                } else {
                    bits<<=16; len+=16;
                    bits|= level & 0xff;
                    if (level < 0) {
                        bits|= 0x8001 + level + 255;
                    } else {
                        bits|= level & 0xffff;
                    }
                }
            }

            uni_ac_vlc_bits[UNI_AC_ENC_INDEX(run, i)]= bits;
            uni_ac_vlc_len [UNI_AC_ENC_INDEX(run, i)]= len;
        }
    }
}
180

F
Fabrice Bellard 已提交
181 182 183
static void put_header(MpegEncContext *s, int header)
{
    align_put_bits(&s->pb);
184 185
    put_bits(&s->pb, 16, header>>16);
    put_bits(&s->pb, 16, header&0xFFFF);
F
Fabrice Bellard 已提交
186 187 188 189 190 191
}

/* put sequence header if needed */
static void mpeg1_encode_sequence_header(MpegEncContext *s)
{
        unsigned int vbv_buffer_size;
F
Fabrice Bellard 已提交
192
        unsigned int fps, v;
193
        int n, i;
194
        uint64_t time_code;
195
        float best_aspect_error= 1E10;
M
Michael Niedermayer 已提交
196
        float aspect_ratio= av_q2d(s->avctx->sample_aspect_ratio);
197
        int constraint_parameter_flag;
198
        
M
Michael Niedermayer 已提交
199
        if(aspect_ratio==0.0) aspect_ratio= 1.0; //pixel aspect 1:1 (VGA)
F
Fabrice Bellard 已提交
200
        
M
cleanup  
Michael Niedermayer 已提交
201
        if (s->current_picture.key_frame) {
F
Fabrice Bellard 已提交
202 203 204 205 206 207 208 209
            /* mpeg1 header repeated every gop */
            put_header(s, SEQ_START_CODE);
            
            /* search closest frame rate */
            {
                int i, dmin, d;
                s->frame_rate_index = 0;
                dmin = 0x7fffffff;
210 211 212 213
                for(i=1;i<14;i++) {
                    if(s->avctx->strict_std_compliance >= 0 && i>=9) break;
                     
                    d = abs(MPEG1_FRAME_RATE_BASE*(int64_t)s->avctx->frame_rate/s->avctx->frame_rate_base - frame_rate_tab[i]);
F
Fabrice Bellard 已提交
214 215 216 217 218 219 220 221 222
                    if (d < dmin) {
                        dmin = d;
                        s->frame_rate_index = i;
                    }
                }
            }
 
            put_bits(&s->pb, 12, s->width);
            put_bits(&s->pb, 12, s->height);
223 224
            
            for(i=1; i<15; i++){
225 226
                float error= aspect_ratio;
                if(s->codec_id == CODEC_ID_MPEG1VIDEO || i <=1)
227
                    error-= 1.0/mpeg1_aspect[i];
228 229 230
                else
                    error-= av_q2d(mpeg2_aspect[i])*s->height/s->width;
             
231 232 233 234 235 236 237 238 239
                error= ABS(error);
                
                if(error < best_aspect_error){
                    best_aspect_error= error;
                    s->aspect_ratio_info= i;
                }
            }
            
            put_bits(&s->pb, 4, s->aspect_ratio_info);
F
Fabrice Bellard 已提交
240
            put_bits(&s->pb, 4, s->frame_rate_index);
241 242 243 244 245 246 247 248
            
            if(s->avctx->rc_max_rate){
                v = (s->avctx->rc_max_rate + 399) / 400;
                if (v > 0x3ffff && s->codec_id == CODEC_ID_MPEG1VIDEO)
                    v = 0x3ffff;
            }else{
                v= 0x3FFFF;
            }
249 250 251 252 253

            if(s->avctx->rc_buffer_size)
                vbv_buffer_size = s->avctx->rc_buffer_size;
            else
                /* VBV calculation: Scaled so that a VCD has the proper VBV size of 40 kilobytes */
254 255 256 257 258
                vbv_buffer_size = (( 20 * s->bit_rate) / (1151929 / 2)) * 8 * 1024;
            vbv_buffer_size= (vbv_buffer_size + 16383) / 16384;

            put_bits(&s->pb, 18, v & 0x3FFFF);
            put_bits(&s->pb, 1, 1); /* marker */
259 260 261 262 263 264 265 266 267 268 269 270
            put_bits(&s->pb, 10, vbv_buffer_size & 0x3FF);

            constraint_parameter_flag= 
                s->width <= 768 && s->height <= 576 && 
                s->mb_width * s->mb_height <= 396 &&
                s->mb_width * s->mb_height * frame_rate_tab[s->frame_rate_index] <= MPEG1_FRAME_RATE_BASE*396*25 &&
                frame_rate_tab[s->frame_rate_index] <= MPEG1_FRAME_RATE_BASE*30 &&
                vbv_buffer_size <= 20 &&
                v <= 1856000/400 &&
                s->codec_id == CODEC_ID_MPEG1VIDEO;
                
            put_bits(&s->pb, 1, constraint_parameter_flag);
271 272 273
            
            ff_write_quant_matrix(&s->pb, s->avctx->intra_matrix);
            ff_write_quant_matrix(&s->pb, s->avctx->inter_matrix);
F
Fabrice Bellard 已提交
274

M
Michael Niedermayer 已提交
275 276 277 278 279 280
            if(s->codec_id == CODEC_ID_MPEG2VIDEO){
                put_header(s, EXT_START_CODE);
                put_bits(&s->pb, 4, 1); //seq ext
                put_bits(&s->pb, 1, 0); //esc
                put_bits(&s->pb, 3, 4); //profile
                put_bits(&s->pb, 4, 8); //level
281
                put_bits(&s->pb, 1, s->progressive_sequence);
M
Michael Niedermayer 已提交
282 283 284 285 286 287 288 289 290 291 292
                put_bits(&s->pb, 2, 1); //chroma format 4:2:0
                put_bits(&s->pb, 2, 0); //horizontal size ext
                put_bits(&s->pb, 2, 0); //vertical size ext
                put_bits(&s->pb, 12, v>>18); //bitrate ext
                put_bits(&s->pb, 1, 1); //marker
                put_bits(&s->pb, 8, vbv_buffer_size >>10); //vbv buffer ext
                put_bits(&s->pb, 1, s->low_delay);
                put_bits(&s->pb, 2, 0); // frame_rate_ext_n
                put_bits(&s->pb, 5, 0); // frame_rate_ext_d
            }
            
F
Fabrice Bellard 已提交
293 294 295 296 297
            put_header(s, GOP_START_CODE);
            put_bits(&s->pb, 1, 0); /* do drop frame */
            /* time code : we must convert from the real frame rate to a
               fake mpeg frame rate in case of low frame rate */
            fps = frame_rate_tab[s->frame_rate_index];
298
            time_code = (int64_t)s->fake_picture_number * MPEG1_FRAME_RATE_BASE;
F
Fabrice Bellard 已提交
299
            s->gop_picture_number = s->fake_picture_number;
300 301
            put_bits(&s->pb, 5, (uint32_t)((time_code / (fps * 3600)) % 24));
            put_bits(&s->pb, 6, (uint32_t)((time_code / (fps * 60)) % 60));
F
Fabrice Bellard 已提交
302
            put_bits(&s->pb, 1, 1);
303
            put_bits(&s->pb, 6, (uint32_t)((time_code / fps) % 60));
304
            put_bits(&s->pb, 6, (uint32_t)((time_code % fps) / MPEG1_FRAME_RATE_BASE));
M
Michael Niedermayer 已提交
305
            put_bits(&s->pb, 1, 0); /* closed gop */
F
Fabrice Bellard 已提交
306 307 308
            put_bits(&s->pb, 1, 0); /* broken link */
        }

309
        if (s->avctx->frame_rate < (24 * s->avctx->frame_rate_base) && s->picture_number > 0) {
F
Fabrice Bellard 已提交
310 311 312
            /* insert empty P pictures to slow down to the desired
               frame rate. Each fake pictures takes about 20 bytes */
            fps = frame_rate_tab[s->frame_rate_index];
313
            n = av_rescale((int64_t)s->picture_number * s->avctx->frame_rate_base, fps, s->avctx->frame_rate) / MPEG1_FRAME_RATE_BASE - 1;
F
Fabrice Bellard 已提交
314 315 316 317 318 319 320 321 322
            while (s->fake_picture_number < n) {
                mpeg1_skip_picture(s, s->fake_picture_number - 
                                   s->gop_picture_number); 
                s->fake_picture_number++;
            }

        }
}

323 324 325 326 327 328 329 330
static inline void encode_mb_skip_run(MpegEncContext *s, int run){
    while (run >= 33) {
        put_bits(&s->pb, 11, 0x008);
        run -= 33;
    }
    put_bits(&s->pb, mbAddrIncrTable[run][1], 
             mbAddrIncrTable[run][0]);
}
F
Fabrice Bellard 已提交
331 332 333 334

/* insert a fake P picture */
static void mpeg1_skip_picture(MpegEncContext *s, int pict_num)
{
M
Michael Niedermayer 已提交
335 336
    assert(s->codec_id == CODEC_ID_MPEG1VIDEO); // mpeg2 can do these repeat things

F
Fabrice Bellard 已提交
337 338 339 340 341 342 343 344 345 346 347 348 349 350 351 352 353 354
    /* mpeg1 picture header */
    put_header(s, PICTURE_START_CODE);
    /* temporal reference */
    put_bits(&s->pb, 10, pict_num & 0x3ff); 
    
    put_bits(&s->pb, 3, P_TYPE);
    put_bits(&s->pb, 16, 0xffff); /* non constant bit rate */
    
    put_bits(&s->pb, 1, 1); /* integer coordinates */
    put_bits(&s->pb, 3, 1); /* forward_f_code */
    
    put_bits(&s->pb, 1, 0); /* extra bit picture */
    
    /* only one slice */
    put_header(s, SLICE_MIN_START_CODE);
    put_bits(&s->pb, 5, 1); /* quantizer scale */
    put_bits(&s->pb, 1, 0); /* slice extra information */
    
355
    encode_mb_skip_run(s, 0);
F
Fabrice Bellard 已提交
356 357 358 359 360 361 362 363 364
    
    /* empty macroblock */
    put_bits(&s->pb, 3, 1); /* motion only */
    
    /* zero motion x & y */
    put_bits(&s->pb, 1, 1); 
    put_bits(&s->pb, 1, 1); 

    /* output a number of empty slice */
365
    encode_mb_skip_run(s, s->mb_width * s->mb_height - 2);
F
Fabrice Bellard 已提交
366 367 368 369 370 371 372 373
    
    /* empty macroblock */
    put_bits(&s->pb, 3, 1); /* motion only */
    
    /* zero motion x & y */
    put_bits(&s->pb, 1, 1); 
    put_bits(&s->pb, 1, 1); 
}
374
#endif //CONFIG_ENCODERS
F
Fabrice Bellard 已提交
375

376
static void common_init(MpegEncContext *s)
F
Fabrice Bellard 已提交
377
{
378 379 380
    s->y_dc_scale_table=
    s->c_dc_scale_table= ff_mpeg1_dc_scale_table;
}
381

M
Michael Niedermayer 已提交
382 383 384 385 386 387 388
void ff_mpeg1_clean_buffers(MpegEncContext *s){
    s->last_dc[0] = 1 << (7 + s->intra_dc_precision);
    s->last_dc[1] = s->last_dc[0];
    s->last_dc[2] = s->last_dc[0];
    memset(s->last_mv, 0, sizeof(s->last_mv));
}

389
#ifdef CONFIG_ENCODERS
390 391 392 393 394 395 396

void ff_mpeg1_encode_slice_header(MpegEncContext *s){
    put_header(s, SLICE_MIN_START_CODE + s->mb_y);
    put_bits(&s->pb, 5, s->qscale); /* quantizer scale */
    put_bits(&s->pb, 1, 0); /* slice extra information */
}

397 398
void mpeg1_encode_picture_header(MpegEncContext *s, int picture_number)
{
F
Fabrice Bellard 已提交
399 400 401 402 403
    mpeg1_encode_sequence_header(s);

    /* mpeg1 picture header */
    put_header(s, PICTURE_START_CODE);
    /* temporal reference */
404 405 406

    // RAL: s->picture_number instead of s->fake_picture_number
    put_bits(&s->pb, 10, (s->picture_number - 
F
Fabrice Bellard 已提交
407
                          s->gop_picture_number) & 0x3ff); 
408
    s->fake_picture_number++;
F
Fabrice Bellard 已提交
409 410
    
    put_bits(&s->pb, 3, s->pict_type);
M
Michael Niedermayer 已提交
411 412 413

    s->vbv_delay_ptr= s->pb.buf + get_bit_count(&s->pb)/8;
    put_bits(&s->pb, 16, 0xFFFF); /* vbv_delay */
F
Fabrice Bellard 已提交
414
    
415 416
    // RAL: Forward f_code also needed for B frames
    if (s->pict_type == P_TYPE || s->pict_type == B_TYPE) {
F
Fabrice Bellard 已提交
417
        put_bits(&s->pb, 1, 0); /* half pel coordinates */
418 419 420 421
        if(s->codec_id == CODEC_ID_MPEG1VIDEO)
            put_bits(&s->pb, 3, s->f_code); /* forward_f_code */
        else
            put_bits(&s->pb, 3, 7); /* forward_f_code */
F
Fabrice Bellard 已提交
422 423
    }
    
424 425 426
    // RAL: Backward f_code necessary for B frames
    if (s->pict_type == B_TYPE) {
        put_bits(&s->pb, 1, 0); /* half pel coordinates */
427 428 429 430
        if(s->codec_id == CODEC_ID_MPEG1VIDEO)
            put_bits(&s->pb, 3, s->b_code); /* backward_f_code */
        else
            put_bits(&s->pb, 3, 7); /* backward_f_code */
M
Michael Niedermayer 已提交
431
    }
432

F
Fabrice Bellard 已提交
433
    put_bits(&s->pb, 1, 0); /* extra bit picture */
434 435

    s->frame_pred_frame_dct = 1;
M
Michael Niedermayer 已提交
436 437 438
    if(s->codec_id == CODEC_ID_MPEG2VIDEO){
        put_header(s, EXT_START_CODE);
        put_bits(&s->pb, 4, 8); //pic ext
439 440 441 442 443 444 445 446 447 448 449 450
        if (s->pict_type == P_TYPE || s->pict_type == B_TYPE) {
            put_bits(&s->pb, 4, s->f_code);
            put_bits(&s->pb, 4, s->f_code);
        }else{
            put_bits(&s->pb, 8, 255);
        }
        if (s->pict_type == B_TYPE) {
            put_bits(&s->pb, 4, s->b_code);
            put_bits(&s->pb, 4, s->b_code);
        }else{
            put_bits(&s->pb, 8, 255);
        }
M
Michael Niedermayer 已提交
451 452
        put_bits(&s->pb, 2, s->intra_dc_precision);
        put_bits(&s->pb, 2, s->picture_structure= PICT_FRAME);
453 454 455 456 457
        if (s->progressive_sequence) {
            put_bits(&s->pb, 1, 0); /* no repeat */
        } else {
            put_bits(&s->pb, 1, s->current_picture_ptr->top_field_first);
        }
458 459 460 461 462
        /* XXX: optimize the generation of this flag with entropy
           measures */
        s->frame_pred_frame_dct = s->progressive_sequence;
        
        put_bits(&s->pb, 1, s->frame_pred_frame_dct);
M
Michael Niedermayer 已提交
463 464 465 466 467 468
        put_bits(&s->pb, 1, s->concealment_motion_vectors);
        put_bits(&s->pb, 1, s->q_scale_type);
        put_bits(&s->pb, 1, s->intra_vlc_format);
        put_bits(&s->pb, 1, s->alternate_scan);
        put_bits(&s->pb, 1, s->repeat_first_field);
        put_bits(&s->pb, 1, s->chroma_420_type=1);
469 470
        s->progressive_frame = s->progressive_sequence;
        put_bits(&s->pb, 1, s->progressive_frame);
M
Michael Niedermayer 已提交
471 472 473
        put_bits(&s->pb, 1, 0); //composite_display_flag
    }
    
474 475
    s->mb_y=0;
    ff_mpeg1_encode_slice_header(s);
F
Fabrice Bellard 已提交
476 477
}

478 479 480 481 482 483 484 485 486 487 488
static inline void put_mb_modes(MpegEncContext *s, int n, int bits, 
                                int has_mv)
{
    put_bits(&s->pb, n, bits);
    if (!s->frame_pred_frame_dct) {
        if (has_mv) 
            put_bits(&s->pb, 2, 2); /* motion_type: frame */
        put_bits(&s->pb, 1, s->interlaced_dct);
    }
}

F
Fabrice Bellard 已提交
489 490 491 492
void mpeg1_encode_mb(MpegEncContext *s,
                     DCTELEM block[6][64],
                     int motion_x, int motion_y)
{
493 494 495 496
    int i, cbp;
    const int mb_x = s->mb_x;
    const int mb_y = s->mb_y;
    const int first_mb= mb_x == s->resync_mb_x && mb_y == s->resync_mb_y;
F
Fabrice Bellard 已提交
497 498 499 500 501 502 503 504

    /* compute cbp */
    cbp = 0;
    for(i=0;i<6;i++) {
        if (s->block_last_index[i] >= 0)
            cbp |= 1 << (5 - i);
    }

M
Michael Niedermayer 已提交
505
    if (cbp == 0 && !first_mb && (mb_x != s->mb_width - 1 || (mb_y != s->mb_height - 1 && s->codec_id == CODEC_ID_MPEG1VIDEO)) && 
506 507 508
        ((s->pict_type == P_TYPE && (motion_x | motion_y) == 0) ||
        (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) |
        ((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))) {
509
        s->mb_skip_run++;
510
        s->qscale -= s->dquant;
M
Michael Niedermayer 已提交
511 512 513
        s->skip_count++;
        s->misc_bits++;
        s->last_bits++;
F
Fabrice Bellard 已提交
514
    } else {
515 516 517 518 519
        if(first_mb){
            assert(s->mb_skip_run == 0);
            encode_mb_skip_run(s, s->mb_x);
        }else{
            encode_mb_skip_run(s, s->mb_skip_run);
F
Fabrice Bellard 已提交
520 521 522
        }
        
        if (s->pict_type == I_TYPE) {
523
            if(s->dquant && cbp){
524
                put_mb_modes(s, 2, 1, 0); /* macroblock_type : macroblock_quant = 1 */
525 526
                put_bits(&s->pb, 5, s->qscale);
            }else{
527
                put_mb_modes(s, 1, 1, 0); /* macroblock_type : macroblock_quant = 0 */
528 529
                s->qscale -= s->dquant;
            }
M
Michael Niedermayer 已提交
530 531
            s->misc_bits+= get_bits_diff(s);
            s->i_count++;
532 533
        } else if (s->mb_intra) {
            if(s->dquant && cbp){
534
                put_mb_modes(s, 6, 0x01, 0);
535 536
                put_bits(&s->pb, 5, s->qscale);
            }else{
537
                put_mb_modes(s, 5, 0x03, 0);
538 539 540 541 542 543 544
                s->qscale -= s->dquant;
            }
            s->misc_bits+= get_bits_diff(s);
            s->i_count++;
            s->last_mv[0][0][0] = 
            s->last_mv[0][0][1] = 0;
        } else if (s->pict_type == P_TYPE) { 
F
Fabrice Bellard 已提交
545 546
                if (cbp != 0) {
                    if (motion_x == 0 && motion_y == 0) {
547
                        if(s->dquant){
548
                            put_mb_modes(s, 5, 1, 0); /* macroblock_pattern & quant */
549 550
                            put_bits(&s->pb, 5, s->qscale);
                        }else{
551
                            put_mb_modes(s, 2, 1, 0); /* macroblock_pattern only */
552
                        }
M
Michael Niedermayer 已提交
553
                        s->misc_bits+= get_bits_diff(s);
F
Fabrice Bellard 已提交
554 555
                        put_bits(&s->pb, mbPatTable[cbp - 1][1], mbPatTable[cbp - 1][0]);
                    } else {
556
                        if(s->dquant){
557
                            put_mb_modes(s, 5, 2, 1); /* motion + cbp */
558 559
                            put_bits(&s->pb, 5, s->qscale);
                        }else{
560
                            put_mb_modes(s, 1, 1, 1); /* motion + cbp */
561
                        }
M
Michael Niedermayer 已提交
562
                        s->misc_bits+= get_bits_diff(s);
563 564
                        mpeg1_encode_motion(s, motion_x - s->last_mv[0][0][0], s->f_code);    // RAL: f_code parameter added
                        mpeg1_encode_motion(s, motion_y - s->last_mv[0][0][1], s->f_code);    // RAL: f_code parameter added
M
Michael Niedermayer 已提交
565
                        s->mv_bits+= get_bits_diff(s);
F
Fabrice Bellard 已提交
566 567 568 569
                        put_bits(&s->pb, mbPatTable[cbp - 1][1], mbPatTable[cbp - 1][0]);
                    }
                } else {
                    put_bits(&s->pb, 3, 1); /* motion only */
570 571
                    if (!s->frame_pred_frame_dct)
                        put_bits(&s->pb, 2, 2); /* motion_type: frame */
572 573
                    mpeg1_encode_motion(s, motion_x - s->last_mv[0][0][0], s->f_code);    // RAL: f_code parameter added
                    mpeg1_encode_motion(s, motion_y - s->last_mv[0][0][1], s->f_code);    // RAL: f_code parameter added
574
                    s->qscale -= s->dquant;
M
Michael Niedermayer 已提交
575
                    s->mv_bits+= get_bits_diff(s);
F
Fabrice Bellard 已提交
576
                }
M
Michael Niedermayer 已提交
577
                s->f_count++;
578 579 580 581 582 583
        } else
            {    // RAL: All the following bloc added for B frames:
                if (cbp != 0)
                    {    // With coded bloc pattern
                    if (s->mv_dir == (MV_DIR_FORWARD | MV_DIR_BACKWARD))
                        {    // Bi-directional motion
584 585
                        if (s->dquant) {
                            put_mb_modes(s, 5, 2, 1);
586
                            put_bits(&s->pb, 5, s->qscale);
587 588 589
                        } else {
                            put_mb_modes(s, 2, 3, 1);
                        }
590 591 592 593 594 595 596 597 598 599 600 601
                        s->misc_bits += get_bits_diff(s);
                        mpeg1_encode_motion(s, s->mv[0][0][0] - s->last_mv[0][0][0], s->f_code);
                        mpeg1_encode_motion(s, s->mv[0][0][1] - s->last_mv[0][0][1], s->f_code);
                        mpeg1_encode_motion(s, s->mv[1][0][0] - s->last_mv[1][0][0], s->b_code);
                        mpeg1_encode_motion(s, s->mv[1][0][1] - s->last_mv[1][0][1], s->b_code);
                        s->b_count++;
                        s->f_count++;
                        s->mv_bits += get_bits_diff(s);
                        put_bits(&s->pb, mbPatTable[cbp - 1][1], mbPatTable[cbp - 1][0]);
                        }
                    else if (s->mv_dir == MV_DIR_BACKWARD)
                        {    // Backward motion
602 603
                        if (s->dquant) {
                            put_mb_modes(s, 6, 2, 1);
604
                            put_bits(&s->pb, 5, s->qscale);
605 606 607
                        } else {
                            put_mb_modes(s, 3, 3, 1);
                        }
608 609 610 611 612 613 614 615 616
                        s->misc_bits += get_bits_diff(s);
                        mpeg1_encode_motion(s, motion_x - s->last_mv[1][0][0], s->b_code); 
                        mpeg1_encode_motion(s, motion_y - s->last_mv[1][0][1], s->b_code); 
                        s->b_count++;
                        s->mv_bits += get_bits_diff(s);
                        put_bits(&s->pb, mbPatTable[cbp - 1][1], mbPatTable[cbp - 1][0]);
                        }
                    else if (s->mv_dir == MV_DIR_FORWARD)
                        {    // Forward motion
617 618
                        if (s->dquant) {
                            put_mb_modes(s, 6, 3, 1);
619
                            put_bits(&s->pb, 5, s->qscale);
620 621 622
                        } else {
                            put_mb_modes(s, 4, 3, 1);
                        }
623 624 625 626 627 628 629 630 631 632 633
                        s->misc_bits += get_bits_diff(s);
                        mpeg1_encode_motion(s, motion_x - s->last_mv[0][0][0], s->f_code); 
                        mpeg1_encode_motion(s, motion_y - s->last_mv[0][0][1], s->f_code); 
                        s->f_count++;
                        s->mv_bits += get_bits_diff(s);
                        put_bits(&s->pb, mbPatTable[cbp - 1][1], mbPatTable[cbp - 1][0]);
                        }
                    }
                else
                    {    // No coded bloc pattern
                    if (s->mv_dir == (MV_DIR_FORWARD | MV_DIR_BACKWARD))
634
                        {    // Bi-directional motion 
635
                        put_bits(&s->pb, 2, 2); /* backward & forward motion */
636 637
                        if (!s->frame_pred_frame_dct)
                            put_bits(&s->pb, 2, 2); /* motion_type: frame */
638 639 640 641 642 643 644 645 646 647
                        mpeg1_encode_motion(s, s->mv[0][0][0] - s->last_mv[0][0][0], s->f_code);
                        mpeg1_encode_motion(s, s->mv[0][0][1] - s->last_mv[0][0][1], s->f_code);
                        mpeg1_encode_motion(s, s->mv[1][0][0] - s->last_mv[1][0][0], s->b_code);
                        mpeg1_encode_motion(s, s->mv[1][0][1] - s->last_mv[1][0][1], s->b_code);
                        s->b_count++;
                        s->f_count++;
                        }
                    else if (s->mv_dir == MV_DIR_BACKWARD)
                        {    // Backward motion
                        put_bits(&s->pb, 3, 2); /* backward motion only */
648 649
                        if (!s->frame_pred_frame_dct)
                            put_bits(&s->pb, 2, 2); /* motion_type: frame */
650 651 652 653 654 655 656
                        mpeg1_encode_motion(s, motion_x - s->last_mv[1][0][0], s->b_code); 
                        mpeg1_encode_motion(s, motion_y - s->last_mv[1][0][1], s->b_code); 
                        s->b_count++;
                        }
                    else if (s->mv_dir == MV_DIR_FORWARD)
                        {    // Forward motion
                        put_bits(&s->pb, 4, 2); /* forward motion only */
657 658
                        if (!s->frame_pred_frame_dct)
                            put_bits(&s->pb, 2, 2); /* motion_type: frame */
659 660 661 662 663 664 665 666
                        mpeg1_encode_motion(s, motion_x - s->last_mv[0][0][0], s->f_code); 
                        mpeg1_encode_motion(s, motion_y - s->last_mv[0][0][1], s->f_code); 
                        s->f_count++;
                        }
                    s->qscale -= s->dquant;
                    s->mv_bits += get_bits_diff(s);
                    }
            // End of bloc from RAL
F
Fabrice Bellard 已提交
667 668 669 670 671 672
            }
        for(i=0;i<6;i++) {
            if (cbp & (1 << (5 - i))) {
                mpeg1_encode_block(s, block[i], i);
            }
        }
673
        s->mb_skip_run = 0;
M
Michael Niedermayer 已提交
674 675 676 677
        if(s->mb_intra)
            s->i_tex_bits+= get_bits_diff(s);
        else
            s->p_tex_bits+= get_bits_diff(s);
F
Fabrice Bellard 已提交
678
    }
679 680 681 682 683 684 685 686 687 688 689 690

    // RAL: By this:
    if (s->mv_dir & MV_DIR_FORWARD)
        {
        s->last_mv[0][0][0]= s->mv[0][0][0];
        s->last_mv[0][0][1]= s->mv[0][0][1];
        }
    if (s->mv_dir & MV_DIR_BACKWARD)
        {
        s->last_mv[1][0][0]= s->mv[1][0][0];
        s->last_mv[1][0][1]= s->mv[1][0][1];
        }
F
Fabrice Bellard 已提交
691 692
}

693 694
// RAL: Parameter added: f_or_b_code
static void mpeg1_encode_motion(MpegEncContext *s, int val, int f_or_b_code)
F
Fabrice Bellard 已提交
695 696 697 698 699 700 701 702 703 704
{
    int code, bit_size, l, m, bits, range, sign;

    if (val == 0) {
        /* zero vector */
        code = 0;
        put_bits(&s->pb,
                 mbMotionVectorTable[0][1], 
                 mbMotionVectorTable[0][0]); 
    } else {
705
        bit_size = f_or_b_code - 1;
F
Fabrice Bellard 已提交
706 707 708 709 710 711 712 713 714 715 716 717 718 719 720 721 722 723 724 725 726 727
        range = 1 << bit_size;
        /* modulo encoding */
        l = 16 * range;
        m = 2 * l;
        if (val < -l) {
            val += m;
        } else if (val >= l) {
            val -= m;
        }

        if (val >= 0) {
            val--;
            code = (val >> bit_size) + 1;
            bits = val & (range - 1);
            sign = 0;
        } else {
            val = -val;
            val--;
            code = (val >> bit_size) + 1;
            bits = val & (range - 1);
            sign = 1;
        }
728 729 730

        assert(code > 0 && code <= 16);

F
Fabrice Bellard 已提交
731 732 733
        put_bits(&s->pb,
                 mbMotionVectorTable[code][1], 
                 mbMotionVectorTable[code][0]); 
734

F
Fabrice Bellard 已提交
735 736 737 738 739 740 741
        put_bits(&s->pb, 1, sign);
        if (bit_size > 0) {
            put_bits(&s->pb, bit_size, bits);
        }
    }
}

742
void ff_mpeg1_encode_init(MpegEncContext *s)
743 744
{
    static int done=0;
745 746 747

    common_init(s);

748 749 750
    if(!done){
        int f_code;
        int mv;
751
	int i;
752 753

        done=1;
754
        init_rl(&rl_mpeg1);
755

756 757 758 759 760
	for(i=0; i<64; i++)
	{
		mpeg1_max_level[0][i]= rl_mpeg1.max_level[0][i];
		mpeg1_index_run[0][i]= rl_mpeg1.index_run[0][i];
	}
761 762
        
        init_uni_ac_vlc(&rl_mpeg1, uni_mpeg1_ac_vlc_bits, uni_mpeg1_ac_vlc_len);
763 764 765 766 767 768 769 770 771 772 773 774 775 776 777 778 779 780 781 782 783

	/* build unified dc encoding tables */
	for(i=-255; i<256; i++)
	{
		int adiff, index;
		int bits, code;
		int diff=i;

		adiff = ABS(diff);
		if(diff<0) diff--;
		index = vlc_dc_table[adiff];

		bits= vlc_dc_lum_bits[index] + index;
		code= (vlc_dc_lum_code[index]<<index) + (diff & ((1 << index) - 1));
		mpeg1_lum_dc_uni[i+255]= bits + (code<<8);
		
		bits= vlc_dc_chroma_bits[index] + index;
		code= (vlc_dc_chroma_code[index]<<index) + (diff & ((1 << index) - 1));
		mpeg1_chr_dc_uni[i+255]= bits + (code<<8);
	}

784 785
        mv_penalty= av_mallocz( sizeof(uint8_t)*(MAX_FCODE+1)*(2*MAX_MV+1) );

786 787 788 789 790 791 792 793 794 795 796 797 798 799 800 801 802 803 804 805 806 807 808 809 810 811 812 813 814 815 816 817 818 819
        for(f_code=1; f_code<=MAX_FCODE; f_code++){
            for(mv=-MAX_MV; mv<=MAX_MV; mv++){
                int len;

                if(mv==0) len= mbMotionVectorTable[0][1];
                else{
                    int val, bit_size, range, code;

                    bit_size = s->f_code - 1;
                    range = 1 << bit_size;

                    val=mv;
                    if (val < 0) 
                        val = -val;
                    val--;
                    code = (val >> bit_size) + 1;
                    if(code<17){
                        len= mbMotionVectorTable[code][1] + 1 + bit_size;
                    }else{
                        len= mbMotionVectorTable[16][1] + 2 + bit_size;
                    }
                }

                mv_penalty[f_code][mv+MAX_MV]= len;
            }
        }
        

        for(f_code=MAX_FCODE; f_code>0; f_code--){
            for(mv=-(8<<f_code); mv<(8<<f_code); mv++){
                fcode_tab[mv+MAX_MV]= f_code;
            }
        }
    }
M
Michael Niedermayer 已提交
820
    s->me.mv_penalty= mv_penalty;
821
    s->fcode_tab= fcode_tab;
M
Michael Niedermayer 已提交
822 823 824 825 826 827 828
    if(s->codec_id == CODEC_ID_MPEG1VIDEO){
        s->min_qcoeff=-255;
        s->max_qcoeff= 255;
    }else{
        s->min_qcoeff=-2047;
        s->max_qcoeff= 2047;
    }
829
    s->intra_ac_vlc_length=
830 831 832
    s->inter_ac_vlc_length=
    s->intra_ac_vlc_last_length=
    s->inter_ac_vlc_last_length= uni_mpeg1_ac_vlc_len;
833
}
834

F
Fabrice Bellard 已提交
835 836 837
static inline void encode_dc(MpegEncContext *s, int diff, int component)
{
    if (component == 0) {
838 839 840 841
        put_bits(
	    &s->pb, 
	    mpeg1_lum_dc_uni[diff+255]&0xFF,
	    mpeg1_lum_dc_uni[diff+255]>>8);
F
Fabrice Bellard 已提交
842
    } else {
843 844 845 846
        put_bits(
            &s->pb, 
	    mpeg1_chr_dc_uni[diff+255]&0xFF,
	    mpeg1_chr_dc_uni[diff+255]>>8);
F
Fabrice Bellard 已提交
847 848 849 850 851 852 853 854 855
    }
}

static void mpeg1_encode_block(MpegEncContext *s, 
                               DCTELEM *block, 
                               int n)
{
    int alevel, level, last_non_zero, dc, diff, i, j, run, last_index, sign;
    int code, component;
856
//    RLTable *rl = &rl_mpeg1;
F
Fabrice Bellard 已提交
857 858 859 860 861 862 863 864 865 866 867

    last_index = s->block_last_index[n];

    /* DC coef */
    if (s->mb_intra) {
        component = (n <= 3 ? 0 : n - 4 + 1);
        dc = block[0]; /* overflow is impossible */
        diff = dc - s->last_dc[component];
        encode_dc(s, diff, component);
        s->last_dc[component] = dc;
        i = 1;
M
Michael Niedermayer 已提交
868 869 870 871 872 873
/*
        if (s->intra_vlc_format)
            rl = &rl_mpeg2;
        else
            rl = &rl_mpeg1;
*/
F
Fabrice Bellard 已提交
874 875 876 877 878
    } else {
        /* encode the first coefficient : needs to be done here because
           it is handled slightly differently */
        level = block[0];
        if (abs(level) == 1) {
879
                code = ((uint32_t)level >> 31); /* the sign bit */
F
Fabrice Bellard 已提交
880 881 882 883 884 885 886 887 888 889 890
                put_bits(&s->pb, 2, code | 0x02);
                i = 1;
        } else {
            i = 0;
            last_non_zero = -1;
            goto next_coef;
        }
    }

    /* now quantify & encode AC coefs */
    last_non_zero = i - 1;
891

F
Fabrice Bellard 已提交
892
    for(;i<=last_index;i++) {
893
        j = s->intra_scantable.permutated[i];
F
Fabrice Bellard 已提交
894 895 896 897 898 899 900 901 902
        level = block[j];
    next_coef:
#if 0
        if (level != 0)
            dprintf("level[%d]=%d\n", i, level);
#endif            
        /* encode using VLC */
        if (level != 0) {
            run = i - last_non_zero - 1;
903 904 905 906 907
            
            alevel= level;
            MASK_ABS(sign, alevel)
            sign&=1;

908
//            code = get_rl_index(rl, 0, run, alevel);
909
            if (alevel <= mpeg1_max_level[0][run]){
910 911 912
                code= mpeg1_index_run[0][run] + alevel - 1;
	    	/* store the vlc & sign at once */
                put_bits(&s->pb, mpeg1_vlc[code][1]+1, (mpeg1_vlc[code][0]<<1) + sign);
F
Fabrice Bellard 已提交
913
            } else {
914 915
		/* escape seems to be pretty rare <5% so i dont optimize it */
                put_bits(&s->pb, mpeg1_vlc[111/*rl->n*/][1], mpeg1_vlc[111/*rl->n*/][0]);
F
Fabrice Bellard 已提交
916 917
                /* escape: only clip in this case */
                put_bits(&s->pb, 6, run);
M
Michael Niedermayer 已提交
918 919 920
                if(s->codec_id == CODEC_ID_MPEG1VIDEO){
                    if (alevel < 128) {
                        put_bits(&s->pb, 8, level & 0xff);
F
Fabrice Bellard 已提交
921
                    } else {
M
Michael Niedermayer 已提交
922 923 924 925 926
                        if (level < 0) {
                            put_bits(&s->pb, 16, 0x8001 + level + 255);
                        } else {
                            put_bits(&s->pb, 16, level & 0xffff);
                        }
F
Fabrice Bellard 已提交
927
                    }
M
Michael Niedermayer 已提交
928 929
                }else{
                    put_bits(&s->pb, 12, level & 0xfff);
F
Fabrice Bellard 已提交
930 931 932 933 934 935 936 937
                }
            }
            last_non_zero = i;
        }
    }
    /* end of block */
    put_bits(&s->pb, 2, 0x2);
}
938
#endif //CONFIG_ENCODERS
F
Fabrice Bellard 已提交
939 940 941 942 943 944 945 946 947 948 949 950

/******************************************/
/* decoding */

static VLC dc_lum_vlc;
static VLC dc_chroma_vlc;
static VLC mv_vlc;
static VLC mbincr_vlc;
static VLC mb_ptype_vlc;
static VLC mb_btype_vlc;
static VLC mb_pat_vlc;

951
static void init_vlcs()
F
Fabrice Bellard 已提交
952 953 954 955
{
    static int done = 0;

    if (!done) {
956
        done = 1;
F
Fabrice Bellard 已提交
957

M
Michael Niedermayer 已提交
958
        init_vlc(&dc_lum_vlc, DC_VLC_BITS, 12, 
F
Fabrice Bellard 已提交
959 960
                 vlc_dc_lum_bits, 1, 1,
                 vlc_dc_lum_code, 2, 2);
M
Michael Niedermayer 已提交
961
        init_vlc(&dc_chroma_vlc,  DC_VLC_BITS, 12, 
F
Fabrice Bellard 已提交
962 963
                 vlc_dc_chroma_bits, 1, 1,
                 vlc_dc_chroma_code, 2, 2);
964
        init_vlc(&mv_vlc, MV_VLC_BITS, 17, 
F
Fabrice Bellard 已提交
965 966
                 &mbMotionVectorTable[0][1], 2, 1,
                 &mbMotionVectorTable[0][0], 2, 1);
967
        init_vlc(&mbincr_vlc, MBINCR_VLC_BITS, 36, 
F
Fabrice Bellard 已提交
968 969
                 &mbAddrIncrTable[0][1], 2, 1,
                 &mbAddrIncrTable[0][0], 2, 1);
970
        init_vlc(&mb_pat_vlc, MB_PAT_VLC_BITS, 63, 
F
Fabrice Bellard 已提交
971 972 973
                 &mbPatTable[0][1], 2, 1,
                 &mbPatTable[0][0], 2, 1);
        
974
        init_vlc(&mb_ptype_vlc, MB_PTYPE_VLC_BITS, 7, 
F
Fabrice Bellard 已提交
975 976
                 &table_mb_ptype[0][1], 2, 1,
                 &table_mb_ptype[0][0], 2, 1);
977
        init_vlc(&mb_btype_vlc, MB_BTYPE_VLC_BITS, 11, 
F
Fabrice Bellard 已提交
978 979 980 981
                 &table_mb_btype[0][1], 2, 1,
                 &table_mb_btype[0][0], 2, 1);
        init_rl(&rl_mpeg1);
        init_rl(&rl_mpeg2);
982 983 984

        init_2d_vlc_rl(&rl_mpeg1);
        init_2d_vlc_rl(&rl_mpeg2);
F
Fabrice Bellard 已提交
985 986 987 988 989
    }
}

static inline int get_dmv(MpegEncContext *s)
{
990 991
    if(get_bits1(&s->gb)) 
        return 1 - (get_bits1(&s->gb) << 1);
F
Fabrice Bellard 已提交
992 993 994 995
    else
        return 0;
}

F
Fabrice Bellard 已提交
996 997
static inline int get_qscale(MpegEncContext *s)
{
998
    int qscale = get_bits(&s->gb, 5);
M
Michael Niedermayer 已提交
999
    if (s->codec_id == CODEC_ID_MPEG2VIDEO) {
F
Fabrice Bellard 已提交
1000
        if (s->q_scale_type) {
1001
            return non_linear_qscale[qscale];
F
Fabrice Bellard 已提交
1002
        } else {
1003
            return qscale << 1;
F
Fabrice Bellard 已提交
1004 1005 1006 1007 1008
        }
    }
    return qscale;
}

F
Fabrice Bellard 已提交
1009 1010 1011 1012 1013 1014 1015 1016 1017
/* motion type (for mpeg2) */
#define MT_FIELD 1
#define MT_FRAME 2
#define MT_16X8  2
#define MT_DMV   3

static int mpeg_decode_mb(MpegEncContext *s,
                          DCTELEM block[6][64])
{
M
Michael Niedermayer 已提交
1018
    int i, j, k, cbp, val, mb_type, motion_type;
F
Fabrice Bellard 已提交
1019 1020 1021
    
    dprintf("decode_mb: x=%d y=%d\n", s->mb_x, s->mb_y);

1022 1023
    assert(s->mb_skiped==0);

1024
    if (s->mb_skip_run-- != 0) {
1025
        if(s->pict_type == I_TYPE){
1026
            av_log(s->avctx, AV_LOG_ERROR, "skiped MB in I frame at %d %d\n", s->mb_x, s->mb_y);
1027 1028 1029
            return -1;
        }
    
F
Fabrice Bellard 已提交
1030 1031 1032 1033 1034 1035 1036 1037 1038 1039
        /* skip mb */
        s->mb_intra = 0;
        for(i=0;i<6;i++)
            s->block_last_index[i] = -1;
        s->mv_type = MV_TYPE_16X16;
        if (s->pict_type == P_TYPE) {
            /* if P type, zero motion vector is implied */
            s->mv_dir = MV_DIR_FORWARD;
            s->mv[0][0][0] = s->mv[0][0][1] = 0;
            s->last_mv[0][0][0] = s->last_mv[0][0][1] = 0;
1040
            s->last_mv[0][1][0] = s->last_mv[0][1][1] = 0;
1041
            s->mb_skiped = 1;
1042
            s->current_picture.mb_type[ s->mb_x + s->mb_y*s->mb_stride ]= MB_TYPE_SKIP | MB_TYPE_L0 | MB_TYPE_16x16;
F
Fabrice Bellard 已提交
1043 1044 1045 1046 1047 1048
        } else {
            /* if B type, reuse previous vectors and directions */
            s->mv[0][0][0] = s->last_mv[0][0][0];
            s->mv[0][0][1] = s->last_mv[0][0][1];
            s->mv[1][0][0] = s->last_mv[1][0][0];
            s->mv[1][0][1] = s->last_mv[1][0][1];
1049

1050 1051 1052 1053
            s->current_picture.mb_type[ s->mb_x + s->mb_y*s->mb_stride ]= 
                s->current_picture.mb_type[ s->mb_x + s->mb_y*s->mb_stride - 1] | MB_TYPE_SKIP;
//            assert(s->current_picture.mb_type[ s->mb_x + s->mb_y*s->mb_stride - 1]&(MB_TYPE_16x16|MB_TYPE_16x8));

1054 1055
            if((s->mv[0][0][0]|s->mv[0][0][1]|s->mv[1][0][0]|s->mv[1][0][1])==0) 
                s->mb_skiped = 1;
F
Fabrice Bellard 已提交
1056
        }
M
Michael Niedermayer 已提交
1057

F
Fabrice Bellard 已提交
1058 1059 1060 1061 1062 1063
        return 0;
    }

    switch(s->pict_type) {
    default:
    case I_TYPE:
1064
        if (get_bits1(&s->gb) == 0) {
M
Michael Niedermayer 已提交
1065
            if (get_bits1(&s->gb) == 0){
1066
                av_log(s->avctx, AV_LOG_ERROR, "invalid mb type in I Frame at %d %d\n", s->mb_x, s->mb_y);
F
Fabrice Bellard 已提交
1067
                return -1;
M
Michael Niedermayer 已提交
1068
            }
1069
            mb_type = MB_TYPE_QUANT | MB_TYPE_INTRA;
F
Fabrice Bellard 已提交
1070
        } else {
1071
            mb_type = MB_TYPE_INTRA;
F
Fabrice Bellard 已提交
1072 1073 1074
        }
        break;
    case P_TYPE:
1075
        mb_type = get_vlc2(&s->gb, mb_ptype_vlc.table, MB_PTYPE_VLC_BITS, 1);
M
Michael Niedermayer 已提交
1076
        if (mb_type < 0){
1077
            av_log(s->avctx, AV_LOG_ERROR, "invalid mb type in P Frame at %d %d\n", s->mb_x, s->mb_y);
F
Fabrice Bellard 已提交
1078
            return -1;
M
Michael Niedermayer 已提交
1079
        }
1080
        mb_type = ptype2mb_type[ mb_type ];
F
Fabrice Bellard 已提交
1081 1082
        break;
    case B_TYPE:
1083
        mb_type = get_vlc2(&s->gb, mb_btype_vlc.table, MB_BTYPE_VLC_BITS, 1);
M
Michael Niedermayer 已提交
1084
        if (mb_type < 0){
1085
            av_log(s->avctx, AV_LOG_ERROR, "invalid mb type in B Frame at %d %d\n", s->mb_x, s->mb_y);
F
Fabrice Bellard 已提交
1086
            return -1;
M
Michael Niedermayer 已提交
1087
        }
1088
        mb_type = btype2mb_type[ mb_type ];
F
Fabrice Bellard 已提交
1089 1090 1091
        break;
    }
    dprintf("mb_type=%x\n", mb_type);
1092 1093 1094 1095 1096 1097 1098
//    motion_type = 0; /* avoid warning */
    if (IS_INTRA(mb_type)) {
        /* compute dct type */
        if (s->picture_structure == PICT_FRAME && //FIXME add a interlaced_dct coded var?
            !s->frame_pred_frame_dct) {
            s->interlaced_dct = get_bits1(&s->gb);
        }
F
Fabrice Bellard 已提交
1099

1100 1101 1102
        if (IS_QUANT(mb_type))
            s->qscale = get_qscale(s);
        
F
Fabrice Bellard 已提交
1103 1104 1105
        if (s->concealment_motion_vectors) {
            /* just parse them */
            if (s->picture_structure != PICT_FRAME) 
1106
                skip_bits1(&s->gb); /* field select */
1107 1108 1109 1110 1111 1112
            
            s->mv[0][0][0]= s->last_mv[0][0][0]= s->last_mv[0][1][0] = 
                mpeg_decode_motion(s, s->mpeg_f_code[0][0], s->last_mv[0][0][0]);
            s->mv[0][0][1]= s->last_mv[0][0][1]= s->last_mv[0][1][1] = 
                mpeg_decode_motion(s, s->mpeg_f_code[0][1], s->last_mv[0][0][1]);

1113
            skip_bits1(&s->gb); /* marker */
1114 1115
        }else
            memset(s->last_mv, 0, sizeof(s->last_mv)); /* reset mv prediction */
F
Fabrice Bellard 已提交
1116
        s->mb_intra = 1;
1117 1118 1119 1120 1121 1122 1123 1124 1125
#ifdef HAVE_XVMC
        //one 1 we memcpy blocks in xvmcvideo
        if(s->avctx->xvmc_acceleration > 1){
            XVMC_pack_pblocks(s,-1);//inter are always full blocks
            if(s->swap_uv){
                exchange_uv(s);
            }
        }
#endif
1126

M
Michael Niedermayer 已提交
1127
        if (s->codec_id == CODEC_ID_MPEG2VIDEO) {
1128
            for(i=0;i<6;i++) {
1129
                if (mpeg2_decode_block_intra(s, s->pblocks[i], i) < 0)
1130 1131 1132 1133
                    return -1;
            }
        } else {
            for(i=0;i<6;i++) {
1134
                if (mpeg1_decode_block_intra(s, s->pblocks[i], i) < 0)
1135 1136 1137
                    return -1;
            }
        }
F
Fabrice Bellard 已提交
1138
    } else {
1139
        if (mb_type & MB_TYPE_ZERO_MV){
1140
            assert(mb_type & MB_TYPE_CBP);
1141 1142 1143 1144 1145 1146 1147 1148 1149 1150 1151 1152 1153 1154 1155 1156 1157 1158 1159 1160 1161 1162 1163 1164 1165 1166 1167 1168 1169 1170

            /* compute dct type */
            if (s->picture_structure == PICT_FRAME && //FIXME add a interlaced_dct coded var?
                !s->frame_pred_frame_dct) {
                s->interlaced_dct = get_bits1(&s->gb);
            }

            if (IS_QUANT(mb_type))
                s->qscale = get_qscale(s);

            s->mv_dir = MV_DIR_FORWARD;
            s->mv_type = MV_TYPE_16X16;
            s->last_mv[0][0][0] = 0;
            s->last_mv[0][0][1] = 0;
            s->last_mv[0][1][0] = 0;
            s->last_mv[0][1][1] = 0;
            s->mv[0][0][0] = 0;
            s->mv[0][0][1] = 0;
        }else{
            assert(mb_type & MB_TYPE_L0L1);
//FIXME decide if MBs in field pictures are MB_TYPE_INTERLACED
            /* get additionnal motion vector type */
            if (s->frame_pred_frame_dct) 
                motion_type = MT_FRAME;
            else{
                motion_type = get_bits(&s->gb, 2);
            }

            /* compute dct type */
            if (s->picture_structure == PICT_FRAME && //FIXME add a interlaced_dct coded var?
1171
                !s->frame_pred_frame_dct && HAS_CBP(mb_type)) {
1172 1173 1174 1175 1176 1177 1178 1179 1180 1181 1182 1183 1184 1185 1186 1187 1188 1189 1190 1191 1192 1193
                s->interlaced_dct = get_bits1(&s->gb);
            }

            if (IS_QUANT(mb_type))
                s->qscale = get_qscale(s);

            /* motion vectors */
            s->mv_dir = 0;
            for(i=0;i<2;i++) {
                if (USES_LIST(mb_type, i)) {
                    s->mv_dir |= (MV_DIR_FORWARD >> i);
                    dprintf("motion_type=%d\n", motion_type);
                    switch(motion_type) {
                    case MT_FRAME: /* or MT_16X8 */
                        if (s->picture_structure == PICT_FRAME) {
                            /* MT_FRAME */
                            mb_type |= MB_TYPE_16x16; 
                            s->mv_type = MV_TYPE_16X16;
                            s->mv[i][0][0]= s->last_mv[i][0][0]= s->last_mv[i][1][0] = 
                                mpeg_decode_motion(s, s->mpeg_f_code[i][0], s->last_mv[i][0][0]);
                            s->mv[i][0][1]= s->last_mv[i][0][1]= s->last_mv[i][1][1] = 
                                mpeg_decode_motion(s, s->mpeg_f_code[i][1], s->last_mv[i][0][1]);
F
Fabrice Bellard 已提交
1194
                            /* full_pel: only for mpeg1 */
1195 1196 1197 1198 1199 1200
                            if (s->full_pel[i]){
                                s->mv[i][0][0] <<= 1;
                                s->mv[i][0][1] <<= 1;
                            }
                        } else {
                            /* MT_16X8 */
M
Michael Niedermayer 已提交
1201
                            mb_type |= MB_TYPE_16x8 | MB_TYPE_INTERLACED; 
1202 1203 1204 1205 1206 1207 1208 1209 1210 1211
                            s->mv_type = MV_TYPE_16X8;
                            for(j=0;j<2;j++) {
                                s->field_select[i][j] = get_bits1(&s->gb);
                                for(k=0;k<2;k++) {
                                    val = mpeg_decode_motion(s, s->mpeg_f_code[i][k],
                                                             s->last_mv[i][j][k]);
                                    s->last_mv[i][j][k] = val;
                                    s->mv[i][j][k] = val;
                                }
                            }
F
Fabrice Bellard 已提交
1212
                        }
1213 1214 1215 1216 1217 1218 1219 1220 1221 1222 1223 1224 1225 1226 1227 1228 1229 1230 1231
                        break;
                    case MT_FIELD:
                        s->mv_type = MV_TYPE_FIELD;
                        if (s->picture_structure == PICT_FRAME) {
                            mb_type |= MB_TYPE_16x8 | MB_TYPE_INTERLACED; 
                            for(j=0;j<2;j++) {
                                s->field_select[i][j] = get_bits1(&s->gb);
                                val = mpeg_decode_motion(s, s->mpeg_f_code[i][0],
                                                         s->last_mv[i][j][0]);
                                s->last_mv[i][j][0] = val;
                                s->mv[i][j][0] = val;
                                dprintf("fmx=%d\n", val);
                                val = mpeg_decode_motion(s, s->mpeg_f_code[i][1],
                                                         s->last_mv[i][j][1] >> 1);
                                s->last_mv[i][j][1] = val << 1;
                                s->mv[i][j][1] = val;
                                dprintf("fmy=%d\n", val);
                            }
                        } else {
M
Michael Niedermayer 已提交
1232
                            mb_type |= MB_TYPE_16x16 | MB_TYPE_INTERLACED; 
1233
                            s->field_select[i][0] = get_bits1(&s->gb);
F
Fabrice Bellard 已提交
1234 1235
                            for(k=0;k<2;k++) {
                                val = mpeg_decode_motion(s, s->mpeg_f_code[i][k],
1236 1237 1238 1239
                                                         s->last_mv[i][0][k]);
                                s->last_mv[i][0][k] = val;
                                s->last_mv[i][1][k] = val;
                                s->mv[i][0][k] = val;
F
Fabrice Bellard 已提交
1240 1241
                            }
                        }
1242 1243 1244 1245 1246 1247 1248 1249 1250 1251 1252 1253 1254 1255
                        break;
                    case MT_DMV:
                        {
                            int dmx, dmy, mx, my, m;

                            mx = mpeg_decode_motion(s, s->mpeg_f_code[i][0], 
                                                    s->last_mv[i][0][0]);
                            s->last_mv[i][0][0] = mx;
                            s->last_mv[i][1][0] = mx;
                            dmx = get_dmv(s);
                            my = mpeg_decode_motion(s, s->mpeg_f_code[i][1], 
                                                    s->last_mv[i][0][1] >> 1);
                            dmy = get_dmv(s);
                            s->mv_type = MV_TYPE_DMV;
1256 1257 1258 1259 1260 1261 1262 1263 1264 1265


                            s->last_mv[i][0][1] = my<<1;
                            s->last_mv[i][1][1] = my<<1;

                            s->mv[i][0][0] = mx;
                            s->mv[i][0][1] = my;
                            s->mv[i][1][0] = mx;//not used
                            s->mv[i][1][1] = my;//not used

1266 1267 1268
                            if (s->picture_structure == PICT_FRAME) {
                                mb_type |= MB_TYPE_16x16 | MB_TYPE_INTERLACED; 

1269
                                //m = 1 + 2 * s->top_field_first;
1270
                                m = s->top_field_first ? 1 : 3;
1271

1272
                                /* top -> top pred */
1273 1274
                                s->mv[i][2][0] = ((mx * m + (mx > 0)) >> 1) + dmx;
                                s->mv[i][2][1] = ((my * m + (my > 0)) >> 1) + dmy - 1;
1275 1276 1277 1278 1279 1280
                                m = 4 - m;
                                s->mv[i][3][0] = ((mx * m + (mx > 0)) >> 1) + dmx;
                                s->mv[i][3][1] = ((my * m + (my > 0)) >> 1) + dmy + 1;
                            } else {
                                mb_type |= MB_TYPE_16x16;

1281 1282 1283 1284 1285 1286
                                s->mv[i][2][0] = ((mx + (mx > 0)) >> 1) + dmx;
                                s->mv[i][2][1] = ((my + (my > 0)) >> 1) + dmy;
                                if(s->picture_structure == PICT_TOP_FIELD)
                                    s->mv[i][2][1]--;
                                else 
                                    s->mv[i][2][1]++;
1287
                            }
F
Fabrice Bellard 已提交
1288
                        }
1289
                        break;
M
Michael Niedermayer 已提交
1290 1291 1292
                    default:
                        av_log(s->avctx, AV_LOG_ERROR, "00 motion_type at %d %d\n", s->mb_x, s->mb_y);
                        return -1;
F
Fabrice Bellard 已提交
1293 1294 1295 1296
                    }
                }
            }
        }
1297 1298
        
        s->mb_intra = 0;
F
Fabrice Bellard 已提交
1299

1300
        if (HAS_CBP(mb_type)) {
1301 1302
            cbp = get_vlc2(&s->gb, mb_pat_vlc.table, MB_PAT_VLC_BITS, 1);
            if (cbp < 0){
1303
                av_log(s->avctx, AV_LOG_ERROR, "invalid cbp at %d %d\n", s->mb_x, s->mb_y);
1304
                return -1;
F
Fabrice Bellard 已提交
1305
            }
1306 1307
            cbp++;

1308 1309 1310 1311 1312 1313 1314 1315 1316 1317
#ifdef HAVE_XVMC
            //on 1 we memcpy blocks in xvmcvideo
            if(s->avctx->xvmc_acceleration > 1){
                XVMC_pack_pblocks(s,cbp);
                if(s->swap_uv){
                    exchange_uv(s);
                }
            }    
#endif

M
Michael Niedermayer 已提交
1318
            if (s->codec_id == CODEC_ID_MPEG2VIDEO) {
1319 1320
                for(i=0;i<6;i++) {
                    if (cbp & 32) {
1321
                        if (mpeg2_decode_block_non_intra(s, s->pblocks[i], i) < 0)
1322 1323 1324 1325 1326
                            return -1;
                    } else {
                        s->block_last_index[i] = -1;
                    }
                    cbp+=cbp;
F
Fabrice Bellard 已提交
1327
                }
1328 1329 1330
            } else {
                for(i=0;i<6;i++) {
                    if (cbp & 32) {
1331
                        if (mpeg1_decode_block_inter(s, s->pblocks[i], i) < 0)
1332 1333 1334 1335 1336
                            return -1;
                    } else {
                        s->block_last_index[i] = -1;
                    }
                    cbp+=cbp;
1337
                }
F
Fabrice Bellard 已提交
1338
            }
1339 1340 1341
        }else{
            for(i=0;i<6;i++)
                s->block_last_index[i] = -1;
F
Fabrice Bellard 已提交
1342 1343
        }
    }
1344 1345 1346

    s->current_picture.mb_type[ s->mb_x + s->mb_y*s->mb_stride ]= mb_type;

F
Fabrice Bellard 已提交
1347 1348 1349 1350 1351 1352
    return 0;
}

/* as h263, but only 17 codes */
static int mpeg_decode_motion(MpegEncContext *s, int fcode, int pred)
{
1353
    int code, sign, val, l, shift;
F
Fabrice Bellard 已提交
1354

1355
    code = get_vlc2(&s->gb, mv_vlc.table, MV_VLC_BITS, 2);
F
Fabrice Bellard 已提交
1356 1357 1358
    if (code == 0) {
        return pred;
    }
1359 1360 1361 1362
    if (code < 0) {
        return 0xffff;
    }

1363
    sign = get_bits1(&s->gb);
F
Fabrice Bellard 已提交
1364
    shift = fcode - 1;
1365 1366 1367
    val = code;
    if (shift) {
        val = (val - 1) << shift;
F
Fabrice Bellard 已提交
1368
        val |= get_bits(&s->gb, shift);
1369 1370
        val++;
    }
F
Fabrice Bellard 已提交
1371 1372 1373 1374 1375
    if (sign)
        val = -val;
    val += pred;
    
    /* modulo decoding */
1376
    l = 1 << (shift+4);
1377
    val = ((val + l)&(l*2-1)) - l;
F
Fabrice Bellard 已提交
1378 1379 1380
    return val;
}

1381
static inline int decode_dc(GetBitContext *gb, int component)
F
Fabrice Bellard 已提交
1382 1383 1384 1385
{
    int code, diff;

    if (component == 0) {
1386
        code = get_vlc2(gb, dc_lum_vlc.table, DC_VLC_BITS, 2);
F
Fabrice Bellard 已提交
1387
    } else {
1388
        code = get_vlc2(gb, dc_chroma_vlc.table, DC_VLC_BITS, 2);
F
Fabrice Bellard 已提交
1389
    }
M
Michael Niedermayer 已提交
1390
    if (code < 0){
1391
        av_log(NULL, AV_LOG_ERROR, "invalid dc code at\n");
F
Fabrice Bellard 已提交
1392
        return 0xffff;
M
Michael Niedermayer 已提交
1393
    }
F
Fabrice Bellard 已提交
1394 1395 1396
    if (code == 0) {
        diff = 0;
    } else {
1397
        diff = get_xbits(gb, code);
F
Fabrice Bellard 已提交
1398 1399 1400 1401
    }
    return diff;
}

1402
static inline int mpeg1_decode_block_intra(MpegEncContext *s, 
F
Fabrice Bellard 已提交
1403 1404 1405 1406
                               DCTELEM *block, 
                               int n)
{
    int level, dc, diff, i, j, run;
1407
    int component;
F
Fabrice Bellard 已提交
1408
    RLTable *rl = &rl_mpeg1;
1409 1410
    uint8_t * const scantable= s->intra_scantable.permutated;
    const uint16_t *quant_matrix= s->intra_matrix;
1411
    const int qscale= s->qscale;
F
Fabrice Bellard 已提交
1412

1413 1414
    /* DC coef */
    component = (n <= 3 ? 0 : n - 4 + 1);
1415
    diff = decode_dc(&s->gb, component);
1416 1417 1418 1419 1420 1421 1422 1423 1424 1425 1426 1427 1428 1429 1430 1431 1432 1433 1434 1435 1436 1437 1438 1439 1440 1441 1442 1443 1444 1445 1446 1447 1448 1449 1450 1451 1452 1453 1454 1455 1456 1457 1458 1459 1460 1461 1462
    if (diff >= 0xffff)
        return -1;
    dc = s->last_dc[component];
    dc += diff;
    s->last_dc[component] = dc;
    block[0] = dc<<3;
    dprintf("dc=%d diff=%d\n", dc, diff);
    i = 0;
    {
        OPEN_READER(re, &s->gb);    
        /* now quantify & encode AC coefs */
        for(;;) {
            UPDATE_CACHE(re, &s->gb);
            GET_RL_VLC(level, run, re, &s->gb, rl->rl_vlc[0], TEX_VLC_BITS, 2);
            
            if(level == 127){
                break;
            } else if(level != 0) {
                i += run;
                j = scantable[i];
                level= (level*qscale*quant_matrix[j])>>3;
                level= (level-1)|1;
                level = (level ^ SHOW_SBITS(re, &s->gb, 1)) - SHOW_SBITS(re, &s->gb, 1);
                LAST_SKIP_BITS(re, &s->gb, 1);
            } else {
                /* escape */
                run = SHOW_UBITS(re, &s->gb, 6)+1; LAST_SKIP_BITS(re, &s->gb, 6);
                UPDATE_CACHE(re, &s->gb);
                level = SHOW_SBITS(re, &s->gb, 8); SKIP_BITS(re, &s->gb, 8);
                if (level == -128) {
                    level = SHOW_UBITS(re, &s->gb, 8) - 256; LAST_SKIP_BITS(re, &s->gb, 8);
                } else if (level == 0) {
                    level = SHOW_UBITS(re, &s->gb, 8)      ; LAST_SKIP_BITS(re, &s->gb, 8);
                }
                i += run;
                j = scantable[i];
                if(level<0){
                    level= -level;
                    level= (level*qscale*quant_matrix[j])>>3;
                    level= (level-1)|1;
                    level= -level;
                }else{
                    level= (level*qscale*quant_matrix[j])>>3;
                    level= (level-1)|1;
                }
            }
            if (i > 63){
1463
                av_log(s->avctx, AV_LOG_ERROR, "ac-tex damaged at %d %d\n", s->mb_x, s->mb_y);
1464 1465 1466 1467 1468 1469 1470 1471 1472 1473 1474 1475 1476 1477 1478 1479 1480
                return -1;
            }

            block[j] = level;
        }
        CLOSE_READER(re, &s->gb);
    }
    s->block_last_index[n] = i;
   return 0;
}

static inline int mpeg1_decode_block_inter(MpegEncContext *s, 
                               DCTELEM *block, 
                               int n)
{
    int level, i, j, run;
    RLTable *rl = &rl_mpeg1;
1481 1482
    uint8_t * const scantable= s->intra_scantable.permutated;
    const uint16_t *quant_matrix= s->inter_matrix;
1483 1484 1485
    const int qscale= s->qscale;

    {
1486 1487
        int v;
        OPEN_READER(re, &s->gb);
1488
        i = -1;
F
Fabrice Bellard 已提交
1489
        /* special case for the first coef. no need to add a second vlc table */
1490 1491
        UPDATE_CACHE(re, &s->gb);
        v= SHOW_UBITS(re, &s->gb, 2);
F
Fabrice Bellard 已提交
1492
        if (v & 2) {
1493
            LAST_SKIP_BITS(re, &s->gb, 2);
1494 1495 1496 1497
            level= (3*qscale*quant_matrix[0])>>4;
            level= (level-1)|1;
            if(v&1)
                level= -level;
1498
            block[0] = level;
1499
            i++;
F
Fabrice Bellard 已提交
1500
        }
1501

1502 1503 1504 1505 1506 1507 1508 1509 1510 1511 1512 1513 1514 1515 1516 1517 1518 1519 1520 1521 1522 1523 1524 1525 1526 1527 1528 1529 1530 1531 1532 1533 1534 1535 1536
        /* now quantify & encode AC coefs */
        for(;;) {
            UPDATE_CACHE(re, &s->gb);
            GET_RL_VLC(level, run, re, &s->gb, rl->rl_vlc[0], TEX_VLC_BITS, 2);
            
            if(level == 127){
                break;
            } else if(level != 0) {
                i += run;
                j = scantable[i];
                level= ((level*2+1)*qscale*quant_matrix[j])>>4;
                level= (level-1)|1;
                level = (level ^ SHOW_SBITS(re, &s->gb, 1)) - SHOW_SBITS(re, &s->gb, 1);
                LAST_SKIP_BITS(re, &s->gb, 1);
            } else {
                /* escape */
                run = SHOW_UBITS(re, &s->gb, 6)+1; LAST_SKIP_BITS(re, &s->gb, 6);
                UPDATE_CACHE(re, &s->gb);
                level = SHOW_SBITS(re, &s->gb, 8); SKIP_BITS(re, &s->gb, 8);
                if (level == -128) {
                    level = SHOW_UBITS(re, &s->gb, 8) - 256; LAST_SKIP_BITS(re, &s->gb, 8);
                } else if (level == 0) {
                    level = SHOW_UBITS(re, &s->gb, 8)      ; LAST_SKIP_BITS(re, &s->gb, 8);
                }
                i += run;
                j = scantable[i];
                if(level<0){
                    level= -level;
                    level= ((level*2+1)*qscale*quant_matrix[j])>>4;
                    level= (level-1)|1;
                    level= -level;
                }else{
                    level= ((level*2+1)*qscale*quant_matrix[j])>>4;
                    level= (level-1)|1;
                }
F
Fabrice Bellard 已提交
1537
            }
1538
            if (i > 63){
1539
                av_log(s->avctx, AV_LOG_ERROR, "ac-tex damaged at %d %d\n", s->mb_x, s->mb_y);
1540 1541 1542 1543
                return -1;
            }

            block[j] = level;
F
Fabrice Bellard 已提交
1544
        }
1545
        CLOSE_READER(re, &s->gb);
F
Fabrice Bellard 已提交
1546
    }
1547
    s->block_last_index[n] = i;
F
Fabrice Bellard 已提交
1548 1549 1550 1551 1552
    return 0;
}

/* Also does unquantization here, since I will never support mpeg2
   encoding */
1553 1554 1555
static inline int mpeg2_decode_block_non_intra(MpegEncContext *s, 
                               DCTELEM *block, 
                               int n)
F
Fabrice Bellard 已提交
1556 1557 1558
{
    int level, i, j, run;
    RLTable *rl = &rl_mpeg1;
1559 1560
    uint8_t * const scantable= s->intra_scantable.permutated;
    const uint16_t *quant_matrix;
1561
    const int qscale= s->qscale;
F
Fabrice Bellard 已提交
1562 1563 1564 1565 1566
    int mismatch;

    mismatch = 1;

    {
1567 1568
        int v;
        OPEN_READER(re, &s->gb);
1569
        i = -1;
1570
        if (n < 4)
1571
            quant_matrix = s->inter_matrix;
F
Fabrice Bellard 已提交
1572
        else
1573
            quant_matrix = s->chroma_inter_matrix;
1574

F
Fabrice Bellard 已提交
1575
        /* special case for the first coef. no need to add a second vlc table */
1576 1577
        UPDATE_CACHE(re, &s->gb);
        v= SHOW_UBITS(re, &s->gb, 2);
F
Fabrice Bellard 已提交
1578
        if (v & 2) {
1579 1580 1581 1582 1583 1584 1585
            LAST_SKIP_BITS(re, &s->gb, 2);
            level= (3*qscale*quant_matrix[0])>>5;
            if(v&1)
                level= -level;
            block[0] = level;
            mismatch ^= level;
            i++;
F
Fabrice Bellard 已提交
1586 1587
        }

1588 1589 1590 1591 1592 1593 1594 1595 1596 1597 1598 1599 1600 1601 1602 1603 1604 1605 1606 1607 1608 1609 1610 1611 1612 1613 1614 1615 1616
        /* now quantify & encode AC coefs */
        for(;;) {
            UPDATE_CACHE(re, &s->gb);
            GET_RL_VLC(level, run, re, &s->gb, rl->rl_vlc[0], TEX_VLC_BITS, 2);
            
            if(level == 127){
                break;
            } else if(level != 0) {
                i += run;
                j = scantable[i];
                level= ((level*2+1)*qscale*quant_matrix[j])>>5;
                level = (level ^ SHOW_SBITS(re, &s->gb, 1)) - SHOW_SBITS(re, &s->gb, 1);
                LAST_SKIP_BITS(re, &s->gb, 1);
            } else {
                /* escape */
                run = SHOW_UBITS(re, &s->gb, 6)+1; LAST_SKIP_BITS(re, &s->gb, 6);
                UPDATE_CACHE(re, &s->gb);
                level = SHOW_SBITS(re, &s->gb, 12); SKIP_BITS(re, &s->gb, 12);

                i += run;
                j = scantable[i];
                if(level<0){
                    level= ((-level*2+1)*qscale*quant_matrix[j])>>5;
                    level= -level;
                }else{
                    level= ((level*2+1)*qscale*quant_matrix[j])>>5;
                }
            }
            if (i > 63){
1617
                av_log(s->avctx, AV_LOG_ERROR, "ac-tex damaged at %d %d\n", s->mb_x, s->mb_y);
1618 1619 1620 1621 1622
                return -1;
            }
            
            mismatch ^= level;
            block[j] = level;
1623
        }
1624
        CLOSE_READER(re, &s->gb);
F
Fabrice Bellard 已提交
1625 1626
    }
    block[63] ^= (mismatch & 1);
1627
    
F
Fabrice Bellard 已提交
1628 1629 1630 1631
    s->block_last_index[n] = i;
    return 0;
}

1632 1633 1634
static inline int mpeg2_decode_block_intra(MpegEncContext *s, 
                               DCTELEM *block, 
                               int n)
F
Fabrice Bellard 已提交
1635 1636
{
    int level, dc, diff, i, j, run;
1637
    int component;
F
Fabrice Bellard 已提交
1638
    RLTable *rl;
1639 1640
    uint8_t * const scantable= s->intra_scantable.permutated;
    const uint16_t *quant_matrix;
1641
    const int qscale= s->qscale;
F
Fabrice Bellard 已提交
1642 1643 1644
    int mismatch;

    /* DC coef */
1645 1646 1647 1648 1649 1650 1651
    if (n < 4){
        quant_matrix = s->intra_matrix;
        component = 0; 
    }else{
        quant_matrix = s->chroma_intra_matrix;
        component = n - 3;
    }
1652
    diff = decode_dc(&s->gb, component);
F
Fabrice Bellard 已提交
1653 1654 1655 1656 1657 1658 1659
    if (diff >= 0xffff)
        return -1;
    dc = s->last_dc[component];
    dc += diff;
    s->last_dc[component] = dc;
    block[0] = dc << (3 - s->intra_dc_precision);
    dprintf("dc=%d\n", block[0]);
1660
    mismatch = block[0] ^ 1;
1661
    i = 0;
F
Fabrice Bellard 已提交
1662 1663 1664 1665
    if (s->intra_vlc_format)
        rl = &rl_mpeg2;
    else
        rl = &rl_mpeg1;
1666

1667 1668 1669 1670 1671 1672 1673 1674 1675 1676 1677 1678 1679 1680 1681 1682 1683 1684 1685 1686 1687 1688 1689 1690 1691 1692 1693 1694 1695 1696
    {
        OPEN_READER(re, &s->gb);    
        /* now quantify & encode AC coefs */
        for(;;) {
            UPDATE_CACHE(re, &s->gb);
            GET_RL_VLC(level, run, re, &s->gb, rl->rl_vlc[0], TEX_VLC_BITS, 2);
            
            if(level == 127){
                break;
            } else if(level != 0) {
                i += run;
                j = scantable[i];
                level= (level*qscale*quant_matrix[j])>>4;
                level = (level ^ SHOW_SBITS(re, &s->gb, 1)) - SHOW_SBITS(re, &s->gb, 1);
                LAST_SKIP_BITS(re, &s->gb, 1);
            } else {
                /* escape */
                run = SHOW_UBITS(re, &s->gb, 6)+1; LAST_SKIP_BITS(re, &s->gb, 6);
                UPDATE_CACHE(re, &s->gb);
                level = SHOW_SBITS(re, &s->gb, 12); SKIP_BITS(re, &s->gb, 12);
                i += run;
                j = scantable[i];
                if(level<0){
                    level= (-level*qscale*quant_matrix[j])>>4;
                    level= -level;
                }else{
                    level= (level*qscale*quant_matrix[j])>>4;
                }
            }
            if (i > 63){
1697
                av_log(s->avctx, AV_LOG_ERROR, "ac-tex damaged at %d %d\n", s->mb_x, s->mb_y);
1698 1699 1700 1701 1702
                return -1;
            }
            
            mismatch^= level;
            block[j] = level;
M
Michael Niedermayer 已提交
1703
        }
1704
        CLOSE_READER(re, &s->gb);
F
Fabrice Bellard 已提交
1705
    }
1706 1707
    block[63]^= mismatch&1;
    
F
Fabrice Bellard 已提交
1708 1709 1710 1711 1712 1713 1714
    s->block_last_index[n] = i;
    return 0;
}

typedef struct Mpeg1Context {
    MpegEncContext mpeg_enc_ctx;
    int mpeg_enc_ctx_allocated; /* true if decoding context allocated */
1715
    int repeat_field; /* true if we must repeat the field */
1716
    AVPanScan pan_scan; /** some temporary storage for the panscan */
F
Fabrice Bellard 已提交
1717 1718 1719 1720 1721
} Mpeg1Context;

static int mpeg_decode_init(AVCodecContext *avctx)
{
    Mpeg1Context *s = avctx->priv_data;
1722
    
M
dr1  
Michael Niedermayer 已提交
1723
    s->mpeg_enc_ctx.flags= avctx->flags;
1724
    common_init(&s->mpeg_enc_ctx);
1725
    init_vlcs();
F
Fabrice Bellard 已提交
1726 1727 1728

    s->mpeg_enc_ctx_allocated = 0;
    s->mpeg_enc_ctx.picture_number = 0;
1729
    s->repeat_field = 0;
1730
    s->mpeg_enc_ctx.codec_id= avctx->codec->id;
F
Fabrice Bellard 已提交
1731 1732 1733 1734 1735
    return 0;
}

/* return the 8 bit start code value and update the search
   state. Return -1 if no start code found */
1736
static int find_start_code(uint8_t **pbuf_ptr, uint8_t *buf_end)
F
Fabrice Bellard 已提交
1737
{
1738
    uint8_t *buf_ptr;
1739
    unsigned int state=0xFFFFFFFF, v;
F
Fabrice Bellard 已提交
1740 1741 1742 1743 1744 1745 1746 1747 1748 1749 1750 1751 1752 1753 1754 1755 1756 1757 1758
    int val;

    buf_ptr = *pbuf_ptr;
    while (buf_ptr < buf_end) {
        v = *buf_ptr++;
        if (state == 0x000001) {
            state = ((state << 8) | v) & 0xffffff;
            val = state;
            goto found;
        }
        state = ((state << 8) | v) & 0xffffff;
    }
    val = -1;
 found:
    *pbuf_ptr = buf_ptr;
    return val;
}

static int mpeg1_decode_picture(AVCodecContext *avctx, 
1759
                                uint8_t *buf, int buf_size)
F
Fabrice Bellard 已提交
1760 1761 1762
{
    Mpeg1Context *s1 = avctx->priv_data;
    MpegEncContext *s = &s1->mpeg_enc_ctx;
M
Michael Niedermayer 已提交
1763
    int ref, f_code, vbv_delay;
F
Fabrice Bellard 已提交
1764

1765
    init_get_bits(&s->gb, buf, buf_size*8);
F
Fabrice Bellard 已提交
1766 1767 1768

    ref = get_bits(&s->gb, 10); /* temporal ref */
    s->pict_type = get_bits(&s->gb, 3);
1769
    dprintf("pict_type=%d number=%d\n", s->pict_type, s->picture_number);
M
Michael Niedermayer 已提交
1770

M
Michael Niedermayer 已提交
1771
    vbv_delay= get_bits(&s->gb, 16);
F
Fabrice Bellard 已提交
1772
    if (s->pict_type == P_TYPE || s->pict_type == B_TYPE) {
1773
        s->full_pel[0] = get_bits1(&s->gb);
F
Fabrice Bellard 已提交
1774 1775 1776 1777 1778 1779 1780
        f_code = get_bits(&s->gb, 3);
        if (f_code == 0)
            return -1;
        s->mpeg_f_code[0][0] = f_code;
        s->mpeg_f_code[0][1] = f_code;
    }
    if (s->pict_type == B_TYPE) {
1781
        s->full_pel[1] = get_bits1(&s->gb);
F
Fabrice Bellard 已提交
1782 1783 1784 1785 1786 1787
        f_code = get_bits(&s->gb, 3);
        if (f_code == 0)
            return -1;
        s->mpeg_f_code[1][0] = f_code;
        s->mpeg_f_code[1][1] = f_code;
    }
M
cleanup  
Michael Niedermayer 已提交
1788 1789
    s->current_picture.pict_type= s->pict_type;
    s->current_picture.key_frame= s->pict_type == I_TYPE;
1790
    
F
Fabrice Bellard 已提交
1791 1792 1793 1794 1795 1796 1797 1798 1799
    s->y_dc_scale = 8;
    s->c_dc_scale = 8;
    s->first_slice = 1;
    return 0;
}

static void mpeg_decode_sequence_extension(MpegEncContext *s)
{
    int horiz_size_ext, vert_size_ext;
M
Michael Niedermayer 已提交
1800
    int bit_rate_ext, vbv_buf_ext;
F
Fabrice Bellard 已提交
1801
    int frame_rate_ext_n, frame_rate_ext_d;
M
Michael Niedermayer 已提交
1802
    int level, profile;
F
Fabrice Bellard 已提交
1803

M
Michael Niedermayer 已提交
1804 1805 1806
    skip_bits(&s->gb, 1); /* profil and level esc*/
    profile= get_bits(&s->gb, 3);
    level= get_bits(&s->gb, 4);
1807
    s->progressive_sequence = get_bits1(&s->gb); /* progressive_sequence */
1808
    skip_bits(&s->gb, 2); /* chroma_format */
F
Fabrice Bellard 已提交
1809 1810 1811 1812 1813 1814
    horiz_size_ext = get_bits(&s->gb, 2);
    vert_size_ext = get_bits(&s->gb, 2);
    s->width |= (horiz_size_ext << 12);
    s->height |= (vert_size_ext << 12);
    bit_rate_ext = get_bits(&s->gb, 12);  /* XXX: handle it */
    s->bit_rate = ((s->bit_rate / 400) | (bit_rate_ext << 12)) * 400;
1815
    skip_bits1(&s->gb); /* marker */
F
Fabrice Bellard 已提交
1816
    vbv_buf_ext = get_bits(&s->gb, 8);
M
Michael Niedermayer 已提交
1817

M
Michael Niedermayer 已提交
1818
    s->low_delay = get_bits1(&s->gb);
M
Michael Niedermayer 已提交
1819 1820
    if(s->flags & CODEC_FLAG_LOW_DELAY) s->low_delay=1;

F
Fabrice Bellard 已提交
1821 1822
    frame_rate_ext_n = get_bits(&s->gb, 2);
    frame_rate_ext_d = get_bits(&s->gb, 5);
1823 1824 1825 1826 1827 1828 1829
    av_reduce(
        &s->avctx->frame_rate, 
        &s->avctx->frame_rate_base, 
        frame_rate_tab[s->frame_rate_index] * (frame_rate_ext_n+1),
        MPEG1_FRAME_RATE_BASE * (frame_rate_ext_d+1),
        1<<30);

F
Fabrice Bellard 已提交
1830
    dprintf("sequence extension\n");
M
Michael Niedermayer 已提交
1831
    s->codec_id= s->avctx->codec_id= CODEC_ID_MPEG2VIDEO;
1832
    s->avctx->sub_id = 2; /* indicates mpeg2 found */
M
Michael Niedermayer 已提交
1833

M
Michael Niedermayer 已提交
1834 1835 1836 1837 1838 1839 1840 1841 1842
    if(s->aspect_ratio_info <= 1)
        s->avctx->sample_aspect_ratio= mpeg2_aspect[s->aspect_ratio_info];
    else{
        s->avctx->sample_aspect_ratio= 
            av_div_q(
                mpeg2_aspect[s->aspect_ratio_info], 
                (AVRational){s->width, s->height}
            );
    }
M
Michael Niedermayer 已提交
1843 1844
    
    if(s->avctx->debug & FF_DEBUG_PICT_INFO)
1845
        av_log(s->avctx, AV_LOG_DEBUG, "profile: %d, level: %d \n", profile, level);
F
Fabrice Bellard 已提交
1846 1847
}

1848 1849 1850 1851 1852 1853 1854 1855 1856 1857 1858 1859 1860 1861 1862 1863 1864 1865 1866 1867
static void mpeg_decode_sequence_display_extension(Mpeg1Context *s1)
{
    MpegEncContext *s= &s1->mpeg_enc_ctx;
    int color_description, w, h;

    skip_bits(&s->gb, 3); /* video format */
    color_description= get_bits1(&s->gb);
    if(color_description){
        skip_bits(&s->gb, 8); /* color primaries */
        skip_bits(&s->gb, 8); /* transfer_characteristics */
        skip_bits(&s->gb, 8); /* matrix_coefficients */
    }
    w= get_bits(&s->gb, 14);
    skip_bits(&s->gb, 1); //marker
    h= get_bits(&s->gb, 14);
    skip_bits(&s->gb, 1); //marker
    
    s1->pan_scan.width= 16*w;
    s1->pan_scan.height=16*h;

M
Michael Niedermayer 已提交
1868 1869 1870 1871 1872 1873
    if(s->aspect_ratio_info > 1)
        s->avctx->sample_aspect_ratio= 
            av_div_q(
                mpeg2_aspect[s->aspect_ratio_info], 
                (AVRational){w, h}
            );
1874 1875
    
    if(s->avctx->debug & FF_DEBUG_PICT_INFO)
1876
        av_log(s->avctx, AV_LOG_DEBUG, "sde w:%d, h:%d\n", w, h);
1877 1878 1879 1880 1881 1882 1883 1884 1885 1886 1887 1888 1889 1890 1891
}

static void mpeg_decode_picture_display_extension(Mpeg1Context *s1)
{
    MpegEncContext *s= &s1->mpeg_enc_ctx;
    int i;

    for(i=0; i<1; i++){ //FIXME count
        s1->pan_scan.position[i][0]= get_sbits(&s->gb, 16);
        skip_bits(&s->gb, 1); //marker
        s1->pan_scan.position[i][1]= get_sbits(&s->gb, 16);
        skip_bits(&s->gb, 1); //marker
    }
   
    if(s->avctx->debug & FF_DEBUG_PICT_INFO)
1892
        av_log(s->avctx, AV_LOG_DEBUG, "pde (%d,%d) (%d,%d) (%d,%d)\n", 
1893 1894 1895 1896 1897 1898
            s1->pan_scan.position[0][0], s1->pan_scan.position[0][1], 
            s1->pan_scan.position[1][0], s1->pan_scan.position[1][1], 
            s1->pan_scan.position[2][0], s1->pan_scan.position[2][1]
        );
}

F
Fabrice Bellard 已提交
1899 1900
static void mpeg_decode_quant_matrix_extension(MpegEncContext *s)
{
1901
    int i, v, j;
F
Fabrice Bellard 已提交
1902

1903 1904
    dprintf("matrix extension\n");

1905
    if (get_bits1(&s->gb)) {
F
Fabrice Bellard 已提交
1906 1907
        for(i=0;i<64;i++) {
            v = get_bits(&s->gb, 8);
1908
            j= s->dsp.idct_permutation[ ff_zigzag_direct[i] ];
1909 1910
            s->intra_matrix[j] = v;
            s->chroma_intra_matrix[j] = v;
F
Fabrice Bellard 已提交
1911 1912
        }
    }
1913
    if (get_bits1(&s->gb)) {
F
Fabrice Bellard 已提交
1914 1915
        for(i=0;i<64;i++) {
            v = get_bits(&s->gb, 8);
1916
            j= s->dsp.idct_permutation[ ff_zigzag_direct[i] ];
1917 1918
            s->inter_matrix[j] = v;
            s->chroma_inter_matrix[j] = v;
F
Fabrice Bellard 已提交
1919 1920
        }
    }
1921
    if (get_bits1(&s->gb)) {
F
Fabrice Bellard 已提交
1922 1923
        for(i=0;i<64;i++) {
            v = get_bits(&s->gb, 8);
1924
            j= s->dsp.idct_permutation[ ff_zigzag_direct[i] ];
1925
            s->chroma_intra_matrix[j] = v;
F
Fabrice Bellard 已提交
1926 1927
        }
    }
1928
    if (get_bits1(&s->gb)) {
F
Fabrice Bellard 已提交
1929 1930
        for(i=0;i<64;i++) {
            v = get_bits(&s->gb, 8);
1931
            j= s->dsp.idct_permutation[ ff_zigzag_direct[i] ];
1932
            s->chroma_inter_matrix[j] = v;
F
Fabrice Bellard 已提交
1933 1934 1935 1936 1937 1938 1939 1940 1941 1942 1943 1944 1945
        }
    }
}

static void mpeg_decode_picture_coding_extension(MpegEncContext *s)
{
    s->full_pel[0] = s->full_pel[1] = 0;
    s->mpeg_f_code[0][0] = get_bits(&s->gb, 4);
    s->mpeg_f_code[0][1] = get_bits(&s->gb, 4);
    s->mpeg_f_code[1][0] = get_bits(&s->gb, 4);
    s->mpeg_f_code[1][1] = get_bits(&s->gb, 4);
    s->intra_dc_precision = get_bits(&s->gb, 2);
    s->picture_structure = get_bits(&s->gb, 2);
1946 1947 1948 1949 1950 1951 1952 1953 1954
    s->top_field_first = get_bits1(&s->gb);
    s->frame_pred_frame_dct = get_bits1(&s->gb);
    s->concealment_motion_vectors = get_bits1(&s->gb);
    s->q_scale_type = get_bits1(&s->gb);
    s->intra_vlc_format = get_bits1(&s->gb);
    s->alternate_scan = get_bits1(&s->gb);
    s->repeat_first_field = get_bits1(&s->gb);
    s->chroma_420_type = get_bits1(&s->gb);
    s->progressive_frame = get_bits1(&s->gb);
1955
    
1956 1957 1958 1959
    if(s->picture_structure == PICT_FRAME)
        s->first_field=0;
    else{
        s->first_field ^= 1;
1960
        memset(s->mbskip_table, 0, s->mb_stride*s->mb_height);
1961 1962
    }
    
1963
    if(s->alternate_scan){
M
Michael Niedermayer 已提交
1964 1965 1966 1967
        ff_init_scantable(s->dsp.idct_permutation, &s->inter_scantable  , ff_alternate_vertical_scan);
        ff_init_scantable(s->dsp.idct_permutation, &s->intra_scantable  , ff_alternate_vertical_scan);
        ff_init_scantable(s->dsp.idct_permutation, &s->intra_h_scantable, ff_alternate_vertical_scan);
        ff_init_scantable(s->dsp.idct_permutation, &s->intra_v_scantable, ff_alternate_vertical_scan);
1968
    }else{
M
Michael Niedermayer 已提交
1969 1970 1971 1972
        ff_init_scantable(s->dsp.idct_permutation, &s->inter_scantable  , ff_zigzag_direct);
        ff_init_scantable(s->dsp.idct_permutation, &s->intra_scantable  , ff_zigzag_direct);
        ff_init_scantable(s->dsp.idct_permutation, &s->intra_h_scantable, ff_alternate_horizontal_scan);
        ff_init_scantable(s->dsp.idct_permutation, &s->intra_v_scantable, ff_alternate_vertical_scan);
1973 1974
    }
    
F
Fabrice Bellard 已提交
1975
    /* composite display not parsed */
1976
    dprintf("intra_dc_precision=%d\n", s->intra_dc_precision);
1977
    dprintf("picture_structure=%d\n", s->picture_structure);
1978 1979
    dprintf("top field first=%d\n", s->top_field_first);
    dprintf("repeat first field=%d\n", s->repeat_first_field);
F
Fabrice Bellard 已提交
1980
    dprintf("conceal=%d\n", s->concealment_motion_vectors);
1981 1982
    dprintf("intra_vlc_format=%d\n", s->intra_vlc_format);
    dprintf("alternate_scan=%d\n", s->alternate_scan);
F
Fabrice Bellard 已提交
1983
    dprintf("frame_pred_frame_dct=%d\n", s->frame_pred_frame_dct);
1984
    dprintf("progressive_frame=%d\n", s->progressive_frame);
F
Fabrice Bellard 已提交
1985 1986 1987
}

static void mpeg_decode_extension(AVCodecContext *avctx, 
1988
                                  uint8_t *buf, int buf_size)
F
Fabrice Bellard 已提交
1989 1990 1991 1992 1993
{
    Mpeg1Context *s1 = avctx->priv_data;
    MpegEncContext *s = &s1->mpeg_enc_ctx;
    int ext_type;

1994
    init_get_bits(&s->gb, buf, buf_size*8);
F
Fabrice Bellard 已提交
1995 1996 1997 1998 1999 2000
    
    ext_type = get_bits(&s->gb, 4);
    switch(ext_type) {
    case 0x1:
        mpeg_decode_sequence_extension(s);
        break;
2001 2002 2003
    case 0x2:
        mpeg_decode_sequence_display_extension(s1);
        break;
F
Fabrice Bellard 已提交
2004 2005 2006
    case 0x3:
        mpeg_decode_quant_matrix_extension(s);
        break;
2007 2008 2009
    case 0x7:
        mpeg_decode_picture_display_extension(s1);
        break;
F
Fabrice Bellard 已提交
2010 2011 2012 2013 2014 2015
    case 0x8:
        mpeg_decode_picture_coding_extension(s);
        break;
    }
}

2016 2017 2018 2019 2020 2021
static void exchange_uv(MpegEncContext *s){
short * tmp;

    tmp = s->pblocks[4];
    s->pblocks[4] = s->pblocks[5];
    s->pblocks[5] = tmp;
M
Michael Niedermayer 已提交
2022 2023
}

M
cleanup  
Michael Niedermayer 已提交
2024 2025 2026 2027 2028 2029 2030 2031 2032 2033
#define DECODE_SLICE_FATAL_ERROR -2
#define DECODE_SLICE_ERROR -1
#define DECODE_SLICE_OK 0

/**
 * decodes a slice.
 * @return DECODE_SLICE_FATAL_ERROR if a non recoverable error occured<br>
 *         DECODE_SLICE_ERROR if the slice is damaged<br>
 *         DECODE_SLICE_OK if this slice is ok<br>
 */
F
Fabrice Bellard 已提交
2034
static int mpeg_decode_slice(AVCodecContext *avctx, 
M
Michael Niedermayer 已提交
2035
                              AVFrame *pict,
F
Fabrice Bellard 已提交
2036
                              int start_code,
2037
                              uint8_t **buf, int buf_size)
F
Fabrice Bellard 已提交
2038 2039 2040 2041
{
    Mpeg1Context *s1 = avctx->priv_data;
    MpegEncContext *s = &s1->mpeg_enc_ctx;
    int ret;
2042
    const int field_pic= s->picture_structure != PICT_FRAME;
F
Fabrice Bellard 已提交
2043

2044 2045 2046
    s->resync_mb_x= s->mb_x = 
    s->resync_mb_y= s->mb_y = -1;
    
F
Fabrice Bellard 已提交
2047
    start_code = (start_code - 1) & 0xff;
M
Michael Niedermayer 已提交
2048
    if (start_code >= s->mb_height){
2049
        av_log(s->avctx, AV_LOG_ERROR, "slice below image (%d >= %d)\n", start_code, s->mb_height);
2050
        return -1;
M
Michael Niedermayer 已提交
2051
    }
2052 2053
    
    ff_mpeg1_clean_buffers(s);
2054
    s->interlaced_dct = 0;
2055
        
F
Fabrice Bellard 已提交
2056
    /* start frame decoding */
M
Michael Niedermayer 已提交
2057 2058
    if (s->first_slice) {
      if(s->first_field || s->picture_structure==PICT_FRAME){
2059
        if(MPV_frame_start(s, avctx) < 0)
M
cleanup  
Michael Niedermayer 已提交
2060
            return DECODE_SLICE_FATAL_ERROR;
2061 2062 2063

        ff_er_frame_start(s);

2064
        /* first check if we must repeat the frame */
F
Fabrice Bellard 已提交
2065
        s->current_picture_ptr->repeat_pict = 0;
2066 2067 2068
        if (s->repeat_first_field) {
            if (s->progressive_sequence) {
                if (s->top_field_first)
F
Fabrice Bellard 已提交
2069
                    s->current_picture_ptr->repeat_pict = 4;
2070
                else
F
Fabrice Bellard 已提交
2071
                    s->current_picture_ptr->repeat_pict = 2;
2072
            } else if (s->progressive_frame) {
F
Fabrice Bellard 已提交
2073
                s->current_picture_ptr->repeat_pict = 1;
2074 2075
            }
        }         
2076 2077

        *s->current_picture_ptr->pan_scan= s1->pan_scan;
M
Michael Niedermayer 已提交
2078 2079
      }else{ //second field
            int i;
2080 2081
            
            if(!s->current_picture_ptr){
2082
                av_log(s->avctx, AV_LOG_ERROR, "first field missing\n");
2083 2084 2085
                return -1;
            }
            
M
Michael Niedermayer 已提交
2086 2087 2088 2089 2090 2091 2092
            for(i=0; i<4; i++){
                s->current_picture.data[i] = s->current_picture_ptr->data[i];
                if(s->picture_structure == PICT_BOTTOM_FIELD){
                    s->current_picture.data[i] += s->current_picture_ptr->linesize[i];
                } 
            }
      }
I
Ivan Kalvachev 已提交
2093 2094 2095 2096 2097 2098 2099
#ifdef HAVE_XVMC
// MPV_frame_start will call this function too,
// but we need to call it on every field
      if(s->avctx->xvmc_acceleration)
         XVMC_field_start(s,avctx);
#endif
    }//fi(s->first_slice)
F
Fabrice Bellard 已提交
2100

2101
    init_get_bits(&s->gb, *buf, buf_size*8);
F
Fabrice Bellard 已提交
2102

F
Fabrice Bellard 已提交
2103
    s->qscale = get_qscale(s);
2104 2105 2106 2107 2108 2109 2110 2111 2112 2113 2114 2115
    if (s->first_slice && (s->first_field || s->picture_structure==PICT_FRAME)) {
        if(s->avctx->debug&FF_DEBUG_PICT_INFO){
             av_log(s->avctx, AV_LOG_DEBUG, "qp:%d fc:%2d%2d%2d%2d %s %s %s %s dc:%d pstruct:%d fdct:%d cmv:%d qtype:%d ivlc:%d rff:%d %s\n", 
                 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],
                 s->pict_type == I_TYPE ? "I" : (s->pict_type == P_TYPE ? "P" : (s->pict_type == B_TYPE ? "B" : "S")), 
                 s->progressive_sequence ? "pro" :"", s->alternate_scan ? "alt" :"", s->top_field_first ? "top" :"", 
                 s->intra_dc_precision, s->picture_structure, s->frame_pred_frame_dct, s->concealment_motion_vectors,
                 s->q_scale_type, s->intra_vlc_format, s->repeat_first_field, s->chroma_420_type ? "420" :"");
        }
    }

    s->first_slice = 0;
2116
    if(s->qscale == 0){
2117
        av_log(s->avctx, AV_LOG_ERROR, "qscale == 0\n");
2118
        return -1;
2119 2120
    }
    
F
Fabrice Bellard 已提交
2121
    /* extra slice info */
2122 2123
    while (get_bits1(&s->gb) != 0) {
        skip_bits(&s->gb, 8);
F
Fabrice Bellard 已提交
2124
    }
2125
    
M
Michael Niedermayer 已提交
2126
    s->mb_x=0;
2127

M
Michael Niedermayer 已提交
2128 2129
    for(;;) {
        int code = get_vlc2(&s->gb, mbincr_vlc.table, MBINCR_VLC_BITS, 2);
2130
        if (code < 0){
2131
            av_log(s->avctx, AV_LOG_ERROR, "first mb_incr damaged\n");
2132
            return -1;
2133
        }
M
Michael Niedermayer 已提交
2134 2135 2136 2137 2138 2139 2140 2141 2142 2143
        if (code >= 33) {
            if (code == 33) {
                s->mb_x += 33;
            }
            /* otherwise, stuffing, nothing to do */
        } else {
            s->mb_x += code;
            break;
        }
    }
2144
    
2145 2146
    s->resync_mb_x= s->mb_x;
    s->resync_mb_y= s->mb_y = start_code;
2147
    s->mb_skip_run= 0;
2148
    ff_init_block_index(s);
M
Michael Niedermayer 已提交
2149

F
Fabrice Bellard 已提交
2150
    for(;;) {
2151 2152 2153 2154 2155 2156
#ifdef HAVE_XVMC
        //one 1 we memcpy blocks in xvmcvideo
        if(s->avctx->xvmc_acceleration > 1)
            XVMC_init_block(s);//set s->block
#endif

2157
	s->dsp.clear_blocks(s->block[0]);
2158

2159
        ret = mpeg_decode_mb(s, s->block);
M
Michael Niedermayer 已提交
2160
        s->chroma_qscale= s->qscale;
2161

F
Fabrice Bellard 已提交
2162 2163 2164
        dprintf("ret=%d\n", ret);
        if (ret < 0)
            return -1;
2165

2166 2167 2168
        if(s->current_picture.motion_val[0] && !s->encoding){ //note motion_val is normally NULL unless we want to extract the MVs
            const int wrap = field_pic ? 2*s->block_wrap[0] : s->block_wrap[0];
            int xy = s->mb_x*2 + 1 + (s->mb_y*2 +1)*wrap;
2169 2170
            int motion_for_top_x, motion_for_top_y, motion_back_top_x, motion_back_top_y;
            int motion_for_bottom_x, motion_for_bottom_y, motion_back_bottom_x, motion_back_bottom_y;
2171 2172 2173
            if(field_pic && !s->first_field)
                xy += wrap/2;

2174
            if (s->mb_intra) {
2175 2176
                motion_for_top_x = motion_for_top_y = motion_back_top_x = motion_back_top_y =
                motion_for_bottom_x = motion_for_bottom_y = motion_back_bottom_x = motion_back_bottom_y = 0;
2177
            }else if (s->mv_type == MV_TYPE_16X16){
2178 2179 2180 2181
                motion_for_top_x = motion_for_bottom_x = s->mv[0][0][0];
                motion_for_top_y = motion_for_bottom_y = s->mv[0][0][1];
                motion_back_top_x = motion_back_bottom_x = s->mv[1][0][0];
                motion_back_top_y = motion_back_bottom_y = s->mv[1][0][1];
2182
            } else /*if ((s->mv_type == MV_TYPE_FIELD) || (s->mv_type == MV_TYPE_16X8))*/ {
2183 2184 2185 2186 2187 2188 2189 2190
                motion_for_top_x = s->mv[0][0][0];
                motion_for_top_y = s->mv[0][0][1];
                motion_for_bottom_x = s->mv[0][1][0];
                motion_for_bottom_y = s->mv[0][1][1];
                motion_back_top_x = s->mv[1][0][0];
                motion_back_top_y = s->mv[1][0][1];
                motion_back_bottom_x = s->mv[1][1][0];
                motion_back_bottom_y = s->mv[1][1][1];
2191 2192
            }

2193 2194 2195 2196 2197 2198 2199 2200
            s->current_picture.motion_val[0][xy][0] = motion_for_top_x;
            s->current_picture.motion_val[0][xy][1] = motion_for_top_y;
            s->current_picture.motion_val[0][xy + 1][0] = motion_for_top_x;
            s->current_picture.motion_val[0][xy + 1][1] = motion_for_top_y;
            s->current_picture.motion_val[0][xy + wrap][0] = motion_for_bottom_x;
            s->current_picture.motion_val[0][xy + wrap][1] = motion_for_bottom_y;
            s->current_picture.motion_val[0][xy + 1 + wrap][0] = motion_for_bottom_x;
            s->current_picture.motion_val[0][xy + 1 + wrap][1] = motion_for_bottom_y;
2201 2202

            if(s->pict_type != B_TYPE){
2203
                motion_back_top_x = motion_back_top_y = motion_back_bottom_x = motion_back_bottom_y = 0;
2204
            }
2205

2206 2207 2208 2209 2210 2211 2212 2213
            s->current_picture.motion_val[1][xy][0] = motion_back_top_x;
            s->current_picture.motion_val[1][xy][1] = motion_back_top_y;
            s->current_picture.motion_val[1][xy + 1][0] = motion_back_top_x;
            s->current_picture.motion_val[1][xy + 1][1] = motion_back_top_y;
            s->current_picture.motion_val[1][xy + wrap][0] = motion_back_bottom_x;
            s->current_picture.motion_val[1][xy + wrap][1] = motion_back_bottom_y;
            s->current_picture.motion_val[1][xy + 1 + wrap][0] = motion_back_bottom_x;
            s->current_picture.motion_val[1][xy + 1 + wrap][1] = motion_back_bottom_y;
2214
        }
2215

2216 2217 2218
        s->dest[0] += 16;
        s->dest[1] += 8;
        s->dest[2] += 8;
M
Michael Niedermayer 已提交
2219

2220 2221
        MPV_decode_mb(s, s->block);
        
M
Michael Niedermayer 已提交
2222
        if (++s->mb_x >= s->mb_width) {
M
Michael Niedermayer 已提交
2223

2224
            ff_draw_horiz_band(s, 16*s->mb_y, 16);
M
Michael Niedermayer 已提交
2225 2226 2227

            s->mb_x = 0;
            s->mb_y++;
2228 2229 2230 2231 2232 2233

            if(s->mb_y<<field_pic >= s->mb_height){
                int left= s->gb.size_in_bits - get_bits_count(&s->gb);

                if(left < 0 || (left && show_bits(&s->gb, FFMIN(left, 23)))
                   || (avctx->error_resilience >= FF_ER_AGGRESSIVE && left>8)){
2234
                    av_log(avctx, AV_LOG_ERROR, "end mismatch left=%d\n", left);
2235 2236 2237 2238
                    return -1;
                }else
                    goto eos;
            }
2239 2240
            
            ff_init_block_index(s);
M
Michael Niedermayer 已提交
2241 2242 2243
        }

        /* skip mb handling */
2244
        if (s->mb_skip_run == -1) {
M
Michael Niedermayer 已提交
2245
            /* read again increment */
2246
            s->mb_skip_run = 0;
M
Michael Niedermayer 已提交
2247 2248
            for(;;) {
                int code = get_vlc2(&s->gb, mbincr_vlc.table, MBINCR_VLC_BITS, 2);
2249
                if (code < 0){
2250
                    av_log(s->avctx, AV_LOG_ERROR, "mb incr damaged\n");
2251
                    return -1;
2252
                }
M
Michael Niedermayer 已提交
2253 2254
                if (code >= 33) {
                    if (code == 33) {
2255
                        s->mb_skip_run += 33;
2256 2257
                    }else if(code == 35){
                        if(s->mb_skip_run != 0 || show_bits(&s->gb, 15) != 0){
2258
                            av_log(s->avctx, AV_LOG_ERROR, "slice mismatch\n");
2259 2260 2261
                            return -1;
                        }
                        goto eos; /* end of slice */
M
Michael Niedermayer 已提交
2262 2263 2264
                    }
                    /* otherwise, stuffing, nothing to do */
                } else {
2265
                    s->mb_skip_run += code;
M
Michael Niedermayer 已提交
2266 2267 2268 2269
                    break;
                }
            }
        }
F
Fabrice Bellard 已提交
2270
    }
2271 2272
eos: // end of slice
    *buf += get_bits_count(&s->gb)/8 - 1;
2273
//printf("y %d %d %d %d\n", s->resync_mb_x, s->resync_mb_y, s->mb_x, s->mb_y);
2274 2275
    return 0;
}
2276

2277 2278 2279 2280 2281 2282 2283 2284 2285
/**
 * handles slice ends.
 * @return 1 if it seems to be the last slice of 
 */
static int slice_end(AVCodecContext *avctx, AVFrame *pict)
{
    Mpeg1Context *s1 = avctx->priv_data;
    MpegEncContext *s = &s1->mpeg_enc_ctx;
       
2286
    if (!s1->mpeg_enc_ctx_allocated || !s->current_picture_ptr)
2287 2288
        return 0;

I
Ivan Kalvachev 已提交
2289 2290 2291 2292
#ifdef HAVE_XVMC
    if(s->avctx->xvmc_acceleration)
        XVMC_field_end(s);
#endif
F
Fabrice Bellard 已提交
2293
    /* end of slice reached */
2294
    if (/*s->mb_y<<field_pic == s->mb_height &&*/ !s->first_field) {
F
Fabrice Bellard 已提交
2295
        /* end of image */
M
cleanup  
Michael Niedermayer 已提交
2296

M
Michael Niedermayer 已提交
2297
        if(s->codec_id == CODEC_ID_MPEG2VIDEO){
M
Michael Niedermayer 已提交
2298 2299 2300 2301
            s->current_picture_ptr->qscale_type= FF_QSCALE_TYPE_MPEG2;
        }else
            s->current_picture_ptr->qscale_type= FF_QSCALE_TYPE_MPEG1;

2302
        ff_er_frame_end(s);
F
Fabrice Bellard 已提交
2303 2304 2305

        MPV_frame_end(s);

M
Michael Niedermayer 已提交
2306
        if (s->pict_type == B_TYPE || s->low_delay) {
M
2x100l  
Michael Niedermayer 已提交
2307
            *pict= *(AVFrame*)s->current_picture_ptr;
2308
            ff_print_debug_info(s, pict);
F
Fabrice Bellard 已提交
2309
        } else {
M
cleanup  
Michael Niedermayer 已提交
2310
            s->picture_number++;
F
Fabrice Bellard 已提交
2311 2312
            /* latency of 1 frame for I and P frames */
            /* XXX: use another variable than picture_number */
2313
            if (s->last_picture_ptr != NULL) {
M
2x100l  
Michael Niedermayer 已提交
2314
                *pict= *(AVFrame*)s->last_picture_ptr;
2315
                 ff_print_debug_info(s, pict);
F
Fabrice Bellard 已提交
2316 2317
            }
        }
M
Michael Niedermayer 已提交
2318

2319
        return 1;
F
Fabrice Bellard 已提交
2320
    } else {
2321
        return 0;
F
Fabrice Bellard 已提交
2322 2323 2324 2325
    }
}

static int mpeg1_decode_sequence(AVCodecContext *avctx, 
2326
                                 uint8_t *buf, int buf_size)
F
Fabrice Bellard 已提交
2327 2328 2329
{
    Mpeg1Context *s1 = avctx->priv_data;
    MpegEncContext *s = &s1->mpeg_enc_ctx;
2330
    int width, height, i, v, j;
M
Michael Niedermayer 已提交
2331
    float aspect;
2332

2333
    init_get_bits(&s->gb, buf, buf_size*8);
F
Fabrice Bellard 已提交
2334 2335 2336

    width = get_bits(&s->gb, 12);
    height = get_bits(&s->gb, 12);
M
Michael Niedermayer 已提交
2337
    s->aspect_ratio_info= get_bits(&s->gb, 4);
2338 2339 2340 2341
    if (s->aspect_ratio_info == 0)
        return -1;
    aspect= 1.0/mpeg1_aspect[s->aspect_ratio_info];
    avctx->sample_aspect_ratio= av_d2q(aspect, 255);
M
Michael Niedermayer 已提交
2342

F
Fabrice Bellard 已提交
2343 2344 2345 2346
    s->frame_rate_index = get_bits(&s->gb, 4);
    if (s->frame_rate_index == 0)
        return -1;
    s->bit_rate = get_bits(&s->gb, 18) * 400;
2347
    if (get_bits1(&s->gb) == 0) /* marker */
F
Fabrice Bellard 已提交
2348 2349 2350 2351 2352 2353 2354 2355 2356 2357 2358 2359 2360
        return -1;
    if (width <= 0 || height <= 0 ||
        (width % 2) != 0 || (height % 2) != 0)
        return -1;
    if (width != s->width ||
        height != s->height) {
        /* start new mpeg1 context decoding */
        s->out_format = FMT_MPEG1;
        if (s1->mpeg_enc_ctx_allocated) {
            MPV_common_end(s);
        }
        s->width = width;
        s->height = height;
M
Michael Niedermayer 已提交
2361
        avctx->has_b_frames= 1;
F
Fabrice Bellard 已提交
2362
        s->avctx = avctx;
F
Fabrice Bellard 已提交
2363 2364
        avctx->width = width;
        avctx->height = height;
2365 2366 2367 2368 2369 2370 2371
        av_reduce(
            &avctx->frame_rate, 
            &avctx->frame_rate_base,
            frame_rate_tab[s->frame_rate_index],
            MPEG1_FRAME_RATE_BASE, //FIXME store in allready reduced form 
            1<<30
            );
F
Fabrice Bellard 已提交
2372 2373
        avctx->bit_rate = s->bit_rate;
        
I
Ivan Kalvachev 已提交
2374 2375 2376
        //get_format() or set_video(width,height,aspect,pix_fmt);
        //until then pix_fmt may be changed right after codec init
        if( avctx->pix_fmt == PIX_FMT_XVMC_MPEG2_IDCT )
2377 2378
            if( avctx->idct_algo == FF_IDCT_AUTO )
                avctx->idct_algo = FF_IDCT_SIMPLE;
I
Ivan Kalvachev 已提交
2379

F
Fabrice Bellard 已提交
2380 2381 2382
        if (MPV_common_init(s) < 0)
            return -1;
        s1->mpeg_enc_ctx_allocated = 1;
2383
        s->swap_uv = 0;//just in case vcr2 and mpeg2 stream have been concatinated
F
Fabrice Bellard 已提交
2384 2385
    }

2386 2387
    skip_bits(&s->gb, 10); /* vbv_buffer_size */
    skip_bits(&s->gb, 1);
F
Fabrice Bellard 已提交
2388 2389

    /* get matrix */
2390
    if (get_bits1(&s->gb)) {
F
Fabrice Bellard 已提交
2391 2392
        for(i=0;i<64;i++) {
            v = get_bits(&s->gb, 8);
2393
            j = s->intra_scantable.permutated[i];
2394 2395
            s->intra_matrix[j] = v;
            s->chroma_intra_matrix[j] = v;
F
Fabrice Bellard 已提交
2396
        }
2397 2398 2399
#ifdef DEBUG
        dprintf("intra matrix present\n");
        for(i=0;i<64;i++)
2400
            dprintf(" %d", s->intra_matrix[s->intra_scantable.permutated[i]]);
2401 2402
        printf("\n");
#endif
F
Fabrice Bellard 已提交
2403 2404
    } else {
        for(i=0;i<64;i++) {
2405
            int j= s->dsp.idct_permutation[i];
2406
            v = ff_mpeg1_default_intra_matrix[i];
2407 2408
            s->intra_matrix[j] = v;
            s->chroma_intra_matrix[j] = v;
F
Fabrice Bellard 已提交
2409 2410
        }
    }
2411
    if (get_bits1(&s->gb)) {
F
Fabrice Bellard 已提交
2412 2413
        for(i=0;i<64;i++) {
            v = get_bits(&s->gb, 8);
2414
            j = s->intra_scantable.permutated[i];
2415 2416
            s->inter_matrix[j] = v;
            s->chroma_inter_matrix[j] = v;
F
Fabrice Bellard 已提交
2417
        }
2418 2419 2420
#ifdef DEBUG
        dprintf("non intra matrix present\n");
        for(i=0;i<64;i++)
2421
            dprintf(" %d", s->inter_matrix[s->intra_scantable.permutated[i]]);
2422 2423
        printf("\n");
#endif
F
Fabrice Bellard 已提交
2424 2425
    } else {
        for(i=0;i<64;i++) {
2426
            int j= s->dsp.idct_permutation[i];
2427
            v = ff_mpeg1_default_non_intra_matrix[i];
2428 2429
            s->inter_matrix[j] = v;
            s->chroma_inter_matrix[j] = v;
F
Fabrice Bellard 已提交
2430 2431 2432 2433 2434 2435 2436 2437
        }
    }

    /* we set mpeg2 parameters so that it emulates mpeg1 */
    s->progressive_sequence = 1;
    s->progressive_frame = 1;
    s->picture_structure = PICT_FRAME;
    s->frame_pred_frame_dct = 1;
M
Michael Niedermayer 已提交
2438
    s->codec_id= s->avctx->codec_id= CODEC_ID_MPEG1VIDEO;
2439
    avctx->sub_id = 1; /* indicates mpeg1 */
F
Fabrice Bellard 已提交
2440
    if(s->flags & CODEC_FLAG_LOW_DELAY) s->low_delay=1;
F
Fabrice Bellard 已提交
2441 2442 2443
    return 0;
}

M
Michael Niedermayer 已提交
2444 2445 2446 2447
static int vcr2_init_sequence(AVCodecContext *avctx)
{
    Mpeg1Context *s1 = avctx->priv_data;
    MpegEncContext *s = &s1->mpeg_enc_ctx;
M
Michael Niedermayer 已提交
2448
    int i, v;
M
Michael Niedermayer 已提交
2449 2450 2451 2452 2453 2454 2455 2456 2457

    /* start new mpeg1 context decoding */
    s->out_format = FMT_MPEG1;
    if (s1->mpeg_enc_ctx_allocated) {
        MPV_common_end(s);
    }
    s->width = avctx->width;
    s->height = avctx->height;
    avctx->has_b_frames= 0; //true?
M
Michael Niedermayer 已提交
2458
    s->low_delay= 1;
M
Michael Niedermayer 已提交
2459
    s->avctx = avctx;
I
Ivan Kalvachev 已提交
2460 2461 2462 2463

    //get_format() or set_video(width,height,aspect,pix_fmt);
    //until then pix_fmt may be changed right after codec init
    if( avctx->pix_fmt == PIX_FMT_XVMC_MPEG2_IDCT )
2464 2465
        if( avctx->idct_algo == FF_IDCT_AUTO )
            avctx->idct_algo = FF_IDCT_SIMPLE;
M
Michael Niedermayer 已提交
2466 2467 2468
    
    if (MPV_common_init(s) < 0)
        return -1;
2469 2470
    exchange_uv(s);//common init reset pblocks, so we swap them here
    s->swap_uv = 1;// in case of xvmc we need to swap uv for each MB 
M
Michael Niedermayer 已提交
2471 2472 2473 2474 2475 2476 2477 2478 2479 2480 2481 2482 2483 2484 2485 2486 2487
    s1->mpeg_enc_ctx_allocated = 1;

    for(i=0;i<64;i++) {
        int j= s->dsp.idct_permutation[i];
        v = ff_mpeg1_default_intra_matrix[i];
        s->intra_matrix[j] = v;
        s->chroma_intra_matrix[j] = v;

        v = ff_mpeg1_default_non_intra_matrix[i];
        s->inter_matrix[j] = v;
        s->chroma_inter_matrix[j] = v;
    }

    s->progressive_sequence = 1;
    s->progressive_frame = 1;
    s->picture_structure = PICT_FRAME;
    s->frame_pred_frame_dct = 1;
M
Michael Niedermayer 已提交
2488
    s->codec_id= s->avctx->codec_id= CODEC_ID_MPEG2VIDEO;
M
Michael Niedermayer 已提交
2489
    avctx->sub_id = 2; /* indicates mpeg2 */
M
Michael Niedermayer 已提交
2490 2491 2492 2493
    return 0;
}


2494 2495 2496 2497 2498 2499 2500 2501 2502 2503 2504 2505 2506 2507 2508 2509 2510 2511 2512 2513 2514 2515 2516 2517 2518 2519 2520 2521 2522
static void mpeg_decode_user_data(AVCodecContext *avctx, 
                                  const uint8_t *buf, int buf_size)
{
    const uint8_t *p;
    int len, flags;
    p = buf;
    len = buf_size;

    /* we parse the DTG active format information */
    if (len >= 5 &&
        p[0] == 'D' && p[1] == 'T' && p[2] == 'G' && p[3] == '1') {
        flags = p[4];
        p += 5;
        len -= 5;
        if (flags & 0x80) {
            /* skip event id */
            if (len < 2)
                return;
            p += 2;
            len -= 2;
        }
        if (flags & 0x40) {
            if (len < 1)
                return;
            avctx->dtg_active_format = p[0] & 0x0f;
        }
    }
}

2523 2524 2525 2526 2527 2528 2529 2530 2531 2532 2533 2534 2535 2536 2537 2538 2539 2540 2541 2542 2543 2544 2545 2546 2547 2548 2549 2550 2551 2552 2553 2554 2555 2556 2557 2558
/**
 * finds the end of the current frame in the bitstream.
 * @return the position of the first byte of the next frame, or -1
 */
static int mpeg1_find_frame_end(MpegEncContext *s, uint8_t *buf, int buf_size){
    ParseContext *pc= &s->parse_context;
    int i;
    uint32_t state;
    
    state= pc->state;
    
    i=0;
    if(!pc->frame_start_found){
        for(i=0; i<buf_size; i++){
            state= (state<<8) | buf[i];
            if(state >= SLICE_MIN_START_CODE && state <= SLICE_MAX_START_CODE){
                i++;
                pc->frame_start_found=1;
                break;
            }
        }
    }
    
    if(pc->frame_start_found){
        for(; i<buf_size; i++){
            state= (state<<8) | buf[i];
            if((state&0xFFFFFF00) == 0x100){
                if(state < SLICE_MIN_START_CODE || state > SLICE_MAX_START_CODE){
                    pc->frame_start_found=0;
                    pc->state=-1; 
                    return i-3;
                }
            }
        }
    }        
    pc->state= state;
2559
    return END_NOT_FOUND;
2560 2561
}

F
Fabrice Bellard 已提交
2562 2563 2564
/* handle buffering and image synchronisation */
static int mpeg_decode_frame(AVCodecContext *avctx, 
                             void *data, int *data_size,
2565
                             uint8_t *buf, int buf_size)
F
Fabrice Bellard 已提交
2566 2567
{
    Mpeg1Context *s = avctx->priv_data;
2568 2569
    uint8_t *buf_end, *buf_ptr;
    int ret, start_code, input_size;
M
Michael Niedermayer 已提交
2570
    AVFrame *picture = data;
2571
    MpegEncContext *s2 = &s->mpeg_enc_ctx;
F
Fabrice Bellard 已提交
2572 2573 2574
    dprintf("fill_buffer\n");

    *data_size = 0;
2575

F
Fabrice Bellard 已提交
2576
    /* special case for last picture */
M
Michael Niedermayer 已提交
2577 2578 2579
    if (buf_size == 0 && s2->low_delay==0 && s2->next_picture_ptr) {
        *picture= *(AVFrame*)s2->next_picture_ptr;
        s2->next_picture_ptr= NULL;
M
cleanup  
Michael Niedermayer 已提交
2580

M
Michael Niedermayer 已提交
2581
        *data_size = sizeof(AVFrame);
F
Fabrice Bellard 已提交
2582 2583 2584
        return 0;
    }

2585
    if(s2->flags&CODEC_FLAG_TRUNCATED){
2586
        int next= mpeg1_find_frame_end(s2, buf, buf_size);
2587 2588 2589 2590 2591
        
        if( ff_combine_frame(s2, next, &buf, &buf_size) < 0 )
            return buf_size;
    }    
    
F
Fabrice Bellard 已提交
2592 2593
    buf_ptr = buf;
    buf_end = buf + buf_size;
2594 2595 2596

#if 0    
    if (s->repeat_field % 2 == 1) { 
2597 2598
        s->repeat_field++;
        //fprintf(stderr,"\nRepeating last frame: %d -> %d! pict: %d %d", avctx->frame_number-1, avctx->frame_number,
2599 2600 2601 2602 2603
        //        s2->picture_number, s->repeat_field);
        if (avctx->flags & CODEC_FLAG_REPEAT_FIELD) {
            *data_size = sizeof(AVPicture);
            goto the_end;
        }
2604
    }
2605
#endif
M
Michael Niedermayer 已提交
2606 2607 2608 2609

    if(s->mpeg_enc_ctx_allocated==0 && avctx->codec_tag == ff_get_fourcc("VCR2"))
        vcr2_init_sequence(avctx);

2610
    for(;;) {
F
Fabrice Bellard 已提交
2611
        /* find start next code */
2612
        start_code = find_start_code(&buf_ptr, buf_end);
M
10l  
Michael Niedermayer 已提交
2613
        if (start_code < 0){
M
Michael Niedermayer 已提交
2614 2615
            if(s2->pict_type != B_TYPE || avctx->hurry_up==0){
                if (slice_end(avctx, picture)) {
F
Fabrice Bellard 已提交
2616
                    if(s2->last_picture_ptr || s2->low_delay) //FIXME merge with the stuff in mpeg_decode_slice
M
Michael Niedermayer 已提交
2617 2618
                        *data_size = sizeof(AVPicture);
                }
2619
            }
M
10l  
Michael Niedermayer 已提交
2620
            return FFMAX(0, buf_ptr - buf - s2->parse_context.last_index);
F
Fabrice Bellard 已提交
2621
        }
2622 2623 2624 2625
        
        input_size = buf_end - buf_ptr;

        if(avctx->debug & FF_DEBUG_STARTCODE){
2626
            av_log(avctx, AV_LOG_DEBUG, "%3X at %d left %d\n", start_code, buf_ptr-buf, input_size);
2627
        }
2628

F
Fabrice Bellard 已提交
2629 2630 2631
                /* prepare data for next start code */
                switch(start_code) {
                case SEQ_START_CODE:
2632
                    mpeg1_decode_sequence(avctx, buf_ptr, 
F
Fabrice Bellard 已提交
2633 2634 2635 2636 2637 2638
                                          input_size);
                    break;
                            
                case PICTURE_START_CODE:
                    /* we have a complete image : we try to decompress it */
                    mpeg1_decode_picture(avctx, 
2639
                                         buf_ptr, input_size);
F
Fabrice Bellard 已提交
2640 2641 2642
                    break;
                case EXT_START_CODE:
                    mpeg_decode_extension(avctx,
2643
                                          buf_ptr, input_size);
2644 2645 2646
                    break;
                case USER_START_CODE:
                    mpeg_decode_user_data(avctx, 
2647
                                          buf_ptr, input_size);
F
Fabrice Bellard 已提交
2648
                    break;
M
Michael Niedermayer 已提交
2649 2650 2651
                case GOP_START_CODE:
                    s2->first_field=0;
                    break;
F
Fabrice Bellard 已提交
2652 2653
                default:
                    if (start_code >= SLICE_MIN_START_CODE &&
2654
                        start_code <= SLICE_MAX_START_CODE) {
M
Michael Niedermayer 已提交
2655 2656
                        
                        /* skip b frames if we dont have reference frames */
M
Michael Niedermayer 已提交
2657
                        if(s2->last_picture_ptr==NULL && s2->pict_type==B_TYPE) break;
M
Michael Niedermayer 已提交
2658 2659 2660 2661
                        /* skip b frames if we are in a hurry */
                        if(avctx->hurry_up && s2->pict_type==B_TYPE) break;
                        /* skip everything if we are in a hurry>=5 */
                        if(avctx->hurry_up>=5) break;
2662
                        
2663 2664
                        if (!s->mpeg_enc_ctx_allocated) break;

F
Fabrice Bellard 已提交
2665
                        ret = mpeg_decode_slice(avctx, picture,
2666
                                                start_code, &buf_ptr, input_size);
2667
                        emms_c();
2668

2669 2670
                        if(ret < 0){
                            if(s2->resync_mb_x>=0 && s2->resync_mb_y>=0)
2671
                                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);
2672 2673 2674
                            if(ret==DECODE_SLICE_FATAL_ERROR) return -1;
                        }else{
                            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);
F
Fabrice Bellard 已提交
2675 2676 2677 2678 2679 2680 2681 2682 2683 2684 2685 2686 2687 2688 2689 2690
                        }
                    }
                    break;
                }
    }
}

static int mpeg_decode_end(AVCodecContext *avctx)
{
    Mpeg1Context *s = avctx->priv_data;

    if (s->mpeg_enc_ctx_allocated)
        MPV_common_end(&s->mpeg_enc_ctx);
    return 0;
}

M
Michael Niedermayer 已提交
2691 2692
AVCodec mpeg1video_decoder = {
    "mpeg1video",
F
Fabrice Bellard 已提交
2693 2694 2695 2696 2697 2698 2699
    CODEC_TYPE_VIDEO,
    CODEC_ID_MPEG1VIDEO,
    sizeof(Mpeg1Context),
    mpeg_decode_init,
    NULL,
    mpeg_decode_end,
    mpeg_decode_frame,
2700
    CODEC_CAP_DRAW_HORIZ_BAND | CODEC_CAP_DR1 | CODEC_CAP_TRUNCATED,
M
Michael Niedermayer 已提交
2701 2702 2703 2704 2705 2706 2707 2708 2709 2710 2711 2712 2713
    .flush= ff_mpeg_flush,
};

AVCodec mpeg2video_decoder = {
    "mpeg2video",
    CODEC_TYPE_VIDEO,
    CODEC_ID_MPEG2VIDEO,
    sizeof(Mpeg1Context),
    mpeg_decode_init,
    NULL,
    mpeg_decode_end,
    mpeg_decode_frame,
    CODEC_CAP_DRAW_HORIZ_BAND | CODEC_CAP_DR1 | CODEC_CAP_TRUNCATED,
I
Ivan Kalvachev 已提交
2714 2715 2716 2717 2718 2719 2720 2721 2722 2723 2724 2725 2726 2727
    .flush= ff_mpeg_flush,
};

//legacy decoder
AVCodec mpegvideo_decoder = {
    "mpegvideo",
    CODEC_TYPE_VIDEO,
    CODEC_ID_MPEG2VIDEO,
    sizeof(Mpeg1Context),
    mpeg_decode_init,
    NULL,
    mpeg_decode_end,
    mpeg_decode_frame,
    CODEC_CAP_DRAW_HORIZ_BAND | CODEC_CAP_DR1 | CODEC_CAP_TRUNCATED,
M
Michael Niedermayer 已提交
2728
    .flush= ff_mpeg_flush,
F
Fabrice Bellard 已提交
2729
};
I
Ivan Kalvachev 已提交
2730 2731 2732 2733 2734 2735 2736

#ifdef HAVE_XVMC
static int mpeg_mc_decode_init(AVCodecContext *avctx){
    Mpeg1Context *s;

    if( !(avctx->slice_flags & SLICE_FLAG_CODED_ORDER) )
        return -1;
2737
    if( !(avctx->slice_flags & SLICE_FLAG_ALLOW_FIELD) ){
I
Ivan Kalvachev 已提交
2738
        dprintf("mpeg12.c: XvMC decoder will work better if SLICE_FLAG_ALLOW_FIELD is set\n");
2739
    }
I
Ivan Kalvachev 已提交
2740 2741 2742 2743
    mpeg_decode_init(avctx);
    s = avctx->priv_data;

    avctx->pix_fmt = PIX_FMT_XVMC_MPEG2_IDCT;
2744
    avctx->xvmc_acceleration = 2;//2 - the blocks are packed!
I
Ivan Kalvachev 已提交
2745 2746 2747 2748 2749 2750 2751 2752 2753 2754 2755 2756 2757 2758 2759 2760 2761

    return 0;
}

AVCodec mpeg_xvmc_decoder = {
    "mpegvideo_xvmc",
    CODEC_TYPE_VIDEO,
    CODEC_ID_MPEG2VIDEO_XVMC,
    sizeof(Mpeg1Context),
    mpeg_mc_decode_init,
    NULL,
    mpeg_decode_end,
    mpeg_decode_frame,
    CODEC_CAP_DRAW_HORIZ_BAND | CODEC_CAP_DR1 | CODEC_CAP_TRUNCATED,
};

#endif
2762 2763 2764 2765 2766

/* this is ugly i know, but the alternative is too make 
   hundreds of vars global and prefix them with ff_mpeg1_
   which is far uglier. */
#include "mdec.c"