Imported Debian version 2.5.3~trusty1
[deb_ffmpeg.git] / ffmpeg / libavcodec / ituh263enc.c
1 /*
2 * ITU H263 bitstream encoder
3 * Copyright (c) 2000,2001 Fabrice Bellard
4 * H263+ support.
5 * Copyright (c) 2001 Juan J. Sierralta P
6 * Copyright (c) 2002-2004 Michael Niedermayer <michaelni@gmx.at>
7 *
8 * This file is part of FFmpeg.
9 *
10 * FFmpeg is free software; you can redistribute it and/or
11 * modify it under the terms of the GNU Lesser General Public
12 * License as published by the Free Software Foundation; either
13 * version 2.1 of the License, or (at your option) any later version.
14 *
15 * FFmpeg is distributed in the hope that it will be useful,
16 * but WITHOUT ANY WARRANTY; without even the implied warranty of
17 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
18 * Lesser General Public License for more details.
19 *
20 * You should have received a copy of the GNU Lesser General Public
21 * License along with FFmpeg; if not, write to the Free Software
22 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
23 */
24
25 /**
26 * @file
27 * h263 bitstream encoder.
28 */
29
30 #include <limits.h>
31
32 #include "libavutil/attributes.h"
33 #include "avcodec.h"
34 #include "mpegvideo.h"
35 #include "h263.h"
36 #include "mathops.h"
37 #include "mpegutils.h"
38 #include "unary.h"
39 #include "flv.h"
40 #include "mpeg4video.h"
41 #include "internal.h"
42
43 /**
44 * Table of number of bits a motion vector component needs.
45 */
46 static uint8_t mv_penalty[MAX_FCODE+1][MAX_MV*2+1];
47
48 /**
49 * Minimal fcode that a motion vector component would need.
50 */
51 static uint8_t fcode_tab[MAX_MV*2+1];
52
53 /**
54 * Minimal fcode that a motion vector component would need in umv.
55 * All entries in this table are 1.
56 */
57 static uint8_t umv_fcode_tab[MAX_MV*2+1];
58
59 //unified encoding tables for run length encoding of coefficients
60 //unified in the sense that the specification specifies the encoding in several steps.
61 static uint8_t uni_h263_intra_aic_rl_len [64*64*2*2];
62 static uint8_t uni_h263_inter_rl_len [64*64*2*2];
63 //#define UNI_MPEG4_ENC_INDEX(last,run,level) ((last)*128 + (run)*256 + (level))
64 //#define UNI_MPEG4_ENC_INDEX(last,run,level) ((last)*128*64 + (run) + (level)*64)
65 #define UNI_MPEG4_ENC_INDEX(last,run,level) ((last)*128*64 + (run)*128 + (level))
66
67 static const uint8_t wrong_run[102] = {
68 1, 2, 3, 5, 4, 10, 9, 8,
69 11, 15, 17, 16, 23, 22, 21, 20,
70 19, 18, 25, 24, 27, 26, 11, 7,
71 6, 1, 2, 13, 2, 2, 2, 2,
72 6, 12, 3, 9, 1, 3, 4, 3,
73 7, 4, 1, 1, 5, 5, 14, 6,
74 1, 7, 1, 8, 1, 1, 1, 1,
75 10, 1, 1, 5, 9, 17, 25, 24,
76 29, 33, 32, 41, 2, 23, 28, 31,
77 3, 22, 30, 4, 27, 40, 8, 26,
78 6, 39, 7, 38, 16, 37, 15, 10,
79 11, 12, 13, 14, 1, 21, 20, 18,
80 19, 2, 1, 34, 35, 36
81 };
82
83 /**
84 * Return the 4 bit value that specifies the given aspect ratio.
85 * This may be one of the standard aspect ratios or it specifies
86 * that the aspect will be stored explicitly later.
87 */
88 av_const int ff_h263_aspect_to_info(AVRational aspect){
89 int i;
90
91 if(aspect.num==0) aspect= (AVRational){1,1};
92
93 for(i=1; i<6; i++){
94 if(av_cmp_q(ff_h263_pixel_aspect[i], aspect) == 0){
95 return i;
96 }
97 }
98
99 return FF_ASPECT_EXTENDED;
100 }
101
102 void ff_h263_encode_picture_header(MpegEncContext * s, int picture_number)
103 {
104 int format, coded_frame_rate, coded_frame_rate_base, i, temp_ref;
105 int best_clock_code=1;
106 int best_divisor=60;
107 int best_error= INT_MAX;
108
109 if(s->h263_plus){
110 for(i=0; i<2; i++){
111 int div, error;
112 div= (s->avctx->time_base.num*1800000LL + 500LL*s->avctx->time_base.den) / ((1000LL+i)*s->avctx->time_base.den);
113 div= av_clip(div, 1, 127);
114 error= FFABS(s->avctx->time_base.num*1800000LL - (1000LL+i)*s->avctx->time_base.den*div);
115 if(error < best_error){
116 best_error= error;
117 best_divisor= div;
118 best_clock_code= i;
119 }
120 }
121 }
122 s->custom_pcf= best_clock_code!=1 || best_divisor!=60;
123 coded_frame_rate= 1800000;
124 coded_frame_rate_base= (1000+best_clock_code)*best_divisor;
125
126 avpriv_align_put_bits(&s->pb);
127
128 /* Update the pointer to last GOB */
129 s->ptr_lastgob = put_bits_ptr(&s->pb);
130 put_bits(&s->pb, 22, 0x20); /* PSC */
131 temp_ref= s->picture_number * (int64_t)coded_frame_rate * s->avctx->time_base.num / //FIXME use timestamp
132 (coded_frame_rate_base * (int64_t)s->avctx->time_base.den);
133 put_sbits(&s->pb, 8, temp_ref); /* TemporalReference */
134
135 put_bits(&s->pb, 1, 1); /* marker */
136 put_bits(&s->pb, 1, 0); /* h263 id */
137 put_bits(&s->pb, 1, 0); /* split screen off */
138 put_bits(&s->pb, 1, 0); /* camera off */
139 put_bits(&s->pb, 1, 0); /* freeze picture release off */
140
141 format = ff_match_2uint16(ff_h263_format, FF_ARRAY_ELEMS(ff_h263_format), s->width, s->height);
142 if (!s->h263_plus) {
143 /* H.263v1 */
144 put_bits(&s->pb, 3, format);
145 put_bits(&s->pb, 1, (s->pict_type == AV_PICTURE_TYPE_P));
146 /* By now UMV IS DISABLED ON H.263v1, since the restrictions
147 of H.263v1 UMV implies to check the predicted MV after
148 calculation of the current MB to see if we're on the limits */
149 put_bits(&s->pb, 1, 0); /* Unrestricted Motion Vector: off */
150 put_bits(&s->pb, 1, 0); /* SAC: off */
151 put_bits(&s->pb, 1, s->obmc); /* Advanced Prediction */
152 put_bits(&s->pb, 1, 0); /* only I/P frames, no PB frame */
153 put_bits(&s->pb, 5, s->qscale);
154 put_bits(&s->pb, 1, 0); /* Continuous Presence Multipoint mode: off */
155 } else {
156 int ufep=1;
157 /* H.263v2 */
158 /* H.263 Plus PTYPE */
159
160 put_bits(&s->pb, 3, 7);
161 put_bits(&s->pb,3,ufep); /* Update Full Extended PTYPE */
162 if (format == 8)
163 put_bits(&s->pb,3,6); /* Custom Source Format */
164 else
165 put_bits(&s->pb, 3, format);
166
167 put_bits(&s->pb,1, s->custom_pcf);
168 put_bits(&s->pb,1, s->umvplus); /* Unrestricted Motion Vector */
169 put_bits(&s->pb,1,0); /* SAC: off */
170 put_bits(&s->pb,1,s->obmc); /* Advanced Prediction Mode */
171 put_bits(&s->pb,1,s->h263_aic); /* Advanced Intra Coding */
172 put_bits(&s->pb,1,s->loop_filter); /* Deblocking Filter */
173 put_bits(&s->pb,1,s->h263_slice_structured); /* Slice Structured */
174 put_bits(&s->pb,1,0); /* Reference Picture Selection: off */
175 put_bits(&s->pb,1,0); /* Independent Segment Decoding: off */
176 put_bits(&s->pb,1,s->alt_inter_vlc); /* Alternative Inter VLC */
177 put_bits(&s->pb,1,s->modified_quant); /* Modified Quantization: */
178 put_bits(&s->pb,1,1); /* "1" to prevent start code emulation */
179 put_bits(&s->pb,3,0); /* Reserved */
180
181 put_bits(&s->pb, 3, s->pict_type == AV_PICTURE_TYPE_P);
182
183 put_bits(&s->pb,1,0); /* Reference Picture Resampling: off */
184 put_bits(&s->pb,1,0); /* Reduced-Resolution Update: off */
185 put_bits(&s->pb,1,s->no_rounding); /* Rounding Type */
186 put_bits(&s->pb,2,0); /* Reserved */
187 put_bits(&s->pb,1,1); /* "1" to prevent start code emulation */
188
189 /* This should be here if PLUSPTYPE */
190 put_bits(&s->pb, 1, 0); /* Continuous Presence Multipoint mode: off */
191
192 if (format == 8) {
193 /* Custom Picture Format (CPFMT) */
194 s->aspect_ratio_info= ff_h263_aspect_to_info(s->avctx->sample_aspect_ratio);
195
196 put_bits(&s->pb,4,s->aspect_ratio_info);
197 put_bits(&s->pb,9,(s->width >> 2) - 1);
198 put_bits(&s->pb,1,1); /* "1" to prevent start code emulation */
199 put_bits(&s->pb,9,(s->height >> 2));
200 if (s->aspect_ratio_info == FF_ASPECT_EXTENDED){
201 put_bits(&s->pb, 8, s->avctx->sample_aspect_ratio.num);
202 put_bits(&s->pb, 8, s->avctx->sample_aspect_ratio.den);
203 }
204 }
205 if(s->custom_pcf){
206 if(ufep){
207 put_bits(&s->pb, 1, best_clock_code);
208 put_bits(&s->pb, 7, best_divisor);
209 }
210 put_sbits(&s->pb, 2, temp_ref>>8);
211 }
212
213 /* Unlimited Unrestricted Motion Vectors Indicator (UUI) */
214 if (s->umvplus)
215 // put_bits(&s->pb,1,1); /* Limited according tables of Annex D */
216 //FIXME check actual requested range
217 put_bits(&s->pb,2,1); /* unlimited */
218 if(s->h263_slice_structured)
219 put_bits(&s->pb,2,0); /* no weird submodes */
220
221 put_bits(&s->pb, 5, s->qscale);
222 }
223
224 put_bits(&s->pb, 1, 0); /* no PEI */
225
226 if(s->h263_slice_structured){
227 put_bits(&s->pb, 1, 1);
228
229 av_assert1(s->mb_x == 0 && s->mb_y == 0);
230 ff_h263_encode_mba(s);
231
232 put_bits(&s->pb, 1, 1);
233 }
234 }
235
236 /**
237 * Encode a group of blocks header.
238 */
239 void ff_h263_encode_gob_header(MpegEncContext * s, int mb_line)
240 {
241 put_bits(&s->pb, 17, 1); /* GBSC */
242
243 if(s->h263_slice_structured){
244 put_bits(&s->pb, 1, 1);
245
246 ff_h263_encode_mba(s);
247
248 if(s->mb_num > 1583)
249 put_bits(&s->pb, 1, 1);
250 put_bits(&s->pb, 5, s->qscale); /* GQUANT */
251 put_bits(&s->pb, 1, 1);
252 put_bits(&s->pb, 2, s->pict_type == AV_PICTURE_TYPE_I); /* GFID */
253 }else{
254 int gob_number= mb_line / s->gob_index;
255
256 put_bits(&s->pb, 5, gob_number); /* GN */
257 put_bits(&s->pb, 2, s->pict_type == AV_PICTURE_TYPE_I); /* GFID */
258 put_bits(&s->pb, 5, s->qscale); /* GQUANT */
259 }
260 }
261
262 /**
263 * modify qscale so that encoding is actually possible in h263 (limit difference to -2..2)
264 */
265 void ff_clean_h263_qscales(MpegEncContext *s){
266 int i;
267 int8_t * const qscale_table = s->current_picture.qscale_table;
268
269 ff_init_qscale_tab(s);
270
271 for(i=1; i<s->mb_num; i++){
272 if(qscale_table[ s->mb_index2xy[i] ] - qscale_table[ s->mb_index2xy[i-1] ] >2)
273 qscale_table[ s->mb_index2xy[i] ]= qscale_table[ s->mb_index2xy[i-1] ]+2;
274 }
275 for(i=s->mb_num-2; i>=0; i--){
276 if(qscale_table[ s->mb_index2xy[i] ] - qscale_table[ s->mb_index2xy[i+1] ] >2)
277 qscale_table[ s->mb_index2xy[i] ]= qscale_table[ s->mb_index2xy[i+1] ]+2;
278 }
279
280 if(s->codec_id != AV_CODEC_ID_H263P){
281 for(i=1; i<s->mb_num; i++){
282 int mb_xy= s->mb_index2xy[i];
283
284 if(qscale_table[mb_xy] != qscale_table[s->mb_index2xy[i-1]] && (s->mb_type[mb_xy]&CANDIDATE_MB_TYPE_INTER4V)){
285 s->mb_type[mb_xy]|= CANDIDATE_MB_TYPE_INTER;
286 }
287 }
288 }
289 }
290
291 static const int dquant_code[5]= {1,0,9,2,3};
292
293 /**
294 * Encode an 8x8 block.
295 * @param block the 8x8 block
296 * @param n block index (0-3 are luma, 4-5 are chroma)
297 */
298 static void h263_encode_block(MpegEncContext * s, int16_t * block, int n)
299 {
300 int level, run, last, i, j, last_index, last_non_zero, sign, slevel, code;
301 RLTable *rl;
302
303 rl = &ff_h263_rl_inter;
304 if (s->mb_intra && !s->h263_aic) {
305 /* DC coef */
306 level = block[0];
307 /* 255 cannot be represented, so we clamp */
308 if (level > 254) {
309 level = 254;
310 block[0] = 254;
311 }
312 /* 0 cannot be represented also */
313 else if (level < 1) {
314 level = 1;
315 block[0] = 1;
316 }
317 if (level == 128) //FIXME check rv10
318 put_bits(&s->pb, 8, 0xff);
319 else
320 put_bits(&s->pb, 8, level);
321 i = 1;
322 } else {
323 i = 0;
324 if (s->h263_aic && s->mb_intra)
325 rl = &ff_rl_intra_aic;
326
327 if(s->alt_inter_vlc && !s->mb_intra){
328 int aic_vlc_bits=0;
329 int inter_vlc_bits=0;
330 int wrong_pos=-1;
331 int aic_code;
332
333 last_index = s->block_last_index[n];
334 last_non_zero = i - 1;
335 for (; i <= last_index; i++) {
336 j = s->intra_scantable.permutated[i];
337 level = block[j];
338 if (level) {
339 run = i - last_non_zero - 1;
340 last = (i == last_index);
341
342 if(level<0) level= -level;
343
344 code = get_rl_index(rl, last, run, level);
345 aic_code = get_rl_index(&ff_rl_intra_aic, last, run, level);
346 inter_vlc_bits += rl->table_vlc[code][1]+1;
347 aic_vlc_bits += ff_rl_intra_aic.table_vlc[aic_code][1]+1;
348
349 if (code == rl->n) {
350 inter_vlc_bits += 1+6+8-1;
351 }
352 if (aic_code == ff_rl_intra_aic.n) {
353 aic_vlc_bits += 1+6+8-1;
354 wrong_pos += run + 1;
355 }else
356 wrong_pos += wrong_run[aic_code];
357 last_non_zero = i;
358 }
359 }
360 i = 0;
361 if(aic_vlc_bits < inter_vlc_bits && wrong_pos > 63)
362 rl = &ff_rl_intra_aic;
363 }
364 }
365
366 /* AC coefs */
367 last_index = s->block_last_index[n];
368 last_non_zero = i - 1;
369 for (; i <= last_index; i++) {
370 j = s->intra_scantable.permutated[i];
371 level = block[j];
372 if (level) {
373 run = i - last_non_zero - 1;
374 last = (i == last_index);
375 sign = 0;
376 slevel = level;
377 if (level < 0) {
378 sign = 1;
379 level = -level;
380 }
381 code = get_rl_index(rl, last, run, level);
382 put_bits(&s->pb, rl->table_vlc[code][1], rl->table_vlc[code][0]);
383 if (code == rl->n) {
384 if(!CONFIG_FLV_ENCODER || s->h263_flv <= 1){
385 put_bits(&s->pb, 1, last);
386 put_bits(&s->pb, 6, run);
387
388 av_assert2(slevel != 0);
389
390 if(level < 128)
391 put_sbits(&s->pb, 8, slevel);
392 else{
393 put_bits(&s->pb, 8, 128);
394 put_sbits(&s->pb, 5, slevel);
395 put_sbits(&s->pb, 6, slevel>>5);
396 }
397 }else{
398 ff_flv2_encode_ac_esc(&s->pb, slevel, level, run, last);
399 }
400 } else {
401 put_bits(&s->pb, 1, sign);
402 }
403 last_non_zero = i;
404 }
405 }
406 }
407
408 /* Encode MV differences on H.263+ with Unrestricted MV mode */
409 static void h263p_encode_umotion(MpegEncContext * s, int val)
410 {
411 short sval = 0;
412 short i = 0;
413 short n_bits = 0;
414 short temp_val;
415 int code = 0;
416 int tcode;
417
418 if ( val == 0)
419 put_bits(&s->pb, 1, 1);
420 else if (val == 1)
421 put_bits(&s->pb, 3, 0);
422 else if (val == -1)
423 put_bits(&s->pb, 3, 2);
424 else {
425
426 sval = ((val < 0) ? (short)(-val):(short)val);
427 temp_val = sval;
428
429 while (temp_val != 0) {
430 temp_val = temp_val >> 1;
431 n_bits++;
432 }
433
434 i = n_bits - 1;
435 while (i > 0) {
436 tcode = (sval & (1 << (i-1))) >> (i-1);
437 tcode = (tcode << 1) | 1;
438 code = (code << 2) | tcode;
439 i--;
440 }
441 code = ((code << 1) | (val < 0)) << 1;
442 put_bits(&s->pb, (2*n_bits)+1, code);
443 }
444 }
445
446 void ff_h263_encode_mb(MpegEncContext * s,
447 int16_t block[6][64],
448 int motion_x, int motion_y)
449 {
450 int cbpc, cbpy, i, cbp, pred_x, pred_y;
451 int16_t pred_dc;
452 int16_t rec_intradc[6];
453 int16_t *dc_ptr[6];
454 const int interleaved_stats= (s->flags&CODEC_FLAG_PASS1);
455
456 if (!s->mb_intra) {
457 /* compute cbp */
458 cbp= get_p_cbp(s, block, motion_x, motion_y);
459
460 if ((cbp | motion_x | motion_y | s->dquant | (s->mv_type - MV_TYPE_16X16)) == 0) {
461 /* skip macroblock */
462 put_bits(&s->pb, 1, 1);
463 if(interleaved_stats){
464 s->misc_bits++;
465 s->last_bits++;
466 }
467 s->skip_count++;
468
469 return;
470 }
471 put_bits(&s->pb, 1, 0); /* mb coded */
472
473 cbpc = cbp & 3;
474 cbpy = cbp >> 2;
475 if(s->alt_inter_vlc==0 || cbpc!=3)
476 cbpy ^= 0xF;
477 if(s->dquant) cbpc+= 8;
478 if(s->mv_type==MV_TYPE_16X16){
479 put_bits(&s->pb,
480 ff_h263_inter_MCBPC_bits[cbpc],
481 ff_h263_inter_MCBPC_code[cbpc]);
482
483 put_bits(&s->pb, ff_h263_cbpy_tab[cbpy][1], ff_h263_cbpy_tab[cbpy][0]);
484 if(s->dquant)
485 put_bits(&s->pb, 2, dquant_code[s->dquant+2]);
486
487 if(interleaved_stats){
488 s->misc_bits+= get_bits_diff(s);
489 }
490
491 /* motion vectors: 16x16 mode */
492 ff_h263_pred_motion(s, 0, 0, &pred_x, &pred_y);
493
494 if (!s->umvplus) {
495 ff_h263_encode_motion_vector(s, motion_x - pred_x,
496 motion_y - pred_y, 1);
497 }
498 else {
499 h263p_encode_umotion(s, motion_x - pred_x);
500 h263p_encode_umotion(s, motion_y - pred_y);
501 if (((motion_x - pred_x) == 1) && ((motion_y - pred_y) == 1))
502 /* To prevent Start Code emulation */
503 put_bits(&s->pb,1,1);
504 }
505 }else{
506 put_bits(&s->pb,
507 ff_h263_inter_MCBPC_bits[cbpc+16],
508 ff_h263_inter_MCBPC_code[cbpc+16]);
509 put_bits(&s->pb, ff_h263_cbpy_tab[cbpy][1], ff_h263_cbpy_tab[cbpy][0]);
510 if(s->dquant)
511 put_bits(&s->pb, 2, dquant_code[s->dquant+2]);
512
513 if(interleaved_stats){
514 s->misc_bits+= get_bits_diff(s);
515 }
516
517 for(i=0; i<4; i++){
518 /* motion vectors: 8x8 mode*/
519 ff_h263_pred_motion(s, i, 0, &pred_x, &pred_y);
520
521 motion_x = s->current_picture.motion_val[0][s->block_index[i]][0];
522 motion_y = s->current_picture.motion_val[0][s->block_index[i]][1];
523 if (!s->umvplus) {
524 ff_h263_encode_motion_vector(s, motion_x - pred_x,
525 motion_y - pred_y, 1);
526 }
527 else {
528 h263p_encode_umotion(s, motion_x - pred_x);
529 h263p_encode_umotion(s, motion_y - pred_y);
530 if (((motion_x - pred_x) == 1) && ((motion_y - pred_y) == 1))
531 /* To prevent Start Code emulation */
532 put_bits(&s->pb,1,1);
533 }
534 }
535 }
536
537 if(interleaved_stats){
538 s->mv_bits+= get_bits_diff(s);
539 }
540 } else {
541 av_assert2(s->mb_intra);
542
543 cbp = 0;
544 if (s->h263_aic) {
545 /* Predict DC */
546 for(i=0; i<6; i++) {
547 int16_t level = block[i][0];
548 int scale;
549
550 if(i<4) scale= s->y_dc_scale;
551 else scale= s->c_dc_scale;
552
553 pred_dc = ff_h263_pred_dc(s, i, &dc_ptr[i]);
554 level -= pred_dc;
555 /* Quant */
556 if (level >= 0)
557 level = (level + (scale>>1))/scale;
558 else
559 level = (level - (scale>>1))/scale;
560
561 if(!s->modified_quant){
562 if (level < -127)
563 level = -127;
564 else if (level > 127)
565 level = 127;
566 }
567
568 block[i][0] = level;
569 /* Reconstruction */
570 rec_intradc[i] = scale*level + pred_dc;
571 /* Oddify */
572 rec_intradc[i] |= 1;
573 //if ((rec_intradc[i] % 2) == 0)
574 // rec_intradc[i]++;
575 /* Clipping */
576 if (rec_intradc[i] < 0)
577 rec_intradc[i] = 0;
578 else if (rec_intradc[i] > 2047)
579 rec_intradc[i] = 2047;
580
581 /* Update AC/DC tables */
582 *dc_ptr[i] = rec_intradc[i];
583 /* AIC can change CBP */
584 if (s->block_last_index[i] > 0 ||
585 (s->block_last_index[i] == 0 && level !=0))
586 cbp |= 1 << (5 - i);
587 }
588 }else{
589 for(i=0; i<6; i++) {
590 /* compute cbp */
591 if (s->block_last_index[i] >= 1)
592 cbp |= 1 << (5 - i);
593 }
594 }
595
596 cbpc = cbp & 3;
597 if (s->pict_type == AV_PICTURE_TYPE_I) {
598 if(s->dquant) cbpc+=4;
599 put_bits(&s->pb,
600 ff_h263_intra_MCBPC_bits[cbpc],
601 ff_h263_intra_MCBPC_code[cbpc]);
602 } else {
603 if(s->dquant) cbpc+=8;
604 put_bits(&s->pb, 1, 0); /* mb coded */
605 put_bits(&s->pb,
606 ff_h263_inter_MCBPC_bits[cbpc + 4],
607 ff_h263_inter_MCBPC_code[cbpc + 4]);
608 }
609 if (s->h263_aic) {
610 /* XXX: currently, we do not try to use ac prediction */
611 put_bits(&s->pb, 1, 0); /* no AC prediction */
612 }
613 cbpy = cbp >> 2;
614 put_bits(&s->pb, ff_h263_cbpy_tab[cbpy][1], ff_h263_cbpy_tab[cbpy][0]);
615 if(s->dquant)
616 put_bits(&s->pb, 2, dquant_code[s->dquant+2]);
617
618 if(interleaved_stats){
619 s->misc_bits+= get_bits_diff(s);
620 }
621 }
622
623 for(i=0; i<6; i++) {
624 /* encode each block */
625 h263_encode_block(s, block[i], i);
626
627 /* Update INTRADC for decoding */
628 if (s->h263_aic && s->mb_intra) {
629 block[i][0] = rec_intradc[i];
630
631 }
632 }
633
634 if(interleaved_stats){
635 if (!s->mb_intra) {
636 s->p_tex_bits+= get_bits_diff(s);
637 s->f_count++;
638 }else{
639 s->i_tex_bits+= get_bits_diff(s);
640 s->i_count++;
641 }
642 }
643 }
644
645 void ff_h263_encode_motion(MpegEncContext * s, int val, int f_code)
646 {
647 int range, bit_size, sign, code, bits;
648
649 if (val == 0) {
650 /* zero vector */
651 code = 0;
652 put_bits(&s->pb, ff_mvtab[code][1], ff_mvtab[code][0]);
653 } else {
654 bit_size = f_code - 1;
655 range = 1 << bit_size;
656 /* modulo encoding */
657 val = sign_extend(val, 6 + bit_size);
658 sign = val>>31;
659 val= (val^sign)-sign;
660 sign&=1;
661
662 val--;
663 code = (val >> bit_size) + 1;
664 bits = val & (range - 1);
665
666 put_bits(&s->pb, ff_mvtab[code][1] + 1, (ff_mvtab[code][0] << 1) | sign);
667 if (bit_size > 0) {
668 put_bits(&s->pb, bit_size, bits);
669 }
670 }
671 }
672
673 static av_cold void init_mv_penalty_and_fcode(MpegEncContext *s)
674 {
675 int f_code;
676 int mv;
677
678 for(f_code=1; f_code<=MAX_FCODE; f_code++){
679 for(mv=-MAX_MV; mv<=MAX_MV; mv++){
680 int len;
681
682 if(mv==0) len= ff_mvtab[0][1];
683 else{
684 int val, bit_size, code;
685
686 bit_size = f_code - 1;
687
688 val=mv;
689 if (val < 0)
690 val = -val;
691 val--;
692 code = (val >> bit_size) + 1;
693 if(code<33){
694 len= ff_mvtab[code][1] + 1 + bit_size;
695 }else{
696 len= ff_mvtab[32][1] + av_log2(code>>5) + 2 + bit_size;
697 }
698 }
699
700 mv_penalty[f_code][mv+MAX_MV]= len;
701 }
702 }
703
704 for(f_code=MAX_FCODE; f_code>0; f_code--){
705 for(mv=-(16<<f_code); mv<(16<<f_code); mv++){
706 fcode_tab[mv+MAX_MV]= f_code;
707 }
708 }
709
710 for(mv=0; mv<MAX_MV*2+1; mv++){
711 umv_fcode_tab[mv]= 1;
712 }
713 }
714
715 static av_cold void init_uni_h263_rl_tab(RLTable *rl, uint32_t *bits_tab,
716 uint8_t *len_tab)
717 {
718 int slevel, run, last;
719
720 av_assert0(MAX_LEVEL >= 64);
721 av_assert0(MAX_RUN >= 63);
722
723 for(slevel=-64; slevel<64; slevel++){
724 if(slevel==0) continue;
725 for(run=0; run<64; run++){
726 for(last=0; last<=1; last++){
727 const int index= UNI_MPEG4_ENC_INDEX(last, run, slevel+64);
728 int level= slevel < 0 ? -slevel : slevel;
729 int sign= slevel < 0 ? 1 : 0;
730 int bits, len, code;
731
732 len_tab[index]= 100;
733
734 /* ESC0 */
735 code= get_rl_index(rl, last, run, level);
736 bits= rl->table_vlc[code][0];
737 len= rl->table_vlc[code][1];
738 bits=bits*2+sign; len++;
739
740 if(code!=rl->n && len < len_tab[index]){
741 if(bits_tab) bits_tab[index]= bits;
742 len_tab [index]= len;
743 }
744 /* ESC */
745 bits= rl->table_vlc[rl->n][0];
746 len = rl->table_vlc[rl->n][1];
747 bits=bits*2+last; len++;
748 bits=bits*64+run; len+=6;
749 bits=bits*256+(level&0xff); len+=8;
750
751 if(len < len_tab[index]){
752 if(bits_tab) bits_tab[index]= bits;
753 len_tab [index]= len;
754 }
755 }
756 }
757 }
758 }
759
760 av_cold void ff_h263_encode_init(MpegEncContext *s)
761 {
762 static int done = 0;
763
764 if (!done) {
765 done = 1;
766
767 ff_init_rl(&ff_h263_rl_inter, ff_h263_static_rl_table_store[0]);
768 ff_init_rl(&ff_rl_intra_aic, ff_h263_static_rl_table_store[1]);
769
770 init_uni_h263_rl_tab(&ff_rl_intra_aic, NULL, uni_h263_intra_aic_rl_len);
771 init_uni_h263_rl_tab(&ff_h263_rl_inter , NULL, uni_h263_inter_rl_len);
772
773 init_mv_penalty_and_fcode(s);
774 }
775 s->me.mv_penalty= mv_penalty; //FIXME exact table for msmpeg4 & h263p
776
777 s->intra_ac_vlc_length =s->inter_ac_vlc_length = uni_h263_inter_rl_len;
778 s->intra_ac_vlc_last_length=s->inter_ac_vlc_last_length= uni_h263_inter_rl_len + 128*64;
779 if(s->h263_aic){
780 s->intra_ac_vlc_length = uni_h263_intra_aic_rl_len;
781 s->intra_ac_vlc_last_length= uni_h263_intra_aic_rl_len + 128*64;
782 }
783 s->ac_esc_length= 7+1+6+8;
784
785 // use fcodes >1 only for mpeg4 & h263 & h263p FIXME
786 switch(s->codec_id){
787 case AV_CODEC_ID_MPEG4:
788 s->fcode_tab= fcode_tab;
789 break;
790 case AV_CODEC_ID_H263P:
791 if(s->umvplus)
792 s->fcode_tab= umv_fcode_tab;
793 if(s->modified_quant){
794 s->min_qcoeff= -2047;
795 s->max_qcoeff= 2047;
796 }else{
797 s->min_qcoeff= -127;
798 s->max_qcoeff= 127;
799 }
800 break;
801 //Note for mpeg4 & h263 the dc-scale table will be set per frame as needed later
802 case AV_CODEC_ID_FLV1:
803 if (s->h263_flv > 1) {
804 s->min_qcoeff= -1023;
805 s->max_qcoeff= 1023;
806 } else {
807 s->min_qcoeff= -127;
808 s->max_qcoeff= 127;
809 }
810 break;
811 default: //nothing needed - default table already set in mpegvideo.c
812 s->min_qcoeff= -127;
813 s->max_qcoeff= 127;
814 }
815 if(s->h263_aic){
816 s->y_dc_scale_table=
817 s->c_dc_scale_table= ff_aic_dc_scale_table;
818 }else{
819 s->y_dc_scale_table=
820 s->c_dc_scale_table= ff_mpeg1_dc_scale_table;
821 }
822 }
823
824 void ff_h263_encode_mba(MpegEncContext *s)
825 {
826 int i, mb_pos;
827
828 for(i=0; i<6; i++){
829 if(s->mb_num-1 <= ff_mba_max[i]) break;
830 }
831 mb_pos= s->mb_x + s->mb_width*s->mb_y;
832 put_bits(&s->pb, ff_mba_length[i], mb_pos);
833 }