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