2 * Copyright (c) 2010 The WebM project authors. All Rights Reserved.
4 * Use of this source code is governed by a BSD-style license
5 * that can be found in the LICENSE file in the root of the source
6 * tree. An additional intellectual property rights grant can be found
7 * in the file PATENTS. All contributing project authors may
8 * be found in the AUTHORS file in the root of the source tree.
13 #include "vpx_config.h"
15 #include "modecosts.h"
16 #include "encodeintra.h"
17 #include "vp8/common/common.h"
18 #include "vp8/common/entropymode.h"
19 #include "pickinter.h"
20 #include "vp8/common/findnearmv.h"
22 #include "vp8/common/reconinter.h"
23 #include "vp8/common/reconintra4x4.h"
24 #include "vp8/common/variance.h"
27 #include "vpx_mem/vpx_mem.h"
28 #if CONFIG_TEMPORAL_DENOISING
29 #include "denoising.h"
32 extern int VP8_UVSSE(MACROBLOCK *x);
35 extern unsigned int cnt_pm;
38 extern const int vp8_ref_frame_order[MAX_MODES];
39 extern const MB_PREDICTION_MODE vp8_mode_order[MAX_MODES];
41 extern int vp8_cost_mv_ref(MB_PREDICTION_MODE m, const int near_mv_ref_ct[4]);
44 int vp8_skip_fractional_mv_step(MACROBLOCK *mb, BLOCK *b, BLOCKD *d,
45 int_mv *bestmv, int_mv *ref_mv,
47 const vp8_variance_fn_ptr_t *vfp,
48 int *mvcost[2], int *distortion,
59 bestmv->as_mv.row <<= 3;
60 bestmv->as_mv.col <<= 3;
65 int vp8_get_inter_mbpred_error(MACROBLOCK *mb,
66 const vp8_variance_fn_ptr_t *vfp,
71 BLOCK *b = &mb->block[0];
72 BLOCKD *d = &mb->e_mbd.block[0];
73 unsigned char *what = (*(b->base_src) + b->src);
74 int what_stride = b->src_stride;
75 int pre_stride = mb->e_mbd.pre.y_stride;
76 unsigned char *in_what = mb->e_mbd.pre.y_buffer + d->offset ;
77 int in_what_stride = pre_stride;
78 int xoffset = this_mv.as_mv.col & 7;
79 int yoffset = this_mv.as_mv.row & 7;
81 in_what += (this_mv.as_mv.row >> 3) * pre_stride + (this_mv.as_mv.col >> 3);
83 if (xoffset | yoffset)
85 return vfp->svf(in_what, in_what_stride, xoffset, yoffset, what, what_stride, sse);
89 return vfp->vf(what, what_stride, in_what, in_what_stride, sse);
95 unsigned int vp8_get4x4sse_cs_c
97 const unsigned char *src_ptr,
99 const unsigned char *ref_ptr,
106 for (r = 0; r < 4; r++)
108 for (c = 0; c < 4; c++)
110 int diff = src_ptr[c] - ref_ptr[c];
111 distortion += diff * diff;
114 src_ptr += source_stride;
115 ref_ptr += recon_stride;
121 static int get_prediction_error(BLOCK *be, BLOCKD *b)
125 sptr = (*(be->base_src) + be->src);
128 return vp8_get4x4sse_cs(sptr, be->src_stride, dptr, 16);
132 static int pick_intra4x4block(
135 B_PREDICTION_MODE *best_mode,
136 const int *mode_costs,
142 BLOCKD *b = &x->e_mbd.block[ib];
143 BLOCK *be = &x->block[ib];
144 int dst_stride = x->e_mbd.dst.y_stride;
145 unsigned char *dst = x->e_mbd.dst.y_buffer + b->offset;
146 B_PREDICTION_MODE mode;
147 int best_rd = INT_MAX;
151 unsigned char *Above = dst - dst_stride;
152 unsigned char *yleft = dst - 1;
153 unsigned char top_left = Above[-1];
155 for (mode = B_DC_PRED; mode <= B_HE_PRED; mode++)
159 rate = mode_costs[mode];
161 vp8_intra4x4_predict(Above, yleft, dst_stride, mode,
162 b->predictor, 16, top_left);
163 distortion = get_prediction_error(be, b);
164 this_rd = RDCOST(x->rdmult, x->rddiv, rate, distortion);
166 if (this_rd < best_rd)
169 *bestdistortion = distortion;
175 b->bmi.as_mode = *best_mode;
176 vp8_encode_intra4x4block(x, ib);
181 static int pick_intra4x4mby_modes
188 MACROBLOCKD *const xd = &mb->e_mbd;
190 int cost = mb->mbmode_cost [xd->frame_type] [B_PRED];
193 const int *bmode_costs;
195 intra_prediction_down_copy(xd, xd->dst.y_buffer - xd->dst.y_stride + 16);
197 bmode_costs = mb->inter_bmode_costs;
199 for (i = 0; i < 16; i++)
201 MODE_INFO *const mic = xd->mode_info_context;
202 const int mis = xd->mode_info_stride;
204 B_PREDICTION_MODE UNINITIALIZED_IS_SAFE(best_mode);
205 int UNINITIALIZED_IS_SAFE(r), UNINITIALIZED_IS_SAFE(d);
207 if (mb->e_mbd.frame_type == KEY_FRAME)
209 const B_PREDICTION_MODE A = above_block_mode(mic, i, mis);
210 const B_PREDICTION_MODE L = left_block_mode(mic, i);
212 bmode_costs = mb->bmode_costs[A][L];
216 pick_intra4x4block(mb, i, &best_mode, bmode_costs, &r, &d);
220 mic->bmi[i].as_mode = best_mode;
222 /* Break out case where we have already exceeded best so far value
225 if (distortion > *best_dist)
233 *best_dist = distortion;
234 error = RDCOST(mb->rdmult, mb->rddiv, cost, distortion);
238 *best_dist = INT_MAX;
245 static void pick_intra_mbuv_mode(MACROBLOCK *mb)
248 MACROBLOCKD *x = &mb->e_mbd;
249 unsigned char *uabove_row = x->dst.u_buffer - x->dst.uv_stride;
250 unsigned char *vabove_row = x->dst.v_buffer - x->dst.uv_stride;
251 unsigned char *usrc_ptr = (mb->block[16].src + *mb->block[16].base_src);
252 unsigned char *vsrc_ptr = (mb->block[20].src + *mb->block[20].base_src);
253 int uvsrc_stride = mb->block[16].src_stride;
254 unsigned char uleft_col[8];
255 unsigned char vleft_col[8];
256 unsigned char utop_left = uabove_row[-1];
257 unsigned char vtop_left = vabove_row[-1];
265 int pred_error[4] = {0, 0, 0, 0}, best_error = INT_MAX;
266 MB_PREDICTION_MODE UNINITIALIZED_IS_SAFE(best_mode);
269 for (i = 0; i < 8; i++)
271 uleft_col[i] = x->dst.u_buffer [i* x->dst.uv_stride -1];
272 vleft_col[i] = x->dst.v_buffer [i* x->dst.uv_stride -1];
275 if (!x->up_available && !x->left_available)
287 for (i = 0; i < 8; i++)
289 Uaverage += uabove_row[i];
290 Vaverage += vabove_row[i];
297 if (x->left_available)
299 for (i = 0; i < 8; i++)
301 Uaverage += uleft_col[i];
302 Vaverage += vleft_col[i];
309 expected_udc = (Uaverage + (1 << (shift - 1))) >> shift;
310 expected_vdc = (Vaverage + (1 << (shift - 1))) >> shift;
314 for (i = 0; i < 8; i++)
316 for (j = 0; j < 8; j++)
319 int predu = uleft_col[i] + uabove_row[j] - utop_left;
320 int predv = vleft_col[i] + vabove_row[j] - vtop_left;
339 diff = u_p - expected_udc;
340 pred_error[DC_PRED] += diff * diff;
341 diff = v_p - expected_vdc;
342 pred_error[DC_PRED] += diff * diff;
345 diff = u_p - uabove_row[j];
346 pred_error[V_PRED] += diff * diff;
347 diff = v_p - vabove_row[j];
348 pred_error[V_PRED] += diff * diff;
351 diff = u_p - uleft_col[i];
352 pred_error[H_PRED] += diff * diff;
353 diff = v_p - vleft_col[i];
354 pred_error[H_PRED] += diff * diff;
358 pred_error[TM_PRED] += diff * diff;
360 pred_error[TM_PRED] += diff * diff;
365 usrc_ptr += uvsrc_stride;
366 vsrc_ptr += uvsrc_stride;
370 usrc_ptr = (mb->block[18].src + *mb->block[18].base_src);
371 vsrc_ptr = (mb->block[22].src + *mb->block[22].base_src);
379 for (i = DC_PRED; i <= TM_PRED; i++)
381 if (best_error > pred_error[i])
383 best_error = pred_error[i];
384 best_mode = (MB_PREDICTION_MODE)i;
389 mb->e_mbd.mode_info_context->mbmi.uv_mode = best_mode;
393 static void update_mvcount(MACROBLOCK *x, int_mv *best_ref_mv)
395 MACROBLOCKD *xd = &x->e_mbd;
396 /* Split MV modes currently not supported when RD is nopt enabled,
397 * therefore, only need to modify MVcount in NEWMV mode. */
398 if (xd->mode_info_context->mbmi.mode == NEWMV)
400 x->MVcount[0][mv_max+((xd->mode_info_context->mbmi.mv.as_mv.row -
401 best_ref_mv->as_mv.row) >> 1)]++;
402 x->MVcount[1][mv_max+((xd->mode_info_context->mbmi.mv.as_mv.col -
403 best_ref_mv->as_mv.col) >> 1)]++;
408 #if CONFIG_MULTI_RES_ENCODING
410 void get_lower_res_motion_info(VP8_COMP *cpi, MACROBLOCKD *xd, int *dissim,
411 int *parent_ref_frame,
412 MB_PREDICTION_MODE *parent_mode,
413 int_mv *parent_ref_mv, int mb_row, int mb_col)
415 LOWER_RES_MB_INFO* store_mode_info
416 = ((LOWER_RES_FRAME_INFO*)cpi->oxcf.mr_low_res_mode_info)->mb_info;
417 unsigned int parent_mb_index;
419 /* Consider different down_sampling_factor. */
421 /* TODO: Removed the loop that supports special down_sampling_factor
422 * such as 2, 4, 8. Will revisit it if needed.
423 * Should also try using a look-up table to see if it helps
425 int parent_mb_row, parent_mb_col;
427 parent_mb_row = mb_row*cpi->oxcf.mr_down_sampling_factor.den
428 /cpi->oxcf.mr_down_sampling_factor.num;
429 parent_mb_col = mb_col*cpi->oxcf.mr_down_sampling_factor.den
430 /cpi->oxcf.mr_down_sampling_factor.num;
431 parent_mb_index = parent_mb_row*cpi->mr_low_res_mb_cols + parent_mb_col;
434 /* Read lower-resolution mode & motion result from memory.*/
435 *parent_ref_frame = store_mode_info[parent_mb_index].ref_frame;
436 *parent_mode = store_mode_info[parent_mb_index].mode;
437 *dissim = store_mode_info[parent_mb_index].dissim;
439 /* For highest-resolution encoder, adjust dissim value. Lower its quality
440 * for good performance. */
441 if (cpi->oxcf.mr_encoder_id == (cpi->oxcf.mr_total_resolutions - 1))
444 if(*parent_ref_frame != INTRA_FRAME)
446 /* Consider different down_sampling_factor.
447 * The result can be rounded to be more precise, but it takes more time.
449 (*parent_ref_mv).as_mv.row = store_mode_info[parent_mb_index].mv.as_mv.row
450 *cpi->oxcf.mr_down_sampling_factor.num
451 /cpi->oxcf.mr_down_sampling_factor.den;
452 (*parent_ref_mv).as_mv.col = store_mode_info[parent_mb_index].mv.as_mv.col
453 *cpi->oxcf.mr_down_sampling_factor.num
454 /cpi->oxcf.mr_down_sampling_factor.den;
456 vp8_clamp_mv2(parent_ref_mv, xd);
461 static void check_for_encode_breakout(unsigned int sse, MACROBLOCK* x)
463 MACROBLOCKD *xd = &x->e_mbd;
465 unsigned int threshold = (xd->block[0].dequant[1]
466 * xd->block[0].dequant[1] >>4);
468 if(threshold < x->encode_breakout)
469 threshold = x->encode_breakout;
471 if (sse < threshold )
473 /* Check u and v to make sure skip is ok */
474 unsigned int sse2 = 0;
478 if (sse2 * 2 < x->encode_breakout)
485 static int evaluate_inter_mode(unsigned int* sse, int rate2, int* distortion2,
486 VP8_COMP *cpi, MACROBLOCK *x, int rd_adj)
488 MB_PREDICTION_MODE this_mode = x->e_mbd.mode_info_context->mbmi.mode;
489 int_mv mv = x->e_mbd.mode_info_context->mbmi.mv;
491 /* Exit early and don't compute the distortion if this macroblock
492 * is marked inactive. */
493 if (cpi->active_map_enabled && x->active_ptr[0] == 0)
501 if((this_mode != NEWMV) ||
502 !(cpi->sf.half_pixel_search) || cpi->common.full_pixel==1)
503 *distortion2 = vp8_get_inter_mbpred_error(x,
504 &cpi->fn_ptr[BLOCK_16X16],
507 this_rd = RDCOST(x->rdmult, x->rddiv, rate2, *distortion2);
509 /* Adjust rd to bias to ZEROMV */
510 if(this_mode == ZEROMV)
512 /* Bias to ZEROMV on LAST_FRAME reference when it is available. */
513 if ((cpi->ref_frame_flags & VP8_LAST_FRAME &
514 cpi->common.refresh_last_frame)
515 && x->e_mbd.mode_info_context->mbmi.ref_frame != LAST_FRAME)
519 this_rd = ((int64_t)this_rd) * rd_adj / 100;
522 check_for_encode_breakout(*sse, x);
526 static void calculate_zeromv_rd_adjustment(VP8_COMP *cpi, MACROBLOCK *x,
529 MODE_INFO *mic = x->e_mbd.mode_info_context;
530 int_mv mv_l, mv_a, mv_al;
531 int local_motion_check = 0;
533 if (cpi->lf_zeromv_pct > 40)
539 if (mic->mbmi.ref_frame != INTRA_FRAME)
540 if( abs(mv_l.as_mv.row) < 8 && abs(mv_l.as_mv.col) < 8)
541 local_motion_check++;
544 mic -= x->e_mbd.mode_info_stride;
545 mv_al = mic->mbmi.mv;
547 if (mic->mbmi.ref_frame != INTRA_FRAME)
548 if( abs(mv_al.as_mv.row) < 8 && abs(mv_al.as_mv.col) < 8)
549 local_motion_check++;
555 if (mic->mbmi.ref_frame != INTRA_FRAME)
556 if( abs(mv_a.as_mv.row) < 8 && abs(mv_a.as_mv.col) < 8)
557 local_motion_check++;
559 if (((!x->e_mbd.mb_to_top_edge || !x->e_mbd.mb_to_left_edge)
560 && local_motion_check >0) || local_motion_check >2 )
562 else if (local_motion_check > 0)
567 void vp8_pick_inter_mode(VP8_COMP *cpi, MACROBLOCK *x, int recon_yoffset,
568 int recon_uvoffset, int *returnrate,
569 int *returndistortion, int *returnintra, int mb_row,
572 BLOCK *b = &x->block[0];
573 BLOCKD *d = &x->e_mbd.block[0];
574 MACROBLOCKD *xd = &x->e_mbd;
575 MB_MODE_INFO best_mbmode;
577 int_mv best_ref_mv_sb[2];
578 int_mv mode_mv_sb[2][MB_MODE_COUNT];
581 MB_PREDICTION_MODE this_mode;
584 int best_rd = INT_MAX;
585 int rd_adjustment = 100;
586 int best_intra_rd = INT_MAX;
591 int bestsme = INT_MAX;
592 int best_mode_index = 0;
593 unsigned int sse = INT_MAX, best_rd_sse = INT_MAX;
594 #if CONFIG_TEMPORAL_DENOISING
595 unsigned int zero_mv_sse = INT_MAX, best_sse = INT_MAX;
598 int sf_improved_mv_pred = cpi->sf.improved_mv_pred;
601 int near_sadidx[8] = {0, 1, 2, 3, 4, 5, 6, 7};
603 /* search range got from mv_pred(). It uses step_param levels. (0-7) */
606 unsigned char *plane[4][3];
607 int ref_frame_map[4];
610 #if CONFIG_MULTI_RES_ENCODING
611 int dissim = INT_MAX;
612 int parent_ref_frame = 0;
613 int parent_ref_valid = cpi->oxcf.mr_encoder_id && cpi->mr_low_res_mv_avail;
614 int_mv parent_ref_mv;
615 MB_PREDICTION_MODE parent_mode = 0;
617 if (parent_ref_valid)
621 get_lower_res_motion_info(cpi, xd, &dissim, &parent_ref_frame,
622 &parent_mode, &parent_ref_mv, mb_row, mb_col);
624 /* TODO(jkoleszar): The references available (ref_frame_flags) to the
625 * lower res encoder should match those available to this encoder, but
626 * there seems to be a situation where this mismatch can happen in the
627 * case of frame dropping and temporal layers. For example,
628 * GOLD being disallowed in ref_frame_flags, but being returned as
631 * In this event, take the conservative approach of disabling the
632 * lower res info for this MB.
635 if (parent_ref_frame == LAST_FRAME)
636 parent_ref_flag = (cpi->ref_frame_flags & VP8_LAST_FRAME);
637 else if (parent_ref_frame == GOLDEN_FRAME)
638 parent_ref_flag = (cpi->ref_frame_flags & VP8_GOLD_FRAME);
639 else if (parent_ref_frame == ALTREF_FRAME)
640 parent_ref_flag = (cpi->ref_frame_flags & VP8_ALTR_FRAME);
642 //assert(!parent_ref_frame || parent_ref_flag);
643 if (parent_ref_frame && !parent_ref_flag)
644 parent_ref_valid = 0;
648 mode_mv = mode_mv_sb[sign_bias];
649 best_ref_mv.as_int = 0;
650 vpx_memset(mode_mv_sb, 0, sizeof(mode_mv_sb));
651 vpx_memset(&best_mbmode, 0, sizeof(best_mbmode));
653 /* Setup search priorities */
654 #if CONFIG_MULTI_RES_ENCODING
655 if (parent_ref_valid && parent_ref_frame && dissim < 8)
657 ref_frame_map[0] = -1;
658 ref_frame_map[1] = parent_ref_frame;
659 ref_frame_map[2] = -1;
660 ref_frame_map[3] = -1;
663 get_reference_search_order(cpi, ref_frame_map);
665 /* Check to see if there is at least 1 valid reference frame that we need
666 * to calculate near_mvs.
668 if (ref_frame_map[1] > 0)
670 sign_bias = vp8_find_near_mvs_bias(&x->e_mbd,
671 x->e_mbd.mode_info_context,
676 cpi->common.ref_frame_sign_bias);
678 mode_mv = mode_mv_sb[sign_bias];
679 best_ref_mv.as_int = best_ref_mv_sb[sign_bias].as_int;
682 get_predictor_pointers(cpi, plane, recon_yoffset, recon_uvoffset);
684 /* Count of the number of MBs tested so far this frame */
685 x->mbs_tested_so_far++;
687 *returnintra = INT_MAX;
690 x->e_mbd.mode_info_context->mbmi.ref_frame = INTRA_FRAME;
692 /* If the frame has big static background and current MB is in low
693 * motion area, its mode decision is biased to ZEROMV mode.
695 calculate_zeromv_rd_adjustment(cpi, x, &rd_adjustment);
697 /* if we encode a new mv this is important
698 * find the best new motion vector
700 for (mode_index = 0; mode_index < MAX_MODES; mode_index++)
703 int this_rd = INT_MAX;
704 int this_ref_frame = ref_frame_map[vp8_ref_frame_order[mode_index]];
706 if (best_rd <= x->rd_threshes[mode_index])
709 if (this_ref_frame < 0)
712 x->e_mbd.mode_info_context->mbmi.ref_frame = this_ref_frame;
714 /* everything but intra */
715 if (x->e_mbd.mode_info_context->mbmi.ref_frame)
717 x->e_mbd.pre.y_buffer = plane[this_ref_frame][0];
718 x->e_mbd.pre.u_buffer = plane[this_ref_frame][1];
719 x->e_mbd.pre.v_buffer = plane[this_ref_frame][2];
721 if (sign_bias != cpi->common.ref_frame_sign_bias[this_ref_frame])
723 sign_bias = cpi->common.ref_frame_sign_bias[this_ref_frame];
724 mode_mv = mode_mv_sb[sign_bias];
725 best_ref_mv.as_int = best_ref_mv_sb[sign_bias].as_int;
728 #if CONFIG_MULTI_RES_ENCODING
729 if (parent_ref_valid)
731 if (vp8_mode_order[mode_index] == NEARESTMV &&
732 mode_mv[NEARESTMV].as_int ==0)
734 if (vp8_mode_order[mode_index] == NEARMV &&
735 mode_mv[NEARMV].as_int ==0)
738 if (vp8_mode_order[mode_index] == NEWMV && parent_mode == ZEROMV
739 && best_ref_mv.as_int==0)
741 else if(vp8_mode_order[mode_index] == NEWMV && dissim==0
742 && best_ref_mv.as_int==parent_ref_mv.as_int)
748 /* Check to see if the testing frequency for this mode is at its max
749 * If so then prevent it from being tested and increase the threshold
751 if (x->mode_test_hit_counts[mode_index] &&
752 (cpi->mode_check_freq[mode_index] > 1))
754 if (x->mbs_tested_so_far <= (cpi->mode_check_freq[mode_index] *
755 x->mode_test_hit_counts[mode_index]))
757 /* Increase the threshold for coding this mode to make it less
758 * likely to be chosen */
759 x->rd_thresh_mult[mode_index] += 4;
761 if (x->rd_thresh_mult[mode_index] > MAX_THRESHMULT)
762 x->rd_thresh_mult[mode_index] = MAX_THRESHMULT;
764 x->rd_threshes[mode_index] =
765 (cpi->rd_baseline_thresh[mode_index] >> 7) *
766 x->rd_thresh_mult[mode_index];
771 /* We have now reached the point where we are going to test the current
772 * mode so increment the counter for the number of times it has been
774 x->mode_test_hit_counts[mode_index] ++;
779 this_mode = vp8_mode_order[mode_index];
781 x->e_mbd.mode_info_context->mbmi.mode = this_mode;
782 x->e_mbd.mode_info_context->mbmi.uv_mode = DC_PRED;
784 /* Work out the cost assosciated with selecting the reference frame */
786 x->ref_frame_cost[x->e_mbd.mode_info_context->mbmi.ref_frame];
789 /* Only consider ZEROMV/ALTREF_FRAME for alt ref frame,
790 * unless ARNR filtering is enabled in which case we want
791 * an unfiltered alternative */
792 if (cpi->is_src_frame_alt_ref && (cpi->oxcf.arnr_max_frames == 0))
794 if (this_mode != ZEROMV ||
795 x->e_mbd.mode_info_context->mbmi.ref_frame != ALTREF_FRAME)
802 /* Pass best so far to pick_intra4x4mby_modes to use as breakout */
803 distortion2 = best_rd_sse;
804 pick_intra4x4mby_modes(x, &rate, &distortion2);
806 if (distortion2 == INT_MAX)
813 distortion2 = vp8_variance16x16(
814 *(b->base_src), b->src_stride,
815 x->e_mbd.predictor, 16, &sse);
816 this_rd = RDCOST(x->rdmult, x->rddiv, rate2, distortion2);
818 if (this_rd < best_intra_rd)
820 best_intra_rd = this_rd;
821 *returnintra = distortion2;
829 /* Split MV modes currently not supported when RD is not enabled. */
836 vp8_build_intra_predictors_mby_s(xd,
837 xd->dst.y_buffer - xd->dst.y_stride,
838 xd->dst.y_buffer - 1,
842 distortion2 = vp8_variance16x16
843 (*(b->base_src), b->src_stride,
844 x->e_mbd.predictor, 16, &sse);
845 rate2 += x->mbmode_cost[x->e_mbd.frame_type][x->e_mbd.mode_info_context->mbmi.mode];
846 this_rd = RDCOST(x->rdmult, x->rddiv, rate2, distortion2);
848 if (this_rd < best_intra_rd)
850 best_intra_rd = this_rd;
851 *returnintra = distortion2;
861 int sadpb = x->sadperbit16;
864 int col_min = ((best_ref_mv.as_mv.col+7)>>3) - MAX_FULL_PEL_VAL;
865 int row_min = ((best_ref_mv.as_mv.row+7)>>3) - MAX_FULL_PEL_VAL;
866 int col_max = (best_ref_mv.as_mv.col>>3)
868 int row_max = (best_ref_mv.as_mv.row>>3)
871 int tmp_col_min = x->mv_col_min;
872 int tmp_col_max = x->mv_col_max;
873 int tmp_row_min = x->mv_row_min;
874 int tmp_row_max = x->mv_row_max;
876 int speed_adjust = (cpi->Speed > 5) ? ((cpi->Speed >= 8)? 3 : 2) : 1;
878 /* Further step/diamond searches as necessary */
879 step_param = cpi->sf.first_step + speed_adjust;
881 #if CONFIG_MULTI_RES_ENCODING
882 /* If lower-res drops this frame, then higher-res encoder does
883 motion search without any previous knowledge. Also, since
884 last frame motion info is not stored, then we can not
885 use improved_mv_pred. */
886 if (cpi->oxcf.mr_encoder_id && !parent_ref_valid)
887 sf_improved_mv_pred = 0;
889 if (parent_ref_valid && parent_ref_frame)
891 /* Use parent MV as predictor. Adjust search range
894 mvp.as_int = parent_ref_mv.as_int;
895 mvp_full.as_mv.col = parent_ref_mv.as_mv.col>>3;
896 mvp_full.as_mv.row = parent_ref_mv.as_mv.row>>3;
898 if(dissim <=32) step_param += 3;
899 else if(dissim <=128) step_param += 2;
900 else step_param += 1;
904 if(sf_improved_mv_pred)
908 vp8_cal_sad(cpi,xd,x, recon_yoffset ,&near_sadidx[0] );
912 vp8_mv_pred(cpi, &x->e_mbd, x->e_mbd.mode_info_context,
913 &mvp,x->e_mbd.mode_info_context->mbmi.ref_frame,
914 cpi->common.ref_frame_sign_bias, &sr,
918 /* adjust search range according to sr from mv prediction */
922 mvp_full.as_mv.col = mvp.as_mv.col>>3;
923 mvp_full.as_mv.row = mvp.as_mv.row>>3;
926 mvp.as_int = best_ref_mv.as_int;
927 mvp_full.as_mv.col = best_ref_mv.as_mv.col>>3;
928 mvp_full.as_mv.row = best_ref_mv.as_mv.row>>3;
932 #if CONFIG_MULTI_RES_ENCODING
933 if (parent_ref_valid && parent_ref_frame && dissim <= 2 &&
934 MAX(abs(best_ref_mv.as_mv.row - parent_ref_mv.as_mv.row),
935 abs(best_ref_mv.as_mv.col - parent_ref_mv.as_mv.col)) <= 4)
937 d->bmi.mv.as_int = mvp_full.as_int;
938 mode_mv[NEWMV].as_int = mvp_full.as_int;
940 cpi->find_fractional_mv_step(x, b, d, &d->bmi.mv, &best_ref_mv,
942 &cpi->fn_ptr[BLOCK_16X16],
948 /* Get intersection of UMV window and valid MV window to
949 * reduce # of checks in diamond search. */
950 if (x->mv_col_min < col_min )
951 x->mv_col_min = col_min;
952 if (x->mv_col_max > col_max )
953 x->mv_col_max = col_max;
954 if (x->mv_row_min < row_min )
955 x->mv_row_min = row_min;
956 if (x->mv_row_max > row_max )
957 x->mv_row_max = row_max;
959 further_steps = (cpi->Speed >= 8)?
960 0: (cpi->sf.max_step_search_steps - 1 - step_param);
962 if (cpi->sf.search_method == HEX)
964 #if CONFIG_MULTI_RES_ENCODING
965 /* TODO: In higher-res pick_inter_mode, step_param is used to
966 * modify hex search range. Here, set step_param to 0 not to
967 * change the behavior in lowest-resolution encoder.
968 * Will improve it later.
970 /* Set step_param to 0 to ensure large-range motion search
971 when encoder drops this frame at lower-resolution.
973 if (!parent_ref_valid)
976 bestsme = vp8_hex_search(x, b, d, &mvp_full, &d->bmi.mv,
978 &cpi->fn_ptr[BLOCK_16X16],
979 x->mvsadcost, x->mvcost, &best_ref_mv);
980 mode_mv[NEWMV].as_int = d->bmi.mv.as_int;
984 bestsme = cpi->diamond_search_sad(x, b, d, &mvp_full,
985 &d->bmi.mv, step_param, sadpb, &num00,
986 &cpi->fn_ptr[BLOCK_16X16],
987 x->mvcost, &best_ref_mv);
988 mode_mv[NEWMV].as_int = d->bmi.mv.as_int;
990 /* Further step/diamond searches as necessary */
994 while (n < further_steps)
1003 cpi->diamond_search_sad(x, b, d, &mvp_full,
1007 &cpi->fn_ptr[BLOCK_16X16],
1008 x->mvcost, &best_ref_mv);
1009 if (thissme < bestsme)
1012 mode_mv[NEWMV].as_int = d->bmi.mv.as_int;
1016 d->bmi.mv.as_int = mode_mv[NEWMV].as_int;
1022 x->mv_col_min = tmp_col_min;
1023 x->mv_col_max = tmp_col_max;
1024 x->mv_row_min = tmp_row_min;
1025 x->mv_row_max = tmp_row_max;
1027 if (bestsme < INT_MAX)
1028 cpi->find_fractional_mv_step(x, b, d, &d->bmi.mv,
1029 &best_ref_mv, x->errorperbit,
1030 &cpi->fn_ptr[BLOCK_16X16],
1035 mode_mv[NEWMV].as_int = d->bmi.mv.as_int;
1038 rate2 += vp8_mv_bit_cost(&mode_mv[NEWMV], &best_ref_mv,
1039 cpi->mb.mvcost, 128);
1045 if (mode_mv[this_mode].as_int == 0)
1050 /* Trap vectors that reach beyond the UMV borders
1051 * Note that ALL New MV, Nearest MV Near MV and Zero MV code drops
1052 * through to this point because of the lack of break statements
1053 * in the previous two cases.
1055 if (((mode_mv[this_mode].as_mv.row >> 3) < x->mv_row_min) ||
1056 ((mode_mv[this_mode].as_mv.row >> 3) > x->mv_row_max) ||
1057 ((mode_mv[this_mode].as_mv.col >> 3) < x->mv_col_min) ||
1058 ((mode_mv[this_mode].as_mv.col >> 3) > x->mv_col_max))
1061 rate2 += vp8_cost_mv_ref(this_mode, mdcounts);
1062 x->e_mbd.mode_info_context->mbmi.mv.as_int =
1063 mode_mv[this_mode].as_int;
1064 this_rd = evaluate_inter_mode(&sse, rate2, &distortion2, cpi, x,
1072 #if CONFIG_TEMPORAL_DENOISING
1073 if (cpi->oxcf.noise_sensitivity)
1076 /* Store for later use by denoiser. */
1077 if (this_mode == ZEROMV && sse < zero_mv_sse )
1080 x->best_zeromv_reference_frame =
1081 x->e_mbd.mode_info_context->mbmi.ref_frame;
1084 /* Store the best NEWMV in x for later use in the denoiser. */
1085 if (x->e_mbd.mode_info_context->mbmi.mode == NEWMV &&
1089 x->best_sse_inter_mode = NEWMV;
1090 x->best_sse_mv = x->e_mbd.mode_info_context->mbmi.mv;
1091 x->need_to_clamp_best_mvs =
1092 x->e_mbd.mode_info_context->mbmi.need_to_clamp_mvs;
1093 x->best_reference_frame =
1094 x->e_mbd.mode_info_context->mbmi.ref_frame;
1099 if (this_rd < best_rd || x->skip)
1101 /* Note index of best mode */
1102 best_mode_index = mode_index;
1104 *returnrate = rate2;
1105 *returndistortion = distortion2;
1108 vpx_memcpy(&best_mbmode, &x->e_mbd.mode_info_context->mbmi,
1109 sizeof(MB_MODE_INFO));
1111 /* Testing this mode gave rise to an improvement in best error
1112 * score. Lower threshold a bit for next time
1114 x->rd_thresh_mult[mode_index] =
1115 (x->rd_thresh_mult[mode_index] >= (MIN_THRESHMULT + 2)) ?
1116 x->rd_thresh_mult[mode_index] - 2 : MIN_THRESHMULT;
1117 x->rd_threshes[mode_index] =
1118 (cpi->rd_baseline_thresh[mode_index] >> 7) *
1119 x->rd_thresh_mult[mode_index];
1122 /* If the mode did not help improve the best error case then raise the
1123 * threshold for testing that mode next time around.
1127 x->rd_thresh_mult[mode_index] += 4;
1129 if (x->rd_thresh_mult[mode_index] > MAX_THRESHMULT)
1130 x->rd_thresh_mult[mode_index] = MAX_THRESHMULT;
1132 x->rd_threshes[mode_index] =
1133 (cpi->rd_baseline_thresh[mode_index] >> 7) *
1134 x->rd_thresh_mult[mode_index];
1141 /* Reduce the activation RD thresholds for the best choice mode */
1142 if ((cpi->rd_baseline_thresh[best_mode_index] > 0) && (cpi->rd_baseline_thresh[best_mode_index] < (INT_MAX >> 2)))
1144 int best_adjustment = (x->rd_thresh_mult[best_mode_index] >> 3);
1146 x->rd_thresh_mult[best_mode_index] =
1147 (x->rd_thresh_mult[best_mode_index]
1148 >= (MIN_THRESHMULT + best_adjustment)) ?
1149 x->rd_thresh_mult[best_mode_index] - best_adjustment :
1151 x->rd_threshes[best_mode_index] =
1152 (cpi->rd_baseline_thresh[best_mode_index] >> 7) *
1153 x->rd_thresh_mult[best_mode_index];
1158 int this_rdbin = (*returndistortion >> 7);
1160 if (this_rdbin >= 1024)
1165 x->error_bins[this_rdbin] ++;
1168 #if CONFIG_TEMPORAL_DENOISING
1169 if (cpi->oxcf.noise_sensitivity)
1171 if (x->best_sse_inter_mode == DC_PRED)
1173 /* No best MV found. */
1174 x->best_sse_inter_mode = best_mbmode.mode;
1175 x->best_sse_mv = best_mbmode.mv;
1176 x->need_to_clamp_best_mvs = best_mbmode.need_to_clamp_mvs;
1177 x->best_reference_frame = best_mbmode.ref_frame;
1178 best_sse = best_rd_sse;
1180 vp8_denoiser_denoise_mb(&cpi->denoiser, x, best_sse, zero_mv_sse,
1181 recon_yoffset, recon_uvoffset);
1184 /* Reevaluate ZEROMV after denoising. */
1185 if (best_mbmode.ref_frame == INTRA_FRAME &&
1186 x->best_zeromv_reference_frame != INTRA_FRAME)
1189 int this_ref_frame = x->best_zeromv_reference_frame;
1190 rate2 = x->ref_frame_cost[this_ref_frame] +
1191 vp8_cost_mv_ref(ZEROMV, mdcounts);
1194 /* set up the proper prediction buffers for the frame */
1195 x->e_mbd.mode_info_context->mbmi.ref_frame = this_ref_frame;
1196 x->e_mbd.pre.y_buffer = plane[this_ref_frame][0];
1197 x->e_mbd.pre.u_buffer = plane[this_ref_frame][1];
1198 x->e_mbd.pre.v_buffer = plane[this_ref_frame][2];
1200 x->e_mbd.mode_info_context->mbmi.mode = ZEROMV;
1201 x->e_mbd.mode_info_context->mbmi.uv_mode = DC_PRED;
1202 x->e_mbd.mode_info_context->mbmi.mv.as_int = 0;
1203 this_rd = evaluate_inter_mode(&sse, rate2, &distortion2, cpi, x,
1206 if (this_rd < best_rd)
1208 vpx_memcpy(&best_mbmode, &x->e_mbd.mode_info_context->mbmi,
1209 sizeof(MB_MODE_INFO));
1216 if (cpi->is_src_frame_alt_ref &&
1217 (best_mbmode.mode != ZEROMV || best_mbmode.ref_frame != ALTREF_FRAME))
1219 x->e_mbd.mode_info_context->mbmi.mode = ZEROMV;
1220 x->e_mbd.mode_info_context->mbmi.ref_frame = ALTREF_FRAME;
1221 x->e_mbd.mode_info_context->mbmi.mv.as_int = 0;
1222 x->e_mbd.mode_info_context->mbmi.uv_mode = DC_PRED;
1223 x->e_mbd.mode_info_context->mbmi.mb_skip_coeff =
1224 (cpi->common.mb_no_coeff_skip);
1225 x->e_mbd.mode_info_context->mbmi.partitioning = 0;
1230 /* set to the best mb mode, this copy can be skip if x->skip since it
1231 * already has the right content */
1233 vpx_memcpy(&x->e_mbd.mode_info_context->mbmi, &best_mbmode,
1234 sizeof(MB_MODE_INFO));
1236 if (best_mbmode.mode <= B_PRED)
1238 /* set mode_info_context->mbmi.uv_mode */
1239 pick_intra_mbuv_mode(x);
1243 != cpi->common.ref_frame_sign_bias[xd->mode_info_context->mbmi.ref_frame])
1244 best_ref_mv.as_int = best_ref_mv_sb[!sign_bias].as_int;
1246 update_mvcount(x, &best_ref_mv);
1250 void vp8_pick_intra_mode(MACROBLOCK *x, int *rate_)
1252 int error4x4, error16x16 = INT_MAX;
1253 int rate, best_rate = 0, distortion, best_sse;
1254 MB_PREDICTION_MODE mode, best_mode = DC_PRED;
1257 BLOCK *b = &x->block[0];
1258 MACROBLOCKD *xd = &x->e_mbd;
1260 xd->mode_info_context->mbmi.ref_frame = INTRA_FRAME;
1262 pick_intra_mbuv_mode(x);
1264 for (mode = DC_PRED; mode <= TM_PRED; mode ++)
1266 xd->mode_info_context->mbmi.mode = mode;
1267 vp8_build_intra_predictors_mby_s(xd,
1268 xd->dst.y_buffer - xd->dst.y_stride,
1269 xd->dst.y_buffer - 1,
1273 distortion = vp8_variance16x16
1274 (*(b->base_src), b->src_stride, xd->predictor, 16, &sse);
1275 rate = x->mbmode_cost[xd->frame_type][mode];
1276 this_rd = RDCOST(x->rdmult, x->rddiv, rate, distortion);
1278 if (error16x16 > this_rd)
1280 error16x16 = this_rd;
1286 xd->mode_info_context->mbmi.mode = best_mode;
1288 error4x4 = pick_intra4x4mby_modes(x, &rate,
1290 if (error4x4 < error16x16)
1292 xd->mode_info_context->mbmi.mode = B_PRED;