seems wmv8 doesnt use intra-inter prediction
[ffmpeg.git] / libavcodec / wmv2.c
1 /*
2  * Copyright (c) 2002 The FFmpeg Project.
3  *
4  * This library is free software; you can redistribute it and/or
5  * modify it under the terms of the GNU Lesser General Public
6  * License as published by the Free Software Foundation; either
7  * version 2 of the License, or (at your option) any later version.
8  *
9  * This library is distributed in the hope that it will be useful,
10  * but WITHOUT ANY WARRANTY; without even the implied warranty of
11  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
12  * Lesser General Public License for more details.
13  *
14  * You should have received a copy of the GNU Lesser General Public
15  * License along with this library; if not, write to the Free Software
16  * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
17  *
18  */
19
20 /**
21  * @file wmv2.c
22  * wmv2 codec.
23  */
24  
25 #include "simple_idct.h"
26  
27 #define SKIP_TYPE_NONE 0
28 #define SKIP_TYPE_MPEG 1
29 #define SKIP_TYPE_ROW  2
30 #define SKIP_TYPE_COL  3
31
32
33 typedef struct Wmv2Context{
34     MpegEncContext s;
35     int j_type_bit;
36     int j_type;
37     int flag3;
38     int flag63;
39     int abt_flag;
40     int abt_type;
41     int abt_type_table[6];
42     int per_mb_abt;
43     int per_block_abt;
44     int mspel_bit;
45     int cbp_table_index;
46     int top_left_mv_flag;
47     int per_mb_rl_bit;
48     int skip_type;
49     int hshift;
50     
51     ScanTable abt_scantable[2];
52     DCTELEM abt_block2[6][64] __align8;
53 }Wmv2Context;
54
55 static void wmv2_common_init(Wmv2Context * w){
56     MpegEncContext * const s= &w->s;
57         
58     ff_init_scantable(s->dsp.idct_permutation, &w->abt_scantable[0], wmv2_scantableA);
59     ff_init_scantable(s->dsp.idct_permutation, &w->abt_scantable[1], wmv2_scantableB);
60 }
61
62 static int encode_ext_header(Wmv2Context *w){
63     MpegEncContext * const s= &w->s;
64     PutBitContext pb;
65     int code;
66         
67     init_put_bits(&pb, s->avctx->extradata, s->avctx->extradata_size);
68
69     put_bits(&pb, 5, s->avctx->frame_rate / s->avctx->frame_rate_base); //yes 29.97 -> 29
70     put_bits(&pb, 11, FFMIN(s->bit_rate/1024, 2047));
71     
72     put_bits(&pb, 1, w->mspel_bit=1);
73     put_bits(&pb, 1, w->flag3=1);
74     put_bits(&pb, 1, w->abt_flag=1);
75     put_bits(&pb, 1, w->j_type_bit=1);
76     put_bits(&pb, 1, w->top_left_mv_flag=0);
77     put_bits(&pb, 1, w->per_mb_rl_bit=1);
78     put_bits(&pb, 3, code=1);
79     
80     flush_put_bits(&pb);
81
82     s->slice_height = s->mb_height / code;
83     
84     return 0;
85 }
86
87 #ifdef CONFIG_ENCODERS
88 static int wmv2_encode_init(AVCodecContext *avctx){
89     Wmv2Context * const w= avctx->priv_data;
90     
91     if(MPV_encode_init(avctx) < 0)
92         return -1;
93     
94     wmv2_common_init(w);
95
96     avctx->extradata_size= 4;
97     avctx->extradata= av_mallocz(avctx->extradata_size + 10);
98     encode_ext_header(w);
99     
100     return 0;
101 }
102
103 static int wmv2_encode_end(AVCodecContext *avctx){
104     
105     if(MPV_encode_end(avctx) < 0)
106         return -1;
107     
108     avctx->extradata_size= 0;
109     av_freep(&avctx->extradata);
110     
111     return 0;
112 }
113
114 int ff_wmv2_encode_picture_header(MpegEncContext * s, int picture_number)
115 {
116     Wmv2Context * const w= (Wmv2Context*)s;
117
118     put_bits(&s->pb, 1, s->pict_type - 1);
119     if(s->pict_type == I_TYPE){
120         put_bits(&s->pb, 7, 0);
121     }
122     put_bits(&s->pb, 5, s->qscale);
123
124     s->dc_table_index = 1;
125     s->mv_table_index = 1; /* only if P frame */
126 //    s->use_skip_mb_code = 1; /* only if P frame */
127     s->per_mb_rl_table = 0;
128     s->mspel= 0;
129     w->per_mb_abt=0;
130     w->abt_type=0;
131     w->j_type=0;
132
133     assert(s->flipflop_rounding);
134
135     if (s->pict_type == I_TYPE) {
136         assert(s->no_rounding==1);
137         if(w->j_type_bit) put_bits(&s->pb, 1, w->j_type);
138         
139         if(w->per_mb_rl_bit) put_bits(&s->pb, 1, s->per_mb_rl_table);
140         
141         if(!s->per_mb_rl_table){
142             code012(&s->pb, s->rl_chroma_table_index);
143             code012(&s->pb, s->rl_table_index);
144         }
145
146         put_bits(&s->pb, 1, s->dc_table_index);
147
148         s->inter_intra_pred= 0;
149     }else{
150         int cbp_index;
151
152         put_bits(&s->pb, 2, SKIP_TYPE_NONE);
153         
154         code012(&s->pb, cbp_index=0);
155         if(s->qscale <= 10){
156             int map[3]= {0,2,1};
157             w->cbp_table_index= map[cbp_index];
158         }else if(s->qscale <= 20){
159             int map[3]= {1,0,2};
160             w->cbp_table_index= map[cbp_index];
161         }else{
162             int map[3]= {2,1,0};
163             w->cbp_table_index= map[cbp_index];
164         }
165
166         if(w->mspel_bit) put_bits(&s->pb, 1, s->mspel);
167     
168         if(w->abt_flag){
169             put_bits(&s->pb, 1, w->per_mb_abt^1);
170             if(!w->per_mb_abt){
171                 code012(&s->pb, w->abt_type);
172             }
173         }
174
175         if(w->per_mb_rl_bit) put_bits(&s->pb, 1, s->per_mb_rl_table);
176         
177         if(!s->per_mb_rl_table){
178             code012(&s->pb, s->rl_table_index);
179             s->rl_chroma_table_index = s->rl_table_index;
180         }
181         put_bits(&s->pb, 1, s->dc_table_index);
182         put_bits(&s->pb, 1, s->mv_table_index);
183     
184         s->inter_intra_pred= 0;//(s->width*s->height < 320*240 && s->bit_rate<=II_BITRATE);
185     }
186     s->esc3_level_length= 0;
187     s->esc3_run_length= 0;
188
189     return 0;
190 }
191
192 // nearly idential to wmv1 but thats just because we dont use the useless M$ crap features
193 // its duplicated here in case someone wants to add support for these carp features
194 void ff_wmv2_encode_mb(MpegEncContext * s, 
195                        DCTELEM block[6][64],
196                        int motion_x, int motion_y)
197 {
198     Wmv2Context * const w= (Wmv2Context*)s;
199     int cbp, coded_cbp, i;
200     int pred_x, pred_y;
201     uint8_t *coded_block;
202
203     handle_slices(s);
204     
205     if (!s->mb_intra) {
206         /* compute cbp */
207         set_stat(ST_INTER_MB);
208         cbp = 0;
209         for (i = 0; i < 6; i++) {
210             if (s->block_last_index[i] >= 0)
211                 cbp |= 1 << (5 - i);
212         }
213         
214         put_bits(&s->pb, 
215                  wmv2_inter_table[w->cbp_table_index][cbp + 64][1], 
216                  wmv2_inter_table[w->cbp_table_index][cbp + 64][0]);
217
218         /* motion vector */
219         h263_pred_motion(s, 0, 0, &pred_x, &pred_y);
220         msmpeg4_encode_motion(s, motion_x - pred_x, 
221                               motion_y - pred_y);
222     } else {
223         /* compute cbp */
224         cbp = 0;
225         coded_cbp = 0;
226         for (i = 0; i < 6; i++) {
227             int val, pred;
228             val = (s->block_last_index[i] >= 1);
229             cbp |= val << (5 - i);
230             if (i < 4) {
231                 /* predict value for close blocks only for luma */
232                 pred = coded_block_pred(s, i, &coded_block);
233                 *coded_block = val;
234                 val = val ^ pred;
235             }
236             coded_cbp |= val << (5 - i);
237         }
238 #if 0
239         if (coded_cbp)
240             printf("cbp=%x %x\n", cbp, coded_cbp);
241 #endif
242
243         if (s->pict_type == I_TYPE) {
244             set_stat(ST_INTRA_MB);
245             put_bits(&s->pb, 
246                      table_mb_intra[coded_cbp][1], table_mb_intra[coded_cbp][0]);
247         } else {
248             put_bits(&s->pb, 
249                      wmv2_inter_table[w->cbp_table_index][cbp][1], 
250                      wmv2_inter_table[w->cbp_table_index][cbp][0]);
251         }
252         set_stat(ST_INTRA_MB);
253         put_bits(&s->pb, 1, 0); /* no AC prediction yet */
254         if(s->inter_intra_pred){
255             s->h263_aic_dir=0;
256             put_bits(&s->pb, table_inter_intra[s->h263_aic_dir][1], table_inter_intra[s->h263_aic_dir][0]);
257         }
258     }
259
260     for (i = 0; i < 6; i++) {
261         msmpeg4_encode_block(s, block[i], i);
262     }
263 }
264 #endif //CONFIG_ENCODERS
265
266 static void parse_mb_skip(Wmv2Context * w){
267     int mb_x, mb_y;
268     MpegEncContext * const s= &w->s;
269     uint32_t * const mb_type= s->current_picture_ptr->mb_type;
270
271     w->skip_type= get_bits(&s->gb, 2);
272     switch(w->skip_type){
273     case SKIP_TYPE_NONE:
274         for(mb_y=0; mb_y<s->mb_height; mb_y++){
275             for(mb_x=0; mb_x<s->mb_width; mb_x++){
276                 mb_type[mb_y*s->mb_stride + mb_x]= MB_TYPE_16x16 | MB_TYPE_L0;
277             }
278         }
279         break;
280     case SKIP_TYPE_MPEG:
281         for(mb_y=0; mb_y<s->mb_height; mb_y++){
282             for(mb_x=0; mb_x<s->mb_width; mb_x++){
283                 mb_type[mb_y*s->mb_stride + mb_x]= (get_bits1(&s->gb) ? MB_TYPE_SKIP : 0) | MB_TYPE_16x16 | MB_TYPE_L0;
284             }
285         }
286         break;
287     case SKIP_TYPE_ROW:
288         for(mb_y=0; mb_y<s->mb_height; mb_y++){
289             if(get_bits1(&s->gb)){
290                 for(mb_x=0; mb_x<s->mb_width; mb_x++){
291                     mb_type[mb_y*s->mb_stride + mb_x]=  MB_TYPE_SKIP | MB_TYPE_16x16 | MB_TYPE_L0;
292                 }
293             }else{
294                 for(mb_x=0; mb_x<s->mb_width; mb_x++){
295                     mb_type[mb_y*s->mb_stride + mb_x]= (get_bits1(&s->gb) ? MB_TYPE_SKIP : 0) | MB_TYPE_16x16 | MB_TYPE_L0;
296                 }
297             }
298         }
299         break;
300     case SKIP_TYPE_COL:
301         for(mb_x=0; mb_x<s->mb_width; mb_x++){
302             if(get_bits1(&s->gb)){
303                 for(mb_y=0; mb_y<s->mb_height; mb_y++){
304                     mb_type[mb_y*s->mb_stride + mb_x]=  MB_TYPE_SKIP | MB_TYPE_16x16 | MB_TYPE_L0;
305                 }
306             }else{
307                 for(mb_y=0; mb_y<s->mb_height; mb_y++){
308                     mb_type[mb_y*s->mb_stride + mb_x]= (get_bits1(&s->gb) ? MB_TYPE_SKIP : 0) | MB_TYPE_16x16 | MB_TYPE_L0;
309                 }
310             }
311         }
312         break;
313     }
314 }
315
316 static int decode_ext_header(Wmv2Context *w){
317     MpegEncContext * const s= &w->s;
318     GetBitContext gb;
319     int fps;
320     int code;
321
322     if(s->avctx->extradata_size<4) return -1;
323     
324     init_get_bits(&gb, s->avctx->extradata, s->avctx->extradata_size*8);
325
326     fps                = get_bits(&gb, 5);
327     s->bit_rate        = get_bits(&gb, 11)*1024;
328     w->mspel_bit       = get_bits1(&gb);
329     w->flag3           = get_bits1(&gb);
330     w->abt_flag        = get_bits1(&gb);
331     w->j_type_bit      = get_bits1(&gb);
332     w->top_left_mv_flag= get_bits1(&gb);
333     w->per_mb_rl_bit   = get_bits1(&gb);
334     code               = get_bits(&gb, 3);
335     
336     if(code==0) return -1;
337
338     s->slice_height = s->mb_height / code;
339
340     if(s->avctx->debug&FF_DEBUG_PICT_INFO){
341         av_log(s->avctx, AV_LOG_DEBUG, "fps:%d, br:%d, qpbit:%d, abt_flag:%d, j_type_bit:%d, tl_mv_flag:%d, mbrl_bit:%d, code:%d, flag3:%d, slices:%d\n", 
342         fps, s->bit_rate, w->mspel_bit, w->abt_flag, w->j_type_bit, w->top_left_mv_flag, w->per_mb_rl_bit, code, w->flag3, 
343         code);
344     }
345     return 0;
346 }
347
348 int ff_wmv2_decode_picture_header(MpegEncContext * s)
349 {
350     Wmv2Context * const w= (Wmv2Context*)s;
351     int code;
352
353 #if 0
354 {
355 int i;
356 for(i=0; i<s->gb.size*8; i++)
357     printf("%d", get_bits1(&s->gb));
358 //    get_bits1(&s->gb);
359 printf("END\n");
360 return -1;
361 }
362 #endif
363     if(s->picture_number==0)
364         decode_ext_header(w);
365
366     s->pict_type = get_bits(&s->gb, 1) + 1;
367     if(s->pict_type == I_TYPE){
368         code = get_bits(&s->gb, 7);
369         av_log(s->avctx, AV_LOG_ERROR, "I7:%X/\n", code);
370     }
371     s->chroma_qscale= s->qscale = get_bits(&s->gb, 5);
372     if(s->qscale < 0)
373        return -1;
374        
375     return 0;
376 }
377
378 int ff_wmv2_decode_secondary_picture_header(MpegEncContext * s)
379 {
380     Wmv2Context * const w= (Wmv2Context*)s;
381
382     if (s->pict_type == I_TYPE) {
383         if(w->j_type_bit) w->j_type= get_bits1(&s->gb);
384         else              w->j_type= 0; //FIXME check
385         
386         if(!w->j_type){
387             if(w->per_mb_rl_bit) s->per_mb_rl_table= get_bits1(&s->gb);
388             else                 s->per_mb_rl_table= 0;
389         
390             if(!s->per_mb_rl_table){
391                 s->rl_chroma_table_index = decode012(&s->gb);
392                 s->rl_table_index = decode012(&s->gb);
393             }
394
395             s->dc_table_index = get_bits1(&s->gb);
396         }
397         s->inter_intra_pred= 0;
398         s->no_rounding = 1;
399         if(s->avctx->debug&FF_DEBUG_PICT_INFO){
400             av_log(s->avctx, AV_LOG_DEBUG, "qscale:%d rlc:%d rl:%d dc:%d mbrl:%d j_type:%d \n", 
401                 s->qscale,
402                 s->rl_chroma_table_index,
403                 s->rl_table_index, 
404                 s->dc_table_index,
405                 s->per_mb_rl_table,
406                 w->j_type);
407         }
408     }else{
409         int cbp_index;
410         w->j_type=0;
411
412         parse_mb_skip(w);
413         cbp_index= decode012(&s->gb);
414         if(s->qscale <= 10){
415             int map[3]= {0,2,1};
416             w->cbp_table_index= map[cbp_index];
417         }else if(s->qscale <= 20){
418             int map[3]= {1,0,2};
419             w->cbp_table_index= map[cbp_index];
420         }else{
421             int map[3]= {2,1,0};
422             w->cbp_table_index= map[cbp_index];
423         }
424
425         if(w->mspel_bit) s->mspel= get_bits1(&s->gb);
426         else             s->mspel= 0; //FIXME check
427     
428         if(w->abt_flag){
429             w->per_mb_abt= get_bits1(&s->gb)^1;
430             if(!w->per_mb_abt){
431                 w->abt_type= decode012(&s->gb);
432             }
433         }
434
435         if(w->per_mb_rl_bit) s->per_mb_rl_table= get_bits1(&s->gb);
436         else                 s->per_mb_rl_table= 0;
437         
438         if(!s->per_mb_rl_table){
439             s->rl_table_index = decode012(&s->gb);
440             s->rl_chroma_table_index = s->rl_table_index;
441         }
442
443         s->dc_table_index = get_bits1(&s->gb);
444         s->mv_table_index = get_bits1(&s->gb);
445     
446         s->inter_intra_pred= 0;//(s->width*s->height < 320*240 && s->bit_rate<=II_BITRATE);
447         s->no_rounding ^= 1;
448         
449         if(s->avctx->debug&FF_DEBUG_PICT_INFO){
450             av_log(s->avctx, AV_LOG_DEBUG, "rl:%d rlc:%d dc:%d mv:%d mbrl:%d qp:%d mspel:%d per_mb_abt:%d abt_type:%d cbp:%d ii:%d\n", 
451                 s->rl_table_index, 
452                 s->rl_chroma_table_index, 
453                 s->dc_table_index,
454                 s->mv_table_index,
455                 s->per_mb_rl_table,
456                 s->qscale,
457                 s->mspel,
458                 w->per_mb_abt,
459                 w->abt_type,
460                 w->cbp_table_index,
461                 s->inter_intra_pred);
462         }
463     }
464     s->esc3_level_length= 0;
465     s->esc3_run_length= 0;
466     
467 s->picture_number++; //FIXME ?
468
469
470 //    if(w->j_type)
471 //        return wmv2_decode_j_picture(w); //FIXME
472
473     if(w->j_type){
474         av_log(s->avctx, AV_LOG_ERROR, "J-type picture isnt supported\n");
475         return -1;
476     }
477
478     return 0;
479 }
480
481 static void ff_wmv2_decode_init(MpegEncContext *s){
482 }
483
484 static inline int wmv2_decode_motion(Wmv2Context *w, int *mx_ptr, int *my_ptr){
485     MpegEncContext * const s= &w->s;
486     int ret;
487    
488     ret= msmpeg4_decode_motion(s, mx_ptr, my_ptr);
489    
490     if(ret<0) return -1;
491    
492     if((((*mx_ptr)|(*my_ptr)) & 1) && s->mspel)
493         w->hshift= get_bits1(&s->gb);
494     else 
495         w->hshift= 0;
496
497 //printf("%d %d  ", *mx_ptr, *my_ptr);
498    
499     return 0;
500 }
501
502 static int16_t *wmv2_pred_motion(Wmv2Context *w, int *px, int *py){
503     MpegEncContext * const s= &w->s;
504     int xy, wrap, diff, type;
505     int16_t *A, *B, *C, *mot_val;
506
507     wrap = s->b8_stride;
508     xy = s->block_index[0];
509
510     mot_val = s->current_picture.motion_val[0][xy];
511
512     A = s->current_picture.motion_val[0][xy - 1];
513     B = s->current_picture.motion_val[0][xy - wrap];
514     C = s->current_picture.motion_val[0][xy + 2 - wrap];
515     
516     diff= FFMAX(ABS(A[0] - B[0]), ABS(A[1] - B[1]));
517     
518     if(s->mb_x && !s->first_slice_line && !s->mspel && w->top_left_mv_flag && diff >= 8)
519         type= get_bits1(&s->gb);
520     else
521         type= 2;
522     
523     if(type == 0){
524         *px= A[0];
525         *py= A[1];
526     }else if(type == 1){
527         *px= B[0];
528         *py= B[1];
529     }else{
530         /* special case for first (slice) line */
531         if (s->first_slice_line) {
532             *px = A[0];
533             *py = A[1];
534         } else {
535             *px = mid_pred(A[0], B[0], C[0]);
536             *py = mid_pred(A[1], B[1], C[1]);
537         }
538     }
539
540     return mot_val;
541 }
542
543 static inline int wmv2_decode_inter_block(Wmv2Context *w, DCTELEM *block, int n, int cbp){
544     MpegEncContext * const s= &w->s;
545     static const int sub_cbp_table[3]= {2,3,1};
546     int sub_cbp;
547
548     if(!cbp){ 
549         s->block_last_index[n] = -1;
550
551         return 0;
552     }
553     
554     if(w->per_block_abt)
555         w->abt_type= decode012(&s->gb);
556 #if 0
557     if(w->per_block_abt)
558         printf("B%d", w->abt_type);
559 #endif
560     w->abt_type_table[n]= w->abt_type;
561
562     if(w->abt_type){
563 //        const uint8_t *scantable= w->abt_scantable[w->abt_type-1].permutated;
564         const uint8_t *scantable= w->abt_scantable[w->abt_type-1].scantable;
565 //        const uint8_t *scantable= w->abt_type-1 ? w->abt_scantable[1].permutated : w->abt_scantable[0].scantable;
566
567         sub_cbp= sub_cbp_table[ decode012(&s->gb) ];
568 //        printf("S%d", sub_cbp);
569
570         if(sub_cbp&1){
571             if (msmpeg4_decode_block(s, block, n, 1, scantable) < 0)
572                 return -1;
573         }
574         
575         if(sub_cbp&2){
576             if (msmpeg4_decode_block(s, w->abt_block2[n], n, 1, scantable) < 0)
577                 return -1;
578         }
579         s->block_last_index[n] = 63;
580
581         return 0;
582     }else{
583         return msmpeg4_decode_block(s, block, n, 1, s->inter_scantable.permutated);
584     }
585 }
586
587 static void wmv2_add_block(Wmv2Context *w, DCTELEM *block1, uint8_t *dst, int stride, int n){
588     MpegEncContext * const s= &w->s;
589
590     switch(w->abt_type_table[n]){
591     case 0:
592         if (s->block_last_index[n] >= 0) {
593             s->dsp.idct_add (dst, stride, block1);
594         }
595         break;
596     case 1:
597         simple_idct84_add(dst           , stride, block1);
598         simple_idct84_add(dst + 4*stride, stride, w->abt_block2[n]);
599         memset(w->abt_block2[n], 0, 64*sizeof(DCTELEM));
600         break;
601     case 2:
602         simple_idct48_add(dst           , stride, block1);
603         simple_idct48_add(dst + 4       , stride, w->abt_block2[n]);
604         memset(w->abt_block2[n], 0, 64*sizeof(DCTELEM));
605         break;
606     default:
607         av_log(s->avctx, AV_LOG_ERROR, "internal error in WMV2 abt\n");
608     }
609 }
610
611 void ff_wmv2_add_mb(MpegEncContext *s, DCTELEM block1[6][64], uint8_t *dest_y, uint8_t *dest_cb, uint8_t *dest_cr){
612     Wmv2Context * const w= (Wmv2Context*)s;
613
614     wmv2_add_block(w, block1[0], dest_y                    , s->linesize, 0);
615     wmv2_add_block(w, block1[1], dest_y + 8                , s->linesize, 1);
616     wmv2_add_block(w, block1[2], dest_y +     8*s->linesize, s->linesize, 2);
617     wmv2_add_block(w, block1[3], dest_y + 8 + 8*s->linesize, s->linesize, 3);
618     
619     if(s->flags&CODEC_FLAG_GRAY) return;
620     
621     wmv2_add_block(w, block1[4], dest_cb                   , s->uvlinesize, 4);
622     wmv2_add_block(w, block1[5], dest_cr                   , s->uvlinesize, 5);
623 }
624
625 void ff_mspel_motion(MpegEncContext *s,
626                                uint8_t *dest_y, uint8_t *dest_cb, uint8_t *dest_cr,
627                                uint8_t **ref_picture, op_pixels_func (*pix_op)[4],
628                                int motion_x, int motion_y, int h)
629 {
630     Wmv2Context * const w= (Wmv2Context*)s;
631     uint8_t *ptr;
632     int dxy, offset, mx, my, src_x, src_y, v_edge_pos, linesize, uvlinesize;
633     int emu=0;
634     
635     dxy = ((motion_y & 1) << 1) | (motion_x & 1);
636     dxy = 2*dxy + w->hshift;
637     src_x = s->mb_x * 16 + (motion_x >> 1);
638     src_y = s->mb_y * 16 + (motion_y >> 1);
639                 
640     /* WARNING: do no forget half pels */
641     v_edge_pos = s->v_edge_pos;
642     src_x = clip(src_x, -16, s->width);
643     src_y = clip(src_y, -16, s->height);
644     linesize   = s->linesize;
645     uvlinesize = s->uvlinesize;
646     ptr = ref_picture[0] + (src_y * linesize) + src_x;
647
648     if(s->flags&CODEC_FLAG_EMU_EDGE){
649         if(src_x<1 || src_y<1 || src_x + 17  >= s->h_edge_pos
650                               || src_y + h+1 >= v_edge_pos){
651             ff_emulated_edge_mc(s->edge_emu_buffer, ptr - 1 - s->linesize, s->linesize, 19, 19, 
652                              src_x-1, src_y-1, s->h_edge_pos, s->v_edge_pos);
653             ptr= s->edge_emu_buffer + 1 + s->linesize;
654             emu=1;
655         }
656     }
657
658     s->dsp.put_mspel_pixels_tab[dxy](dest_y             , ptr             , linesize);
659     s->dsp.put_mspel_pixels_tab[dxy](dest_y+8           , ptr+8           , linesize);
660     s->dsp.put_mspel_pixels_tab[dxy](dest_y  +8*linesize, ptr  +8*linesize, linesize);
661     s->dsp.put_mspel_pixels_tab[dxy](dest_y+8+8*linesize, ptr+8+8*linesize, linesize);
662
663     if(s->flags&CODEC_FLAG_GRAY) return;
664
665     if (s->out_format == FMT_H263) {
666         dxy = 0;
667         if ((motion_x & 3) != 0)
668             dxy |= 1;
669         if ((motion_y & 3) != 0)
670             dxy |= 2;
671         mx = motion_x >> 2;
672         my = motion_y >> 2;
673     } else {
674         mx = motion_x / 2;
675         my = motion_y / 2;
676         dxy = ((my & 1) << 1) | (mx & 1);
677         mx >>= 1;
678         my >>= 1;
679     }
680     
681     src_x = s->mb_x * 8 + mx;
682     src_y = s->mb_y * 8 + my;
683     src_x = clip(src_x, -8, s->width >> 1);
684     if (src_x == (s->width >> 1))
685         dxy &= ~1;
686     src_y = clip(src_y, -8, s->height >> 1);
687     if (src_y == (s->height >> 1))
688         dxy &= ~2;
689     offset = (src_y * uvlinesize) + src_x;
690     ptr = ref_picture[1] + offset;
691     if(emu){
692         ff_emulated_edge_mc(s->edge_emu_buffer, ptr, s->uvlinesize, 9, 9, 
693                          src_x, src_y, s->h_edge_pos>>1, s->v_edge_pos>>1);
694         ptr= s->edge_emu_buffer;
695     }
696     pix_op[1][dxy](dest_cb, ptr, uvlinesize, h >> 1);
697
698     ptr = ref_picture[2] + offset;
699     if(emu){
700         ff_emulated_edge_mc(s->edge_emu_buffer, ptr, s->uvlinesize, 9, 9, 
701                          src_x, src_y, s->h_edge_pos>>1, s->v_edge_pos>>1);
702         ptr= s->edge_emu_buffer;
703     }
704     pix_op[1][dxy](dest_cr, ptr, uvlinesize, h >> 1);
705 }
706
707
708 static int wmv2_decode_mb(MpegEncContext *s, DCTELEM block[6][64])
709 {
710     Wmv2Context * const w= (Wmv2Context*)s;
711     int cbp, code, i;
712     uint8_t *coded_val;
713
714     if(w->j_type) return 0;
715     
716     if (s->pict_type == P_TYPE) {
717         if(IS_SKIP(s->current_picture.mb_type[s->mb_y * s->mb_stride + s->mb_x])){
718             /* skip mb */
719             s->mb_intra = 0;
720             for(i=0;i<6;i++)
721                 s->block_last_index[i] = -1;
722             s->mv_dir = MV_DIR_FORWARD;
723             s->mv_type = MV_TYPE_16X16;
724             s->mv[0][0][0] = 0;
725             s->mv[0][0][1] = 0;
726             s->mb_skiped = 1;
727             w->hshift=0;
728             return 0;
729         }
730
731         code = get_vlc2(&s->gb, mb_non_intra_vlc[w->cbp_table_index].table, MB_NON_INTRA_VLC_BITS, 3);
732         if (code < 0)
733             return -1;
734         s->mb_intra = (~code & 0x40) >> 6;
735             
736         cbp = code & 0x3f;
737     } else {
738         s->mb_intra = 1;
739         code = get_vlc2(&s->gb, mb_intra_vlc.table, MB_INTRA_VLC_BITS, 2);
740         if (code < 0){
741             av_log(s->avctx, AV_LOG_ERROR, "II-cbp illegal at %d %d\n", s->mb_x, s->mb_y);
742             return -1;
743         }
744         /* predict coded block pattern */
745         cbp = 0;
746         for(i=0;i<6;i++) {
747             int val = ((code >> (5 - i)) & 1);
748             if (i < 4) {
749                 int pred = coded_block_pred(s, i, &coded_val);
750                 val = val ^ pred;
751                 *coded_val = val;
752             }
753             cbp |= val << (5 - i);
754         }
755     }
756
757     if (!s->mb_intra) {
758         int mx, my;
759 //printf("P at %d %d\n", s->mb_x, s->mb_y);
760         wmv2_pred_motion(w, &mx, &my);
761         
762         if(cbp){
763             if(s->per_mb_rl_table){
764                 s->rl_table_index = decode012(&s->gb);
765                 s->rl_chroma_table_index = s->rl_table_index;
766             }
767
768             if(w->abt_flag && w->per_mb_abt){
769                 w->per_block_abt= get_bits1(&s->gb);
770                 if(!w->per_block_abt)
771                     w->abt_type= decode012(&s->gb);
772             }else
773                 w->per_block_abt=0;
774         }
775         
776         if (wmv2_decode_motion(w, &mx, &my) < 0)
777             return -1;
778
779         s->mv_dir = MV_DIR_FORWARD;
780         s->mv_type = MV_TYPE_16X16;
781         s->mv[0][0][0] = mx;
782         s->mv[0][0][1] = my;
783
784         for (i = 0; i < 6; i++) {
785             if (wmv2_decode_inter_block(w, block[i], i, (cbp >> (5 - i)) & 1) < 0)
786             {
787                 av_log(s->avctx, AV_LOG_ERROR, "\nerror while decoding inter block: %d x %d (%d)\n", s->mb_x, s->mb_y, i);
788                 return -1;
789             }
790         }    
791     } else {
792 //if(s->pict_type==P_TYPE)
793 //   printf("%d%d ", s->inter_intra_pred, cbp);
794 //printf("I at %d %d %d %06X\n", s->mb_x, s->mb_y, ((cbp&3)? 1 : 0) +((cbp&0x3C)? 2 : 0), show_bits(&s->gb, 24));
795         s->ac_pred = get_bits1(&s->gb);
796         if(s->inter_intra_pred){
797             s->h263_aic_dir= get_vlc2(&s->gb, inter_intra_vlc.table, INTER_INTRA_VLC_BITS, 1);
798 //            printf("%d%d %d %d/", s->ac_pred, s->h263_aic_dir, s->mb_x, s->mb_y);
799         }
800         if(s->per_mb_rl_table && cbp){
801             s->rl_table_index = decode012(&s->gb);
802             s->rl_chroma_table_index = s->rl_table_index;
803         }
804     
805         for (i = 0; i < 6; i++) {
806             if (msmpeg4_decode_block(s, block[i], i, (cbp >> (5 - i)) & 1, NULL) < 0)
807             {
808                 av_log(s->avctx, AV_LOG_ERROR, "\nerror while decoding intra block: %d x %d (%d)\n", s->mb_x, s->mb_y, i);
809                 return -1;
810             }
811         }    
812     }
813
814     return 0;
815 }
816
817 static int wmv2_decode_init(AVCodecContext *avctx){
818     Wmv2Context * const w= avctx->priv_data;
819     
820     if(ff_h263_decode_init(avctx) < 0)
821         return -1;
822     
823     wmv2_common_init(w);
824     
825     return 0;
826 }
827
828 AVCodec wmv2_decoder = {
829     "wmv2",
830     CODEC_TYPE_VIDEO,
831     CODEC_ID_WMV2,
832     sizeof(Wmv2Context),
833     wmv2_decode_init,
834     NULL,
835     ff_h263_decode_end,
836     ff_h263_decode_frame,
837     CODEC_CAP_DRAW_HORIZ_BAND | CODEC_CAP_DR1,
838 };
839
840 #ifdef CONFIG_ENCODERS
841 AVCodec wmv2_encoder = {
842     "wmv2",
843     CODEC_TYPE_VIDEO,
844     CODEC_ID_WMV2,
845     sizeof(Wmv2Context),
846     wmv2_encode_init,
847     MPV_encode_picture,
848     MPV_encode_end,
849 };
850 #endif