h263.c 157 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
//#undef NDEBUG
//#include <assert.h>

Michael Niedermayer's avatar
Michael Niedermayer committed
37 38 39
//rounded divison & shift
#define RSHIFT(a,b) ((a) > 0 ? ((a) + (1<<((b)-1)))>>(b) : ((a) + (1<<((b)-1))-1)>>(b))

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

46 47 48 49 50 51 52 53 54
#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

Fabrice Bellard's avatar
Fabrice Bellard committed
55 56
static void h263_encode_block(MpegEncContext * s, DCTELEM * block,
			      int n);
57
static void h263_encode_motion(MpegEncContext * s, int val, int fcode);
58
static void h263p_encode_umotion(MpegEncContext * s, int val);
59
static inline void mpeg4_encode_block(MpegEncContext * s, DCTELEM * block,
60 61
			       int n, int dc, UINT8 *scan_table, 
                               PutBitContext *dc_pb, PutBitContext *ac_pb);
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

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

78
static UINT16 mv_penalty[MAX_FCODE+1][MAX_MV*2+1];
79
static UINT8 fcode_tab[MAX_MV*2+1];
80
static UINT8 umv_fcode_tab[MAX_MV*2+1];
81

82 83
static UINT16 uni_DCtab_lum  [512][2];
static UINT16 uni_DCtab_chrom[512][2];
84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99
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
*/
100

Fabrice Bellard's avatar
Fabrice Bellard committed
101 102 103 104 105
int h263_get_picture_format(int width, int height)
{
    int format;

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

void h263_encode_picture_header(MpegEncContext * s, int picture_number)
{
122
    int format;
Fabrice Bellard's avatar
Fabrice Bellard committed
123 124

    align_put_bits(&s->pb);
125 126

    /* Update the pointer to last GOB */
127
    s->ptr_lastgob = pbBufPtr(&s->pb);
128 129 130
    s->gob_number = 0;

    put_bits(&s->pb, 22, 0x20); /* PSC */
131
    put_bits(&s->pb, 8, (((INT64)s->picture_number * 30 * FRAME_RATE_BASE) / 
Fabrice Bellard's avatar
Fabrice Bellard committed
132 133 134 135 136 137 138
                         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 */
139 140
    
    format = h263_get_picture_format(s->width, s->height);
Fabrice Bellard's avatar
Fabrice Bellard committed
141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158
    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 */
159 160 161 162 163
        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
164
        put_bits(&s->pb,1,0); /* Custom PCF: off */
165 166
        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
167 168
        put_bits(&s->pb,1,0); /* SAC: off */
        put_bits(&s->pb,1,0); /* Advanced Prediction Mode: off */
169
        put_bits(&s->pb,1,s->h263_aic); /* Advanced Intra Coding */
Fabrice Bellard's avatar
Fabrice Bellard committed
170 171 172 173 174 175 176 177 178 179 180 181 182
        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 */
183 184 185 186 187
        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
188 189 190 191 192 193
        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 */
		
194 195
		if (format == 7) {
            /* Custom Picture Format (CPFMT) */
Fabrice Bellard's avatar
Fabrice Bellard committed
196
		
197
	    if (s->aspect_ratio_info)
198
        	put_bits(&s->pb,4,s->aspect_ratio_info);
199
	    else
200
        	put_bits(&s->pb,4,2); /* Aspect ratio: CIF 12:11 (4:3) picture */
201 202 203
            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));
204 205 206 207 208
	    if (s->aspect_ratio_info == FF_ASPECT_EXTENDED)
	    {
		put_bits(&s->pb, 8, s->aspected_width);
		put_bits(&s->pb, 8, s->aspected_height);
	    }
209 210
        }
        
Fabrice Bellard's avatar
Fabrice Bellard committed
211
        /* Unlimited Unrestricted Motion Vectors Indicator (UUI) */
212
        if (s->umvplus)
Fabrice Bellard's avatar
Fabrice Bellard committed
213 214 215 216 217
            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 */
218 219 220 221 222 223 224 225

    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
226 227
}

228 229 230 231 232 233 234
int h263_encode_gob_header(MpegEncContext * s, int mb_line)
{
    int pdif=0;
    
    /* Check to see if we need to put a new GBSC */
    /* for RTP packetization                    */
    if (s->rtp_mode) {
235
        pdif = pbBufPtr(&s->pb) - s->ptr_lastgob;
236 237 238
        if (pdif >= s->rtp_payload_size) {
            /* Bad luck, packet must be cut before */
            align_put_bits(&s->pb);
239 240 241
            flush_put_bits(&s->pb);
            /* Call the RTP callback to send the last GOB */
            if (s->rtp_callback) {
242
                pdif = pbBufPtr(&s->pb) - s->ptr_lastgob;
243 244
                s->rtp_callback(s->ptr_lastgob, pdif, s->gob_number);
            }
245
            s->ptr_lastgob = pbBufPtr(&s->pb);
246
            put_bits(&s->pb, 17, 1); /* GBSC */
247
            s->gob_number = mb_line / s->gob_index;
248
            put_bits(&s->pb, 5, s->gob_number); /* GN */
249
            put_bits(&s->pb, 2, s->pict_type == I_TYPE); /* GFID */
250
            put_bits(&s->pb, 5, s->qscale); /* GQUANT */
251
            //fprintf(stderr,"\nGOB: %2d size: %d", s->gob_number - 1, pdif);
252 253 254 255
            return pdif;
       } else if (pdif + s->mb_line_avgsize >= s->rtp_payload_size) {
           /* Cut the packet before we can't */
           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) {
259
               pdif = pbBufPtr(&s->pb) - s->ptr_lastgob;
260 261
               s->rtp_callback(s->ptr_lastgob, pdif, s->gob_number);
           }
262
           s->ptr_lastgob = pbBufPtr(&s->pb);
263
           put_bits(&s->pb, 17, 1); /* GBSC */
264
           s->gob_number = mb_line / s->gob_index;
265
           put_bits(&s->pb, 5, s->gob_number); /* GN */
266
           put_bits(&s->pb, 2, s->pict_type == I_TYPE); /* GFID */
267
           put_bits(&s->pb, 5, s->qscale); /* GQUANT */
268
           //fprintf(stderr,"\nGOB: %2d size: %d", s->gob_number - 1, pdif);
269 270 271 272 273
           return pdif;
       }
   }
   return 0;
}
274 275 276 277 278 279 280 281 282

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

    return score0 > score1 ? 1 : 0;    
}

337 338 339
/**
 * modify qscale so that encoding is acually possible in h263 (limit difference to -2..2)
 */
340
void ff_clean_h263_qscales(MpegEncContext *s){
341 342 343 344 345 346 347 348 349 350
    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;
    }
351 352
}

353 354 355
/**
 * modify mb_type & qscale so that encoding is acually possible in mpeg4
 */
356 357 358 359
void ff_clean_mpeg4_qscales(MpegEncContext *s){
    int i;
    
    ff_clean_h263_qscales(s);
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 392 393 394 395
    
    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;
            }
        }
    }
}

396 397 398 399
void mpeg4_encode_mb(MpegEncContext * s,
		    DCTELEM block[6][64],
		    int motion_x, int motion_y)
{
400
    int cbpc, cbpy, i, pred_x, pred_y;
401
    int bits;
402 403 404 405
    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;
406
    const int dquant_code[5]= {1,0,9,2,3};
407
    
408 409 410
    //    printf("**mb x=%d y=%d\n", s->mb_x, s->mb_y);
    if (!s->mb_intra) {
        /* compute cbp */
411
        int cbp = 0;
412
        for (i = 0; i < 6; i++) {
413 414
            if (s->block_last_index[i] >= 0)
                cbp |= 1 << (5 - i);
415
        }
416 417 418 419 420 421 422 423 424 425 426

        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;
            }
427 428 429 430
            
            assert(s->dquant>=-2 && s->dquant<=2);
            assert((s->dquant&1)==0);
            assert(mb_type>=0);
431 432

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

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

                if(interleaved_stats){
                    s->misc_bits++;
                    s->last_bits++;
                }
454 455 456 457 458 459 460 461
                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);
            
462 463 464 465 466 467 468
            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;
469 470 471 472 473 474 475
            
            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
            }
476

477 478 479 480 481
            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
482

483 484 485 486 487
            switch(mb_type)
            {
            case 0: /* direct */
                h263_encode_motion(s, motion_x, 1);
                h263_encode_motion(s, motion_y, 1);                
488 489
                s->b_count++;
                s->f_count++;
490 491 492 493 494 495 496 497 498 499
                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];
500 501
                s->b_count++;
                s->f_count++;
502 503 504 505 506 507
                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;
508
                s->b_count++;
509 510 511 512 513 514
                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;
515
                s->f_count++;
516
                break;
517 518
            default:
                printf("unknown mb type\n");
519 520
                return;
            }
521 522 523 524 525 526

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

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

            if(interleaved_stats){
                bits= get_bit_count(&s->pb);
                s->p_tex_bits+= bits - s->last_bits;
                s->last_bits=bits;
537 538
            }
        }else{ /* s->pict_type==B_TYPE */
539
            if ((cbp | motion_x | motion_y | s->dquant) == 0 && s->mv_type==MV_TYPE_16X16) {
Michael Niedermayer's avatar
Michael Niedermayer committed
540 541 542 543
                /* 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;
544 545 546 547 548 549 550 551 552 553 554
                    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
555 556
                    s->mb_skiped=1;
                    for(i=0; i<s->max_b_frames; i++){
Michael Niedermayer's avatar
Michael Niedermayer committed
557 558 559 560 561 562 563
                        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;
                        diff= pix_abs16x16(p_pic, b_pic, s->linesize);
564
                        if(diff>s->qscale*70){ //FIXME check that 70 is optimal
Michael Niedermayer's avatar
Michael Niedermayer committed
565 566 567 568 569 570 571 572 573 574
                            s->mb_skiped=0;
                            break;
                        }
                    }
                }else
                    s->mb_skiped=1; 

                if(s->mb_skiped==1){
                    /* skip macroblock */
                    put_bits(&s->pb, 1, 1);
575 576 577 578 579

                    if(interleaved_stats){
                        s->misc_bits++;
                        s->last_bits++;
                    }
Michael Niedermayer's avatar
Michael Niedermayer committed
580 581 582
                    s->skip_count++;
                    return;
                }
583
            }
584

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

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

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

611 612 613 614 615 616 617 618 619 620 621 622
                /* 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;
623
                put_bits(pb2, cbpy_tab[cbpy][1], cbpy_tab[cbpy][0]);
624

625 626 627 628 629
                if(interleaved_stats){
                    bits= get_bit_count(&s->pb);
                    s->misc_bits+= bits - s->last_bits;
                    s->last_bits=bits;
                }
630 631 632 633 634 635 636 637

                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
638
            }
639 640 641 642 643 644

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

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

            if(interleaved_stats){
                bits= get_bit_count(&s->pb);
                s->p_tex_bits+= bits - s->last_bits;
                s->last_bits=bits;
655
            }
656
            s->f_count++;
657 658
        }
    } else {
659
        int cbp;
660 661 662 663 664 665 666 667 668
        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;

669
            dc_diff[i]= level - ff_mpeg4_pred_dc(s, i, &dc_ptr, &dir[i]);
670 671 672 673 674 675 676 677 678 679 680 681 682 683 684
            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]);
685 686
                if (dir[i]==0) st = s->intra_v_scantable.permutated; /* left */
                else           st = s->intra_h_scantable.permutated; /* top */
687 688 689 690 691 692 693 694 695

                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++)
696
                scan_table[i]= s->intra_scantable.permutated;
697 698 699 700 701 702 703 704 705 706 707
        }

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

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

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

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

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

747 748 749 750 751 752
        /* 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
753
                ac_val = s->ac_val[0][0] + s->block_index[i] * 16;
754 755 756

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

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

803
        /* motion vectors: 16x16 mode only now */
804
        h263_pred_motion(s, 0, &pred_x, &pred_y);
805
      
806 807 808 809 810 811 812 813
        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))
814
                /* To prevent Start Code emulation */
815 816
                put_bits(&s->pb,1,1);
        }
817 818 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
    } 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
862

863 864
        cbpc = cbp & 3;
        if (s->pict_type == I_TYPE) {
865
            if(s->dquant) cbpc+=4;
866 867 868 869
            put_bits(&s->pb,
                intra_MCBPC_bits[cbpc],
                intra_MCBPC_code[cbpc]);
        } else {
870
            if(s->dquant) cbpc+=8;
871 872 873 874 875 876 877 878 879 880 881
            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]);
882 883
        if(s->dquant)
            put_bits(&s->pb, 2, dquant_code[s->dquant+2]);
Fabrice Bellard's avatar
Fabrice Bellard committed
884 885
    }

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

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
918
    }
919 920 921 922 923 924 925
    /* B C
     * A X 
     */
    a = dc_val[(x - 1) + (y) * wrap];
    c = dc_val[(x) + (y - 1) * wrap];
    
    /* No prediction outside GOB boundary */
926
    if (s->first_slice_line && ((n < 2) || (n > 3)))
927 928 929 930 931 932 933 934 935 936 937 938 939 940
        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
941 942
}

943

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

    /* 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];