Statistics
| Branch: | Revision:

ffmpeg / libavcodec / mpeg12.c @ da9c9637

History | View | Annotate | Download (107 KB)

1 de6d9b64 Fabrice Bellard
/*
2 983e3246 Michael Niedermayer
 * MPEG1 codec / MPEG2 decoder
3 ff4ec49e Fabrice Bellard
 * Copyright (c) 2000,2001 Fabrice Bellard.
4 8f2ab833 Michael Niedermayer
 * Copyright (c) 2002-2004 Michael Niedermayer <michaelni@gmx.at> 
5 de6d9b64 Fabrice Bellard
 *
6 ff4ec49e Fabrice Bellard
 * This library is free software; you can redistribute it and/or
7
 * modify it under the terms of the GNU Lesser General Public
8
 * License as published by the Free Software Foundation; either
9
 * version 2 of the License, or (at your option) any later version.
10 de6d9b64 Fabrice Bellard
 *
11 ff4ec49e Fabrice Bellard
 * This library is distributed in the hope that it will be useful,
12 de6d9b64 Fabrice Bellard
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
13 ff4ec49e Fabrice Bellard
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
14
 * Lesser General Public License for more details.
15 de6d9b64 Fabrice Bellard
 *
16 ff4ec49e Fabrice Bellard
 * You should have received a copy of the GNU Lesser General Public
17
 * License along with this library; if not, write to the Free Software
18
 * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
19 de6d9b64 Fabrice Bellard
 */
20 983e3246 Michael Niedermayer
 
21
/**
22
 * @file mpeg12.c
23 029911d1 Michael Niedermayer
 * MPEG1/2 codec
24 983e3246 Michael Niedermayer
 */
25
 
26 a74127c0 Fabrice Bellard
//#define DEBUG
27 de6d9b64 Fabrice Bellard
#include "avcodec.h"
28
#include "dsputil.h"
29
#include "mpegvideo.h"
30
31
#include "mpeg12data.h"
32
33 bb198e19 Michael Niedermayer
//#undef NDEBUG
34
//#include <assert.h>
35
36 d87c0267 Michael Niedermayer
37 de6d9b64 Fabrice Bellard
/* Start codes. */
38
#define SEQ_END_CODE                0x000001b7
39
#define SEQ_START_CODE                0x000001b3
40
#define GOP_START_CODE                0x000001b8
41
#define PICTURE_START_CODE        0x00000100
42
#define SLICE_MIN_START_CODE        0x00000101
43
#define SLICE_MAX_START_CODE        0x000001af
44
#define EXT_START_CODE                0x000001b5
45
#define USER_START_CODE                0x000001b2
46
47 4f68b084 Michael Niedermayer
#define DC_VLC_BITS 9
48
#define MV_VLC_BITS 9
49
#define MBINCR_VLC_BITS 9
50
#define MB_PAT_VLC_BITS 9
51
#define MB_PTYPE_VLC_BITS 6
52
#define MB_BTYPE_VLC_BITS 6
53
#define TEX_VLC_BITS 9
54
55 764ef400 Mike Melanson
#ifdef CONFIG_ENCODERS
56 de6d9b64 Fabrice Bellard
static void mpeg1_encode_block(MpegEncContext *s, 
57
                         DCTELEM *block, 
58
                         int component);
59 63b15e55 Michael Niedermayer
static void mpeg1_encode_motion(MpegEncContext *s, int val, int f_or_b_code);    // RAL: f_code parameter added
60 764ef400 Mike Melanson
#endif //CONFIG_ENCODERS
61 a0201736 Michael Niedermayer
static inline int mpeg1_decode_block_inter(MpegEncContext *s, 
62
                              DCTELEM *block, 
63
                              int n);
64
static inline int mpeg1_decode_block_intra(MpegEncContext *s, 
65 de6d9b64 Fabrice Bellard
                              DCTELEM *block, 
66
                              int n);
67 628b210f Michael Niedermayer
static inline int mpeg1_fast_decode_block_inter(MpegEncContext *s, DCTELEM *block, int n);
68 3729c912 Michael Niedermayer
static inline int mpeg2_decode_block_non_intra(MpegEncContext *s, 
69 de6d9b64 Fabrice Bellard
                                        DCTELEM *block, 
70
                                        int n);
71 3729c912 Michael Niedermayer
static inline int mpeg2_decode_block_intra(MpegEncContext *s, 
72 de6d9b64 Fabrice Bellard
                                    DCTELEM *block, 
73
                                    int n);
74 6fc5b059 Michael Niedermayer
static inline int mpeg2_fast_decode_block_non_intra(MpegEncContext *s, DCTELEM *block, int n);
75 de6d9b64 Fabrice Bellard
static int mpeg_decode_motion(MpegEncContext *s, int fcode, int pred);
76 57489bc5 Ivan Kalvachev
static void exchange_uv(MpegEncContext *s);
77 de6d9b64 Fabrice Bellard
78 2e7b4c84 Ivan Kalvachev
#ifdef HAVE_XVMC
79
extern int XVMC_field_start(MpegEncContext *s, AVCodecContext *avctx);
80
extern int XVMC_field_end(MpegEncContext *s);
81 a579db0c Ivan Kalvachev
extern void XVMC_pack_pblocks(MpegEncContext *s,int cbp);
82 57489bc5 Ivan Kalvachev
extern void XVMC_init_block(MpegEncContext *s);//set s->block
83 2e7b4c84 Ivan Kalvachev
#endif
84
85 530d5740 Ivan Kalvachev
const enum PixelFormat pixfmt_yuv_420[]= {PIX_FMT_YUV420P,-1};
86
const enum PixelFormat pixfmt_yuv_422[]= {PIX_FMT_YUV422P,-1};
87
const enum PixelFormat pixfmt_yuv_444[]= {PIX_FMT_YUV444P,-1};
88
const enum PixelFormat pixfmt_xvmc_mpg2_420[] = {
89
                                           PIX_FMT_XVMC_MPEG2_IDCT,
90
                                           PIX_FMT_XVMC_MPEG2_MC,
91
                                           -1};
92 c442d75c Michael Niedermayer
#ifdef CONFIG_ENCODERS
93 30952237 Michael Niedermayer
static uint8_t (*mv_penalty)[MAX_MV*2+1]= NULL;
94 0c1a9eda Zdenek Kabelac
static uint8_t fcode_tab[MAX_MV*2+1];
95 11ce8834 Michael Niedermayer
96 c442d75c Michael Niedermayer
static uint32_t uni_mpeg1_ac_vlc_bits[64*64*2];
97
static uint8_t  uni_mpeg1_ac_vlc_len [64*64*2];
98 2a250222 Michael Niedermayer
99
/* simple include everything table for dc, first byte is bits number next 3 are code*/
100
static uint32_t mpeg1_lum_dc_uni[512];
101
static uint32_t mpeg1_chr_dc_uni[512];
102
103
static uint8_t mpeg1_index_run[2][64];
104
static int8_t mpeg1_max_level[2][64];
105 764ef400 Mike Melanson
#endif //CONFIG_ENCODERS
106 c442d75c Michael Niedermayer
107 4f68b084 Michael Niedermayer
static void init_2d_vlc_rl(RLTable *rl)
108
{
109 07787186 Dominik Mierzejewski
    int i;
110 4f68b084 Michael Niedermayer
    
111
    init_vlc(&rl->vlc, TEX_VLC_BITS, rl->n + 2, 
112
             &rl->table_vlc[0][1], 4, 2,
113
             &rl->table_vlc[0][0], 4, 2);
114
115
    
116
    rl->rl_vlc[0]= av_malloc(rl->vlc.table_size*sizeof(RL_VLC_ELEM));
117
    for(i=0; i<rl->vlc.table_size; i++){
118
        int code= rl->vlc.table[i][0];
119
        int len = rl->vlc.table[i][1];
120
        int level, run;
121
    
122
        if(len==0){ // illegal code
123
            run= 65;
124
            level= MAX_LEVEL;
125
        }else if(len<0){ //more bits needed
126
            run= 0;
127
            level= code;
128
        }else{
129
            if(code==rl->n){ //esc
130
                run= 65;
131
                level= 0;
132
            }else if(code==rl->n+1){ //eob
133 a0201736 Michael Niedermayer
                run= 0;
134
                level= 127;
135 4f68b084 Michael Niedermayer
            }else{
136
                run=   rl->table_run  [code] + 1;
137
                level= rl->table_level[code];
138
            }
139
        }
140
        rl->rl_vlc[0][i].len= len;
141
        rl->rl_vlc[0][i].level= level;
142
        rl->rl_vlc[0][i].run= run;
143
    }
144
}
145
146 2a250222 Michael Niedermayer
#ifdef CONFIG_ENCODERS
147 c442d75c Michael Niedermayer
static void init_uni_ac_vlc(RLTable *rl, uint32_t *uni_ac_vlc_bits, uint8_t *uni_ac_vlc_len){
148
    int i;
149
150
    for(i=0; i<128; i++){
151
        int level= i-64;
152
        int run;
153
        for(run=0; run<64; run++){
154
            int len, bits, code;
155
            
156
            int alevel= ABS(level);
157
            int sign= (level>>31)&1;
158
159
            if (alevel > rl->max_level[0][run])
160
                code= 111; /*rl->n*/
161
            else
162
                code= rl->index_run[0][run] + alevel - 1;
163
164
            if (code < 111 /* rl->n */) {
165
                    /* store the vlc & sign at once */
166
                len=   mpeg1_vlc[code][1]+1;
167
                bits= (mpeg1_vlc[code][0]<<1) + sign;
168
            } else {
169
                len=  mpeg1_vlc[111/*rl->n*/][1]+6;
170
                bits= mpeg1_vlc[111/*rl->n*/][0]<<6;
171
172
                bits|= run;
173
                if (alevel < 128) {
174
                    bits<<=8; len+=8;
175
                    bits|= level & 0xff;
176
                } else {
177
                    bits<<=16; len+=16;
178
                    bits|= level & 0xff;
179
                    if (level < 0) {
180
                        bits|= 0x8001 + level + 255;
181
                    } else {
182
                        bits|= level & 0xffff;
183
                    }
184
                }
185
            }
186
187
            uni_ac_vlc_bits[UNI_AC_ENC_INDEX(run, i)]= bits;
188
            uni_ac_vlc_len [UNI_AC_ENC_INDEX(run, i)]= len;
189
        }
190
    }
191
}
192 4f68b084 Michael Niedermayer
193 2c492e94 Michael Niedermayer
194
static int find_frame_rate_index(MpegEncContext *s){
195
    int i;
196
    int64_t dmin= INT64_MAX;
197
    int64_t d;
198
199
    for(i=1;i<14;i++) {
200 8d52ec7e Michael Niedermayer
        int64_t n0= 1001LL/frame_rate_tab[i].den*frame_rate_tab[i].num*s->avctx->frame_rate_base;
201
        int64_t n1= 1001LL*s->avctx->frame_rate;
202 2c492e94 Michael Niedermayer
        if(s->avctx->strict_std_compliance >= 0 && i>=9) break;
203 8d52ec7e Michael Niedermayer
204
        d = ABS(n0 - n1);
205 2c492e94 Michael Niedermayer
        if(d < dmin){
206
            dmin=d;
207
            s->frame_rate_index= i;
208
        }
209
    }
210
    if(dmin)
211
        return -1;
212
    else
213
        return 0;
214
}
215
216
static int encode_init(AVCodecContext *avctx)
217
{
218
    MpegEncContext *s = avctx->priv_data;
219
220
    if(MPV_encode_init(avctx) < 0)
221
        return -1;
222
223
    if(find_frame_rate_index(s) < 0){
224
        if(s->strict_std_compliance >=0){
225
            av_log(avctx, AV_LOG_ERROR, "MPEG1/2 doesnt support %d/%d fps\n", avctx->frame_rate, avctx->frame_rate_base);
226
            return -1;
227
        }else{
228
            av_log(avctx, AV_LOG_INFO, "MPEG1/2 doesnt support %d/%d fps, there may be AV sync issues\n", avctx->frame_rate, avctx->frame_rate_base);
229
        }
230
    }
231
    
232
    return 0;
233
}
234
235 de6d9b64 Fabrice Bellard
static void put_header(MpegEncContext *s, int header)
236
{
237
    align_put_bits(&s->pb);
238 80feb2a2 Michael Niedermayer
    put_bits(&s->pb, 16, header>>16);
239
    put_bits(&s->pb, 16, header&0xFFFF);
240 de6d9b64 Fabrice Bellard
}
241
242
/* put sequence header if needed */
243
static void mpeg1_encode_sequence_header(MpegEncContext *s)
244
{
245
        unsigned int vbv_buffer_size;
246 d78647e8 Fabrice Bellard
        unsigned int fps, v;
247 7906085f Falk Hüffner
        int i;
248 0c1a9eda Zdenek Kabelac
        uint64_t time_code;
249 426b8061 Michael Niedermayer
        float best_aspect_error= 1E10;
250 5ff85f1d Michael Niedermayer
        float aspect_ratio= av_q2d(s->avctx->sample_aspect_ratio);
251 98b0858e Michael Niedermayer
        int constraint_parameter_flag;
252 426b8061 Michael Niedermayer
        
253 5ff85f1d Michael Niedermayer
        if(aspect_ratio==0.0) aspect_ratio= 1.0; //pixel aspect 1:1 (VGA)
254 de6d9b64 Fabrice Bellard
        
255 1e491e29 Michael Niedermayer
        if (s->current_picture.key_frame) {
256 8d52ec7e Michael Niedermayer
            AVRational framerate= frame_rate_tab[s->frame_rate_index];
257
258 de6d9b64 Fabrice Bellard
            /* mpeg1 header repeated every gop */
259
            put_header(s, SEQ_START_CODE);
260
 
261
            put_bits(&s->pb, 12, s->width);
262
            put_bits(&s->pb, 12, s->height);
263 426b8061 Michael Niedermayer
            
264
            for(i=1; i<15; i++){
265 5c9e4723 Michael Niedermayer
                float error= aspect_ratio;
266
                if(s->codec_id == CODEC_ID_MPEG1VIDEO || i <=1)
267 f259747b Michael Niedermayer
                    error-= 1.0/mpeg1_aspect[i];
268 5c9e4723 Michael Niedermayer
                else
269
                    error-= av_q2d(mpeg2_aspect[i])*s->height/s->width;
270
             
271 426b8061 Michael Niedermayer
                error= ABS(error);
272
                
273
                if(error < best_aspect_error){
274
                    best_aspect_error= error;
275
                    s->aspect_ratio_info= i;
276
                }
277
            }
278
            
279
            put_bits(&s->pb, 4, s->aspect_ratio_info);
280 de6d9b64 Fabrice Bellard
            put_bits(&s->pb, 4, s->frame_rate_index);
281 12dccd4e Michael Niedermayer
            
282
            if(s->avctx->rc_max_rate){
283
                v = (s->avctx->rc_max_rate + 399) / 400;
284
                if (v > 0x3ffff && s->codec_id == CODEC_ID_MPEG1VIDEO)
285
                    v = 0x3ffff;
286
            }else{
287
                v= 0x3FFFF;
288
            }
289 ecfd40b3 Henry Mason
290
            if(s->avctx->rc_buffer_size)
291
                vbv_buffer_size = s->avctx->rc_buffer_size;
292
            else
293
                /* VBV calculation: Scaled so that a VCD has the proper VBV size of 40 kilobytes */
294 12dccd4e Michael Niedermayer
                vbv_buffer_size = (( 20 * s->bit_rate) / (1151929 / 2)) * 8 * 1024;
295
            vbv_buffer_size= (vbv_buffer_size + 16383) / 16384;
296
297
            put_bits(&s->pb, 18, v & 0x3FFFF);
298
            put_bits(&s->pb, 1, 1); /* marker */
299 98b0858e Michael Niedermayer
            put_bits(&s->pb, 10, vbv_buffer_size & 0x3FF);
300
301
            constraint_parameter_flag= 
302
                s->width <= 768 && s->height <= 576 && 
303
                s->mb_width * s->mb_height <= 396 &&
304 8d52ec7e Michael Niedermayer
                s->mb_width * s->mb_height * framerate.num <= framerate.den*396*25 &&
305
                framerate.num <= framerate.den*30 &&
306 98b0858e Michael Niedermayer
                vbv_buffer_size <= 20 &&
307
                v <= 1856000/400 &&
308
                s->codec_id == CODEC_ID_MPEG1VIDEO;
309
                
310
            put_bits(&s->pb, 1, constraint_parameter_flag);
311 d6eb3c50 Michael Niedermayer
            
312
            ff_write_quant_matrix(&s->pb, s->avctx->intra_matrix);
313
            ff_write_quant_matrix(&s->pb, s->avctx->inter_matrix);
314 de6d9b64 Fabrice Bellard
315 029911d1 Michael Niedermayer
            if(s->codec_id == CODEC_ID_MPEG2VIDEO){
316
                put_header(s, EXT_START_CODE);
317
                put_bits(&s->pb, 4, 1); //seq ext
318
                put_bits(&s->pb, 1, 0); //esc
319 baced9f5 Michael Niedermayer
                
320
                if(s->avctx->profile == FF_PROFILE_UNKNOWN){
321
                    put_bits(&s->pb, 3, 4); //profile
322
                }else{
323
                    put_bits(&s->pb, 3, s->avctx->profile); //profile
324
                }
325
326
                if(s->avctx->level == FF_LEVEL_UNKNOWN){
327
                    put_bits(&s->pb, 4, 8); //level
328
                }else{
329
                    put_bits(&s->pb, 4, s->avctx->level); //level
330
                }
331
332 3ab5b8cc Fabrice Bellard
                put_bits(&s->pb, 1, s->progressive_sequence);
333 029911d1 Michael Niedermayer
                put_bits(&s->pb, 2, 1); //chroma format 4:2:0
334
                put_bits(&s->pb, 2, 0); //horizontal size ext
335
                put_bits(&s->pb, 2, 0); //vertical size ext
336
                put_bits(&s->pb, 12, v>>18); //bitrate ext
337
                put_bits(&s->pb, 1, 1); //marker
338
                put_bits(&s->pb, 8, vbv_buffer_size >>10); //vbv buffer ext
339
                put_bits(&s->pb, 1, s->low_delay);
340
                put_bits(&s->pb, 2, 0); // frame_rate_ext_n
341
                put_bits(&s->pb, 5, 0); // frame_rate_ext_d
342
            }
343
            
344 de6d9b64 Fabrice Bellard
            put_header(s, GOP_START_CODE);
345
            put_bits(&s->pb, 1, 0); /* do drop frame */
346
            /* time code : we must convert from the real frame rate to a
347
               fake mpeg frame rate in case of low frame rate */
348 8d52ec7e Michael Niedermayer
            fps = (framerate.num + framerate.den/2)/ framerate.den;
349 2c492e94 Michael Niedermayer
            time_code = s->current_picture_ptr->coded_picture_number;
350
351
            s->gop_picture_number = time_code;
352 0c1a9eda Zdenek Kabelac
            put_bits(&s->pb, 5, (uint32_t)((time_code / (fps * 3600)) % 24));
353
            put_bits(&s->pb, 6, (uint32_t)((time_code / (fps * 60)) % 60));
354 de6d9b64 Fabrice Bellard
            put_bits(&s->pb, 1, 1);
355 0c1a9eda Zdenek Kabelac
            put_bits(&s->pb, 6, (uint32_t)((time_code / fps) % 60));
356 b1b8f9fc Michael Niedermayer
            put_bits(&s->pb, 6, (uint32_t)((time_code % fps)));
357 303e50e6 Michael Niedermayer
            put_bits(&s->pb, 1, !!(s->flags & CODEC_FLAG_CLOSED_GOP));
358 de6d9b64 Fabrice Bellard
            put_bits(&s->pb, 1, 0); /* broken link */
359
        }
360
}
361
362 9b8709d1 Michael Niedermayer
static inline void encode_mb_skip_run(MpegEncContext *s, int run){
363
    while (run >= 33) {
364
        put_bits(&s->pb, 11, 0x008);
365
        run -= 33;
366
    }
367
    put_bits(&s->pb, mbAddrIncrTable[run][1], 
368
             mbAddrIncrTable[run][0]);
369
}
370 764ef400 Mike Melanson
#endif //CONFIG_ENCODERS
371 de6d9b64 Fabrice Bellard
372 8f8402e4 Michael Niedermayer
static void common_init(MpegEncContext *s)
373 de6d9b64 Fabrice Bellard
{
374 a6b9ffbf Ivan Kalvachev
375 8f8402e4 Michael Niedermayer
    s->y_dc_scale_table=
376 bf266e19 Michael Niedermayer
    s->c_dc_scale_table= mpeg2_dc_scale_table[s->intra_dc_precision];
377 8fe2c60e Michael Niedermayer
378 8f8402e4 Michael Niedermayer
}
379 f004ca1c Michael Niedermayer
380 2a250222 Michael Niedermayer
void ff_mpeg1_clean_buffers(MpegEncContext *s){
381
    s->last_dc[0] = 1 << (7 + s->intra_dc_precision);
382
    s->last_dc[1] = s->last_dc[0];
383
    s->last_dc[2] = s->last_dc[0];
384
    memset(s->last_mv, 0, sizeof(s->last_mv));
385
}
386
387 7604246d Wolfgang Hesseler
#ifdef CONFIG_ENCODERS
388 9b8709d1 Michael Niedermayer
389
void ff_mpeg1_encode_slice_header(MpegEncContext *s){
390
    put_header(s, SLICE_MIN_START_CODE + s->mb_y);
391
    put_bits(&s->pb, 5, s->qscale); /* quantizer scale */
392
    put_bits(&s->pb, 1, 0); /* slice extra information */
393
}
394
395 8f8402e4 Michael Niedermayer
void mpeg1_encode_picture_header(MpegEncContext *s, int picture_number)
396
{
397 de6d9b64 Fabrice Bellard
    mpeg1_encode_sequence_header(s);
398
399
    /* mpeg1 picture header */
400
    put_header(s, PICTURE_START_CODE);
401
    /* temporal reference */
402 63b15e55 Michael Niedermayer
403
    // RAL: s->picture_number instead of s->fake_picture_number
404
    put_bits(&s->pb, 10, (s->picture_number - 
405 de6d9b64 Fabrice Bellard
                          s->gop_picture_number) & 0x3ff); 
406
    put_bits(&s->pb, 3, s->pict_type);
407 d60a8f85 Michael Niedermayer
408 fe455f33 Alex Beregszaszi
    s->vbv_delay_ptr= s->pb.buf + put_bits_count(&s->pb)/8;
409 d60a8f85 Michael Niedermayer
    put_bits(&s->pb, 16, 0xFFFF); /* vbv_delay */
410 de6d9b64 Fabrice Bellard
    
411 63b15e55 Michael Niedermayer
    // RAL: Forward f_code also needed for B frames
412
    if (s->pict_type == P_TYPE || s->pict_type == B_TYPE) {
413 de6d9b64 Fabrice Bellard
        put_bits(&s->pb, 1, 0); /* half pel coordinates */
414 cb231c48 Michael Niedermayer
        if(s->codec_id == CODEC_ID_MPEG1VIDEO)
415
            put_bits(&s->pb, 3, s->f_code); /* forward_f_code */
416
        else
417
            put_bits(&s->pb, 3, 7); /* forward_f_code */
418 de6d9b64 Fabrice Bellard
    }
419
    
420 63b15e55 Michael Niedermayer
    // RAL: Backward f_code necessary for B frames
421
    if (s->pict_type == B_TYPE) {
422
        put_bits(&s->pb, 1, 0); /* half pel coordinates */
423 cb231c48 Michael Niedermayer
        if(s->codec_id == CODEC_ID_MPEG1VIDEO)
424
            put_bits(&s->pb, 3, s->b_code); /* backward_f_code */
425
        else
426
            put_bits(&s->pb, 3, 7); /* backward_f_code */
427 029911d1 Michael Niedermayer
    }
428 63b15e55 Michael Niedermayer
429 de6d9b64 Fabrice Bellard
    put_bits(&s->pb, 1, 0); /* extra bit picture */
430 3ab5b8cc Fabrice Bellard
431
    s->frame_pred_frame_dct = 1;
432 029911d1 Michael Niedermayer
    if(s->codec_id == CODEC_ID_MPEG2VIDEO){
433
        put_header(s, EXT_START_CODE);
434
        put_bits(&s->pb, 4, 8); //pic ext
435 cb231c48 Michael Niedermayer
        if (s->pict_type == P_TYPE || s->pict_type == B_TYPE) {
436
            put_bits(&s->pb, 4, s->f_code);
437
            put_bits(&s->pb, 4, s->f_code);
438
        }else{
439
            put_bits(&s->pb, 8, 255);
440
        }
441
        if (s->pict_type == B_TYPE) {
442
            put_bits(&s->pb, 4, s->b_code);
443
            put_bits(&s->pb, 4, s->b_code);
444
        }else{
445
            put_bits(&s->pb, 8, 255);
446
        }
447 029911d1 Michael Niedermayer
        put_bits(&s->pb, 2, s->intra_dc_precision);
448 9c3d33d6 Michael Niedermayer
        
449
        assert(s->picture_structure == PICT_FRAME);
450
        put_bits(&s->pb, 2, s->picture_structure);
451 7a0f9d7e Fabrice Bellard
        if (s->progressive_sequence) {
452
            put_bits(&s->pb, 1, 0); /* no repeat */
453
        } else {
454
            put_bits(&s->pb, 1, s->current_picture_ptr->top_field_first);
455
        }
456 3ab5b8cc Fabrice Bellard
        /* XXX: optimize the generation of this flag with entropy
457
           measures */
458
        s->frame_pred_frame_dct = s->progressive_sequence;
459
        
460
        put_bits(&s->pb, 1, s->frame_pred_frame_dct);
461 029911d1 Michael Niedermayer
        put_bits(&s->pb, 1, s->concealment_motion_vectors);
462
        put_bits(&s->pb, 1, s->q_scale_type);
463
        put_bits(&s->pb, 1, s->intra_vlc_format);
464
        put_bits(&s->pb, 1, s->alternate_scan);
465
        put_bits(&s->pb, 1, s->repeat_first_field);
466
        put_bits(&s->pb, 1, s->chroma_420_type=1);
467 3ab5b8cc Fabrice Bellard
        s->progressive_frame = s->progressive_sequence;
468
        put_bits(&s->pb, 1, s->progressive_frame);
469 029911d1 Michael Niedermayer
        put_bits(&s->pb, 1, 0); //composite_display_flag
470
    }
471 baaf3f46 Michael Niedermayer
    if(s->flags & CODEC_FLAG_SVCD_SCAN_OFFSET){
472
        int i;
473
474
        put_header(s, USER_START_CODE);
475
        for(i=0; i<sizeof(svcd_scan_offset_placeholder); i++){
476
            put_bits(&s->pb, 8, svcd_scan_offset_placeholder[i]);
477
        }
478
    }
479 029911d1 Michael Niedermayer
    
480 9b8709d1 Michael Niedermayer
    s->mb_y=0;
481
    ff_mpeg1_encode_slice_header(s);
482 de6d9b64 Fabrice Bellard
}
483
484 3ab5b8cc Fabrice Bellard
static inline void put_mb_modes(MpegEncContext *s, int n, int bits, 
485 bb198e19 Michael Niedermayer
                                int has_mv, int field_motion)
486 3ab5b8cc Fabrice Bellard
{
487
    put_bits(&s->pb, n, bits);
488
    if (!s->frame_pred_frame_dct) {
489
        if (has_mv) 
490 bb198e19 Michael Niedermayer
            put_bits(&s->pb, 2, 2 - field_motion); /* motion_type: frame/field */
491 3ab5b8cc Fabrice Bellard
        put_bits(&s->pb, 1, s->interlaced_dct);
492
    }
493
}
494
495 de6d9b64 Fabrice Bellard
void mpeg1_encode_mb(MpegEncContext *s,
496
                     DCTELEM block[6][64],
497
                     int motion_x, int motion_y)
498
{
499 9b8709d1 Michael Niedermayer
    int i, cbp;
500
    const int mb_x = s->mb_x;
501
    const int mb_y = s->mb_y;
502
    const int first_mb= mb_x == s->resync_mb_x && mb_y == s->resync_mb_y;
503 de6d9b64 Fabrice Bellard
504
    /* compute cbp */
505
    cbp = 0;
506
    for(i=0;i<6;i++) {
507
        if (s->block_last_index[i] >= 0)
508
            cbp |= 1 << (5 - i);
509
    }
510 bb198e19 Michael Niedermayer
    
511 86807410 Michael Niedermayer
    if (cbp == 0 && !first_mb && s->mv_type == MV_TYPE_16X16 &&
512
        (mb_x != s->mb_width - 1 || (mb_y != s->mb_height - 1 && s->codec_id == CODEC_ID_MPEG1VIDEO)) && 
513
        ((s->pict_type == P_TYPE && (motion_x | motion_y) == 0) ||
514 63b15e55 Michael Niedermayer
        (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) |
515
        ((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))) {
516 9b8709d1 Michael Niedermayer
        s->mb_skip_run++;
517 d87c0267 Michael Niedermayer
        s->qscale -= s->dquant;
518 49733979 Michael Niedermayer
        s->skip_count++;
519
        s->misc_bits++;
520
        s->last_bits++;
521 bb198e19 Michael Niedermayer
        if(s->pict_type == P_TYPE){
522
            s->last_mv[0][1][0]= s->last_mv[0][0][0]= 
523
            s->last_mv[0][1][1]= s->last_mv[0][0][1]= 0;
524
        }
525 de6d9b64 Fabrice Bellard
    } else {
526 9b8709d1 Michael Niedermayer
        if(first_mb){
527
            assert(s->mb_skip_run == 0);
528
            encode_mb_skip_run(s, s->mb_x);
529
        }else{
530
            encode_mb_skip_run(s, s->mb_skip_run);
531 de6d9b64 Fabrice Bellard
        }
532
        
533
        if (s->pict_type == I_TYPE) {
534 d87c0267 Michael Niedermayer
            if(s->dquant && cbp){
535 bb198e19 Michael Niedermayer
                put_mb_modes(s, 2, 1, 0, 0); /* macroblock_type : macroblock_quant = 1 */
536 d87c0267 Michael Niedermayer
                put_bits(&s->pb, 5, s->qscale);
537
            }else{
538 bb198e19 Michael Niedermayer
                put_mb_modes(s, 1, 1, 0, 0); /* macroblock_type : macroblock_quant = 0 */
539 d87c0267 Michael Niedermayer
                s->qscale -= s->dquant;
540
            }
541 49733979 Michael Niedermayer
            s->misc_bits+= get_bits_diff(s);
542
            s->i_count++;
543 63b15e55 Michael Niedermayer
        } else if (s->mb_intra) {
544
            if(s->dquant && cbp){
545 bb198e19 Michael Niedermayer
                put_mb_modes(s, 6, 0x01, 0, 0);
546 63b15e55 Michael Niedermayer
                put_bits(&s->pb, 5, s->qscale);
547
            }else{
548 bb198e19 Michael Niedermayer
                put_mb_modes(s, 5, 0x03, 0, 0);
549 63b15e55 Michael Niedermayer
                s->qscale -= s->dquant;
550
            }
551
            s->misc_bits+= get_bits_diff(s);
552
            s->i_count++;
553 bb198e19 Michael Niedermayer
            memset(s->last_mv, 0, sizeof(s->last_mv));
554 63b15e55 Michael Niedermayer
        } else if (s->pict_type == P_TYPE) { 
555 bb198e19 Michael Niedermayer
            if(s->mv_type == MV_TYPE_16X16){
556 de6d9b64 Fabrice Bellard
                if (cbp != 0) {
557 bb198e19 Michael Niedermayer
                    if ((motion_x|motion_y) == 0) {
558 d87c0267 Michael Niedermayer
                        if(s->dquant){
559 bb198e19 Michael Niedermayer
                            put_mb_modes(s, 5, 1, 0, 0); /* macroblock_pattern & quant */
560 d87c0267 Michael Niedermayer
                            put_bits(&s->pb, 5, s->qscale);
561
                        }else{
562 bb198e19 Michael Niedermayer
                            put_mb_modes(s, 2, 1, 0, 0); /* macroblock_pattern only */
563 d87c0267 Michael Niedermayer
                        }
564 49733979 Michael Niedermayer
                        s->misc_bits+= get_bits_diff(s);
565 de6d9b64 Fabrice Bellard
                    } else {
566 d87c0267 Michael Niedermayer
                        if(s->dquant){
567 bb198e19 Michael Niedermayer
                            put_mb_modes(s, 5, 2, 1, 0); /* motion + cbp */
568 d87c0267 Michael Niedermayer
                            put_bits(&s->pb, 5, s->qscale);
569
                        }else{
570 bb198e19 Michael Niedermayer
                            put_mb_modes(s, 1, 1, 1, 0); /* motion + cbp */
571 d87c0267 Michael Niedermayer
                        }
572 49733979 Michael Niedermayer
                        s->misc_bits+= get_bits_diff(s);
573 63b15e55 Michael Niedermayer
                        mpeg1_encode_motion(s, motion_x - s->last_mv[0][0][0], s->f_code);    // RAL: f_code parameter added
574
                        mpeg1_encode_motion(s, motion_y - s->last_mv[0][0][1], s->f_code);    // RAL: f_code parameter added
575 49733979 Michael Niedermayer
                        s->mv_bits+= get_bits_diff(s);
576 de6d9b64 Fabrice Bellard
                    }
577
                } else {
578
                    put_bits(&s->pb, 3, 1); /* motion only */
579 3ab5b8cc Fabrice Bellard
                    if (!s->frame_pred_frame_dct)
580
                        put_bits(&s->pb, 2, 2); /* motion_type: frame */
581 bb198e19 Michael Niedermayer
                    s->misc_bits+= get_bits_diff(s);
582 63b15e55 Michael Niedermayer
                    mpeg1_encode_motion(s, motion_x - s->last_mv[0][0][0], s->f_code);    // RAL: f_code parameter added
583
                    mpeg1_encode_motion(s, motion_y - s->last_mv[0][0][1], s->f_code);    // RAL: f_code parameter added
584 d87c0267 Michael Niedermayer
                    s->qscale -= s->dquant;
585 49733979 Michael Niedermayer
                    s->mv_bits+= get_bits_diff(s);
586 de6d9b64 Fabrice Bellard
                }
587 bb198e19 Michael Niedermayer
                s->last_mv[0][1][0]= s->last_mv[0][0][0]= motion_x;
588
                s->last_mv[0][1][1]= s->last_mv[0][0][1]= motion_y;
589
            }else{
590
                assert(!s->frame_pred_frame_dct && s->mv_type == MV_TYPE_FIELD);
591
592
                if (cbp) {
593
                    if(s->dquant){
594
                        put_mb_modes(s, 5, 2, 1, 1); /* motion + cbp */
595
                        put_bits(&s->pb, 5, s->qscale);
596
                    }else{
597
                        put_mb_modes(s, 1, 1, 1, 1); /* motion + cbp */
598 63b15e55 Michael Niedermayer
                    }
599 bb198e19 Michael Niedermayer
                } else {
600
                    put_bits(&s->pb, 3, 1); /* motion only */
601
                    put_bits(&s->pb, 2, 1); /* motion_type: field */
602 63b15e55 Michael Niedermayer
                    s->qscale -= s->dquant;
603 bb198e19 Michael Niedermayer
                }
604
                s->misc_bits+= get_bits_diff(s);
605
                for(i=0; i<2; i++){
606
                    put_bits(&s->pb, 1, s->field_select[0][i]);
607
                    mpeg1_encode_motion(s, s->mv[0][i][0] -  s->last_mv[0][i][0]    , s->f_code);
608
                    mpeg1_encode_motion(s, s->mv[0][i][1] - (s->last_mv[0][i][1]>>1), s->f_code);
609
                    s->last_mv[0][i][0]=   s->mv[0][i][0];
610
                    s->last_mv[0][i][1]= 2*s->mv[0][i][1];
611
                }
612
                s->mv_bits+= get_bits_diff(s);
613
            }
614
            if(cbp)
615 461cd5bf Ivan Kalvachev
                put_bits(&s->pb, mbPatTable[cbp][1], mbPatTable[cbp][0]);
616 bb198e19 Michael Niedermayer
            s->f_count++;
617
        } else{  
618
            static const int mb_type_len[4]={0,3,4,2}; //bak,for,bi
619
620
            if(s->mv_type == MV_TYPE_16X16){
621
                if (cbp){    // With coded bloc pattern
622
                    if (s->dquant) {
623
                        if(s->mv_dir == MV_DIR_FORWARD)
624
                            put_mb_modes(s, 6, 3, 1, 0);
625
                        else
626
                            put_mb_modes(s, mb_type_len[s->mv_dir]+3, 2, 1, 0);
627
                        put_bits(&s->pb, 5, s->qscale);
628
                    } else {
629
                        put_mb_modes(s, mb_type_len[s->mv_dir], 3, 1, 0);
630
                    }
631
                }else{    // No coded bloc pattern
632
                    put_bits(&s->pb, mb_type_len[s->mv_dir], 2);
633
                    if (!s->frame_pred_frame_dct)
634
                        put_bits(&s->pb, 2, 2); /* motion_type: frame */
635
                    s->qscale -= s->dquant;
636
                }
637
                s->misc_bits += get_bits_diff(s);
638
                if (s->mv_dir&MV_DIR_FORWARD){
639
                    mpeg1_encode_motion(s, s->mv[0][0][0] - s->last_mv[0][0][0], s->f_code); 
640
                    mpeg1_encode_motion(s, s->mv[0][0][1] - s->last_mv[0][0][1], s->f_code); 
641
                    s->last_mv[0][0][0]=s->last_mv[0][1][0]= s->mv[0][0][0];
642
                    s->last_mv[0][0][1]=s->last_mv[0][1][1]= s->mv[0][0][1];
643
                    s->f_count++;
644
                }
645
                if (s->mv_dir&MV_DIR_BACKWARD){
646
                    mpeg1_encode_motion(s, s->mv[1][0][0] - s->last_mv[1][0][0], s->b_code); 
647
                    mpeg1_encode_motion(s, s->mv[1][0][1] - s->last_mv[1][0][1], s->b_code); 
648
                    s->last_mv[1][0][0]=s->last_mv[1][1][0]= s->mv[1][0][0];
649
                    s->last_mv[1][0][1]=s->last_mv[1][1][1]= s->mv[1][0][1];
650
                    s->b_count++;
651
                }
652
            }else{
653
                assert(s->mv_type == MV_TYPE_FIELD);
654
                assert(!s->frame_pred_frame_dct);
655
                if (cbp){    // With coded bloc pattern
656
                    if (s->dquant) {
657
                        if(s->mv_dir == MV_DIR_FORWARD)
658
                            put_mb_modes(s, 6, 3, 1, 1);
659
                        else
660
                            put_mb_modes(s, mb_type_len[s->mv_dir]+3, 2, 1, 1);
661
                        put_bits(&s->pb, 5, s->qscale);
662
                    } else {
663
                        put_mb_modes(s, mb_type_len[s->mv_dir], 3, 1, 1);
664 63b15e55 Michael Niedermayer
                    }
665 bb198e19 Michael Niedermayer
                }else{    // No coded bloc pattern
666
                    put_bits(&s->pb, mb_type_len[s->mv_dir], 2);
667
                    put_bits(&s->pb, 2, 1); /* motion_type: field */
668
                    s->qscale -= s->dquant;
669
                }
670
                s->misc_bits += get_bits_diff(s);
671
                if (s->mv_dir&MV_DIR_FORWARD){
672
                    for(i=0; i<2; i++){
673
                        put_bits(&s->pb, 1, s->field_select[0][i]);
674
                        mpeg1_encode_motion(s, s->mv[0][i][0] -  s->last_mv[0][i][0]    , s->f_code);
675
                        mpeg1_encode_motion(s, s->mv[0][i][1] - (s->last_mv[0][i][1]>>1), s->f_code);
676
                        s->last_mv[0][i][0]=   s->mv[0][i][0];
677
                        s->last_mv[0][i][1]= 2*s->mv[0][i][1];
678
                    }
679
                    s->f_count++;
680
                }
681
                if (s->mv_dir&MV_DIR_BACKWARD){
682
                    for(i=0; i<2; i++){
683
                        put_bits(&s->pb, 1, s->field_select[1][i]);
684
                        mpeg1_encode_motion(s, s->mv[1][i][0] -  s->last_mv[1][i][0]    , s->b_code);
685
                        mpeg1_encode_motion(s, s->mv[1][i][1] - (s->last_mv[1][i][1]>>1), s->b_code);
686
                        s->last_mv[1][i][0]=   s->mv[1][i][0];
687
                        s->last_mv[1][i][1]= 2*s->mv[1][i][1];
688
                    }
689
                    s->b_count++;
690
                }
691 de6d9b64 Fabrice Bellard
            }
692 bb198e19 Michael Niedermayer
            s->mv_bits += get_bits_diff(s);
693
            if(cbp)
694 461cd5bf Ivan Kalvachev
                put_bits(&s->pb, mbPatTable[cbp][1], mbPatTable[cbp][0]);
695 bb198e19 Michael Niedermayer
        }
696 de6d9b64 Fabrice Bellard
        for(i=0;i<6;i++) {
697
            if (cbp & (1 << (5 - i))) {
698
                mpeg1_encode_block(s, block[i], i);
699
            }
700
        }
701 9b8709d1 Michael Niedermayer
        s->mb_skip_run = 0;
702 49733979 Michael Niedermayer
        if(s->mb_intra)
703
            s->i_tex_bits+= get_bits_diff(s);
704
        else
705
            s->p_tex_bits+= get_bits_diff(s);
706 de6d9b64 Fabrice Bellard
    }
707
}
708
709 63b15e55 Michael Niedermayer
// RAL: Parameter added: f_or_b_code
710
static void mpeg1_encode_motion(MpegEncContext *s, int val, int f_or_b_code)
711 de6d9b64 Fabrice Bellard
{
712 0ff93477 Michael Niedermayer
    int code, bit_size, l, bits, range, sign;
713 de6d9b64 Fabrice Bellard
714
    if (val == 0) {
715
        /* zero vector */
716
        code = 0;
717
        put_bits(&s->pb,
718
                 mbMotionVectorTable[0][1], 
719
                 mbMotionVectorTable[0][0]); 
720
    } else {
721 63b15e55 Michael Niedermayer
        bit_size = f_or_b_code - 1;
722 de6d9b64 Fabrice Bellard
        range = 1 << bit_size;
723
        /* modulo encoding */
724 0ff93477 Michael Niedermayer
        l= INT_BIT - 5 - bit_size;
725
        val= (val<<l)>>l;
726 de6d9b64 Fabrice Bellard
727
        if (val >= 0) {
728
            val--;
729
            code = (val >> bit_size) + 1;
730
            bits = val & (range - 1);
731
            sign = 0;
732
        } else {
733
            val = -val;
734
            val--;
735
            code = (val >> bit_size) + 1;
736
            bits = val & (range - 1);
737
            sign = 1;
738
        }
739 63b15e55 Michael Niedermayer
740
        assert(code > 0 && code <= 16);
741
742 de6d9b64 Fabrice Bellard
        put_bits(&s->pb,
743
                 mbMotionVectorTable[code][1], 
744
                 mbMotionVectorTable[code][0]); 
745 63b15e55 Michael Niedermayer
746 de6d9b64 Fabrice Bellard
        put_bits(&s->pb, 1, sign);
747
        if (bit_size > 0) {
748
            put_bits(&s->pb, bit_size, bits);
749
        }
750
    }
751
}
752
753 8f8402e4 Michael Niedermayer
void ff_mpeg1_encode_init(MpegEncContext *s)
754 11ce8834 Michael Niedermayer
{
755
    static int done=0;
756 8f8402e4 Michael Niedermayer
757
    common_init(s);
758
759 11ce8834 Michael Niedermayer
    if(!done){
760
        int f_code;
761
        int mv;
762 8f8402e4 Michael Niedermayer
        int i;
763 11ce8834 Michael Niedermayer
764
        done=1;
765 8f8402e4 Michael Niedermayer
        init_rl(&rl_mpeg1);
766 c442d75c Michael Niedermayer
767 8f8402e4 Michael Niedermayer
        for(i=0; i<64; i++)
768
        {
769
                mpeg1_max_level[0][i]= rl_mpeg1.max_level[0][i];
770
                mpeg1_index_run[0][i]= rl_mpeg1.index_run[0][i];
771
        }
772 c442d75c Michael Niedermayer
        
773
        init_uni_ac_vlc(&rl_mpeg1, uni_mpeg1_ac_vlc_bits, uni_mpeg1_ac_vlc_len);
774 8f8402e4 Michael Niedermayer
775
        /* build unified dc encoding tables */
776
        for(i=-255; i<256; i++)
777
        {
778
                int adiff, index;
779
                int bits, code;
780
                int diff=i;
781
782
                adiff = ABS(diff);
783
                if(diff<0) diff--;
784 8aedc51d Michael Niedermayer
                index = av_log2(2*adiff);
785 8f8402e4 Michael Niedermayer
786
                bits= vlc_dc_lum_bits[index] + index;
787
                code= (vlc_dc_lum_code[index]<<index) + (diff & ((1 << index) - 1));
788
                mpeg1_lum_dc_uni[i+255]= bits + (code<<8);
789
                
790
                bits= vlc_dc_chroma_bits[index] + index;
791
                code= (vlc_dc_chroma_code[index]<<index) + (diff & ((1 << index) - 1));
792
                mpeg1_chr_dc_uni[i+255]= bits + (code<<8);
793
        }
794
795 30952237 Michael Niedermayer
        mv_penalty= av_mallocz( sizeof(uint8_t)*(MAX_FCODE+1)*(2*MAX_MV+1) );
796
797 11ce8834 Michael Niedermayer
        for(f_code=1; f_code<=MAX_FCODE; f_code++){
798
            for(mv=-MAX_MV; mv<=MAX_MV; mv++){
799
                int len;
800
801
                if(mv==0) len= mbMotionVectorTable[0][1];
802
                else{
803
                    int val, bit_size, range, code;
804
805 7813d6a3 Michael Niedermayer
                    bit_size = f_code - 1;
806 11ce8834 Michael Niedermayer
                    range = 1 << bit_size;
807
808
                    val=mv;
809
                    if (val < 0) 
810
                        val = -val;
811
                    val--;
812
                    code = (val >> bit_size) + 1;
813
                    if(code<17){
814
                        len= mbMotionVectorTable[code][1] + 1 + bit_size;
815
                    }else{
816
                        len= mbMotionVectorTable[16][1] + 2 + bit_size;
817
                    }
818
                }
819
820
                mv_penalty[f_code][mv+MAX_MV]= len;
821
            }
822
        }
823
        
824
825
        for(f_code=MAX_FCODE; f_code>0; f_code--){
826
            for(mv=-(8<<f_code); mv<(8<<f_code); mv++){
827
                fcode_tab[mv+MAX_MV]= f_code;
828
            }
829
        }
830
    }
831 1457ab52 Michael Niedermayer
    s->me.mv_penalty= mv_penalty;
832 11ce8834 Michael Niedermayer
    s->fcode_tab= fcode_tab;
833 029911d1 Michael Niedermayer
    if(s->codec_id == CODEC_ID_MPEG1VIDEO){
834
        s->min_qcoeff=-255;
835
        s->max_qcoeff= 255;
836
    }else{
837
        s->min_qcoeff=-2047;
838
        s->max_qcoeff= 2047;
839
    }
840 c442d75c Michael Niedermayer
    s->intra_ac_vlc_length=
841 e5021fff Michael Niedermayer
    s->inter_ac_vlc_length=
842
    s->intra_ac_vlc_last_length=
843
    s->inter_ac_vlc_last_length= uni_mpeg1_ac_vlc_len;
844 11ce8834 Michael Niedermayer
}
845 8f8402e4 Michael Niedermayer
846 de6d9b64 Fabrice Bellard
static inline void encode_dc(MpegEncContext *s, int diff, int component)
847
{
848 bf266e19 Michael Niedermayer
  if(((unsigned) (diff+255)) >= 511){
849
        int index;
850
851
        if(diff<0){
852
            index= av_log2_16bit(-2*diff);
853
            diff--;
854
        }else{
855
            index= av_log2_16bit(2*diff);
856
        }
857
        if (component == 0) {
858
            put_bits(
859
                &s->pb, 
860
                vlc_dc_lum_bits[index] + index,
861
                (vlc_dc_lum_code[index]<<index) + (diff & ((1 << index) - 1)));
862
        }else{
863
            put_bits(
864
                &s->pb, 
865
                vlc_dc_chroma_bits[index] + index,
866
                (vlc_dc_chroma_code[index]<<index) + (diff & ((1 << index) - 1)));
867
        }
868
  }else{
869 de6d9b64 Fabrice Bellard
    if (component == 0) {
870 f004ca1c Michael Niedermayer
        put_bits(
871
            &s->pb, 
872
            mpeg1_lum_dc_uni[diff+255]&0xFF,
873
            mpeg1_lum_dc_uni[diff+255]>>8);
874 de6d9b64 Fabrice Bellard
    } else {
875 f004ca1c Michael Niedermayer
        put_bits(
876
            &s->pb, 
877
            mpeg1_chr_dc_uni[diff+255]&0xFF,
878
            mpeg1_chr_dc_uni[diff+255]>>8);
879 de6d9b64 Fabrice Bellard
    }
880 bf266e19 Michael Niedermayer
  }
881 de6d9b64 Fabrice Bellard
}
882
883
static void mpeg1_encode_block(MpegEncContext *s, 
884
                               DCTELEM *block, 
885
                               int n)
886
{
887
    int alevel, level, last_non_zero, dc, diff, i, j, run, last_index, sign;
888
    int code, component;
889 80feb2a2 Michael Niedermayer
//    RLTable *rl = &rl_mpeg1;
890 de6d9b64 Fabrice Bellard
891
    last_index = s->block_last_index[n];
892
893
    /* DC coef */
894
    if (s->mb_intra) {
895
        component = (n <= 3 ? 0 : n - 4 + 1);
896
        dc = block[0]; /* overflow is impossible */
897
        diff = dc - s->last_dc[component];
898
        encode_dc(s, diff, component);
899
        s->last_dc[component] = dc;
900
        i = 1;
901 029911d1 Michael Niedermayer
/*
902
        if (s->intra_vlc_format)
903
            rl = &rl_mpeg2;
904
        else
905
            rl = &rl_mpeg1;
906
*/
907 de6d9b64 Fabrice Bellard
    } else {
908
        /* encode the first coefficient : needs to be done here because
909
           it is handled slightly differently */
910
        level = block[0];
911
        if (abs(level) == 1) {
912 0c1a9eda Zdenek Kabelac
                code = ((uint32_t)level >> 31); /* the sign bit */
913 de6d9b64 Fabrice Bellard
                put_bits(&s->pb, 2, code | 0x02);
914
                i = 1;
915
        } else {
916
            i = 0;
917
            last_non_zero = -1;
918
            goto next_coef;
919
        }
920
    }
921
922
    /* now quantify & encode AC coefs */
923
    last_non_zero = i - 1;
924 80feb2a2 Michael Niedermayer
925 de6d9b64 Fabrice Bellard
    for(;i<=last_index;i++) {
926 2ad1516a Michael Niedermayer
        j = s->intra_scantable.permutated[i];
927 de6d9b64 Fabrice Bellard
        level = block[j];
928
    next_coef:
929
#if 0
930
        if (level != 0)
931
            dprintf("level[%d]=%d\n", i, level);
932
#endif            
933
        /* encode using VLC */
934
        if (level != 0) {
935
            run = i - last_non_zero - 1;
936 2ad1516a Michael Niedermayer
            
937
            alevel= level;
938
            MASK_ABS(sign, alevel)
939
            sign&=1;
940
941 80feb2a2 Michael Niedermayer
//            code = get_rl_index(rl, 0, run, alevel);
942 c442d75c Michael Niedermayer
            if (alevel <= mpeg1_max_level[0][run]){
943 80feb2a2 Michael Niedermayer
                code= mpeg1_index_run[0][run] + alevel - 1;
944
                    /* store the vlc & sign at once */
945
                put_bits(&s->pb, mpeg1_vlc[code][1]+1, (mpeg1_vlc[code][0]<<1) + sign);
946 de6d9b64 Fabrice Bellard
            } else {
947 80feb2a2 Michael Niedermayer
                /* escape seems to be pretty rare <5% so i dont optimize it */
948
                put_bits(&s->pb, mpeg1_vlc[111/*rl->n*/][1], mpeg1_vlc[111/*rl->n*/][0]);
949 de6d9b64 Fabrice Bellard
                /* escape: only clip in this case */
950
                put_bits(&s->pb, 6, run);
951 029911d1 Michael Niedermayer
                if(s->codec_id == CODEC_ID_MPEG1VIDEO){
952
                    if (alevel < 128) {
953
                        put_bits(&s->pb, 8, level & 0xff);
954 de6d9b64 Fabrice Bellard
                    } else {
955 029911d1 Michael Niedermayer
                        if (level < 0) {
956
                            put_bits(&s->pb, 16, 0x8001 + level + 255);
957
                        } else {
958
                            put_bits(&s->pb, 16, level & 0xffff);
959
                        }
960 de6d9b64 Fabrice Bellard
                    }
961 029911d1 Michael Niedermayer
                }else{
962
                    put_bits(&s->pb, 12, level & 0xfff);
963 de6d9b64 Fabrice Bellard
                }
964
            }
965
            last_non_zero = i;
966
        }
967
    }
968
    /* end of block */
969
    put_bits(&s->pb, 2, 0x2);
970
}
971 7604246d Wolfgang Hesseler
#endif //CONFIG_ENCODERS
972 de6d9b64 Fabrice Bellard
973
/******************************************/
974
/* decoding */
975
976
static VLC dc_lum_vlc;
977
static VLC dc_chroma_vlc;
978
static VLC mv_vlc;
979
static VLC mbincr_vlc;
980
static VLC mb_ptype_vlc;
981
static VLC mb_btype_vlc;
982
static VLC mb_pat_vlc;
983
984 7906085f Falk Hüffner
static void init_vlcs(void)
985 de6d9b64 Fabrice Bellard
{
986
    static int done = 0;
987
988
    if (!done) {
989 915bbac6 Fabrice Bellard
        done = 1;
990 de6d9b64 Fabrice Bellard
991 9ac7ecd6 Michael Niedermayer
        init_vlc(&dc_lum_vlc, DC_VLC_BITS, 12, 
992 de6d9b64 Fabrice Bellard
                 vlc_dc_lum_bits, 1, 1,
993
                 vlc_dc_lum_code, 2, 2);
994 9ac7ecd6 Michael Niedermayer
        init_vlc(&dc_chroma_vlc,  DC_VLC_BITS, 12, 
995 de6d9b64 Fabrice Bellard
                 vlc_dc_chroma_bits, 1, 1,
996
                 vlc_dc_chroma_code, 2, 2);
997 8ed2ddb2 Michael Niedermayer
        init_vlc(&mv_vlc, MV_VLC_BITS, 17, 
998 de6d9b64 Fabrice Bellard
                 &mbMotionVectorTable[0][1], 2, 1,
999
                 &mbMotionVectorTable[0][0], 2, 1);
1000 9c00c3af Michael Niedermayer
        init_vlc(&mbincr_vlc, MBINCR_VLC_BITS, 36, 
1001 de6d9b64 Fabrice Bellard
                 &mbAddrIncrTable[0][1], 2, 1,
1002
                 &mbAddrIncrTable[0][0], 2, 1);
1003 461cd5bf Ivan Kalvachev
        init_vlc(&mb_pat_vlc, MB_PAT_VLC_BITS, 64,
1004 de6d9b64 Fabrice Bellard
                 &mbPatTable[0][1], 2, 1,
1005
                 &mbPatTable[0][0], 2, 1);
1006
        
1007 7bc9090a Michael Niedermayer
        init_vlc(&mb_ptype_vlc, MB_PTYPE_VLC_BITS, 7, 
1008 de6d9b64 Fabrice Bellard
                 &table_mb_ptype[0][1], 2, 1,
1009
                 &table_mb_ptype[0][0], 2, 1);
1010 7bc9090a Michael Niedermayer
        init_vlc(&mb_btype_vlc, MB_BTYPE_VLC_BITS, 11, 
1011 de6d9b64 Fabrice Bellard
                 &table_mb_btype[0][1], 2, 1,
1012
                 &table_mb_btype[0][0], 2, 1);
1013
        init_rl(&rl_mpeg1);
1014
        init_rl(&rl_mpeg2);
1015 4f68b084 Michael Niedermayer
1016
        init_2d_vlc_rl(&rl_mpeg1);
1017
        init_2d_vlc_rl(&rl_mpeg2);
1018 de6d9b64 Fabrice Bellard
    }
1019
}
1020
1021
static inline int get_dmv(MpegEncContext *s)
1022
{
1023 612476ef Arpi
    if(get_bits1(&s->gb)) 
1024
        return 1 - (get_bits1(&s->gb) << 1);
1025 de6d9b64 Fabrice Bellard
    else
1026
        return 0;
1027
}
1028
1029 0ee50938 Fabrice Bellard
static inline int get_qscale(MpegEncContext *s)
1030
{
1031 7062fad6 BERO
    int qscale = get_bits(&s->gb, 5);
1032 51929fd3 Michael Niedermayer
    if (s->q_scale_type) {
1033
        return non_linear_qscale[qscale];
1034
    } else {
1035
        return qscale << 1;
1036 0ee50938 Fabrice Bellard
    }
1037
}
1038
1039 de6d9b64 Fabrice Bellard
/* motion type (for mpeg2) */
1040
#define MT_FIELD 1
1041
#define MT_FRAME 2
1042
#define MT_16X8  2
1043
#define MT_DMV   3
1044
1045
static int mpeg_decode_mb(MpegEncContext *s,
1046 5e5c247a Ivan Kalvachev
                          DCTELEM block[12][64])
1047 de6d9b64 Fabrice Bellard
{
1048 d2975f8d Michael Niedermayer
    int i, j, k, cbp, val, mb_type, motion_type;
1049 715731a3 Ivan Kalvachev
    const int mb_block_count = 4 + (1<< s->chroma_format)
1050
1051 de6d9b64 Fabrice Bellard
    dprintf("decode_mb: x=%d y=%d\n", s->mb_x, s->mb_y);
1052
1053 f943e138 Michael Niedermayer
    assert(s->mb_skiped==0);
1054
1055 9b8709d1 Michael Niedermayer
    if (s->mb_skip_run-- != 0) {
1056 7bc9090a Michael Niedermayer
        if(s->pict_type == I_TYPE){
1057 9b879566 Michel Bardiaux
            av_log(s->avctx, AV_LOG_ERROR, "skiped MB in I frame at %d %d\n", s->mb_x, s->mb_y);
1058 7bc9090a Michael Niedermayer
            return -1;
1059
        }
1060
    
1061 de6d9b64 Fabrice Bellard
        /* skip mb */
1062
        s->mb_intra = 0;
1063 5e5c247a Ivan Kalvachev
        for(i=0;i<12;i++)
1064 de6d9b64 Fabrice Bellard
            s->block_last_index[i] = -1;
1065 4c9544d0 Michael Niedermayer
        if(s->picture_structure == PICT_FRAME)
1066
            s->mv_type = MV_TYPE_16X16;
1067
        else
1068
            s->mv_type = MV_TYPE_FIELD;
1069 de6d9b64 Fabrice Bellard
        if (s->pict_type == P_TYPE) {
1070
            /* if P type, zero motion vector is implied */
1071
            s->mv_dir = MV_DIR_FORWARD;
1072
            s->mv[0][0][0] = s->mv[0][0][1] = 0;
1073
            s->last_mv[0][0][0] = s->last_mv[0][0][1] = 0;
1074 25ed7f92 Fabrice Bellard
            s->last_mv[0][1][0] = s->last_mv[0][1][1] = 0;
1075 4c9544d0 Michael Niedermayer
            s->field_select[0][0]= s->picture_structure - 1;
1076 f943e138 Michael Niedermayer
            s->mb_skiped = 1;
1077 7bc9090a Michael Niedermayer
            s->current_picture.mb_type[ s->mb_x + s->mb_y*s->mb_stride ]= MB_TYPE_SKIP | MB_TYPE_L0 | MB_TYPE_16x16;
1078 de6d9b64 Fabrice Bellard
        } else {
1079
            /* if B type, reuse previous vectors and directions */
1080
            s->mv[0][0][0] = s->last_mv[0][0][0];
1081
            s->mv[0][0][1] = s->last_mv[0][0][1];
1082
            s->mv[1][0][0] = s->last_mv[1][0][0];
1083
            s->mv[1][0][1] = s->last_mv[1][0][1];
1084 f943e138 Michael Niedermayer
1085 7bc9090a Michael Niedermayer
            s->current_picture.mb_type[ s->mb_x + s->mb_y*s->mb_stride ]= 
1086
                s->current_picture.mb_type[ s->mb_x + s->mb_y*s->mb_stride - 1] | MB_TYPE_SKIP;
1087
//            assert(s->current_picture.mb_type[ s->mb_x + s->mb_y*s->mb_stride - 1]&(MB_TYPE_16x16|MB_TYPE_16x8));
1088
1089 f943e138 Michael Niedermayer
            if((s->mv[0][0][0]|s->mv[0][0][1]|s->mv[1][0][0]|s->mv[1][0][1])==0) 
1090
                s->mb_skiped = 1;
1091 de6d9b64 Fabrice Bellard
        }
1092 59b571c1 Michael Niedermayer
1093 de6d9b64 Fabrice Bellard
        return 0;
1094
    }
1095
1096
    switch(s->pict_type) {
1097
    default:
1098
    case I_TYPE:
1099 612476ef Arpi
        if (get_bits1(&s->gb) == 0) {
1100 e94bc100 Michael Niedermayer
            if (get_bits1(&s->gb) == 0){
1101 9b879566 Michel Bardiaux
                av_log(s->avctx, AV_LOG_ERROR, "invalid mb type in I Frame at %d %d\n", s->mb_x, s->mb_y);
1102 de6d9b64 Fabrice Bellard
                return -1;
1103 e94bc100 Michael Niedermayer
            }
1104 7bc9090a Michael Niedermayer
            mb_type = MB_TYPE_QUANT | MB_TYPE_INTRA;
1105 de6d9b64 Fabrice Bellard
        } else {
1106 7bc9090a Michael Niedermayer
            mb_type = MB_TYPE_INTRA;
1107 de6d9b64 Fabrice Bellard
        }
1108
        break;
1109
    case P_TYPE:
1110 8ed2ddb2 Michael Niedermayer
        mb_type = get_vlc2(&s->gb, mb_ptype_vlc.table, MB_PTYPE_VLC_BITS, 1);
1111 9ac7ecd6 Michael Niedermayer
        if (mb_type < 0){
1112 9b879566 Michel Bardiaux
            av_log(s->avctx, AV_LOG_ERROR, "invalid mb type in P Frame at %d %d\n", s->mb_x, s->mb_y);
1113 de6d9b64 Fabrice Bellard
            return -1;
1114 9ac7ecd6 Michael Niedermayer
        }
1115 7bc9090a Michael Niedermayer
        mb_type = ptype2mb_type[ mb_type ];
1116 de6d9b64 Fabrice Bellard
        break;
1117
    case B_TYPE:
1118 8ed2ddb2 Michael Niedermayer
        mb_type = get_vlc2(&s->gb, mb_btype_vlc.table, MB_BTYPE_VLC_BITS, 1);
1119 9ac7ecd6 Michael Niedermayer
        if (mb_type < 0){
1120 9b879566 Michel Bardiaux
            av_log(s->avctx, AV_LOG_ERROR, "invalid mb type in B Frame at %d %d\n", s->mb_x, s->mb_y);
1121 de6d9b64 Fabrice Bellard
            return -1;
1122 9ac7ecd6 Michael Niedermayer
        }
1123 7bc9090a Michael Niedermayer
        mb_type = btype2mb_type[ mb_type ];
1124 de6d9b64 Fabrice Bellard
        break;
1125
    }
1126
    dprintf("mb_type=%x\n", mb_type);
1127 7bc9090a Michael Niedermayer
//    motion_type = 0; /* avoid warning */
1128
    if (IS_INTRA(mb_type)) {
1129
        /* compute dct type */
1130
        if (s->picture_structure == PICT_FRAME && //FIXME add a interlaced_dct coded var?
1131
            !s->frame_pred_frame_dct) {
1132
            s->interlaced_dct = get_bits1(&s->gb);
1133
        }
1134 de6d9b64 Fabrice Bellard
1135 7bc9090a Michael Niedermayer
        if (IS_QUANT(mb_type))
1136
            s->qscale = get_qscale(s);
1137
        
1138 de6d9b64 Fabrice Bellard
        if (s->concealment_motion_vectors) {
1139
            /* just parse them */
1140
            if (s->picture_structure != PICT_FRAME) 
1141 612476ef Arpi
                skip_bits1(&s->gb); /* field select */
1142 daab3296 Michael Niedermayer
            
1143
            s->mv[0][0][0]= s->last_mv[0][0][0]= s->last_mv[0][1][0] = 
1144
                mpeg_decode_motion(s, s->mpeg_f_code[0][0], s->last_mv[0][0][0]);
1145
            s->mv[0][0][1]= s->last_mv[0][0][1]= s->last_mv[0][1][1] = 
1146
                mpeg_decode_motion(s, s->mpeg_f_code[0][1], s->last_mv[0][0][1]);
1147
1148 7bc9090a Michael Niedermayer
            skip_bits1(&s->gb); /* marker */
1149 daab3296 Michael Niedermayer
        }else
1150
            memset(s->last_mv, 0, sizeof(s->last_mv)); /* reset mv prediction */
1151 de6d9b64 Fabrice Bellard
        s->mb_intra = 1;
1152 a579db0c Ivan Kalvachev
#ifdef HAVE_XVMC
1153
        //one 1 we memcpy blocks in xvmcvideo
1154
        if(s->avctx->xvmc_acceleration > 1){
1155
            XVMC_pack_pblocks(s,-1);//inter are always full blocks
1156
            if(s->swap_uv){
1157
                exchange_uv(s);
1158
            }
1159
        }
1160
#endif
1161 7bc9090a Michael Niedermayer
1162 029911d1 Michael Niedermayer
        if (s->codec_id == CODEC_ID_MPEG2VIDEO) {
1163 715731a3 Ivan Kalvachev
            for(i=0;i<mb_block_count;i++) {
1164 a579db0c Ivan Kalvachev
                if (mpeg2_decode_block_intra(s, s->pblocks[i], i) < 0)
1165 7bc9090a Michael Niedermayer
                    return -1;
1166
            }
1167
        } else {
1168
            for(i=0;i<6;i++) {
1169 a579db0c Ivan Kalvachev
                if (mpeg1_decode_block_intra(s, s->pblocks[i], i) < 0)
1170 7bc9090a Michael Niedermayer
                    return -1;
1171
            }
1172
        }
1173 de6d9b64 Fabrice Bellard
    } else {
1174 7bc9090a Michael Niedermayer
        if (mb_type & MB_TYPE_ZERO_MV){
1175 b40cd4e0 Michael Niedermayer
            assert(mb_type & MB_TYPE_CBP);
1176 7bc9090a Michael Niedermayer
1177
            /* compute dct type */
1178
            if (s->picture_structure == PICT_FRAME && //FIXME add a interlaced_dct coded var?
1179
                !s->frame_pred_frame_dct) {
1180
                s->interlaced_dct = get_bits1(&s->gb);
1181
            }
1182
1183
            if (IS_QUANT(mb_type))
1184
                s->qscale = get_qscale(s);
1185
1186
            s->mv_dir = MV_DIR_FORWARD;
1187 4c9544d0 Michael Niedermayer
            if(s->picture_structure == PICT_FRAME)
1188
                s->mv_type = MV_TYPE_16X16;
1189
            else{
1190
                s->mv_type = MV_TYPE_FIELD;
1191
                mb_type |= MB_TYPE_INTERLACED;
1192
                s->field_select[0][0]= s->picture_structure - 1;
1193
            }
1194 7bc9090a Michael Niedermayer
            s->last_mv[0][0][0] = 0;
1195
            s->last_mv[0][0][1] = 0;
1196
            s->last_mv[0][1][0] = 0;
1197
            s->last_mv[0][1][1] = 0;
1198
            s->mv[0][0][0] = 0;
1199
            s->mv[0][0][1] = 0;
1200
        }else{
1201
            assert(mb_type & MB_TYPE_L0L1);
1202
//FIXME decide if MBs in field pictures are MB_TYPE_INTERLACED
1203
            /* get additionnal motion vector type */
1204
            if (s->frame_pred_frame_dct) 
1205
                motion_type = MT_FRAME;
1206
            else{
1207
                motion_type = get_bits(&s->gb, 2);
1208
            }
1209
1210
            /* compute dct type */
1211
            if (s->picture_structure == PICT_FRAME && //FIXME add a interlaced_dct coded var?
1212 b40cd4e0 Michael Niedermayer
                !s->frame_pred_frame_dct && HAS_CBP(mb_type)) {
1213 7bc9090a Michael Niedermayer
                s->interlaced_dct = get_bits1(&s->gb);
1214
            }
1215
1216
            if (IS_QUANT(mb_type))
1217
                s->qscale = get_qscale(s);
1218
1219
            /* motion vectors */
1220
            s->mv_dir = 0;
1221
            for(i=0;i<2;i++) {
1222
                if (USES_LIST(mb_type, i)) {
1223
                    s->mv_dir |= (MV_DIR_FORWARD >> i);
1224
                    dprintf("motion_type=%d\n", motion_type);
1225
                    switch(motion_type) {
1226
                    case MT_FRAME: /* or MT_16X8 */
1227
                        if (s->picture_structure == PICT_FRAME) {
1228
                            /* MT_FRAME */
1229
                            mb_type |= MB_TYPE_16x16; 
1230
                            s->mv_type = MV_TYPE_16X16;
1231
                            s->mv[i][0][0]= s->last_mv[i][0][0]= s->last_mv[i][1][0] = 
1232
                                mpeg_decode_motion(s, s->mpeg_f_code[i][0], s->last_mv[i][0][0]);
1233
                            s->mv[i][0][1]= s->last_mv[i][0][1]= s->last_mv[i][1][1] = 
1234
                                mpeg_decode_motion(s, s->mpeg_f_code[i][1], s->last_mv[i][0][1]);
1235 de6d9b64 Fabrice Bellard
                            /* full_pel: only for mpeg1 */
1236 7bc9090a Michael Niedermayer
                            if (s->full_pel[i]){
1237
                                s->mv[i][0][0] <<= 1;
1238
                                s->mv[i][0][1] <<= 1;
1239
                            }
1240
                        } else {
1241
                            /* MT_16X8 */
1242 26f548a7 Michael Niedermayer
                            mb_type |= MB_TYPE_16x8 | MB_TYPE_INTERLACED; 
1243 7bc9090a Michael Niedermayer
                            s->mv_type = MV_TYPE_16X8;
1244
                            for(j=0;j<2;j++) {
1245
                                s->field_select[i][j] = get_bits1(&s->gb);
1246
                                for(k=0;k<2;k++) {
1247
                                    val = mpeg_decode_motion(s, s->mpeg_f_code[i][k],
1248
                                                             s->last_mv[i][j][k]);
1249
                                    s->last_mv[i][j][k] = val;
1250
                                    s->mv[i][j][k] = val;
1251
                                }
1252
                            }
1253 de6d9b64 Fabrice Bellard
                        }
1254 7bc9090a Michael Niedermayer
                        break;
1255
                    case MT_FIELD:
1256
                        s->mv_type = MV_TYPE_FIELD;
1257
                        if (s->picture_structure == PICT_FRAME) {
1258
                            mb_type |= MB_TYPE_16x8 | MB_TYPE_INTERLACED; 
1259
                            for(j=0;j<2;j++) {
1260
                                s->field_select[i][j] = get_bits1(&s->gb);
1261
                                val = mpeg_decode_motion(s, s->mpeg_f_code[i][0],
1262
                                                         s->last_mv[i][j][0]);
1263
                                s->last_mv[i][j][0] = val;
1264
                                s->mv[i][j][0] = val;
1265
                                dprintf("fmx=%d\n", val);
1266
                                val = mpeg_decode_motion(s, s->mpeg_f_code[i][1],
1267
                                                         s->last_mv[i][j][1] >> 1);
1268
                                s->last_mv[i][j][1] = val << 1;
1269
                                s->mv[i][j][1] = val;
1270
                                dprintf("fmy=%d\n", val);
1271
                            }
1272
                        } else {
1273 26f548a7 Michael Niedermayer
                            mb_type |= MB_TYPE_16x16 | MB_TYPE_INTERLACED; 
1274 7bc9090a Michael Niedermayer
                            s->field_select[i][0] = get_bits1(&s->gb);
1275 de6d9b64 Fabrice Bellard
                            for(k=0;k<2;k++) {
1276
                                val = mpeg_decode_motion(s, s->mpeg_f_code[i][k],
1277 7bc9090a Michael Niedermayer
                                                         s->last_mv[i][0][k]);
1278
                                s->last_mv[i][0][k] = val;
1279
                                s->last_mv[i][1][k] = val;
1280
                                s->mv[i][0][k] = val;
1281 de6d9b64 Fabrice Bellard
                            }
1282
                        }
1283 7bc9090a Michael Niedermayer
                        break;
1284
                    case MT_DMV:
1285
                        {
1286
                            int dmx, dmy, mx, my, m;
1287
1288
                            mx = mpeg_decode_motion(s, s->mpeg_f_code[i][0], 
1289
                                                    s->last_mv[i][0][0]);
1290
                            s->last_mv[i][0][0] = mx;
1291
                            s->last_mv[i][1][0] = mx;
1292
                            dmx = get_dmv(s);
1293
                            my = mpeg_decode_motion(s, s->mpeg_f_code[i][1], 
1294
                                                    s->last_mv[i][0][1] >> 1);
1295
                            dmy = get_dmv(s);
1296
                            s->mv_type = MV_TYPE_DMV;
1297 1dff7d56 Ivan Kalvachev
1298
1299
                            s->last_mv[i][0][1] = my<<1;
1300
                            s->last_mv[i][1][1] = my<<1;
1301
1302
                            s->mv[i][0][0] = mx;
1303
                            s->mv[i][0][1] = my;
1304
                            s->mv[i][1][0] = mx;//not used
1305
                            s->mv[i][1][1] = my;//not used
1306
1307 7bc9090a Michael Niedermayer
                            if (s->picture_structure == PICT_FRAME) {
1308
                                mb_type |= MB_TYPE_16x16 | MB_TYPE_INTERLACED; 
1309
1310 1dff7d56 Ivan Kalvachev
                                //m = 1 + 2 * s->top_field_first;
1311 7bc9090a Michael Niedermayer
                                m = s->top_field_first ? 1 : 3;
1312 1dff7d56 Ivan Kalvachev
1313 7bc9090a Michael Niedermayer
                                /* top -> top pred */
1314 1dff7d56 Ivan Kalvachev
                                s->mv[i][2][0] = ((mx * m + (mx > 0)) >> 1) + dmx;
1315
                                s->mv[i][2][1] = ((my * m + (my > 0)) >> 1) + dmy - 1;
1316 7bc9090a Michael Niedermayer
                                m = 4 - m;
1317
                                s->mv[i][3][0] = ((mx * m + (mx > 0)) >> 1) + dmx;
1318
                                s->mv[i][3][1] = ((my * m + (my > 0)) >> 1) + dmy + 1;
1319
                            } else {
1320
                                mb_type |= MB_TYPE_16x16;
1321
1322 1dff7d56 Ivan Kalvachev
                                s->mv[i][2][0] = ((mx + (mx > 0)) >> 1) + dmx;
1323
                                s->mv[i][2][1] = ((my + (my > 0)) >> 1) + dmy;
1324
                                if(s->picture_structure == PICT_TOP_FIELD)
1325
                                    s->mv[i][2][1]--;
1326
                                else 
1327
                                    s->mv[i][2][1]++;
1328 7bc9090a Michael Niedermayer
                            }
1329 de6d9b64 Fabrice Bellard
                        }
1330 7bc9090a Michael Niedermayer
                        break;
1331 26f548a7 Michael Niedermayer
                    default:
1332
                        av_log(s->avctx, AV_LOG_ERROR, "00 motion_type at %d %d\n", s->mb_x, s->mb_y);
1333
                        return -1;
1334 de6d9b64 Fabrice Bellard
                    }
1335
                }
1336
            }
1337
        }
1338 7bc9090a Michael Niedermayer
        
1339
        s->mb_intra = 0;
1340 de6d9b64 Fabrice Bellard
1341 b40cd4e0 Michael Niedermayer
        if (HAS_CBP(mb_type)) {
1342 7bc9090a Michael Niedermayer
            cbp = get_vlc2(&s->gb, mb_pat_vlc.table, MB_PAT_VLC_BITS, 1);
1343 494c56d3 Ivan Kalvachev
            if (cbp < 0 || ((cbp == 0) && (s->chroma_format < 2)) ){
1344 9b879566 Michel Bardiaux
                av_log(s->avctx, AV_LOG_ERROR, "invalid cbp at %d %d\n", s->mb_x, s->mb_y);
1345 7bc9090a Michael Niedermayer
                return -1;
1346 de6d9b64 Fabrice Bellard
            }
1347 715731a3 Ivan Kalvachev
            if(mb_block_count > 6){
1348
                 cbp<<= mb_block_count-6;
1349
                 cbp |= get_bits(&s->gb, mb_block_count-6);
1350 5e5c247a Ivan Kalvachev
            }
1351 7bc9090a Michael Niedermayer
1352 a579db0c Ivan Kalvachev
#ifdef HAVE_XVMC
1353
            //on 1 we memcpy blocks in xvmcvideo
1354
            if(s->avctx->xvmc_acceleration > 1){
1355
                XVMC_pack_pblocks(s,cbp);
1356
                if(s->swap_uv){
1357
                    exchange_uv(s);
1358
                }
1359
            }    
1360
#endif
1361
1362 029911d1 Michael Niedermayer
            if (s->codec_id == CODEC_ID_MPEG2VIDEO) {
1363 6fc5b059 Michael Niedermayer
                if(s->flags2 & CODEC_FLAG2_FAST){
1364
                    for(i=0;i<6;i++) {
1365
                        if(cbp & 32) {
1366
                            mpeg2_fast_decode_block_non_intra(s, s->pblocks[i], i);
1367
                        } else {
1368
                            s->block_last_index[i] = -1;
1369
                        }
1370
                        cbp+=cbp;
1371
                    }
1372
                }else{
1373
                    cbp<<= 12-mb_block_count;
1374
    
1375
                    for(i=0;i<mb_block_count;i++) {
1376
                        if ( cbp & (1<<11) ) {
1377
                            if (mpeg2_decode_block_non_intra(s, s->pblocks[i], i) < 0)
1378
                                return -1;
1379
                        } else {
1380
                            s->block_last_index[i] = -1;
1381
                        }
1382
                        cbp+=cbp;
1383 7bc9090a Michael Niedermayer
                    }
1384 de6d9b64 Fabrice Bellard
                }
1385 7bc9090a Michael Niedermayer
            } else {
1386 628b210f Michael Niedermayer
                if(s->flags2 & CODEC_FLAG2_FAST){
1387
                    for(i=0;i<6;i++) {
1388
                        if (cbp & 32) {
1389
                            mpeg1_fast_decode_block_inter(s, s->pblocks[i], i);
1390
                        } else {
1391
                            s->block_last_index[i] = -1;
1392
                        }
1393
                        cbp+=cbp;
1394
                    }
1395
                }else{
1396
                    for(i=0;i<6;i++) {
1397
                        if (cbp & 32) {
1398
                            if (mpeg1_decode_block_inter(s, s->pblocks[i], i) < 0)
1399
                                return -1;
1400
                        } else {
1401
                            s->block_last_index[i] = -1;
1402
                        }
1403
                        cbp+=cbp;
1404 7bc9090a Michael Niedermayer
                    }
1405 a0201736 Michael Niedermayer
                }
1406 de6d9b64 Fabrice Bellard
            }
1407 7bc9090a Michael Niedermayer
        }else{
1408
            for(i=0;i<6;i++)
1409
                s->block_last_index[i] = -1;
1410 de6d9b64 Fabrice Bellard
        }
1411
    }
1412 7bc9090a Michael Niedermayer
1413
    s->current_picture.mb_type[ s->mb_x + s->mb_y*s->mb_stride ]= mb_type;
1414
1415 de6d9b64 Fabrice Bellard
    return 0;
1416
}
1417
1418
/* as h263, but only 17 codes */
1419
static int mpeg_decode_motion(MpegEncContext *s, int fcode, int pred)
1420
{
1421 05858889 BERO
    int code, sign, val, l, shift;
1422 de6d9b64 Fabrice Bellard
1423 8ed2ddb2 Michael Niedermayer
    code = get_vlc2(&s->gb, mv_vlc.table, MV_VLC_BITS, 2);
1424 de6d9b64 Fabrice Bellard
    if (code == 0) {
1425
        return pred;
1426
    }
1427 7bc9090a Michael Niedermayer
    if (code < 0) {
1428
        return 0xffff;
1429
    }
1430
1431 612476ef Arpi
    sign = get_bits1(&s->gb);
1432 de6d9b64 Fabrice Bellard
    shift = fcode - 1;
1433 05858889 BERO
    val = code;
1434
    if (shift) {
1435
        val = (val - 1) << shift;
1436 de6d9b64 Fabrice Bellard
        val |= get_bits(&s->gb, shift);
1437 05858889 BERO
        val++;
1438
    }
1439 de6d9b64 Fabrice Bellard
    if (sign)
1440
        val = -val;
1441
    val += pred;
1442
    
1443
    /* modulo decoding */
1444 0ff93477 Michael Niedermayer
    l= INT_BIT - 5 - shift;
1445
    val = (val<<l)>>l;
1446 de6d9b64 Fabrice Bellard
    return val;
1447
}
1448
1449 c3bf0288 Michael Niedermayer
static inline int decode_dc(GetBitContext *gb, int component)
1450 de6d9b64 Fabrice Bellard
{
1451
    int code, diff;
1452
1453
    if (component == 0) {
1454 c3bf0288 Michael Niedermayer
        code = get_vlc2(gb, dc_lum_vlc.table, DC_VLC_BITS, 2);
1455 de6d9b64 Fabrice Bellard
    } else {
1456 c3bf0288 Michael Niedermayer
        code = get_vlc2(gb, dc_chroma_vlc.table, DC_VLC_BITS, 2);
1457 de6d9b64 Fabrice Bellard
    }
1458 9ac7ecd6 Michael Niedermayer
    if (code < 0){
1459 9b879566 Michel Bardiaux
        av_log(NULL, AV_LOG_ERROR, "invalid dc code at\n");
1460 de6d9b64 Fabrice Bellard
        return 0xffff;
1461 9ac7ecd6 Michael Niedermayer
    }
1462 de6d9b64 Fabrice Bellard
    if (code == 0) {
1463
        diff = 0;
1464
    } else {
1465 c3bf0288 Michael Niedermayer
        diff = get_xbits(gb, code);
1466 de6d9b64 Fabrice Bellard
    }
1467
    return diff;
1468
}
1469
1470 a0201736 Michael Niedermayer
static inline int mpeg1_decode_block_intra(MpegEncContext *s, 
1471 de6d9b64 Fabrice Bellard
                               DCTELEM *block, 
1472
                               int n)
1473
{
1474
    int level, dc, diff, i, j, run;
1475 a0201736 Michael Niedermayer
    int component;
1476 de6d9b64 Fabrice Bellard
    RLTable *rl = &rl_mpeg1;
1477 0c1a9eda Zdenek Kabelac
    uint8_t * const scantable= s->intra_scantable.permutated;
1478
    const uint16_t *quant_matrix= s->intra_matrix;
1479 a0201736 Michael Niedermayer
    const int qscale= s->qscale;
1480 de6d9b64 Fabrice Bellard
1481 a0201736 Michael Niedermayer
    /* DC coef */
1482
    component = (n <= 3 ? 0 : n - 4 + 1);
1483 c3bf0288 Michael Niedermayer
    diff = decode_dc(&s->gb, component);
1484 a0201736 Michael Niedermayer
    if (diff >= 0xffff)
1485
        return -1;
1486
    dc = s->last_dc[component];
1487
    dc += diff;
1488
    s->last_dc[component] = dc;
1489
    block[0] = dc<<3;
1490
    dprintf("dc=%d diff=%d\n", dc, diff);
1491
    i = 0;
1492
    {
1493
        OPEN_READER(re, &s->gb);    
1494
        /* now quantify & encode AC coefs */
1495
        for(;;) {
1496
            UPDATE_CACHE(re, &s->gb);
1497
            GET_RL_VLC(level, run, re, &s->gb, rl->rl_vlc[0], TEX_VLC_BITS, 2);
1498
            
1499
            if(level == 127){
1500
                break;
1501
            } else if(level != 0) {
1502
                i += run;
1503
                j = scantable[i];
1504 51929fd3 Michael Niedermayer
                level= (level*qscale*quant_matrix[j])>>4;
1505 a0201736 Michael Niedermayer
                level= (level-1)|1;
1506
                level = (level ^ SHOW_SBITS(re, &s->gb, 1)) - SHOW_SBITS(re, &s->gb, 1);
1507
                LAST_SKIP_BITS(re, &s->gb, 1);
1508
            } else {
1509
                /* escape */
1510
                run = SHOW_UBITS(re, &s->gb, 6)+1; LAST_SKIP_BITS(re, &s->gb, 6);
1511
                UPDATE_CACHE(re, &s->gb);
1512
                level = SHOW_SBITS(re, &s->gb, 8); SKIP_BITS(re, &s->gb, 8);
1513
                if (level == -128) {
1514
                    level = SHOW_UBITS(re, &s->gb, 8) - 256; LAST_SKIP_BITS(re, &s->gb, 8);
1515
                } else if (level == 0) {
1516
                    level = SHOW_UBITS(re, &s->gb, 8)      ; LAST_SKIP_BITS(re, &s->gb, 8);
1517
                }
1518
                i += run;
1519
                j = scantable[i];
1520
                if(level<0){
1521
                    level= -level;
1522 51929fd3 Michael Niedermayer
                    level= (level*qscale*quant_matrix[j])>>4;
1523 a0201736 Michael Niedermayer
                    level= (level-1)|1;
1524
                    level= -level;
1525
                }else{
1526 51929fd3 Michael Niedermayer
                    level= (level*qscale*quant_matrix[j])>>4;
1527 a0201736 Michael Niedermayer
                    level= (level-1)|1;
1528
                }
1529
            }
1530
            if (i > 63){
1531 9b879566 Michel Bardiaux
                av_log(s->avctx, AV_LOG_ERROR, "ac-tex damaged at %d %d\n", s->mb_x, s->mb_y);
1532 a0201736 Michael Niedermayer
                return -1;
1533
            }
1534
1535
            block[j] = level;
1536
        }
1537
        CLOSE_READER(re, &s->gb);
1538
    }
1539
    s->block_last_index[n] = i;
1540
   return 0;
1541
}
1542
1543
static inline int mpeg1_decode_block_inter(MpegEncContext *s, 
1544
                               DCTELEM *block, 
1545
                               int n)
1546
{
1547
    int level, i, j, run;
1548
    RLTable *rl = &rl_mpeg1;
1549 0c1a9eda Zdenek Kabelac
    uint8_t * const scantable= s->intra_scantable.permutated;
1550
    const uint16_t *quant_matrix= s->inter_matrix;
1551 a0201736 Michael Niedermayer
    const int qscale= s->qscale;
1552
1553
    {
1554 8db1a1dd Michael Niedermayer
        int v;
1555
        OPEN_READER(re, &s->gb);
1556 a0201736 Michael Niedermayer
        i = -1;
1557 de6d9b64 Fabrice Bellard
        /* special case for the first coef. no need to add a second vlc table */
1558 8db1a1dd Michael Niedermayer
        UPDATE_CACHE(re, &s->gb);
1559
        v= SHOW_UBITS(re, &s->gb, 2);
1560 de6d9b64 Fabrice Bellard
        if (v & 2) {
1561 3729c912 Michael Niedermayer
            LAST_SKIP_BITS(re, &s->gb, 2);
1562 51929fd3 Michael Niedermayer
            level= (3*qscale*quant_matrix[0])>>5;
1563 a0201736 Michael Niedermayer
            level= (level-1)|1;
1564
            if(v&1)
1565
                level= -level;
1566 3729c912 Michael Niedermayer
            block[0] = level;
1567 a0201736 Michael Niedermayer
            i++;
1568 de6d9b64 Fabrice Bellard
        }
1569 3729c912 Michael Niedermayer
1570 a0201736 Michael Niedermayer
        /* now quantify & encode AC coefs */
1571
        for(;;) {
1572
            UPDATE_CACHE(re, &s->gb);
1573
            GET_RL_VLC(level, run, re, &s->gb, rl->rl_vlc[0], TEX_VLC_BITS, 2);
1574
            
1575
            if(level == 127){
1576
                break;
1577
            } else if(level != 0) {
1578
                i += run;
1579
                j = scantable[i];
1580 51929fd3 Michael Niedermayer
                level= ((level*2+1)*qscale*quant_matrix[j])>>5;
1581 a0201736 Michael Niedermayer
                level= (level-1)|1;
1582
                level = (level ^ SHOW_SBITS(re, &s->gb, 1)) - SHOW_SBITS(re, &s->gb, 1);
1583
                LAST_SKIP_BITS(re, &s->gb, 1);
1584
            } else {
1585
                /* escape */
1586
                run = SHOW_UBITS(re, &s->gb, 6)+1; LAST_SKIP_BITS(re, &s->gb, 6);
1587
                UPDATE_CACHE(re, &s->gb);
1588
                level = SHOW_SBITS(re, &s->gb, 8); SKIP_BITS(re, &s->gb, 8);
1589
                if (level == -128) {
1590
                    level = SHOW_UBITS(re, &s->gb, 8) - 256; LAST_SKIP_BITS(re, &s->gb, 8);
1591
                } else if (level == 0) {
1592
                    level = SHOW_UBITS(re, &s->gb, 8)      ; LAST_SKIP_BITS(re, &s->gb, 8);
1593
                }
1594
                i += run;
1595
                j = scantable[i];
1596
                if(level<0){
1597
                    level= -level;
1598 51929fd3 Michael Niedermayer
                    level= ((level*2+1)*qscale*quant_matrix[j])>>5;
1599 a0201736 Michael Niedermayer
                    level= (level-1)|1;
1600
                    level= -level;
1601
                }else{
1602 51929fd3 Michael Niedermayer
                    level= ((level*2+1)*qscale*quant_matrix[j])>>5;
1603 a0201736 Michael Niedermayer
                    level= (level-1)|1;
1604
                }
1605 de6d9b64 Fabrice Bellard
            }
1606 a0201736 Michael Niedermayer
            if (i > 63){
1607 9b879566 Michel Bardiaux
                av_log(s->avctx, AV_LOG_ERROR, "ac-tex damaged at %d %d\n", s->mb_x, s->mb_y);
1608 a0201736 Michael Niedermayer
                return -1;
1609
            }
1610
1611
            block[j] = level;
1612 de6d9b64 Fabrice Bellard
        }
1613 a0201736 Michael Niedermayer
        CLOSE_READER(re, &s->gb);
1614 de6d9b64 Fabrice Bellard
    }
1615 a0201736 Michael Niedermayer
    s->block_last_index[n] = i;
1616 de6d9b64 Fabrice Bellard
    return 0;
1617
}
1618
1619 628b210f Michael Niedermayer
static inline int mpeg1_fast_decode_block_inter(MpegEncContext *s, DCTELEM *block, int n)
1620
{
1621
    int level, i, j, run;
1622
    RLTable *rl = &rl_mpeg1;
1623
    uint8_t * const scantable= s->intra_scantable.permutated;
1624
    const int qscale= s->qscale;
1625
1626
    {
1627
        int v;
1628
        OPEN_READER(re, &s->gb);
1629
        i = -1;
1630
        /* special case for the first coef. no need to add a second vlc table */
1631
        UPDATE_CACHE(re, &s->gb);
1632
        v= SHOW_UBITS(re, &s->gb, 2);
1633
        if (v & 2) {
1634
            LAST_SKIP_BITS(re, &s->gb, 2);
1635 3cb01872 Michael Niedermayer
            level= (3*qscale)>>1;
1636 628b210f Michael Niedermayer
            level= (level-1)|1;
1637
            if(v&1)
1638
                level= -level;
1639
            block[0] = level;
1640
            i++;
1641
        }
1642
1643
        /* now quantify & encode AC coefs */
1644
        for(;;) {
1645
            UPDATE_CACHE(re, &s->gb);
1646
            GET_RL_VLC(level, run, re, &s->gb, rl->rl_vlc[0], TEX_VLC_BITS, 2);
1647
            
1648
            if(level == 127){
1649
                break;
1650
            } else if(level != 0) {
1651
                i += run;
1652
                j = scantable[i];
1653
                level= ((level*2+1)*qscale)>>1;
1654
                level= (level-1)|1;
1655
                level = (level ^ SHOW_SBITS(re, &s->gb, 1)) - SHOW_SBITS(re, &s->gb, 1);
1656
                LAST_SKIP_BITS(re, &s->gb, 1);
1657
            } else {
1658
                /* escape */
1659
                run = SHOW_UBITS(re, &s->gb, 6)+1; LAST_SKIP_BITS(re, &s->gb, 6);
1660
                UPDATE_CACHE(re, &s->gb);
1661
                level = SHOW_SBITS(re, &s->gb, 8); SKIP_BITS(re, &s->gb, 8);
1662
                if (level == -128) {
1663
                    level = SHOW_UBITS(re, &s->gb, 8) - 256; LAST_SKIP_BITS(re, &s->gb, 8);
1664
                } else if (level == 0) {
1665
                    level = SHOW_UBITS(re, &s->gb, 8)      ; LAST_SKIP_BITS(re, &s->gb, 8);
1666
                }
1667
                i += run;
1668
                j = scantable[i];
1669
                if(level<0){
1670
                    level= -level;
1671
                    level= ((level*2+1)*qscale)>>1;
1672
                    level= (level-1)|1;
1673
                    level= -level;
1674
                }else{
1675
                    level= ((level*2+1)*qscale)>>1;
1676
                    level= (level-1)|1;
1677
                }
1678
            }
1679
1680
            block[j] = level;
1681
        }
1682
        CLOSE_READER(re, &s->gb);
1683
    }
1684
    s->block_last_index[n] = i;
1685
    return 0;
1686
}
1687
1688
1689 3729c912 Michael Niedermayer
static inline int mpeg2_decode_block_non_intra(MpegEncContext *s, 
1690
                               DCTELEM *block, 
1691
                               int n)
1692 de6d9b64 Fabrice Bellard
{
1693
    int level, i, j, run;
1694
    RLTable *rl = &rl_mpeg1;
1695 0c1a9eda Zdenek Kabelac
    uint8_t * const scantable= s->intra_scantable.permutated;
1696
    const uint16_t *quant_matrix;
1697 3729c912 Michael Niedermayer
    const int qscale= s->qscale;
1698 de6d9b64 Fabrice Bellard
    int mismatch;
1699
1700
    mismatch = 1;
1701
1702
    {
1703 8db1a1dd Michael Niedermayer
        int v;
1704
        OPEN_READER(re, &s->gb);
1705 3729c912 Michael Niedermayer
        i = -1;
1706 8db1a1dd Michael Niedermayer
        if (n < 4)
1707 3729c912 Michael Niedermayer
            quant_matrix = s->inter_matrix;
1708 de6d9b64 Fabrice Bellard
        else
1709 3729c912 Michael Niedermayer
            quant_matrix = s->chroma_inter_matrix;
1710 8db1a1dd Michael Niedermayer
1711 de6d9b64 Fabrice Bellard
        /* special case for the first coef. no need to add a second vlc table */
1712 8db1a1dd Michael Niedermayer
        UPDATE_CACHE(re, &s->gb);
1713
        v= SHOW_UBITS(re, &s->gb, 2);
1714 de6d9b64 Fabrice Bellard
        if (v & 2) {
1715 3729c912 Michael Niedermayer
            LAST_SKIP_BITS(re, &s->gb, 2);
1716
            level= (3*qscale*quant_matrix[0])>>5;
1717
            if(v&1)
1718
                level= -level;
1719
            block[0] = level;
1720
            mismatch ^= level;
1721
            i++;
1722 de6d9b64 Fabrice Bellard
        }
1723
1724 3729c912 Michael Niedermayer
        /* now quantify & encode AC coefs */
1725
        for(;;) {
1726
            UPDATE_CACHE(re, &s->gb);
1727
            GET_RL_VLC(level, run, re, &s->gb, rl->rl_vlc[0], TEX_VLC_BITS, 2);
1728
            
1729
            if(level == 127){
1730
                break;
1731
            } else if(level != 0) {
1732
                i += run;
1733
                j = scantable[i];
1734
                level= ((level*2+1)*qscale*quant_matrix[j])>>5;
1735
                level = (level ^ SHOW_SBITS(re, &s->gb, 1)) - SHOW_SBITS(re, &s->gb, 1);
1736
                LAST_SKIP_BITS(re, &s->gb, 1);
1737
            } else {
1738
                /* escape */
1739
                run = SHOW_UBITS(re, &s->gb, 6)+1; LAST_SKIP_BITS(re, &s->gb, 6);
1740
                UPDATE_CACHE(re, &s->gb);
1741
                level = SHOW_SBITS(re, &s->gb, 12); SKIP_BITS(re, &s->gb, 12);
1742
1743
                i += run;
1744
                j = scantable[i];
1745
                if(level<0){
1746
                    level= ((-level*2+1)*qscale*quant_matrix[j])>>5;
1747
                    level= -level;
1748
                }else{
1749
                    level= ((level*2+1)*qscale*quant_matrix[j])>>5;
1750
                }
1751
            }
1752
            if (i > 63){
1753 9b879566 Michel Bardiaux
                av_log(s->avctx, AV_LOG_ERROR, "ac-tex damaged at %d %d\n", s->mb_x, s->mb_y);
1754 3729c912 Michael Niedermayer
                return -1;
1755
            }
1756
            
1757
            mismatch ^= level;
1758
            block[j] = level;
1759 25ed7f92 Fabrice Bellard
        }
1760 3729c912 Michael Niedermayer
        CLOSE_READER(re, &s->gb);
1761 de6d9b64 Fabrice Bellard
    }
1762
    block[63] ^= (mismatch & 1);
1763 3729c912 Michael Niedermayer
    
1764 de6d9b64 Fabrice Bellard
    s->block_last_index[n] = i;
1765
    return 0;
1766
}
1767
1768 6fc5b059 Michael Niedermayer
static inline int mpeg2_fast_decode_block_non_intra(MpegEncContext *s, 
1769
                               DCTELEM *block, 
1770
                               int n)
1771
{
1772
    int level, i, j, run;
1773
    RLTable *rl = &rl_mpeg1;
1774
    uint8_t * const scantable= s->intra_scantable.permutated;
1775
    const int qscale= s->qscale;
1776
    int v;
1777
    OPEN_READER(re, &s->gb);
1778
    i = -1;
1779
1780
    /* special case for the first coef. no need to add a second vlc table */
1781
    UPDATE_CACHE(re, &s->gb);
1782
    v= SHOW_UBITS(re, &s->gb, 2);
1783
    if (v & 2) {
1784
        LAST_SKIP_BITS(re, &s->gb, 2);
1785
        level= (3*qscale)>>1;
1786
        if(v&1)
1787
            level= -level;
1788
        block[0] = level;
1789
        i++;
1790
    }
1791
1792
    /* now quantify & encode AC coefs */
1793
    for(;;) {
1794
        UPDATE_CACHE(re, &s->gb);
1795
        GET_RL_VLC(level, run, re, &s->gb, rl->rl_vlc[0], TEX_VLC_BITS, 2);
1796
        
1797
        if(level == 127){
1798
            break;
1799
        } else if(level != 0) {
1800
            i += run;
1801
            j = scantable[i];
1802
            level= ((level*2+1)*qscale)>>1;
1803
            level = (level ^ SHOW_SBITS(re, &s->gb, 1)) - SHOW_SBITS(re, &s->gb, 1);
1804
            LAST_SKIP_BITS(re, &s->gb, 1);
1805
        } else {
1806
            /* escape */
1807
            run = SHOW_UBITS(re, &s->gb, 6)+1; LAST_SKIP_BITS(re, &s->gb, 6);
1808
            UPDATE_CACHE(re, &s->gb);
1809
            level = SHOW_SBITS(re, &s->gb, 12); SKIP_BITS(re, &s->gb, 12);
1810
1811
            i += run;
1812
            j = scantable[i];
1813
            if(level<0){
1814
                level= ((-level*2+1)*qscale)>>1;
1815
                level= -level;
1816
            }else{
1817
                level= ((level*2+1)*qscale)>>1;
1818
            }
1819
        }
1820
        
1821
        block[j] = level;
1822
    }
1823
    CLOSE_READER(re, &s->gb);
1824
    s->block_last_index[n] = i;
1825
    return 0;
1826
}
1827
1828
1829 3729c912 Michael Niedermayer
static inline int mpeg2_decode_block_intra(MpegEncContext *s, 
1830
                               DCTELEM *block, 
1831
                               int n)
1832 de6d9b64 Fabrice Bellard
{
1833
    int level, dc, diff, i, j, run;
1834 3729c912 Michael Niedermayer
    int component;
1835 de6d9b64 Fabrice Bellard
    RLTable *rl;
1836 0c1a9eda Zdenek Kabelac
    uint8_t * const scantable= s->intra_scantable.permutated;
1837
    const uint16_t *quant_matrix;
1838 3729c912 Michael Niedermayer
    const int qscale= s->qscale;
1839 de6d9b64 Fabrice Bellard
    int mismatch;
1840
1841
    /* DC coef */
1842 3729c912 Michael Niedermayer
    if (n < 4){
1843
        quant_matrix = s->intra_matrix;
1844
        component = 0; 
1845
    }else{
1846
        quant_matrix = s->chroma_intra_matrix;
1847 5e5c247a Ivan Kalvachev
        component = (n&1) + 1;
1848 3729c912 Michael Niedermayer
    }
1849 c3bf0288 Michael Niedermayer
    diff = decode_dc(&s->gb, component);
1850 de6d9b64 Fabrice Bellard
    if (diff >= 0xffff)
1851
        return -1;
1852
    dc = s->last_dc[component];
1853
    dc += diff;
1854
    s->last_dc[component] = dc;
1855
    block[0] = dc << (3 - s->intra_dc_precision);
1856
    dprintf("dc=%d\n", block[0]);
1857 d753173a Fabrice Bellard
    mismatch = block[0] ^ 1;
1858 3729c912 Michael Niedermayer
    i = 0;
1859 de6d9b64 Fabrice Bellard
    if (s->intra_vlc_format)
1860
        rl = &rl_mpeg2;
1861
    else
1862
        rl = &rl_mpeg1;
1863 25ed7f92 Fabrice Bellard
1864 3729c912 Michael Niedermayer
    {
1865
        OPEN_READER(re, &s->gb);    
1866
        /* now quantify & encode AC coefs */
1867
        for(;;) {
1868
            UPDATE_CACHE(re, &s->gb);
1869
            GET_RL_VLC(level, run, re, &s->gb, rl->rl_vlc[0], TEX_VLC_BITS, 2);
1870
            
1871
            if(level == 127){
1872
                break;
1873
            } else if(level != 0) {
1874
                i += run;
1875
                j = scantable[i];
1876
                level= (level*qscale*quant_matrix[j])>>4;
1877
                level = (level ^ SHOW_SBITS(re, &s->gb, 1)) - SHOW_SBITS(re, &s->gb, 1);
1878
                LAST_SKIP_BITS(re, &s->gb, 1);
1879
            } else {
1880
                /* escape */
1881
                run = SHOW_UBITS(re, &s->gb, 6)+1; LAST_SKIP_BITS(re, &s->gb, 6);
1882
                UPDATE_CACHE(re, &s->gb);
1883
                level = SHOW_SBITS(re, &s->gb, 12); SKIP_BITS(re, &s->gb, 12);
1884
                i += run;
1885
                j = scantable[i];
1886
                if(level<0){
1887
                    level= (-level*qscale*quant_matrix[j])>>4;
1888
                    level= -level;
1889
                }else{
1890
                    level= (level*qscale*quant_matrix[j])>>4;
1891
                }
1892
            }
1893
            if (i > 63){
1894 9b879566 Michel Bardiaux
                av_log(s->avctx, AV_LOG_ERROR, "ac-tex damaged at %d %d\n", s->mb_x, s->mb_y);
1895 3729c912 Michael Niedermayer
                return -1;
1896
            }
1897
            
1898
            mismatch^= level;
1899
            block[j] = level;
1900 9ac7ecd6 Michael Niedermayer
        }
1901 3729c912 Michael Niedermayer
        CLOSE_READER(re, &s->gb);
1902 de6d9b64 Fabrice Bellard
    }
1903 3729c912 Michael Niedermayer
    block[63]^= mismatch&1;
1904
    
1905 de6d9b64 Fabrice Bellard
    s->block_last_index[n] = i;
1906
    return 0;
1907
}
1908
1909
typedef struct Mpeg1Context {
1910
    MpegEncContext mpeg_enc_ctx;
1911
    int mpeg_enc_ctx_allocated; /* true if decoding context allocated */
1912 1cb0edb4 Juanjo
    int repeat_field; /* true if we must repeat the field */
1913 fa384dcc Michael Niedermayer
    AVPanScan pan_scan; /** some temporary storage for the panscan */
1914 c62c07d3 Michael Niedermayer
    int slice_count;
1915 a6b9ffbf Ivan Kalvachev
    int swap_uv;//indicate VCR2
1916
    int save_aspect_info;
1917 c32b14bb Michael Niedermayer
    AVRational frame_rate_ext;       ///< MPEG-2 specific framerate modificator
1918 a6b9ffbf Ivan Kalvachev
1919 de6d9b64 Fabrice Bellard
} Mpeg1Context;
1920
1921
static int mpeg_decode_init(AVCodecContext *avctx)
1922
{
1923
    Mpeg1Context *s = avctx->priv_data;
1924 3edcacde Michael Niedermayer
    MpegEncContext *s2 = &s->mpeg_enc_ctx;
1925 d1700ead Ivan Kalvachev
    int i;
1926 3edcacde Michael Niedermayer
    
1927 d1700ead Ivan Kalvachev
    //we need some parmutation to store
1928
    //matrixes, until MPV_common_init()
1929
    //set the real permutatuon 
1930
    for(i=0;i<64;i++)
1931
       s2->dsp.idct_permutation[i]=i;
1932
1933 3edcacde Michael Niedermayer
    MPV_decode_defaults(s2);
1934 8f8402e4 Michael Niedermayer
    
1935 303aebf9 Michael Niedermayer
    s->mpeg_enc_ctx.avctx= avctx;
1936 94aec31f Michael Niedermayer
    s->mpeg_enc_ctx.flags= avctx->flags;
1937 303e50e6 Michael Niedermayer
    s->mpeg_enc_ctx.flags2= avctx->flags2;
1938 8f8402e4 Michael Niedermayer
    common_init(&s->mpeg_enc_ctx);
1939 c3bf0288 Michael Niedermayer
    init_vlcs();
1940 de6d9b64 Fabrice Bellard
1941
    s->mpeg_enc_ctx_allocated = 0;
1942
    s->mpeg_enc_ctx.picture_number = 0;
1943 1cb0edb4 Juanjo
    s->repeat_field = 0;
1944 d7e9533a Michael Niedermayer
    s->mpeg_enc_ctx.codec_id= avctx->codec->id;
1945 de6d9b64 Fabrice Bellard
    return 0;
1946
}
1947
1948 a6b9ffbf Ivan Kalvachev
static void quant_matrix_rebuild(uint16_t *matrix, const uint8_t *old_perm, 
1949
                                     const uint8_t *new_perm){
1950 c32b14bb Michael Niedermayer
    uint16_t temp_matrix[64];
1951
    int i;
1952 a6b9ffbf Ivan Kalvachev
1953
    memcpy(temp_matrix,matrix,64*sizeof(uint16_t));
1954
    
1955
    for(i=0;i<64;i++){
1956
        matrix[new_perm[i]] = temp_matrix[old_perm[i]];
1957
    }      
1958
}
1959
1960
//Call this function when we know all parameters
1961
//it may be called in different places for mpeg1 and mpeg2
1962
static int mpeg_decode_postinit(AVCodecContext *avctx){
1963 c32b14bb Michael Niedermayer
    Mpeg1Context *s1 = avctx->priv_data;
1964
    MpegEncContext *s = &s1->mpeg_enc_ctx;
1965
    uint8_t old_permutation[64];
1966 a6b9ffbf Ivan Kalvachev
1967
    if (
1968
            (s1->mpeg_enc_ctx_allocated == 0)|| 
1969
        avctx->width  != s->width ||
1970
        avctx->height != s->height||
1971
//      s1->save_aspect_info != avctx->aspect_ratio_info||
1972
        0)
1973
    {
1974
    
1975
        if (s1->mpeg_enc_ctx_allocated) {
1976
            MPV_common_end(s);
1977
        }
1978
1979
        if( (s->width == 0 )||(s->height == 0))
1980
            return -2;
1981
1982
        avctx->width = s->width;
1983
        avctx->height = s->height;
1984
        avctx->bit_rate = s->bit_rate;
1985
        s1->save_aspect_info = s->aspect_ratio_info;
1986
1987
     //low_delay may be forced, in this case we will have B frames
1988
     //that behave like P frames
1989
        avctx->has_b_frames = !(s->low_delay);
1990
1991
        if(avctx->sub_id==1){//s->codec_id==avctx->codec_id==CODEC_ID
1992
            //mpeg1 fps
1993
            avctx->frame_rate     = frame_rate_tab[s->frame_rate_index].num;
1994
            avctx->frame_rate_base= frame_rate_tab[s->frame_rate_index].den;
1995
            //mpeg1 aspect
1996
            avctx->sample_aspect_ratio= av_d2q(
1997
                    1.0/mpeg1_aspect[s->aspect_ratio_info], 255);
1998
1999
        }else{//mpeg2
2000
        //mpeg2 fps
2001
            av_reduce(
2002
                &s->avctx->frame_rate, 
2003
                &s->avctx->frame_rate_base, 
2004 c32b14bb Michael Niedermayer
                frame_rate_tab[s->frame_rate_index].num * s1->frame_rate_ext.num,
2005
                frame_rate_tab[s->frame_rate_index].den * s1->frame_rate_ext.den,
2006 a6b9ffbf Ivan Kalvachev
                1<<30);
2007
        //mpeg2 aspect
2008
            if(s->aspect_ratio_info > 1){
2009 6fe1a1a5 Ivan Kalvachev
                if( (s1->pan_scan.width == 0 )||(s1->pan_scan.height == 0) ){
2010 a6b9ffbf Ivan Kalvachev
                    s->avctx->sample_aspect_ratio= 
2011 5ac47683 Ivan Kalvachev
                        av_div_q(
2012
                         mpeg2_aspect[s->aspect_ratio_info], 
2013
                         (AVRational){s->width, s->height}
2014
                         );
2015 a6b9ffbf Ivan Kalvachev
                }else{
2016
                    s->avctx->sample_aspect_ratio= 
2017
                        av_div_q(
2018
                         mpeg2_aspect[s->aspect_ratio_info], 
2019 6fe1a1a5 Ivan Kalvachev
                         (AVRational){s1->pan_scan.width, s1->pan_scan.height}
2020 a6b9ffbf Ivan Kalvachev
                        );
2021
                }
2022
            }else{
2023
                s->avctx->sample_aspect_ratio= 
2024 5ac47683 Ivan Kalvachev
                    mpeg2_aspect[s->aspect_ratio_info];
2025 a6b9ffbf Ivan Kalvachev
            }
2026
        }//mpeg2
2027
2028
        if(avctx->xvmc_acceleration){
2029
            avctx->pix_fmt = avctx->get_format(avctx,pixfmt_xvmc_mpg2_420);
2030
        }else{
2031
            if(s->chroma_format <  2){
2032
                avctx->pix_fmt = avctx->get_format(avctx,pixfmt_yuv_420);
2033
            }else
2034
            if(s->chroma_format == 2){
2035
                avctx->pix_fmt = avctx->get_format(avctx,pixfmt_yuv_422);
2036
            }else
2037
            if(s->chroma_format >  2){
2038
                avctx->pix_fmt = avctx->get_format(avctx,pixfmt_yuv_444);
2039
            }
2040
        }
2041
        //until then pix_fmt may be changed right after codec init
2042
        if( avctx->pix_fmt == PIX_FMT_XVMC_MPEG2_IDCT )
2043
            if( avctx->idct_algo == FF_IDCT_AUTO )
2044
                avctx->idct_algo = FF_IDCT_SIMPLE;
2045
2046
        //quantization matrixes may need reordering 
2047
        //if dct permutation is changed
2048
        memcpy(old_permutation,s->dsp.idct_permutation,64*sizeof(uint8_t));
2049
2050
        if (MPV_common_init(s) < 0)
2051
            return -2;
2052
2053
        quant_matrix_rebuild(s->intra_matrix,       old_permutation,s->dsp.idct_permutation);
2054
        quant_matrix_rebuild(s->inter_matrix,       old_permutation,s->dsp.idct_permutation);
2055
        quant_matrix_rebuild(s->chroma_intra_matrix,old_permutation,s->dsp.idct_permutation);
2056
        quant_matrix_rebuild(s->chroma_inter_matrix,old_permutation,s->dsp.idct_permutation);
2057
2058
        s1->mpeg_enc_ctx_allocated = 1;
2059
    }
2060
    return 0;
2061
}
2062
2063 de6d9b64 Fabrice Bellard
/* return the 8 bit start code value and update the search
2064
   state. Return -1 if no start code found */
2065 5f0f7713 Zdenek Kabelac
static int find_start_code(const uint8_t **pbuf_ptr, const uint8_t *buf_end)
2066 de6d9b64 Fabrice Bellard
{
2067 c7516a00 Michael Niedermayer
    const uint8_t *buf_ptr= *pbuf_ptr;
2068
2069
    buf_ptr++; //gurantees that -1 is within the array
2070
    buf_end -= 2; // gurantees that +2 is within the array
2071 de6d9b64 Fabrice Bellard
2072
    while (buf_ptr < buf_end) {
2073 c7516a00 Michael Niedermayer
        if(*buf_ptr==0){
2074
            while(buf_ptr < buf_end && buf_ptr[1]==0)
2075
                buf_ptr++;
2076
2077
            if(buf_ptr[-1] == 0 && buf_ptr[1] == 1){
2078
                *pbuf_ptr = buf_ptr+3;
2079
                return buf_ptr[2] + 0x100;
2080
            }
2081 de6d9b64 Fabrice Bellard
        }
2082 c7516a00 Michael Niedermayer
        buf_ptr += 2;
2083 de6d9b64 Fabrice Bellard
    }
2084 c7516a00 Michael Niedermayer
    buf_end += 2; //undo the hack above
2085
    
2086
    *pbuf_ptr = buf_end;
2087
    return -1;
2088 de6d9b64 Fabrice Bellard
}
2089
2090
static int mpeg1_decode_picture(AVCodecContext *avctx, 
2091 5f0f7713 Zdenek Kabelac
                                const uint8_t *buf, int buf_size)
2092 de6d9b64 Fabrice Bellard
{
2093
    Mpeg1Context *s1 = avctx->priv_data;
2094
    MpegEncContext *s = &s1->mpeg_enc_ctx;
2095 d60a8f85 Michael Niedermayer
    int ref, f_code, vbv_delay;
2096 de6d9b64 Fabrice Bellard
2097 a6b9ffbf Ivan Kalvachev
    if(mpeg_decode_postinit(s->avctx) < 0) 
2098
       return -2;
2099
2100 68f593b4 Michael Niedermayer
    init_get_bits(&s->gb, buf, buf_size*8);
2101 de6d9b64 Fabrice Bellard
2102
    ref = get_bits(&s->gb, 10); /* temporal ref */
2103
    s->pict_type = get_bits(&s->gb, 3);
2104 58bfafbe Michael Niedermayer
2105 d60a8f85 Michael Niedermayer
    vbv_delay= get_bits(&s->gb, 16);
2106 de6d9b64 Fabrice Bellard
    if (s->pict_type == P_TYPE || s->pict_type == B_TYPE) {
2107 612476ef Arpi
        s->full_pel[0] = get_bits1(&s->gb);
2108 de6d9b64 Fabrice Bellard
        f_code = get_bits(&s->gb, 3);
2109
        if (f_code == 0)
2110
            return -1;
2111
        s->mpeg_f_code[0][0] = f_code;
2112
        s->mpeg_f_code[0][1] = f_code;
2113
    }
2114
    if (s->pict_type == B_TYPE) {
2115 612476ef Arpi
        s->full_pel[1] = get_bits1(&s->gb);
2116 de6d9b64 Fabrice Bellard
        f_code = get_bits(&s->gb, 3);
2117
        if (f_code == 0)
2118
            return -1;
2119
        s->mpeg_f_code[1][0] = f_code;
2120
        s->mpeg_f_code[1][1] = f_code;
2121
    }
2122 1e491e29 Michael Niedermayer
    s->current_picture.pict_type= s->pict_type;
2123
    s->current_picture.key_frame= s->pict_type == I_TYPE;
2124 d9cb5429 Michael Niedermayer
    
2125 b1b8f9fc Michael Niedermayer
//    if(avctx->debug & FF_DEBUG_PICT_INFO)
2126
//        av_log(avctx, AV_LOG_DEBUG, "vbv_delay %d, ref %d\n", vbv_delay, ref);
2127
    
2128 de6d9b64 Fabrice Bellard
    s->y_dc_scale = 8;
2129
    s->c_dc_scale = 8;
2130
    s->first_slice = 1;
2131
    return 0;
2132
}
2133
2134 c32b14bb Michael Niedermayer
static void mpeg_decode_sequence_extension(Mpeg1Context *s1)
2135 de6d9b64 Fabrice Bellard
{
2136 c32b14bb Michael Niedermayer
    MpegEncContext *s= &s1->mpeg_enc_ctx;
2137 de6d9b64 Fabrice Bellard
    int horiz_size_ext, vert_size_ext;
2138 baaf3f46 Michael Niedermayer
    int bit_rate_ext;
2139 de6d9b64 Fabrice Bellard
2140 029911d1 Michael Niedermayer
    skip_bits(&s->gb, 1); /* profil and level esc*/
2141 baced9f5 Michael Niedermayer
    s->avctx->profile= get_bits(&s->gb, 3);
2142
    s->avctx->level= get_bits(&s->gb, 4);
2143 1cb0edb4 Juanjo
    s->progressive_sequence = get_bits1(&s->gb); /* progressive_sequence */
2144 5e5c247a Ivan Kalvachev
    s->chroma_format = get_bits(&s->gb, 2); /* chroma_format 1=420, 2=422, 3=444 */
2145 de6d9b64 Fabrice Bellard
    horiz_size_ext = get_bits(&s->gb, 2);
2146
    vert_size_ext = get_bits(&s->gb, 2);
2147
    s->width |= (horiz_size_ext << 12);
2148
    s->height |= (vert_size_ext << 12);
2149
    bit_rate_ext = get_bits(&s->gb, 12);  /* XXX: handle it */
2150 5e5c247a Ivan Kalvachev
    s->bit_rate += (bit_rate_ext << 12) * 400;
2151 612476ef Arpi
    skip_bits1(&s->gb); /* marker */
2152 0368c72d Michael Niedermayer
    s->avctx->rc_buffer_size += get_bits(&s->gb, 8)*1024*16<<10;
2153 e0560448 Michael Niedermayer
2154 945f15b7 Michael Niedermayer
    s->low_delay = get_bits1(&s->gb);
2155 e0560448 Michael Niedermayer
    if(s->flags & CODEC_FLAG_LOW_DELAY) s->low_delay=1;
2156
2157 c32b14bb Michael Niedermayer
    s1->frame_rate_ext.num = get_bits(&s->gb, 2)+1;
2158
    s1->frame_rate_ext.den = get_bits(&s->gb, 5)+1;
2159 14bea432 Michael Niedermayer
2160 de6d9b64 Fabrice Bellard
    dprintf("sequence extension\n");
2161 029911d1 Michael Niedermayer
    s->codec_id= s->avctx->codec_id= CODEC_ID_MPEG2VIDEO;
2162 fb4a4a56 Fabrice Bellard
    s->avctx->sub_id = 2; /* indicates mpeg2 found */
2163 945f15b7 Michael Niedermayer
2164 029911d1 Michael Niedermayer
    if(s->avctx->debug & FF_DEBUG_PICT_INFO)
2165 0368c72d Michael Niedermayer
        av_log(s->avctx, AV_LOG_DEBUG, "profile: %d, level: %d vbv buffer: %d, bitrate:%d\n", 
2166 baced9f5 Michael Niedermayer
               s->avctx->profile, s->avctx->level, s->avctx->rc_buffer_size, s->bit_rate);
2167 a6b9ffbf Ivan Kalvachev
2168 de6d9b64 Fabrice Bellard
}
2169
2170 fa384dcc Michael Niedermayer
static void mpeg_decode_sequence_display_extension(Mpeg1Context *s1)
2171
{
2172
    MpegEncContext *s= &s1->mpeg_enc_ctx;
2173
    int color_description, w, h;
2174
2175
    skip_bits(&s->gb, 3); /* video format */
2176
    color_description= get_bits1(&s->gb);
2177
    if(color_description){
2178
        skip_bits(&s->gb, 8); /* color primaries */
2179
        skip_bits(&s->gb, 8); /* transfer_characteristics */
2180
        skip_bits(&s->gb, 8); /* matrix_coefficients */
2181
    }
2182
    w= get_bits(&s->gb, 14);
2183
    skip_bits(&s->gb, 1); //marker
2184
    h= get_bits(&s->gb, 14);
2185
    skip_bits(&s->gb, 1); //marker
2186
    
2187
    s1->pan_scan.width= 16*w;
2188
    s1->pan_scan.height=16*h;
2189 6fe1a1a5 Ivan Kalvachev
        
2190 fa384dcc Michael Niedermayer
    if(s->avctx->debug & FF_DEBUG_PICT_INFO)
2191 9b879566 Michel Bardiaux
        av_log(s->avctx, AV_LOG_DEBUG, "sde w:%d, h:%d\n", w, h);
2192 fa384dcc Michael Niedermayer
}
2193
2194
static void mpeg_decode_picture_display_extension(Mpeg1Context *s1)
2195
{
2196
    MpegEncContext *s= &s1->mpeg_enc_ctx;
2197 6fe1a1a5 Ivan Kalvachev
    int i,nofco;
2198
2199
    nofco = 1;
2200
    if(s->progressive_sequence){
2201
        if(s->repeat_first_field){
2202
            nofco++;
2203
            if(s->top_field_first)
2204
                nofco++;        
2205
        }
2206
    }else{
2207
        if(s->picture_structure == PICT_FRAME){
2208
            nofco++;
2209
            if(s->repeat_first_field)
2210
                nofco++;
2211
        }
2212
    }
2213
    for(i=0; i<nofco; i++){
2214 fa384dcc Michael Niedermayer
        s1->pan_scan.position[i][0]= get_sbits(&s->gb, 16);
2215
        skip_bits(&s->gb, 1); //marker
2216
        s1->pan_scan.position[i][1]= get_sbits(&s->gb, 16);
2217
        skip_bits(&s->gb, 1); //marker
2218
    }
2219
   
2220
    if(s->avctx->debug & FF_DEBUG_PICT_INFO)
2221 9b879566 Michel Bardiaux
        av_log(s->avctx, AV_LOG_DEBUG, "pde (%d,%d) (%d,%d) (%d,%d)\n", 
2222 fa384dcc Michael Niedermayer
            s1->pan_scan.position[0][0], s1->pan_scan.position[0][1], 
2223
            s1->pan_scan.position[1][0], s1->pan_scan.position[1][1], 
2224
            s1->pan_scan.position[2][0], s1->pan_scan.position[2][1]
2225
        );
2226
}
2227
2228 de6d9b64 Fabrice Bellard
static void mpeg_decode_quant_matrix_extension(MpegEncContext *s)
2229
{
2230 60832448 Fabrice Bellard
    int i, v, j;
2231 de6d9b64 Fabrice Bellard
2232 25ed7f92 Fabrice Bellard
    dprintf("matrix extension\n");
2233
2234 612476ef Arpi
    if (get_bits1(&s->gb)) {
2235 de6d9b64 Fabrice Bellard
        for(i=0;i<64;i++) {
2236
            v = get_bits(&s->gb, 8);
2237 b0368839 Michael Niedermayer
            j= s->dsp.idct_permutation[ ff_zigzag_direct[i] ];
2238 60832448 Fabrice Bellard
            s->intra_matrix[j] = v;
2239
            s->chroma_intra_matrix[j] = v;
2240 de6d9b64 Fabrice Bellard
        }
2241
    }
2242 612476ef Arpi
    if (get_bits1(&s->gb)) {
2243 de6d9b64 Fabrice Bellard
        for(i=0;i<64;i++) {
2244
            v = get_bits(&s->gb, 8);
2245 b0368839 Michael Niedermayer
            j= s->dsp.idct_permutation[ ff_zigzag_direct[i] ];
2246 d7e9533a Michael Niedermayer
            s->inter_matrix[j] = v;
2247
            s->chroma_inter_matrix[j] = v;
2248 de6d9b64 Fabrice Bellard
        }
2249
    }
2250 612476ef Arpi
    if (get_bits1(&s->gb)) {
2251 de6d9b64 Fabrice Bellard
        for(i=0;i<64;i++) {
2252
            v = get_bits(&s->gb, 8);
2253 b0368839 Michael Niedermayer
            j= s->dsp.idct_permutation[ ff_zigzag_direct[i] ];
2254 60832448 Fabrice Bellard
            s->chroma_intra_matrix[j] = v;
2255 de6d9b64 Fabrice Bellard
        }
2256
    }
2257 612476ef Arpi
    if (get_bits1(&s->gb)) {
2258 de6d9b64 Fabrice Bellard
        for(i=0;i<64;i++) {
2259
            v = get_bits(&s->gb, 8);
2260 b0368839 Michael Niedermayer
            j= s->dsp.idct_permutation[ ff_zigzag_direct[i] ];
2261 d7e9533a Michael Niedermayer
            s->chroma_inter_matrix[j] = v;
2262 de6d9b64 Fabrice Bellard
        }
2263
    }
2264
}
2265
2266
static void mpeg_decode_picture_coding_extension(MpegEncContext *s)
2267
{
2268
    s->full_pel[0] = s->full_pel[1] = 0;
2269
    s->mpeg_f_code[0][0] = get_bits(&s->gb, 4);
2270
    s->mpeg_f_code[0][1] = get_bits(&s->gb, 4);
2271
    s->mpeg_f_code[1][0] = get_bits(&s->gb, 4);
2272
    s->mpeg_f_code[1][1] = get_bits(&s->gb, 4);
2273
    s->intra_dc_precision = get_bits(&s->gb, 2);
2274
    s->picture_structure = get_bits(&s->gb, 2);
2275 612476ef Arpi
    s->top_field_first = get_bits1(&s->gb);
2276
    s->frame_pred_frame_dct = get_bits1(&s->gb);
2277
    s->concealment_motion_vectors = get_bits1(&s->gb);
2278
    s->q_scale_type = get_bits1(&s->gb);
2279
    s->intra_vlc_format = get_bits1(&s->gb);
2280
    s->alternate_scan = get_bits1(&s->gb);
2281
    s->repeat_first_field = get_bits1(&s->gb);
2282
    s->chroma_420_type = get_bits1(&s->gb);
2283
    s->progressive_frame = get_bits1(&s->gb);
2284 bb198e19 Michael Niedermayer
2285 dfb476cb Michael Niedermayer
    if(s->picture_structure == PICT_FRAME)
2286
        s->first_field=0;
2287
    else{
2288
        s->first_field ^= 1;
2289 7bc9090a Michael Niedermayer
        memset(s->mbskip_table, 0, s->mb_stride*s->mb_height);
2290 dfb476cb Michael Niedermayer
    }
2291
    
2292 acf44abb Michael Niedermayer
    if(s->alternate_scan){
2293 3d2e8cce Michael Niedermayer
        ff_init_scantable(s->dsp.idct_permutation, &s->inter_scantable  , ff_alternate_vertical_scan);
2294
        ff_init_scantable(s->dsp.idct_permutation, &s->intra_scantable  , ff_alternate_vertical_scan);
2295 acf44abb Michael Niedermayer
    }else{
2296 3d2e8cce Michael Niedermayer
        ff_init_scantable(s->dsp.idct_permutation, &s->inter_scantable  , ff_zigzag_direct);
2297
        ff_init_scantable(s->dsp.idct_permutation, &s->intra_scantable  , ff_zigzag_direct);
2298 acf44abb Michael Niedermayer
    }
2299
    
2300 de6d9b64 Fabrice Bellard
    /* composite display not parsed */
2301 1cb0edb4 Juanjo
    dprintf("intra_dc_precision=%d\n", s->intra_dc_precision);
2302 25ed7f92 Fabrice Bellard
    dprintf("picture_structure=%d\n", s->picture_structure);
2303 e0a3d744 Juanjo
    dprintf("top field first=%d\n", s->top_field_first);
2304
    dprintf("repeat first field=%d\n", s->repeat_first_field);
2305 de6d9b64 Fabrice Bellard
    dprintf("conceal=%d\n", s->concealment_motion_vectors);
2306 25ed7f92 Fabrice Bellard
    dprintf("intra_vlc_format=%d\n", s->intra_vlc_format);
2307
    dprintf("alternate_scan=%d\n", s->alternate_scan);
2308 de6d9b64 Fabrice Bellard
    dprintf("frame_pred_frame_dct=%d\n", s->frame_pred_frame_dct);
2309 1cb0edb4 Juanjo
    dprintf("progressive_frame=%d\n", s->progressive_frame);
2310 de6d9b64 Fabrice Bellard
}
2311
2312
static void mpeg_decode_extension(AVCodecContext *avctx, 
2313 5f0f7713 Zdenek Kabelac
                                  const uint8_t *buf, int buf_size)
2314 de6d9b64 Fabrice Bellard
{
2315
    Mpeg1Context *s1 = avctx->priv_data;
2316
    MpegEncContext *s = &s1->mpeg_enc_ctx;
2317
    int ext_type;
2318
2319 68f593b4 Michael Niedermayer
    init_get_bits(&s->gb, buf, buf_size*8);
2320 de6d9b64 Fabrice Bellard
    
2321
    ext_type = get_bits(&s->gb, 4);
2322
    switch(ext_type) {
2323
    case 0x1:
2324 c32b14bb Michael Niedermayer
        mpeg_decode_sequence_extension(s1);
2325 de6d9b64 Fabrice Bellard
        break;
2326 fa384dcc Michael Niedermayer
    case 0x2:
2327
        mpeg_decode_sequence_display_extension(s1);
2328
        break;
2329 de6d9b64 Fabrice Bellard
    case 0x3:
2330
        mpeg_decode_quant_matrix_extension(s);
2331
        break;
2332 fa384dcc Michael Niedermayer
    case 0x7:
2333
        mpeg_decode_picture_display_extension(s1);
2334
        break;
2335 de6d9b64 Fabrice Bellard
    case 0x8:
2336
        mpeg_decode_picture_coding_extension(s);
2337
        break;
2338
    }
2339
}
2340
2341 a579db0c Ivan Kalvachev
static void exchange_uv(MpegEncContext *s){
2342 c32b14bb Michael Niedermayer
    short * tmp = s->pblocks[4];
2343 a579db0c Ivan Kalvachev
    s->pblocks[4] = s->pblocks[5];
2344
    s->pblocks[5] = tmp;
2345 ff862be5 Michael Niedermayer
}
2346
2347 c62c07d3 Michael Niedermayer
static int mpeg_field_start(MpegEncContext *s){
2348
    AVCodecContext *avctx= s->avctx;
2349
    Mpeg1Context *s1 = (Mpeg1Context*)s;
2350 de6d9b64 Fabrice Bellard
2351
    /* start frame decoding */
2352 c62c07d3 Michael Niedermayer
    if(s->first_field || s->picture_structure==PICT_FRAME){
2353 d6db1c9c Michael Niedermayer
        if(MPV_frame_start(s, avctx) < 0)
2354 c62c07d3 Michael Niedermayer
            return -1;
2355 7bc9090a Michael Niedermayer
2356
        ff_er_frame_start(s);
2357
2358 2ec23b6d Michael Niedermayer
        /* first check if we must repeat the frame */
2359 9ee2c20e Fabrice Bellard
        s->current_picture_ptr->repeat_pict = 0;
2360 2ec23b6d Michael Niedermayer
        if (s->repeat_first_field) {
2361
            if (s->progressive_sequence) {
2362
                if (s->top_field_first)
2363 9ee2c20e Fabrice Bellard
                    s->current_picture_ptr->repeat_pict = 4;
2364 2ec23b6d Michael Niedermayer
                else
2365 9ee2c20e Fabrice Bellard
                    s->current_picture_ptr->repeat_pict = 2;
2366 2ec23b6d Michael Niedermayer
            } else if (s->progressive_frame) {
2367 9ee2c20e Fabrice Bellard
                s->current_picture_ptr->repeat_pict = 1;
2368 2ec23b6d Michael Niedermayer
            }
2369
        }         
2370 fa384dcc Michael Niedermayer
2371
        *s->current_picture_ptr->pan_scan= s1->pan_scan;
2372 c62c07d3 Michael Niedermayer
    }else{ //second field
2373 b536d0aa Michael Niedermayer
            int i;
2374 37b787f1 Michael Niedermayer
            
2375
            if(!s->current_picture_ptr){
2376 9b879566 Michel Bardiaux
                av_log(s->avctx, AV_LOG_ERROR, "first field missing\n");
2377 37b787f1 Michael Niedermayer
                return -1;
2378
            }
2379
            
2380 b536d0aa Michael Niedermayer
            for(i=0; i<4; i++){
2381
                s->current_picture.data[i] = s->current_picture_ptr->data[i];
2382
                if(s->picture_structure == PICT_BOTTOM_FIELD){
2383
                    s->current_picture.data[i] += s->current_picture_ptr->linesize[i];
2384
                } 
2385
            }
2386 c62c07d3 Michael Niedermayer
    }
2387 2e7b4c84 Ivan Kalvachev
#ifdef HAVE_XVMC
2388
// MPV_frame_start will call this function too,
2389
// but we need to call it on every field
2390 c62c07d3 Michael Niedermayer
    if(s->avctx->xvmc_acceleration)
2391 2e7b4c84 Ivan Kalvachev
         XVMC_field_start(s,avctx);
2392
#endif
2393 de6d9b64 Fabrice Bellard
2394 c62c07d3 Michael Niedermayer
    return 0;
2395
}
2396
2397
#define DECODE_SLICE_ERROR -1
2398
#define DECODE_SLICE_OK 0
2399
2400
/**
2401
 * decodes a slice. MpegEncContext.mb_y must be set to the MB row from the startcode
2402
 * @return DECODE_SLICE_ERROR if the slice is damaged<br>
2403
 *         DECODE_SLICE_OK if this slice is ok<br>
2404
 */
2405
static int mpeg_decode_slice(Mpeg1Context *s1, int mb_y,
2406 5f0f7713 Zdenek Kabelac
                             const uint8_t **buf, int buf_size)
2407 c62c07d3 Michael Niedermayer
{
2408
    MpegEncContext *s = &s1->mpeg_enc_ctx;
2409
    AVCodecContext *avctx= s->avctx;
2410
    int ret;
2411
    const int field_pic= s->picture_structure != PICT_FRAME;
2412 178fcca8 Michael Niedermayer
    const int lowres= s->avctx->lowres;
2413 c62c07d3 Michael Niedermayer
2414
    s->resync_mb_x=
2415
    s->resync_mb_y= -1;
2416
2417 25ef43bb Michael Niedermayer
    if (mb_y<<field_pic >= s->mb_height){
2418
        av_log(s->avctx, AV_LOG_ERROR, "slice below image (%d >= %d)\n", mb_y, s->mb_height);
2419 c62c07d3 Michael Niedermayer
        return -1;
2420
    }
2421
    
2422 80097bbf Michael Niedermayer
    init_get_bits(&s->gb, *buf, buf_size*8);
2423 de6d9b64 Fabrice Bellard
2424 c62c07d3 Michael Niedermayer
    ff_mpeg1_clean_buffers(s);
2425
    s->interlaced_dct = 0;
2426
2427 0ee50938 Fabrice Bellard
    s->qscale = get_qscale(s);
2428 a4337a51 Michael Niedermayer
2429 7bc9090a Michael Niedermayer
    if(s->qscale == 0){
2430 9b879566 Michel Bardiaux
        av_log(s->avctx, AV_LOG_ERROR, "qscale == 0\n");
2431 db6e7795 Michael Niedermayer
        return -1;
2432 7bc9090a Michael Niedermayer
    }
2433
    
2434 de6d9b64 Fabrice Bellard
    /* extra slice info */
2435 612476ef Arpi
    while (get_bits1(&s->gb) != 0) {
2436
        skip_bits(&s->gb, 8);
2437 de6d9b64 Fabrice Bellard
    }
2438 7bc9090a Michael Niedermayer
    
2439 ce5b7c5e Michael Niedermayer
    s->mb_x=0;
2440 7bc9090a Michael Niedermayer
2441 ce5b7c5e Michael Niedermayer
    for(;;) {
2442
        int code = get_vlc2(&s->gb, mbincr_vlc.table, MBINCR_VLC_BITS, 2);
2443 7bc9090a Michael Niedermayer
        if (code < 0){
2444 9b879566 Michel Bardiaux
            av_log(s->avctx, AV_LOG_ERROR, "first mb_incr damaged\n");
2445 db6e7795 Michael Niedermayer
            return -1;
2446 7bc9090a Michael Niedermayer
        }
2447 ce5b7c5e Michael Niedermayer
        if (code >= 33) {
2448
            if (code == 33) {
2449
                s->mb_x += 33;
2450
            }
2451
            /* otherwise, stuffing, nothing to do */
2452
        } else {
2453
            s->mb_x += code;
2454
            break;
2455
        }
2456
    }
2457 c62c07d3 Michael Niedermayer
2458 7bc9090a Michael Niedermayer
    s->resync_mb_x= s->mb_x;
2459 c62c07d3 Michael Niedermayer
    s->resync_mb_y= s->mb_y= mb_y;
2460 9b8709d1 Michael Niedermayer
    s->mb_skip_run= 0;
2461 7d1c3fc1 Michael Niedermayer
    ff_init_block_index(s);
2462 ce5b7c5e Michael Niedermayer
2463 c62c07d3 Michael Niedermayer
    if (s->mb_y==0 && s->mb_x==0 && (s->first_field || s->picture_structure==PICT_FRAME)) {
2464
        if(s->avctx->debug&FF_DEBUG_PICT_INFO){
2465
             av_log(s->avctx, AV_LOG_DEBUG, "qp:%d fc:%2d%2d%2d%2d %s %s %s %s %s dc:%d pstruct:%d fdct:%d cmv:%d qtype:%d ivlc:%d rff:%d %s\n", 
2466
                 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],
2467
                 s->pict_type == I_TYPE ? "I" : (s->pict_type == P_TYPE ? "P" : (s->pict_type == B_TYPE ? "B" : "S")), 
2468
                 s->progressive_sequence ? "ps" :"", s->progressive_frame ? "pf" : "", s->alternate_scan ? "alt" :"", s->top_field_first ? "top" :"", 
2469
                 s->intra_dc_precision, s->picture_structure, s->frame_pred_frame_dct, s->concealment_motion_vectors,
2470
                 s->q_scale_type, s->intra_vlc_format, s->repeat_first_field, s->chroma_420_type ? "420" :"");
2471
        }
2472
    }    
2473
    
2474 de6d9b64 Fabrice Bellard
    for(;;) {
2475 a579db0c Ivan Kalvachev
#ifdef HAVE_XVMC
2476
        //one 1 we memcpy blocks in xvmcvideo
2477
        if(s->avctx->xvmc_acceleration > 1)
2478
            XVMC_init_block(s);//set s->block
2479
#endif
2480
2481 eb4b3dd3 Zdenek Kabelac
        s->dsp.clear_blocks(s->block[0]);
2482 ffdff4d7 Ivan Kalvachev
        if(!s->chroma_y_shift){
2483
            s->dsp.clear_blocks(s->block[6]);
2484
        }
2485 b7ec19d3 Fabrice Bellard
        ret = mpeg_decode_mb(s, s->block);
2486 332f9ac4 Michael Niedermayer
        s->chroma_qscale= s->qscale;
2487 7bc9090a Michael Niedermayer
2488 de6d9b64 Fabrice Bellard
        dprintf("ret=%d\n", ret);
2489
        if (ret < 0)
2490
            return -1;
2491 8d7ec294 Wolfgang Hesseler
2492 31b1ec5d Michael Niedermayer
        if(s->current_picture.motion_val[0] && !s->encoding){ //note motion_val is normally NULL unless we want to extract the MVs
2493 137c8468 Michael Niedermayer
            const int wrap = field_pic ? 2*s->b8_stride : s->b8_stride;
2494
            int xy = s->mb_x*2 + s->mb_y*2*wrap;
2495 4c9544d0 Michael Niedermayer
            int motion_x, motion_y, dir, i;
2496 31b1ec5d Michael Niedermayer
            if(field_pic && !s->first_field)
2497
                xy += wrap/2;
2498
2499 50c93f74 Michael Niedermayer
            for(i=0; i<2; i++){
2500
                for(dir=0; dir<2; dir++){
2501 4c9544d0 Michael Niedermayer
                    if (s->mb_intra || (dir==1 && s->pict_type != B_TYPE)) {
2502
                        motion_x = motion_y = 0;
2503 0c352cad Michael Niedermayer
                    }else if (s->mv_type == MV_TYPE_16X16 || (s->mv_type == MV_TYPE_FIELD && field_pic)){
2504 4c9544d0 Michael Niedermayer
                        motion_x = s->mv[dir][0][0];
2505
                        motion_y = s->mv[dir][0][1];
2506
                    } else /*if ((s->mv_type == MV_TYPE_FIELD) || (s->mv_type == MV_TYPE_16X8))*/ {
2507
                        motion_x = s->mv[dir][i][0];
2508
                        motion_y = s->mv[dir][i][1];
2509
                    }
2510 7c4f71c4 Michael Niedermayer
2511 4c9544d0 Michael Niedermayer
                    s->current_picture.motion_val[dir][xy    ][0] = motion_x;
2512
                    s->current_picture.motion_val[dir][xy    ][1] = motion_y;
2513
                    s->current_picture.motion_val[dir][xy + 1][0] = motion_x;
2514
                    s->current_picture.motion_val[dir][xy + 1][1] = motion_y;
2515 7c4f71c4 Michael Niedermayer
                    s->current_picture.ref_index [dir][xy    ]=
2516
                    s->current_picture.ref_index [dir][xy + 1]= s->field_select[dir][i];
2517 4c9544d0 Michael Niedermayer
                }
2518 50c93f74 Michael Niedermayer
                xy += wrap;
2519 7bc9090a Michael Niedermayer
            }
2520
        }
2521 8d7ec294 Wolfgang Hesseler
2522 178fcca8 Michael Niedermayer
        s->dest[0] += 16 >> lowres;
2523
        s->dest[1] += 16 >> (s->chroma_x_shift + lowres);
2524
        s->dest[2] += 16 >> (s->chroma_x_shift + lowres);
2525 ce5b7c5e Michael Niedermayer
2526 7d1c3fc1 Michael Niedermayer
        MPV_decode_mb(s, s->block);
2527
        
2528 ce5b7c5e Michael Niedermayer
        if (++s->mb_x >= s->mb_width) {
2529 ff862be5 Michael Niedermayer
2530 3bb07d61 Michael Niedermayer
            ff_draw_horiz_band(s, 16*s->mb_y, 16);
2531 ce5b7c5e Michael Niedermayer
2532
            s->mb_x = 0;
2533
            s->mb_y++;
2534 cf713bb8 Michael Niedermayer
2535
            if(s->mb_y<<field_pic >= s->mb_height){
2536
                int left= s->gb.size_in_bits - get_bits_count(&s->gb);
2537
2538
                if(left < 0 || (left && show_bits(&s->gb, FFMIN(left, 23)))
2539
                   || (avctx->error_resilience >= FF_ER_AGGRESSIVE && left>8)){
2540 9b879566 Michel Bardiaux
                    av_log(avctx, AV_LOG_ERROR, "end mismatch left=%d\n", left);
2541 cf713bb8 Michael Niedermayer
                    return -1;
2542
                }else
2543
                    goto eos;
2544
            }
2545 7d1c3fc1 Michael Niedermayer
            
2546
            ff_init_block_index(s);
2547 ce5b7c5e Michael Niedermayer
        }
2548
2549
        /* skip mb handling */
2550 9b8709d1 Michael Niedermayer
        if (s->mb_skip_run == -1) {
2551 ce5b7c5e Michael Niedermayer
            /* read again increment */
2552 9b8709d1 Michael Niedermayer
            s->mb_skip_run = 0;
2553 ce5b7c5e Michael Niedermayer
            for(;;) {
2554
                int code = get_vlc2(&s->gb, mbincr_vlc.table, MBINCR_VLC_BITS, 2);
2555 7bc9090a Michael Niedermayer
                if (code < 0){
2556 9b879566 Michel Bardiaux
                    av_log(s->avctx, AV_LOG_ERROR, "mb incr damaged\n");
2557 9c00c3af Michael Niedermayer
                    return -1;
2558 7bc9090a Michael Niedermayer
                }
2559 ce5b7c5e Michael Niedermayer
                if (code >= 33) {
2560
                    if (code == 33) {
2561 9b8709d1 Michael Niedermayer
                        s->mb_skip_run += 33;
2562 9c00c3af Michael Niedermayer
                    }else if(code == 35){
2563
                        if(s->mb_skip_run != 0 || show_bits(&s->gb, 15) != 0){
2564 9b879566 Michel Bardiaux
                            av_log(s->avctx, AV_LOG_ERROR, "slice mismatch\n");
2565 9c00c3af Michael Niedermayer
                            return -1;
2566
                        }
2567
                        goto eos; /* end of slice */
2568 ce5b7c5e Michael Niedermayer
                    }
2569
                    /* otherwise, stuffing, nothing to do */
2570
                } else {
2571 9b8709d1 Michael Niedermayer
                    s->mb_skip_run += code;
2572 ce5b7c5e Michael Niedermayer
                    break;
2573
                }
2574
            }
2575
        }
2576 de6d9b64 Fabrice Bellard
    }
2577 db6e7795 Michael Niedermayer
eos: // end of slice
2578
    *buf += get_bits_count(&s->gb)/8 - 1;
2579 7bc9090a Michael Niedermayer
//printf("y %d %d %d %d\n", s->resync_mb_x, s->resync_mb_y, s->mb_x, s->mb_y);
2580 db6e7795 Michael Niedermayer
    return 0;
2581
}
2582 7bc9090a Michael Niedermayer
2583 c62c07d3 Michael Niedermayer
static int slice_decode_thread(AVCodecContext *c, void *arg){
2584
    MpegEncContext *s= arg;
2585 49a37411 Zdenek Kabelac
    const uint8_t *buf= s->gb.buffer;
2586 c62c07d3 Michael Niedermayer
    int mb_y= s->start_mb_y;
2587
2588
    s->error_count= 3*(s->end_mb_y - s->start_mb_y)*s->mb_width;
2589
2590
    for(;;){
2591
        int start_code, ret;
2592
2593
        ret= mpeg_decode_slice((Mpeg1Context*)s, mb_y, &buf, s->gb.buffer_end - buf);
2594
        emms_c();
2595
//av_log(c, AV_LOG_DEBUG, "ret:%d resync:%d/%d mb:%d/%d ts:%d/%d ec:%d\n", 
2596
//ret, s->resync_mb_x, s->resync_mb_y, s->mb_x, s->mb_y, s->start_mb_y, s->end_mb_y, s->error_count);
2597
        if(ret < 0){
2598
            if(s->resync_mb_x>=0 && s->resync_mb_y>=0)
2599
                ff_er_add_slice(s, s->resync_mb_x, s->resync_mb_y, s->mb_x, s->mb_y, AC_ERROR|DC_ERROR|MV_ERROR);
2600
        }else{
2601
            ff_er_add_slice(s, s->resync_mb_x, s->resync_mb_y, s->mb_x-1, s->mb_y, AC_END|DC_END|MV_END);
2602
        }
2603
        
2604
        if(s->mb_y == s->end_mb_y)
2605
            return 0;
2606
        
2607
        start_code = find_start_code(&buf, s->gb.buffer_end);
2608
        mb_y= start_code - SLICE_MIN_START_CODE;
2609
        if(mb_y < 0 || mb_y >= s->end_mb_y)
2610
            return -1;
2611
    }
2612
    
2613
    return 0; //not reached
2614
}
2615
2616 db6e7795 Michael Niedermayer
/**
2617
 * handles slice ends.
2618
 * @return 1 if it seems to be the last slice of 
2619
 */
2620
static int slice_end(AVCodecContext *avctx, AVFrame *pict)
2621
{
2622
    Mpeg1Context *s1 = avctx->priv_data;
2623
    MpegEncContext *s = &s1->mpeg_enc_ctx;
2624
       
2625 36b58e85 Roman Shaposhnik
    if (!s1->mpeg_enc_ctx_allocated || !s->current_picture_ptr)
2626 855127bc Fabrice Bellard
        return 0;
2627
2628 2e7b4c84 Ivan Kalvachev
#ifdef HAVE_XVMC
2629
    if(s->avctx->xvmc_acceleration)
2630
        XVMC_field_end(s);
2631
#endif
2632 de6d9b64 Fabrice Bellard
    /* end of slice reached */
2633 db6e7795 Michael Niedermayer
    if (/*s->mb_y<<field_pic == s->mb_height &&*/ !s->first_field) {
2634 de6d9b64 Fabrice Bellard
        /* end of image */
2635 1e491e29 Michael Niedermayer
2636 51929fd3 Michael Niedermayer
        s->current_picture_ptr->qscale_type= FF_QSCALE_TYPE_MPEG2;
2637 0426af31 Michael Niedermayer
2638 7bc9090a Michael Niedermayer
        ff_er_frame_end(s);
2639 de6d9b64 Fabrice Bellard
2640
        MPV_frame_end(s);
2641
2642 4e00e76b Michael Niedermayer
        if (s->pict_type == B_TYPE || s->low_delay) {
2643 8e937a4a Michael Niedermayer
            *pict= *(AVFrame*)s->current_picture_ptr;
2644 0c9bbaec Wolfgang Hesseler
            ff_print_debug_info(s, pict);
2645 de6d9b64 Fabrice Bellard
        } else {
2646 1e491e29 Michael Niedermayer
            s->picture_number++;
2647 de6d9b64 Fabrice Bellard
            /* latency of 1 frame for I and P frames */
2648
            /* XXX: use another variable than picture_number */
2649 80097bbf Michael Niedermayer
            if (s->last_picture_ptr != NULL) {
2650 8e937a4a Michael Niedermayer
                *pict= *(AVFrame*)s->last_picture_ptr;
2651 0c9bbaec Wolfgang Hesseler
                 ff_print_debug_info(s, pict);
2652 de6d9b64 Fabrice Bellard
            }
2653
        }
2654 ff862be5 Michael Niedermayer
2655 db6e7795 Michael Niedermayer
        return 1;
2656 de6d9b64 Fabrice Bellard
    } else {
2657 db6e7795 Michael Niedermayer
        return 0;
2658 de6d9b64 Fabrice Bellard
    }
2659
}
2660
2661
static int mpeg1_decode_sequence(AVCodecContext *avctx, 
2662 5f0f7713 Zdenek Kabelac
                                 const uint8_t *buf, int buf_size)
2663 de6d9b64 Fabrice Bellard
{
2664
    Mpeg1Context *s1 = avctx->priv_data;
2665
    MpegEncContext *s = &s1->mpeg_enc_ctx;
2666 5ac47683 Ivan Kalvachev
    int width,height;
2667 a6b9ffbf Ivan Kalvachev
    int i, v, j;
2668 fb4a4a56 Fabrice Bellard
2669 68f593b4 Michael Niedermayer
    init_get_bits(&s->gb, buf, buf_size*8);
2670 de6d9b64 Fabrice Bellard
2671 5ac47683 Ivan Kalvachev
    width = get_bits(&s->gb, 12);
2672
    height = get_bits(&s->gb, 12);
2673
    if (width <= 0 || height <= 0 ||
2674
        (width % 2) != 0 || (height % 2) != 0)
2675
        return -1;
2676 945f15b7 Michael Niedermayer
    s->aspect_ratio_info= get_bits(&s->gb, 4);
2677 32e7b91a Fabrice Bellard
    if (s->aspect_ratio_info == 0)
2678
        return -1;
2679 de6d9b64 Fabrice Bellard
    s->frame_rate_index = get_bits(&s->gb, 4);
2680 8d52ec7e Michael Niedermayer
    if (s->frame_rate_index == 0 || s->frame_rate_index > 13)
2681 de6d9b64 Fabrice Bellard
        return -1;
2682
    s->bit_rate = get_bits(&s->gb, 18) * 400;
2683 612476ef Arpi
    if (get_bits1(&s->gb) == 0) /* marker */
2684 de6d9b64 Fabrice Bellard
        return -1;
2685 5ac47683 Ivan Kalvachev
    s->width = width;
2686
    s->height = height;
2687 de6d9b64 Fabrice Bellard
2688 0368c72d Michael Niedermayer
    s->avctx->rc_buffer_size= get_bits(&s->gb, 10) * 1024*16;
2689 612476ef Arpi
    skip_bits(&s->gb, 1);
2690 de6d9b64 Fabrice Bellard
2691
    /* get matrix */
2692 612476ef Arpi
    if (get_bits1(&s->gb)) {
2693 de6d9b64 Fabrice Bellard
        for(i=0;i<64;i++) {
2694
            v = get_bits(&s->gb, 8);
2695 e6dc9c6f Michael Niedermayer
            if(v==0){
2696
                av_log(s->avctx, AV_LOG_ERROR, "intra matrix damaged\n");
2697
                return -1;
2698
            }
2699 a6b9ffbf Ivan Kalvachev
            j = s->dsp.idct_permutation[ ff_zigzag_direct[i] ];
2700 60832448 Fabrice Bellard
            s->intra_matrix[j] = v;
2701
            s->chroma_intra_matrix[j] = v;
2702 de6d9b64 Fabrice Bellard
        }
2703 25ed7f92 Fabrice Bellard
#ifdef DEBUG
2704
        dprintf("intra matrix present\n");
2705
        for(i=0;i<64;i++)
2706 a6b9ffbf Ivan Kalvachev
            dprintf(" %d", s->intra_matrix[s->dsp.idct_permutation[i]);
2707 25ed7f92 Fabrice Bellard
        printf("\n");
2708
#endif
2709 de6d9b64 Fabrice Bellard
    } else {
2710
        for(i=0;i<64;i++) {
2711 a6b9ffbf Ivan Kalvachev
            j = s->dsp.idct_permutation[i];
2712 adc09b2e Marko Kreen
            v = ff_mpeg1_default_intra_matrix[i];
2713 2ad1516a Michael Niedermayer
            s->intra_matrix[j] = v;
2714
            s->chroma_intra_matrix[j] = v;
2715 de6d9b64 Fabrice Bellard
        }
2716
    }
2717 612476ef Arpi
    if (get_bits1(&s->gb)) {
2718 de6d9b64 Fabrice Bellard
        for(i=0;i<64;i++) {
2719
            v = get_bits(&s->gb, 8);
2720 e6dc9c6f Michael Niedermayer
            if(v==0){
2721
                av_log(s->avctx, AV_LOG_ERROR, "inter matrix damaged\n");
2722
                return -1;
2723
            }
2724 a6b9ffbf Ivan Kalvachev
            j = s->dsp.idct_permutation[ ff_zigzag_direct[i] ];
2725 d7e9533a Michael Niedermayer
            s->inter_matrix[j] = v;
2726
            s->chroma_inter_matrix[j] = v;
2727 de6d9b64 Fabrice Bellard
        }
2728 25ed7f92 Fabrice Bellard
#ifdef DEBUG
2729
        dprintf("non intra matrix present\n");
2730
        for(i=0;i<64;i++)
2731 a6b9ffbf Ivan Kalvachev
            dprintf(" %d", s->inter_matrix[s->dsp.idct_permutation[i]);
2732 25ed7f92 Fabrice Bellard
        printf("\n");
2733
#endif
2734 de6d9b64 Fabrice Bellard
    } else {
2735
        for(i=0;i<64;i++) {
2736 b0368839 Michael Niedermayer
            int j= s->dsp.idct_permutation[i];
2737 adc09b2e Marko Kreen
            v = ff_mpeg1_default_non_intra_matrix[i];
2738 2ad1516a Michael Niedermayer
            s->inter_matrix[j] = v;
2739
            s->chroma_inter_matrix[j] = v;
2740 de6d9b64 Fabrice Bellard
        }
2741
    }
2742 e6dc9c6f Michael Niedermayer
    
2743
    if(show_bits(&s->gb, 23) != 0){
2744
        av_log(s->avctx, AV_LOG_ERROR, "sequence header damaged\n");
2745
        return -1;
2746
    }
2747 de6d9b64 Fabrice Bellard
2748
    /* we set mpeg2 parameters so that it emulates mpeg1 */
2749
    s->progressive_sequence = 1;
2750
    s->progressive_frame = 1;
2751
    s->picture_structure = PICT_FRAME;
2752
    s->frame_pred_frame_dct = 1;
2753 4eaad9c0 Ivan Kalvachev
    s->chroma_format = 1;
2754 029911d1 Michael Niedermayer
    s->codec_id= s->avctx->codec_id= CODEC_ID_MPEG1VIDEO;
2755 fb4a4a56 Fabrice Bellard
    avctx->sub_id = 1; /* indicates mpeg1 */
2756 a6b9ffbf Ivan Kalvachev
    s->out_format = FMT_MPEG1;
2757
    s->swap_uv = 0;//AFAIK VCR2 don't have SEQ_HEADER
2758 248a189a Fabrice Bellard
    if(s->flags & CODEC_FLAG_LOW_DELAY) s->low_delay=1;
2759 0368c72d Michael Niedermayer
    
2760<