h263.c 156 KB
Newer Older
Fabrice Bellard's avatar
Fabrice Bellard committed
1 2
/*
 * H263/MPEG4 backend for ffmpeg encoder and decoder
Fabrice Bellard's avatar
Fabrice Bellard committed
3
 * Copyright (c) 2000,2001 Fabrice Bellard.
4
 * H263+ support.
Fabrice Bellard's avatar
Fabrice Bellard committed
5 6
 * Copyright (c) 2001 Juan J. Sierralta P.
 *
Fabrice Bellard's avatar
Fabrice Bellard committed
7 8 9 10
 * This library is free software; you can redistribute it and/or
 * modify it under the terms of the GNU Lesser General Public
 * License as published by the Free Software Foundation; either
 * version 2 of the License, or (at your option) any later version.
Fabrice Bellard's avatar
Fabrice Bellard committed
11
 *
Fabrice Bellard's avatar
Fabrice Bellard committed
12
 * This library is distributed in the hope that it will be useful,
Fabrice Bellard's avatar
Fabrice Bellard committed
13
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
Fabrice Bellard's avatar
Fabrice Bellard committed
14 15
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
 * Lesser General Public License for more details.
Fabrice Bellard's avatar
Fabrice Bellard committed
16
 *
Fabrice Bellard's avatar
Fabrice Bellard committed
17 18 19
 * You should have received a copy of the GNU Lesser General Public
 * License along with this library; if not, write to the Free Software
 * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
20
 *
21 22 23
 * ac prediction encoding, b-frame support, error resilience, optimizations,
 * qpel decoding, gmc decoding, interlaced decoding, 
 * by Michael Niedermayer <michaelni@gmx.at>
Fabrice Bellard's avatar
Fabrice Bellard committed
24
 */
25 26
 
//#define DEBUG
Fabrice Bellard's avatar
Fabrice Bellard committed
27 28 29 30 31 32 33
#include "common.h"
#include "dsputil.h"
#include "avcodec.h"
#include "mpegvideo.h"
#include "h263data.h"
#include "mpeg4data.h"

34 35 36 37
//#undef NDEBUG
//#include <assert.h>

#if 1
Felix Bünemann's avatar
Felix Bünemann committed
38
#define PRINT_MB_TYPE(a) {}
39 40 41
#else
#define PRINT_MB_TYPE(a) printf(a)
#endif
Michael Niedermayer's avatar
Michael Niedermayer committed
42

43 44 45 46 47 48 49 50 51
#define INTRA_MCBPC_VLC_BITS 6
#define INTER_MCBPC_VLC_BITS 6
#define CBPY_VLC_BITS 6
#define MV_VLC_BITS 9
#define DC_VLC_BITS 9
#define SPRITE_TRAJ_VLC_BITS 6
#define MB_TYPE_B_VLC_BITS 4
#define TEX_VLC_BITS 9

52
#ifdef CONFIG_ENCODERS
Fabrice Bellard's avatar
Fabrice Bellard committed
53 54
static void h263_encode_block(MpegEncContext * s, DCTELEM * block,
			      int n);
55
static void h263_encode_motion(MpegEncContext * s, int val, int fcode);
56
static void h263p_encode_umotion(MpegEncContext * s, int val);
57
static inline void mpeg4_encode_block(MpegEncContext * s, DCTELEM * block,
58 59
			       int n, int dc, UINT8 *scan_table, 
                               PutBitContext *dc_pb, PutBitContext *ac_pb);
60 61
#endif

Michael Niedermayer's avatar
Michael Niedermayer committed
62
static int h263_decode_motion(MpegEncContext * s, int pred, int fcode);
63
static int h263p_decode_umotion(MpegEncContext * s, int pred);
Fabrice Bellard's avatar
Fabrice Bellard committed
64 65
static int h263_decode_block(MpegEncContext * s, DCTELEM * block,
                             int n, int coded);
66 67
static inline int mpeg4_decode_dc(MpegEncContext * s, int n, int *dir_ptr);
static inline int mpeg4_decode_block(MpegEncContext * s, DCTELEM * block,
Michael Niedermayer's avatar
Michael Niedermayer committed
68
                              int n, int coded, int intra);
69
static int h263_pred_dc(MpegEncContext * s, int n, UINT16 **dc_val_ptr);
70 71
static void mpeg4_inv_pred_ac(MpegEncContext * s, INT16 *block, int n,
                              int dir);
72
static void mpeg4_decode_sprite_trajectory(MpegEncContext * s);
73 74
static inline int ff_mpeg4_pred_dc(MpegEncContext * s, int n, UINT16 **dc_val_ptr, int *dir_ptr);

75
extern UINT32 inverse[256];
Fabrice Bellard's avatar
Fabrice Bellard committed
76

77 78 79 80 81
static UINT16 uni_DCtab_lum  [512][2];
static UINT16 uni_DCtab_chrom[512][2];

#ifdef CONFIG_ENCODERS
static UINT16 (*mv_penalty)[MAX_MV*2+1]= NULL;
82
static UINT8 fcode_tab[MAX_MV*2+1];
83
static UINT8 umv_fcode_tab[MAX_MV*2+1];
84

85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100
static UINT32 uni_mpeg4_intra_rl_bits[64*64*2*2];
static UINT8  uni_mpeg4_intra_rl_len [64*64*2*2];
static UINT32 uni_mpeg4_inter_rl_bits[64*64*2*2];
static UINT8  uni_mpeg4_inter_rl_len [64*64*2*2];
#define UNI_MPEG4_ENC_INDEX(last,run,level) ((last)*128 + (run)*256 + (level))
//#define UNI_MPEG4_ENC_INDEX(last,run,level) ((last)*128*64 + (run) + (level)*64)

/* mpeg4
inter
max level: 24/6
max run: 53/63

intra
max level: 53/16
max run: 29/41
*/
101 102
#endif

103

Fabrice Bellard's avatar
Fabrice Bellard committed
104 105 106 107 108
int h263_get_picture_format(int width, int height)
{
    int format;

    if (width == 128 && height == 96)
109
        format = 1;
Fabrice Bellard's avatar
Fabrice Bellard committed
110
    else if (width == 176 && height == 144)
111
        format = 2;
Fabrice Bellard's avatar
Fabrice Bellard committed
112
    else if (width == 352 && height == 288)
113
        format = 3;
Fabrice Bellard's avatar
Fabrice Bellard committed
114
    else if (width == 704 && height == 576)
115
        format = 4;
Fabrice Bellard's avatar
Fabrice Bellard committed
116
    else if (width == 1408 && height == 1152)
117
        format = 5;
Fabrice Bellard's avatar
Fabrice Bellard committed
118 119 120 121 122
    else
        format = 7;
    return format;
}

Michael Niedermayer's avatar
Michael Niedermayer committed
123 124 125 126 127
static void float_aspect_to_info(MpegEncContext * s, float aspect){
    int i;

    aspect*= s->height/(double)s->width;
//printf("%f\n", aspect);
Michael Niedermayer's avatar
Michael Niedermayer committed
128
    
Michael Niedermayer's avatar
Michael Niedermayer committed
129 130
    if(aspect==0) aspect= 1.0;

Michael Niedermayer's avatar
Michael Niedermayer committed
131
    ff_float2fraction(&s->aspected_width, &s->aspected_height, aspect, 255);
Michael Niedermayer's avatar
Michael Niedermayer committed
132 133 134 135 136 137 138 139

//printf("%d %d\n", s->aspected_width, s->aspected_height);
    for(i=1; i<6; i++){
        if(s->aspected_width == pixel_aspect[i][0] && s->aspected_height== pixel_aspect[i][1]){
            s->aspect_ratio_info=i;
            return;
        }
    }
Michael Niedermayer's avatar
Michael Niedermayer committed
140
    
Michael Niedermayer's avatar
Michael Niedermayer committed
141
    s->aspect_ratio_info= FF_ASPECT_EXTENDED;
Michael Niedermayer's avatar
Michael Niedermayer committed
142 143
}

Fabrice Bellard's avatar
Fabrice Bellard committed
144 145
void h263_encode_picture_header(MpegEncContext * s, int picture_number)
{
146
    int format;
Fabrice Bellard's avatar
Fabrice Bellard committed
147 148

    align_put_bits(&s->pb);
149 150

    /* Update the pointer to last GOB */
151
    s->ptr_lastgob = pbBufPtr(&s->pb);
152 153 154
    s->gob_number = 0;

    put_bits(&s->pb, 22, 0x20); /* PSC */
155
    put_bits(&s->pb, 8, (((INT64)s->picture_number * 30 * FRAME_RATE_BASE) / 
Fabrice Bellard's avatar
Fabrice Bellard committed
156 157 158 159 160 161 162
                         s->frame_rate) & 0xff);

    put_bits(&s->pb, 1, 1);	/* marker */
    put_bits(&s->pb, 1, 0);	/* h263 id */
    put_bits(&s->pb, 1, 0);	/* split screen off */
    put_bits(&s->pb, 1, 0);	/* camera  off */
    put_bits(&s->pb, 1, 0);	/* freeze picture release off */
163 164
    
    format = h263_get_picture_format(s->width, s->height);
Fabrice Bellard's avatar
Fabrice Bellard committed
165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182
    if (!s->h263_plus) {
        /* H.263v1 */
        put_bits(&s->pb, 3, format);
        put_bits(&s->pb, 1, (s->pict_type == P_TYPE));
        /* By now UMV IS DISABLED ON H.263v1, since the restrictions
        of H.263v1 UMV implies to check the predicted MV after
        calculation of the current MB to see if we're on the limits */
        put_bits(&s->pb, 1, 0);	/* unrestricted motion vector: off */
        put_bits(&s->pb, 1, 0);	/* SAC: off */
        put_bits(&s->pb, 1, 0);	/* advanced prediction mode: off */
        put_bits(&s->pb, 1, 0);	/* not PB frame */
        put_bits(&s->pb, 5, s->qscale);
        put_bits(&s->pb, 1, 0);	/* Continuous Presence Multipoint mode: off */
    } else {
        /* H.263v2 */
        /* H.263 Plus PTYPE */
        put_bits(&s->pb, 3, 7);
        put_bits(&s->pb,3,1); /* Update Full Extended PTYPE */
183 184 185 186 187
        if (format == 7)
            put_bits(&s->pb,3,6); /* Custom Source Format */
        else
            put_bits(&s->pb, 3, format);
            
Fabrice Bellard's avatar
Fabrice Bellard committed
188
        put_bits(&s->pb,1,0); /* Custom PCF: off */
189 190
        s->umvplus = (s->pict_type == P_TYPE) && s->unrestricted_mv;
        put_bits(&s->pb, 1, s->umvplus); /* Unrestricted Motion Vector */
Fabrice Bellard's avatar
Fabrice Bellard committed
191 192
        put_bits(&s->pb,1,0); /* SAC: off */
        put_bits(&s->pb,1,0); /* Advanced Prediction Mode: off */
193
        put_bits(&s->pb,1,s->h263_aic); /* Advanced Intra Coding */
Fabrice Bellard's avatar
Fabrice Bellard committed
194 195 196 197 198 199 200 201 202 203 204 205 206
        put_bits(&s->pb,1,0); /* Deblocking Filter: off */
        put_bits(&s->pb,1,0); /* Slice Structured: off */
        put_bits(&s->pb,1,0); /* Reference Picture Selection: off */
        put_bits(&s->pb,1,0); /* Independent Segment Decoding: off */
        put_bits(&s->pb,1,0); /* Alternative Inter VLC: off */
        put_bits(&s->pb,1,0); /* Modified Quantization: off */
        put_bits(&s->pb,1,1); /* "1" to prevent start code emulation */
        put_bits(&s->pb,3,0); /* Reserved */
		
        put_bits(&s->pb, 3, s->pict_type == P_TYPE);
		
        put_bits(&s->pb,1,0); /* Reference Picture Resampling: off */
        put_bits(&s->pb,1,0); /* Reduced-Resolution Update: off */
207 208 209 210 211
        if (s->pict_type == I_TYPE)
            s->no_rounding = 0;
        else
            s->no_rounding ^= 1;
        put_bits(&s->pb,1,s->no_rounding); /* Rounding Type */
Fabrice Bellard's avatar
Fabrice Bellard committed
212 213 214 215 216 217
        put_bits(&s->pb,2,0); /* Reserved */
        put_bits(&s->pb,1,1); /* "1" to prevent start code emulation */
		
        /* This should be here if PLUSPTYPE */
        put_bits(&s->pb, 1, 0);	/* Continuous Presence Multipoint mode: off */
		
218 219
		if (format == 7) {
            /* Custom Picture Format (CPFMT) */
Michael Niedermayer's avatar
Michael Niedermayer committed
220
            float_aspect_to_info(s, s->avctx->aspect_ratio);
Michael Niedermayer's avatar
Michael Niedermayer committed
221 222

            put_bits(&s->pb,4,s->aspect_ratio_info);
223 224 225
            put_bits(&s->pb,9,(s->width >> 2) - 1);
            put_bits(&s->pb,1,1); /* "1" to prevent start code emulation */
            put_bits(&s->pb,9,(s->height >> 2));
226 227 228 229 230
	    if (s->aspect_ratio_info == FF_ASPECT_EXTENDED)
	    {
		put_bits(&s->pb, 8, s->aspected_width);
		put_bits(&s->pb, 8, s->aspected_height);
	    }
231 232
        }
        
Fabrice Bellard's avatar
Fabrice Bellard committed
233
        /* Unlimited Unrestricted Motion Vectors Indicator (UUI) */
234
        if (s->umvplus)
Fabrice Bellard's avatar
Fabrice Bellard committed
235 236 237 238 239
            put_bits(&s->pb,1,1); /* Limited according tables of Annex D */
        put_bits(&s->pb, 5, s->qscale);
    }

    put_bits(&s->pb, 1, 0);	/* no PEI */
240 241 242 243 244 245 246 247

    if(s->h263_aic){
         s->y_dc_scale_table= 
         s->c_dc_scale_table= h263_aic_dc_scale_table;
    }else{
        s->y_dc_scale_table=
        s->c_dc_scale_table= ff_mpeg1_dc_scale_table;
    }
Fabrice Bellard's avatar
Fabrice Bellard committed
248 249
}

Michael Niedermayer's avatar
Michael Niedermayer committed
250 251 252
/**
 * Encodes a group of blocks header.
 */
253 254 255
int h263_encode_gob_header(MpegEncContext * s, int mb_line)
{
           align_put_bits(&s->pb);
256 257 258
           flush_put_bits(&s->pb);
           /* Call the RTP callback to send the last GOB */
           if (s->rtp_callback) {
Michael Niedermayer's avatar
Michael Niedermayer committed
259
               int pdif = pbBufPtr(&s->pb) - s->ptr_lastgob;
260 261
               s->rtp_callback(s->ptr_lastgob, pdif, s->gob_number);
           }
262
           put_bits(&s->pb, 17, 1); /* GBSC */
263
           s->gob_number = mb_line / s->gob_index;
264
           put_bits(&s->pb, 5, s->gob_number); /* GN */
265
           put_bits(&s->pb, 2, s->pict_type == I_TYPE); /* GFID */
266
           put_bits(&s->pb, 5, s->qscale); /* GQUANT */
267
           //fprintf(stderr,"\nGOB: %2d size: %d", s->gob_number - 1, pdif);
Michael Niedermayer's avatar
Michael Niedermayer committed
268
    return 0;
269
}
270 271 272 273 274 275 276 277 278

static inline int decide_ac_pred(MpegEncContext * s, DCTELEM block[6][64], int dir[6])
{
    int score0=0, score1=0;
    int i, n;

    for(n=0; n<6; n++){
        INT16 *ac_val, *ac_val1;

Michael Niedermayer's avatar
Michael Niedermayer committed
279
        ac_val = s->ac_val[0][0] + s->block_index[n] * 16;
280 281
        ac_val1= ac_val;
        if(dir[n]){
282 283
            const int xy= s->mb_x + s->mb_y*s->mb_width - s->mb_width;
            /* top prediction */
Michael Niedermayer's avatar
Michael Niedermayer committed
284
            ac_val-= s->block_wrap[n]*16;
285 286 287
            if(s->mb_y==0 || s->qscale == s->qscale_table[xy] || n==2 || n==3){
                /* same qscale */
                for(i=1; i<8; i++){
288
                    const int level= block[n][s->idct_permutation[i   ]];
289 290
                    score0+= ABS(level);
                    score1+= ABS(level - ac_val[i+8]);
291
                    ac_val1[i  ]=    block[n][s->idct_permutation[i<<3]];
292 293 294 295 296
                    ac_val1[i+8]= level;
                }
            }else{
                /* different qscale, we must rescale */
                for(i=1; i<8; i++){
297
                    const int level= block[n][s->idct_permutation[i   ]];
298 299
                    score0+= ABS(level);
                    score1+= ABS(level - ROUNDED_DIV(ac_val[i + 8]*s->qscale_table[xy], s->qscale));
300
                    ac_val1[i  ]=    block[n][s->idct_permutation[i<<3]];
301 302
                    ac_val1[i+8]= level;
                }
303 304
            }
        }else{
305 306
            const int xy= s->mb_x-1 + s->mb_y*s->mb_width;
            /* left prediction */
307
            ac_val-= 16;
308 309 310
            if(s->mb_x==0 || s->qscale == s->qscale_table[xy] || n==1 || n==3){
                /* same qscale */
                for(i=1; i<8; i++){
311
                    const int level= block[n][s->idct_permutation[i<<3]];
312 313 314
                    score0+= ABS(level);
                    score1+= ABS(level - ac_val[i]);
                    ac_val1[i  ]= level;
315
                    ac_val1[i+8]=    block[n][s->idct_permutation[i   ]];
316 317 318 319
                }
            }else{
                /* different qscale, we must rescale */
                for(i=1; i<8; i++){
320
                    const int level= block[n][s->idct_permutation[i<<3]];
321 322 323
                    score0+= ABS(level);
                    score1+= ABS(level - ROUNDED_DIV(ac_val[i]*s->qscale_table[xy], s->qscale));
                    ac_val1[i  ]= level;
324
                    ac_val1[i+8]=    block[n][s->idct_permutation[i   ]];
325
                }
326 327 328 329 330 331 332
            }
        }
    }

    return score0 > score1 ? 1 : 0;    
}

333 334 335
/**
 * modify qscale so that encoding is acually possible in h263 (limit difference to -2..2)
 */
336
void ff_clean_h263_qscales(MpegEncContext *s){
337 338 339 340 341 342 343 344 345 346
    int i;
    
    for(i=1; i<s->mb_num; i++){
        if(s->qscale_table[i] - s->qscale_table[i-1] >2)
            s->qscale_table[i]= s->qscale_table[i-1]+2;
    }
    for(i=s->mb_num-2; i>=0; i--){
        if(s->qscale_table[i] - s->qscale_table[i+1] >2)
            s->qscale_table[i]= s->qscale_table[i+1]+2;
    }
347 348
}

349 350 351
/**
 * modify mb_type & qscale so that encoding is acually possible in mpeg4
 */
352 353 354 355
void ff_clean_mpeg4_qscales(MpegEncContext *s){
    int i;
    
    ff_clean_h263_qscales(s);
356 357 358 359 360 361 362 363 364 365 366 367 368 369 370 371 372 373 374 375 376 377 378 379 380 381 382 383 384 385 386 387 388 389 390 391
    
    for(i=1; i<s->mb_num; i++){
        if(s->qscale_table[i] != s->qscale_table[i-1] && (s->mb_type[i]&MB_TYPE_INTER4V)){
            s->mb_type[i]&= ~MB_TYPE_INTER4V;
            s->mb_type[i]|= MB_TYPE_INTER;
        }
    }

    if(s->pict_type== B_TYPE){
        int odd=0;
        /* ok, come on, this isnt funny anymore, theres more code for handling this mpeg4 mess than
           for the actual adaptive quantization */
        
        for(i=0; i<s->mb_num; i++){
            odd += s->qscale_table[i]&1;
        }
        
        if(2*odd > s->mb_num) odd=1;
        else                  odd=0;
        
        for(i=0; i<s->mb_num; i++){
            if((s->qscale_table[i]&1) != odd)
                s->qscale_table[i]++;
            if(s->qscale_table[i] > 31)
                s->qscale_table[i]= 31;
        }            
    
        for(i=1; i<s->mb_num; i++){
            if(s->qscale_table[i] != s->qscale_table[i-1] && (s->mb_type[i]&MB_TYPE_DIRECT)){
                s->mb_type[i]&= ~MB_TYPE_DIRECT;
                s->mb_type[i]|= MB_TYPE_BIDIR;
            }
        }
    }
}

392
#ifdef CONFIG_ENCODERS
393 394 395 396
void mpeg4_encode_mb(MpegEncContext * s,
		    DCTELEM block[6][64],
		    int motion_x, int motion_y)
{
397
    int cbpc, cbpy, i, pred_x, pred_y;
398
    int bits;
399 400 401 402
    PutBitContext * const pb2    = s->data_partitioning                         ? &s->pb2    : &s->pb;
    PutBitContext * const tex_pb = s->data_partitioning && s->pict_type!=B_TYPE ? &s->tex_pb : &s->pb;
    PutBitContext * const dc_pb  = s->data_partitioning && s->pict_type!=I_TYPE ? &s->pb2    : &s->pb;
    const int interleaved_stats= (s->flags&CODEC_FLAG_PASS1) && !s->data_partitioning ? 1 : 0;
403
    const int dquant_code[5]= {1,0,9,2,3};
404
    
405 406 407
    //    printf("**mb x=%d y=%d\n", s->mb_x, s->mb_y);
    if (!s->mb_intra) {
        /* compute cbp */
408
        int cbp = 0;
409
        for (i = 0; i < 6; i++) {
410 411
            if (s->block_last_index[i] >= 0)
                cbp |= 1 << (5 - i);
412
        }
413 414 415 416 417 418 419 420 421 422 423

        if(s->pict_type==B_TYPE){
            static const int mb_type_table[8]= {-1, 2, 3, 1,-1,-1,-1, 0}; /* convert from mv_dir to type */
            int mb_type=  mb_type_table[s->mv_dir];
            
            if(s->mb_x==0){
                s->last_mv[0][0][0]= 
                s->last_mv[0][0][1]= 
                s->last_mv[1][0][0]= 
                s->last_mv[1][0][1]= 0;
            }
424 425 426 427
            
            assert(s->dquant>=-2 && s->dquant<=2);
            assert((s->dquant&1)==0);
            assert(mb_type>=0);
428 429

            /* nothing to do if this MB was skiped in the next P Frame */
430
            if(s->mbskip_table[s->mb_y * s->mb_width + s->mb_x]){ //FIXME avoid DCT & ...
431 432 433 434 435
                s->skip_count++;
                s->mv[0][0][0]= 
                s->mv[0][0][1]= 
                s->mv[1][0][0]= 
                s->mv[1][0][1]= 0;
436
                s->mv_dir= MV_DIR_FORWARD; //doesnt matter
437
                s->qscale -= s->dquant;
438 439 440 441 442
                return;
            }

            if ((cbp | motion_x | motion_y | mb_type) ==0) {
                /* direct MB with MV={0,0} */
443 444
                assert(s->dquant==0);
                
445
                put_bits(&s->pb, 1, 1); /* mb not coded modb1=1 */
446 447 448 449 450

                if(interleaved_stats){
                    s->misc_bits++;
                    s->last_bits++;
                }
451 452 453 454 455 456 457 458
                s->skip_count++;
                return;
            }
            put_bits(&s->pb, 1, 0);	/* mb coded modb1=0 */
            put_bits(&s->pb, 1, cbp ? 0 : 1); /* modb2 */ //FIXME merge
            put_bits(&s->pb, mb_type+1, 1); // this table is so simple that we dont need it :)
            if(cbp) put_bits(&s->pb, 6, cbp);
            
459 460 461 462 463 464 465
            if(cbp && mb_type){
                if(s->dquant)
                    put_bits(&s->pb, 2, (s->dquant>>2)+3);
                else
                    put_bits(&s->pb, 1, 0);
            }else
                s->qscale -= s->dquant;
466 467 468 469 470 471 472
            
            if(!s->progressive_sequence){
                if(cbp)
                    put_bits(&s->pb, 1, s->interlaced_dct);
                if(mb_type) // not diect mode
                    put_bits(&s->pb, 1, 0); // no interlaced ME yet
            }
473

474 475 476 477 478
            if(interleaved_stats){
                bits= get_bit_count(&s->pb);
                s->misc_bits+= bits - s->last_bits;
                s->last_bits=bits;
            }
Michael Niedermayer's avatar
Michael Niedermayer committed
479

480 481 482 483 484
            switch(mb_type)
            {
            case 0: /* direct */
                h263_encode_motion(s, motion_x, 1);
                h263_encode_motion(s, motion_y, 1);                
485 486
                s->b_count++;
                s->f_count++;
487 488 489 490 491 492 493 494 495 496
                break;
            case 1: /* bidir */
                h263_encode_motion(s, s->mv[0][0][0] - s->last_mv[0][0][0], s->f_code);
                h263_encode_motion(s, s->mv[0][0][1] - s->last_mv[0][0][1], s->f_code);
                h263_encode_motion(s, s->mv[1][0][0] - s->last_mv[1][0][0], s->b_code);
                h263_encode_motion(s, s->mv[1][0][1] - s->last_mv[1][0][1], s->b_code);
                s->last_mv[0][0][0]= s->mv[0][0][0];
                s->last_mv[0][0][1]= s->mv[0][0][1];
                s->last_mv[1][0][0]= s->mv[1][0][0];
                s->last_mv[1][0][1]= s->mv[1][0][1];
497 498
                s->b_count++;
                s->f_count++;
499 500 501 502 503 504
                break;
            case 2: /* backward */
                h263_encode_motion(s, motion_x - s->last_mv[1][0][0], s->b_code);
                h263_encode_motion(s, motion_y - s->last_mv[1][0][1], s->b_code);
                s->last_mv[1][0][0]= motion_x;
                s->last_mv[1][0][1]= motion_y;
505
                s->b_count++;
506 507 508 509 510 511
                break;
            case 3: /* forward */
                h263_encode_motion(s, motion_x - s->last_mv[0][0][0], s->f_code);
                h263_encode_motion(s, motion_y - s->last_mv[0][0][1], s->f_code);
                s->last_mv[0][0][0]= motion_x;
                s->last_mv[0][0][1]= motion_y;
512
                s->f_count++;
513
                break;
514 515
            default:
                printf("unknown mb type\n");
516 517
                return;
            }
518 519 520 521 522 523

            if(interleaved_stats){
                bits= get_bit_count(&s->pb);
                s->mv_bits+= bits - s->last_bits;
                s->last_bits=bits;
            }
524

525 526
            /* encode each block */
            for (i = 0; i < 6; i++) {
527
                mpeg4_encode_block(s, block[i], i, 0, s->intra_scantable.permutated, NULL, &s->pb);
528 529 530 531 532 533
            }

            if(interleaved_stats){
                bits= get_bit_count(&s->pb);
                s->p_tex_bits+= bits - s->last_bits;
                s->last_bits=bits;
534 535
            }
        }else{ /* s->pict_type==B_TYPE */
536
            if ((cbp | motion_x | motion_y | s->dquant) == 0 && s->mv_type==MV_TYPE_16X16) {
Michael Niedermayer's avatar
Michael Niedermayer committed
537 538 539 540
                /* check if the B frames can skip it too, as we must skip it if we skip here 
                   why didnt they just compress the skip-mb bits instead of reusing them ?! */
                if(s->max_b_frames>0){
                    int i;
541 542 543 544 545 546 547 548 549 550 551
                    int x,y, offset;
                    uint8_t *p_pic;

                    x= s->mb_x*16;
                    y= s->mb_y*16;
                    if(x+16 > s->width)  x= s->width-16;
                    if(y+16 > s->height) y= s->height-16;

                    offset= x + y*s->linesize;
                    p_pic= s->new_picture[0] + offset;
                    
Michael Niedermayer's avatar
Michael Niedermayer committed
552 553
                    s->mb_skiped=1;
                    for(i=0; i<s->max_b_frames; i++){
Michael Niedermayer's avatar
Michael Niedermayer committed
554 555 556 557 558 559
                        uint8_t *b_pic;
                        int diff;

                        if(s->coded_order[i+1].pict_type!=B_TYPE) break;

                        b_pic= s->coded_order[i+1].picture[0] + offset;
560
			diff= s->dsp.pix_abs16x16(p_pic, b_pic, s->linesize);
561
                        if(diff>s->qscale*70){ //FIXME check that 70 is optimal
Michael Niedermayer's avatar
Michael Niedermayer committed
562 563 564 565 566 567 568 569 570 571
                            s->mb_skiped=0;
                            break;
                        }
                    }
                }else
                    s->mb_skiped=1; 

                if(s->mb_skiped==1){
                    /* skip macroblock */
                    put_bits(&s->pb, 1, 1);
572 573 574 575 576

                    if(interleaved_stats){
                        s->misc_bits++;
                        s->last_bits++;
                    }
Michael Niedermayer's avatar
Michael Niedermayer committed
577 578 579
                    s->skip_count++;
                    return;
                }
580
            }
581

582 583 584
            put_bits(&s->pb, 1, 0);	/* mb coded */
            if(s->mv_type==MV_TYPE_16X16){
                cbpc = cbp & 3;
585
                if(s->dquant) cbpc+= 8;
586 587 588
                put_bits(&s->pb,
                        inter_MCBPC_bits[cbpc],
                        inter_MCBPC_code[cbpc]);
589

590 591
                cbpy = cbp >> 2;
                cbpy ^= 0xf;
592
                put_bits(pb2, cbpy_tab[cbpy][1], cbpy_tab[cbpy][0]);
593 594
                if(s->dquant)
                    put_bits(pb2, 2, dquant_code[s->dquant+2]);
595 596 597 598 599 600

                if(!s->progressive_sequence){
                    if(cbp)
                        put_bits(pb2, 1, s->interlaced_dct);
                    put_bits(pb2, 1, 0); // no interlaced ME yet
                }
601
                    
602 603 604 605 606
                if(interleaved_stats){
                    bits= get_bit_count(&s->pb);
                    s->misc_bits+= bits - s->last_bits;
                    s->last_bits=bits;
                }
607

608 609 610 611 612 613 614 615 616 617 618 619
                /* motion vectors: 16x16 mode */
                h263_pred_motion(s, 0, &pred_x, &pred_y);
            
                h263_encode_motion(s, motion_x - pred_x, s->f_code);
                h263_encode_motion(s, motion_y - pred_y, s->f_code);
            }else{
                cbpc = (cbp & 3)+16;
                put_bits(&s->pb,
                        inter_MCBPC_bits[cbpc],
                        inter_MCBPC_code[cbpc]);
                cbpy = cbp >> 2;
                cbpy ^= 0xf;
620
                put_bits(pb2, cbpy_tab[cbpy][1], cbpy_tab[cbpy][0]);
621

622 623 624 625 626
                if(!s->progressive_sequence){
                    if(cbp)
                        put_bits(pb2, 1, s->interlaced_dct);
                }
    
627 628 629 630 631
                if(interleaved_stats){
                    bits= get_bit_count(&s->pb);
                    s->misc_bits+= bits - s->last_bits;
                    s->last_bits=bits;
                }
632 633 634 635 636 637 638 639

                for(i=0; i<4; i++){
                    /* motion vectors: 8x8 mode*/
                    h263_pred_motion(s, i, &pred_x, &pred_y);

                    h263_encode_motion(s, s->motion_val[ s->block_index[i] ][0] - pred_x, s->f_code);
                    h263_encode_motion(s, s->motion_val[ s->block_index[i] ][1] - pred_y, s->f_code);
                }
Michael Niedermayer's avatar
Michael Niedermayer committed
640
            }
641 642 643 644 645 646

            if(interleaved_stats){ 
                bits= get_bit_count(&s->pb);
                s->mv_bits+= bits - s->last_bits;
                s->last_bits=bits;
            }
647

648 649
            /* encode each block */
            for (i = 0; i < 6; i++) {
650
                mpeg4_encode_block(s, block[i], i, 0, s->intra_scantable.permutated, NULL, tex_pb);
651 652 653 654 655 656
            }

            if(interleaved_stats){
                bits= get_bit_count(&s->pb);
                s->p_tex_bits+= bits - s->last_bits;
                s->last_bits=bits;
657
            }
658
            s->f_count++;
659 660
        }
    } else {
661
        int cbp;
662 663 664 665 666 667 668 669 670
        int dc_diff[6];   //dc values with the dc prediction subtracted 
        int dir[6];  //prediction direction
        int zigzag_last_index[6];
        UINT8 *scan_table[6];

        for(i=0; i<6; i++){
            const int level= block[i][0];
            UINT16 *dc_ptr;

671
            dc_diff[i]= level - ff_mpeg4_pred_dc(s, i, &dc_ptr, &dir[i]);
672 673 674 675 676 677 678 679 680 681 682 683 684 685 686
            if (i < 4) {
                *dc_ptr = level * s->y_dc_scale;
            } else {
                *dc_ptr = level * s->c_dc_scale;
            }
        }

        s->ac_pred= decide_ac_pred(s, block, dir);

        if(s->ac_pred){
            for(i=0; i<6; i++){
                UINT8 *st;
                int last_index;

                mpeg4_inv_pred_ac(s, block[i], i, dir[i]);
687 688
                if (dir[i]==0) st = s->intra_v_scantable.permutated; /* left */
                else           st = s->intra_h_scantable.permutated; /* top */
689 690 691 692 693 694 695 696 697

                for(last_index=63; last_index>=0; last_index--) //FIXME optimize
                    if(block[i][st[last_index]]) break;
                zigzag_last_index[i]= s->block_last_index[i];
                s->block_last_index[i]= last_index;
                scan_table[i]= st;
            }
        }else{
            for(i=0; i<6; i++)
698
                scan_table[i]= s->intra_scantable.permutated;
699 700 701 702 703 704 705 706 707 708 709
        }

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

        cbpc = cbp & 3;
        if (s->pict_type == I_TYPE) {
710
            if(s->dquant) cbpc+=4;
711 712 713 714
            put_bits(&s->pb,
                intra_MCBPC_bits[cbpc],
                intra_MCBPC_code[cbpc]);
        } else {
715
            if(s->dquant) cbpc+=8;
716 717 718 719 720
            put_bits(&s->pb, 1, 0);	/* mb coded */
            put_bits(&s->pb,
                inter_MCBPC_bits[cbpc + 4],
                inter_MCBPC_code[cbpc + 4]);
        }
721
        put_bits(pb2, 1, s->ac_pred);
722
        cbpy = cbp >> 2;
723
        put_bits(pb2, cbpy_tab[cbpy][1], cbpy_tab[cbpy][0]);
724 725
        if(s->dquant)
            put_bits(dc_pb, 2, dquant_code[s->dquant+2]);
726

727 728 729 730
        if(!s->progressive_sequence){
            put_bits(dc_pb, 1, s->interlaced_dct);
        }

731 732 733 734 735
        if(interleaved_stats){
            bits= get_bit_count(&s->pb);
            s->misc_bits+= bits - s->last_bits;
            s->last_bits=bits;
        }
736

737 738
        /* encode each block */
        for (i = 0; i < 6; i++) {
739
            mpeg4_encode_block(s, block[i], i, dc_diff[i], scan_table[i], dc_pb, tex_pb);
740 741
        }

742 743 744 745 746
        if(interleaved_stats){
            bits= get_bit_count(&s->pb);
            s->i_tex_bits+= bits - s->last_bits;
            s->last_bits=bits;
        }
747 748
        s->i_count++;

749 750 751 752 753 754
        /* restore ac coeffs & last_index stuff if we messed them up with the prediction */
        if(s->ac_pred){
            for(i=0; i<6; i++){
                int j;    
                INT16 *ac_val;

Michael Niedermayer's avatar
Michael Niedermayer committed
755
                ac_val = s->ac_val[0][0] + s->block_index[i] * 16;
756 757 758

                if(dir[i]){
                    for(j=1; j<8; j++) 
759
                        block[i][s->idct_permutation[j   ]]= ac_val[j+8];
760 761
                }else{
                    for(j=1; j<8; j++) 
762
                        block[i][s->idct_permutation[j<<3]]= ac_val[j  ];
763 764 765 766 767 768 769
                }
                s->block_last_index[i]= zigzag_last_index[i];
            }
        }
    }
}

Fabrice Bellard's avatar
Fabrice Bellard committed
770 771 772 773 774
void h263_encode_mb(MpegEncContext * s,
		    DCTELEM block[6][64],
		    int motion_x, int motion_y)
{
    int cbpc, cbpy, i, cbp, pred_x, pred_y;
775 776 777
    INT16 pred_dc;
    INT16 rec_intradc[6];
    UINT16 *dc_ptr[6];
778
    const int dquant_code[5]= {1,0,9,2,3};
779 780
           
    //printf("**mb x=%d y=%d\n", s->mb_x, s->mb_y);
781
    if (!s->mb_intra) {
782
        /* compute cbp */
783 784 785 786 787
        cbp = 0;
        for (i = 0; i < 6; i++) {
            if (s->block_last_index[i] >= 0)
                cbp |= 1 << (5 - i);
        }
788
        if ((cbp | motion_x | motion_y | s->dquant) == 0) {
789 790 791 792
            /* skip macroblock */
            put_bits(&s->pb, 1, 1);
            return;
        }
793
        put_bits(&s->pb, 1, 0);	/* mb coded */
794
        cbpc = cbp & 3;
795
        if(s->dquant) cbpc+= 8;
796
        put_bits(&s->pb,
797 798
		    inter_MCBPC_bits[cbpc],
		    inter_MCBPC_code[cbpc]);
799 800 801
        cbpy = cbp >> 2;
        cbpy ^= 0xf;
        put_bits(&s->pb, cbpy_tab[cbpy][1], cbpy_tab[cbpy][0]);
802 803
        if(s->dquant)
            put_bits(&s->pb, 2, dquant_code[s->dquant+2]);
804

805
        /* motion vectors: 16x16 mode only now */
806
        h263_pred_motion(s, 0, &pred_x, &pred_y);
807
      
808 809 810 811 812 813 814 815
        if (!s->umvplus) {  
            h263_encode_motion(s, motion_x - pred_x, s->f_code);
            h263_encode_motion(s, motion_y - pred_y, s->f_code);
        }
        else {
            h263p_encode_umotion(s, motion_x - pred_x);
            h263p_encode_umotion(s, motion_y - pred_y);
            if (((motion_x - pred_x) == 1) && ((motion_y - pred_y) == 1))
816
                /* To prevent Start Code emulation */
817 818
                put_bits(&s->pb,1,1);
        }
819 820 821 822 823 824 825 826 827 828 829 830 831 832 833 834 835 836 837 838 839 840 841 842 843 844 845 846 847 848 849 850 851 852 853 854 855 856 857 858 859 860 861 862 863
    } else {
        int li = s->h263_aic ? 0 : 1;
        
        cbp = 0;
        for(i=0; i<6; i++) {
            /* Predict DC */
            if (s->h263_aic && s->mb_intra) {
                INT16 level = block[i][0];
            
                pred_dc = h263_pred_dc(s, i, &dc_ptr[i]);
                level -= pred_dc;
                /* Quant */
                if (level < 0)
                    level = (level + (s->qscale >> 1))/(s->y_dc_scale);
                else
                    level = (level - (s->qscale >> 1))/(s->y_dc_scale);
                    
                /* AIC can change CBP */
                if (level == 0 && s->block_last_index[i] == 0)
                    s->block_last_index[i] = -1;
                else if (level < -127)
                    level = -127;
                else if (level > 127)
                    level = 127;
                
                block[i][0] = level;
                /* Reconstruction */ 
                rec_intradc[i] = (s->y_dc_scale*level) + pred_dc;
                /* Oddify */
                rec_intradc[i] |= 1;
                //if ((rec_intradc[i] % 2) == 0)
                //    rec_intradc[i]++;
                /* Clipping */
                if (rec_intradc[i] < 0)
                    rec_intradc[i] = 0;
                else if (rec_intradc[i] > 2047)
                    rec_intradc[i] = 2047;
                                
                /* Update AC/DC tables */
                *dc_ptr[i] = rec_intradc[i];
            }
            /* compute cbp */
            if (s->block_last_index[i] >= li)
                cbp |= 1 << (5 - i);
        }
Fabrice Bellard's avatar
Fabrice Bellard committed
864

865 866
        cbpc = cbp & 3;
        if (s->pict_type == I_TYPE) {
867
            if(s->dquant) cbpc+=4;
868 869 870 871
            put_bits(&s->pb,
                intra_MCBPC_bits[cbpc],
                intra_MCBPC_code[cbpc]);
        } else {
872
            if(s->dquant) cbpc+=8;
873 874 875 876 877 878 879 880 881 882 883
            put_bits(&s->pb, 1, 0);	/* mb coded */
            put_bits(&s->pb,
                inter_MCBPC_bits[cbpc + 4],
                inter_MCBPC_code[cbpc + 4]);
        }
        if (s->h263_aic) {
            /* XXX: currently, we do not try to use ac prediction */
            put_bits(&s->pb, 1, 0);	/* no AC prediction */
        }
        cbpy = cbp >> 2;
        put_bits(&s->pb, cbpy_tab[cbpy][1], cbpy_tab[cbpy][0]);
884 885
        if(s->dquant)
            put_bits(&s->pb, 2, dquant_code[s->dquant+2]);
Fabrice Bellard's avatar
Fabrice Bellard committed
886 887
    }

888 889
    for(i=0; i<6; i++) {
        /* encode each block */
890
        h263_encode_block(s, block[i], i);
891 892 893 894 895 896 897 898
    
        /* Update INTRADC for decoding */
        if (s->h263_aic && s->mb_intra) {
            block[i][0] = rec_intradc[i];
            
        }
    }
}
899
#endif
900 901 902 903 904 905 906 907 908 909 910 911 912 913 914 915 916 917 918 919 920

static int h263_pred_dc(MpegEncContext * s, int n, UINT16 **dc_val_ptr)
{
    int x, y, wrap, a, c, pred_dc, scale;
    INT16 *dc_val, *ac_val;

    /* find prediction */
    if (n < 4) {
        x = 2 * s->mb_x + 1 + (n & 1);
        y = 2 * s->mb_y + 1 + ((n & 2) >> 1);
        wrap = s->mb_width * 2 + 2;
        dc_val = s->dc_val[0];
        ac_val = s->ac_val[0][0];
        scale = s->y_dc_scale;
    } else {
        x = s->mb_x + 1;
        y = s->mb_y + 1;
        wrap = s->mb_width + 2;
        dc_val = s->dc_val[n - 4 + 1];
        ac_val = s->ac_val[n - 4 + 1][0];
        scale = s->c_dc_scale;
Fabrice Bellard's avatar
Fabrice Bellard committed
921
    }
922 923 924 925 926 927 928
    /* B C
     * A X 
     */
    a = dc_val[(x - 1) + (y) * wrap];
    c = dc_val[(x) + (y - 1) * wrap];
    
    /* No prediction outside GOB boundary */
929
    if (s->first_slice_line && ((n < 2) || (n > 3)))
930 931 932 933 934 935 936 937 938 939 940 941 942 943
        c = 1024;
    pred_dc = 1024;
    /* just DC prediction */
    if (a != 1024 && c != 1024)
        pred_dc = (a + c) >> 1;
    else if (a != 1024)
        pred_dc = a;
    else
        pred_dc = c;
    
    /* we assume pred is positive */
    //pred_dc = (pred_dc + (scale >> 1)) / scale;
    *dc_val_ptr = &dc_val[x + y * wrap];
    return pred_dc;
Fabrice Bellard's avatar
Fabrice Bellard committed
944 945
}

946

947
void h263_pred_acdc(MpegEncContext * s, INT16 *block, int n)
948
{
949 950
    int x, y, wrap, a, c, pred_dc, scale, i;
    INT16 *dc_val, *ac_val, *ac_val1;
951 952 953 954 955 956 957

    /* find prediction */
    if (n < 4) {
        x = 2 * s->mb_x + 1 + (n & 1);
        y = 2 * s->mb_y + 1 + ((n & 2) >> 1);
        wrap = s->mb_width * 2 + 2;
        dc_val = s->dc_val[0];
958
        ac_val = s->ac_val[0][0];
959 960 961