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