vp9_pickmode.c 15.2 KB
Newer Older
1 2 3 4 5 6 7 8 9 10 11
/*
 *  Copyright (c) 2014 The WebM project authors. All Rights Reserved.
 *
 *  Use of this source code is governed by a BSD-style license
 *  that can be found in the LICENSE file in the root of the source
 *  tree. An additional intellectual property rights grant can be found
 *  in the file PATENTS.  All contributing project authors may
 *  be found in the AUTHORS file in the root of the source tree.
 */

#include <assert.h>
12 13 14
#include <limits.h>
#include <math.h>
#include <stdio.h>
15

16 17 18 19 20 21
#include "./vp9_rtcd.h"

#include "vpx_mem/vpx_mem.h"

#include "vp9/common/vp9_common.h"
#include "vp9/common/vp9_mvref_common.h"
22 23
#include "vp9/common/vp9_reconinter.h"
#include "vp9/common/vp9_reconintra.h"
24

Dmitry Kovalev's avatar
Dmitry Kovalev committed
25
#include "vp9/encoder/vp9_encoder.h"
26
#include "vp9/encoder/vp9_ratectrl.h"
27
#include "vp9/encoder/vp9_rdopt.h"
28

Yunqing Wang's avatar
Yunqing Wang committed
29
static void full_pixel_motion_search(VP9_COMP *cpi, MACROBLOCK *x,
30
                                    BLOCK_SIZE bsize, int mi_row, int mi_col,
31
                                    int_mv *tmp_mv, int *rate_mv) {
32
  MACROBLOCKD *xd = &x->e_mbd;
33
  MB_MODE_INFO *mbmi = &xd->mi[0]->mbmi;
34
  struct buf_2d backup_yv12[MAX_MB_PLANE] = {{0}};
Deb Mukherjee's avatar
Deb Mukherjee committed
35
  int step_param;
36 37 38
  int sadpb = x->sadperbit16;
  MV mvp_full;
  int ref = mbmi->ref_frame[0];
39
  const MV ref_mv = mbmi->ref_mvs[ref][0].as_mv;
40 41 42 43 44 45 46
  int i;

  int tmp_col_min = x->mv_col_min;
  int tmp_col_max = x->mv_col_max;
  int tmp_row_min = x->mv_row_min;
  int tmp_row_max = x->mv_row_max;

47 48
  const YV12_BUFFER_CONFIG *scaled_ref_frame = vp9_get_scaled_ref_frame(cpi,
                                                                        ref);
49 50 51 52 53 54 55 56
  if (scaled_ref_frame) {
    int i;
    // Swap out the reference frame for a version that's been scaled to
    // match the resolution of the current frame, allowing the existing
    // motion search code to be used without additional modifications.
    for (i = 0; i < MAX_MB_PLANE; i++)
      backup_yv12[i] = xd->plane[i].pre[0];

57
    vp9_setup_pre_planes(xd, 0, scaled_ref_frame, mi_row, mi_col, NULL);
58 59
  }

60
  vp9_set_mv_search_range(x, &ref_mv);
61 62 63 64 65

  // TODO(jingning) exploiting adaptive motion search control in non-RD
  // mode decision too.
  step_param = 6;

Jim Bankoski's avatar
Jim Bankoski committed
66
  for (i = LAST_FRAME; i <= LAST_FRAME && cpi->common.show_frame; ++i) {
67 68 69 70 71 72 73 74
    if ((x->pred_mv_sad[ref] >> 3) > x->pred_mv_sad[i]) {
      tmp_mv->as_int = INVALID_MV;

      if (scaled_ref_frame) {
        int i;
        for (i = 0; i < MAX_MB_PLANE; i++)
          xd->plane[i].pre[0] = backup_yv12[i];
      }
Yunqing Wang's avatar
Yunqing Wang committed
75
      return;
76 77
    }
  }
Yaowu Xu's avatar
Yaowu Xu committed
78 79 80 81 82
  assert(x->mv_best_ref_index[ref] <= 2);
  if (x->mv_best_ref_index[ref] < 2)
    mvp_full = mbmi->ref_mvs[ref][x->mv_best_ref_index[ref]].as_mv;
  else
    mvp_full = x->pred_mv[ref].as_mv;
83 84 85 86

  mvp_full.col >>= 3;
  mvp_full.row >>= 3;

87 88 89
  full_pixel_search(cpi, x, bsize, &mvp_full, step_param, sadpb, &ref_mv,
                    &tmp_mv->as_mv, INT_MAX, 0);

90 91 92 93 94 95 96 97 98 99
  x->mv_col_min = tmp_col_min;
  x->mv_col_max = tmp_col_max;
  x->mv_row_min = tmp_row_min;
  x->mv_row_max = tmp_row_max;

  if (scaled_ref_frame) {
    int i;
    for (i = 0; i < MAX_MB_PLANE; i++)
      xd->plane[i].pre[0] = backup_yv12[i];
  }
100 101 102 103 104 105

  // calculate the bit cost on motion vector
  mvp_full.row = tmp_mv->as_mv.row * 8;
  mvp_full.col = tmp_mv->as_mv.col * 8;
  *rate_mv = vp9_mv_bit_cost(&mvp_full, &ref_mv,
                             x->nmvjointcost, x->mvcost, MV_COST_WEIGHT);
106
}
107

108 109
static void sub_pixel_motion_search(VP9_COMP *cpi, MACROBLOCK *x,
                                    BLOCK_SIZE bsize, int mi_row, int mi_col,
110
                                    MV *tmp_mv) {
111
  MACROBLOCKD *xd = &x->e_mbd;
112
  MB_MODE_INFO *mbmi = &xd->mi[0]->mbmi;
113 114
  struct buf_2d backup_yv12[MAX_MB_PLANE] = {{0}};
  int ref = mbmi->ref_frame[0];
Dmitry Kovalev's avatar
Dmitry Kovalev committed
115
  MV ref_mv = mbmi->ref_mvs[ref][0].as_mv;
116
  int dis;
117

118 119 120 121 122 123 124 125 126 127
  const YV12_BUFFER_CONFIG *scaled_ref_frame = vp9_get_scaled_ref_frame(cpi,
                                                                        ref);
  if (scaled_ref_frame) {
    int i;
    // Swap out the reference frame for a version that's been scaled to
    // match the resolution of the current frame, allowing the existing
    // motion search code to be used without additional modifications.
    for (i = 0; i < MAX_MB_PLANE; i++)
      backup_yv12[i] = xd->plane[i].pre[0];

128
    vp9_setup_pre_planes(xd, 0, scaled_ref_frame, mi_row, mi_col, NULL);
129 130
  }

Dmitry Kovalev's avatar
Dmitry Kovalev committed
131
  cpi->find_fractional_mv_step(x, tmp_mv, &ref_mv,
132 133 134 135 136 137 138 139 140 141 142 143 144
                               cpi->common.allow_high_precision_mv,
                               x->errorperbit,
                               &cpi->fn_ptr[bsize],
                               cpi->sf.subpel_force_stop,
                               cpi->sf.subpel_iters_per_step,
                               x->nmvjointcost, x->mvcost,
                               &dis, &x->pred_sse[ref]);

  if (scaled_ref_frame) {
    int i;
    for (i = 0; i < MAX_MB_PLANE; i++)
      xd->plane[i].pre[0] = backup_yv12[i];
  }
145 146

  x->pred_mv[ref].as_mv = *tmp_mv;
147 148
}

149 150 151 152 153 154 155
static void model_rd_for_sb_y(VP9_COMP *cpi, BLOCK_SIZE bsize,
                              MACROBLOCK *x, MACROBLOCKD *xd,
                              int *out_rate_sum, int64_t *out_dist_sum) {
  // Note our transform coeffs are 8 times an orthogonal transform.
  // Hence quantizer step is also 8 times. To get effective quantizer
  // we need to divide by 8 before sending to modeling function.
  unsigned int sse;
156 157
  int rate;
  int64_t dist;
158 159 160 161

  struct macroblock_plane *const p = &x->plane[0];
  struct macroblockd_plane *const pd = &xd->plane[0];

162 163
  unsigned int var = cpi->fn_ptr[bsize].vf(p->src.buf, p->src.stride,
                                           pd->dst.buf, pd->dst.stride, &sse);
164

165 166 167 168 169 170 171
  // TODO(jingning) This is a temporary solution to account for frames with
  // light changes. Need to customize the rate-distortion modeling for non-RD
  // mode decision.
  if ((sse >> 3) > var)
    sse = var;

  vp9_model_rd_from_var_lapndz(var + sse, 1 << num_pels_log2_lookup[bsize],
172 173
                               pd->dequant[1] >> 3, &rate, &dist);
  *out_rate_sum = rate;
174
  *out_dist_sum = dist << 3;
175 176
}

177 178 179 180 181 182 183
// TODO(jingning) placeholder for inter-frame non-RD mode decision.
// this needs various further optimizations. to be continued..
int64_t vp9_pick_inter_mode(VP9_COMP *cpi, MACROBLOCK *x,
                            const TileInfo *const tile,
                            int mi_row, int mi_col,
                            int *returnrate,
                            int64_t *returndistortion,
Jim Bankoski's avatar
Jim Bankoski committed
184
                            BLOCK_SIZE bsize) {
185
  MACROBLOCKD *xd = &x->e_mbd;
186
  MB_MODE_INFO *mbmi = &xd->mi[0]->mbmi;
187 188
  struct macroblock_plane *const p = &x->plane[0];
  struct macroblockd_plane *const pd = &xd->plane[0];
189
  PREDICTION_MODE this_mode, best_mode = ZEROMV;
190
  MV_REFERENCE_FRAME ref_frame, best_ref_frame = LAST_FRAME;
191
  INTERP_FILTER best_pred_filter = EIGHTTAP;
192 193 194 195 196
  int_mv frame_mv[MB_MODE_COUNT][MAX_REF_FRAMES];
  struct buf_2d yv12_mb[4][MAX_MB_PLANE];
  static const int flag_list[4] = { 0, VP9_LAST_FLAG, VP9_GOLD_FLAG,
                                    VP9_ALT_FLAG };
  int64_t best_rd = INT64_MAX;
197
  int64_t this_rd = INT64_MAX;
198

199 200 201
  int rate = INT_MAX;
  int64_t dist = INT64_MAX;

202 203 204 205 206 207 208
  VP9_COMMON *cm = &cpi->common;
  int intra_cost_penalty = 20 * vp9_dc_quant(cm->base_qindex, cm->y_dc_delta_q);

  const int64_t inter_mode_thresh = RDCOST(x->rdmult, x->rddiv,
                                           intra_cost_penalty, 0);
  const int64_t intra_mode_cost = 50;

209
  unsigned char segment_id = mbmi->segment_id;
Dmitry Kovalev's avatar
Dmitry Kovalev committed
210 211
  const int *const rd_threshes = cpi->rd.threshes[segment_id][bsize];
  const int *const rd_thresh_freq_fact = cpi->rd.thresh_freq_fact[bsize];
212
  // Mode index conversion form THR_MODES to PREDICTION_MODE for a ref frame.
213
  int mode_idx[MB_MODE_COUNT] = {0};
214
  INTERP_FILTER filter_ref = SWITCHABLE;
215 216 217
  int bsl = mi_width_log2_lookup[bsize];
  int pred_filter_search = (((mi_row + mi_col) >> bsl) +
                            cpi->sf.chessboard_index) & 0x01;
218

219 220 221
  x->skip_encode = cpi->sf.skip_encode_frame && x->q_index < QIDX_SKIP_THRESH;

  x->skip = 0;
222
  if (!x->in_active_map)
223 224 225
    x->skip = 1;
  // initialize mode decisions
  *returnrate = INT_MAX;
226
  *returndistortion = INT64_MAX;
227 228 229 230 231 232
  vpx_memset(mbmi, 0, sizeof(MB_MODE_INFO));
  mbmi->sb_type = bsize;
  mbmi->ref_frame[0] = NONE;
  mbmi->ref_frame[1] = NONE;
  mbmi->tx_size = MIN(max_txsize_lookup[bsize],
                      tx_mode_to_biggest_tx_size[cpi->common.tx_mode]);
233 234
  mbmi->interp_filter = cpi->common.interp_filter == SWITCHABLE ?
                        EIGHTTAP : cpi->common.interp_filter;
235
  mbmi->skip = 0;
236
  mbmi->segment_id = segment_id;
237

Jim Bankoski's avatar
Jim Bankoski committed
238
  for (ref_frame = LAST_FRAME; ref_frame <= LAST_FRAME ; ++ref_frame) {
239 240
    x->pred_mv_sad[ref_frame] = INT_MAX;
    if (cpi->ref_frame_flags & flag_list[ref_frame]) {
241
      vp9_setup_buffer_inter(cpi, x, tile,
242
                             ref_frame, bsize, mi_row, mi_col,
243 244 245 246 247 248
                             frame_mv[NEARESTMV], frame_mv[NEARMV], yv12_mb);
    }
    frame_mv[NEWMV][ref_frame].as_int = INVALID_MV;
    frame_mv[ZEROMV][ref_frame].as_int = 0;
  }

249 250 251 252 253
  if (xd->up_available)
    filter_ref = xd->mi[-xd->mi_stride]->mbmi.interp_filter;
  else if (xd->left_available)
    filter_ref = xd->mi[-1]->mbmi.interp_filter;

Jim Bankoski's avatar
Jim Bankoski committed
254
  for (ref_frame = LAST_FRAME; ref_frame <= LAST_FRAME ; ++ref_frame) {
255 256 257 258 259 260 261 262 263
    if (!(cpi->ref_frame_flags & flag_list[ref_frame]))
      continue;

    // Select prediction reference frames.
    xd->plane[0].pre[0] = yv12_mb[ref_frame][0];

    clamp_mv2(&frame_mv[NEARESTMV][ref_frame].as_mv, xd);
    clamp_mv2(&frame_mv[NEARMV][ref_frame].as_mv, xd);

264 265
    mbmi->ref_frame[0] = ref_frame;

266 267 268 269 270 271 272 273
    // Set conversion index for LAST_FRAME.
    if (ref_frame == LAST_FRAME) {
      mode_idx[NEARESTMV] = THR_NEARESTMV;   // LAST_FRAME, NEARESTMV
      mode_idx[NEARMV] = THR_NEARMV;         // LAST_FRAME, NEARMV
      mode_idx[ZEROMV] = THR_ZEROMV;         // LAST_FRAME, ZEROMV
      mode_idx[NEWMV] = THR_NEWMV;           // LAST_FRAME, NEWMV
    }

274
    for (this_mode = NEARESTMV; this_mode <= NEWMV; ++this_mode) {
275 276
      int rate_mv = 0;

277 278 279
      if (cpi->sf.disable_inter_mode_mask[bsize] &
          (1 << INTER_OFFSET(this_mode)))
        continue;
Jim Bankoski's avatar
Jim Bankoski committed
280

281 282 283 284 285
      if (best_rd < ((int64_t)rd_threshes[mode_idx[this_mode]] *
          rd_thresh_freq_fact[this_mode] >> 5) ||
          rd_threshes[mode_idx[this_mode]] == INT_MAX)
        continue;

Jim Bankoski's avatar
Jim Bankoski committed
286
      if (this_mode == NEWMV) {
287
        int rate_mode = 0;
288
        if (this_rd < (int64_t)(1 << num_pels_log2_lookup[bsize]))
289 290
          continue;

291
        full_pixel_motion_search(cpi, x, bsize, mi_row, mi_col,
292
                                 &frame_mv[NEWMV][ref_frame], &rate_mv);
Jim Bankoski's avatar
Jim Bankoski committed
293 294 295

        if (frame_mv[NEWMV][ref_frame].as_int == INVALID_MV)
          continue;
296

297 298
        rate_mode = cpi->inter_mode_cost[mbmi->mode_context[ref_frame]]
                                        [INTER_OFFSET(this_mode)];
299 300 301
        if (RDCOST(x->rdmult, x->rddiv, rate_mv + rate_mode, 0) > best_rd)
          continue;

302
        sub_pixel_motion_search(cpi, x, bsize, mi_row, mi_col,
303
                                &frame_mv[NEWMV][ref_frame].as_mv);
Jim Bankoski's avatar
Jim Bankoski committed
304 305
      }

306 307 308 309 310
      if (this_mode != NEARESTMV)
        if (frame_mv[this_mode][ref_frame].as_int ==
            frame_mv[NEARESTMV][ref_frame].as_int)
          continue;

311 312
      mbmi->mode = this_mode;
      mbmi->mv[0].as_int = frame_mv[this_mode][ref_frame].as_int;
313

314 315 316 317
      // Search for the best prediction filter type, when the resulting
      // motion vector is at sub-pixel accuracy level for luma component, i.e.,
      // the last three bits are all zeros.
      if ((this_mode == NEWMV || filter_ref == SWITCHABLE) &&
318
          pred_filter_search &&
319 320 321 322 323 324 325 326 327 328 329 330 331
          ((mbmi->mv[0].as_mv.row & 0x07) != 0 ||
           (mbmi->mv[0].as_mv.col & 0x07) != 0)) {
        int64_t tmp_rdcost1 = INT64_MAX;
        int64_t tmp_rdcost2 = INT64_MAX;
        int64_t tmp_rdcost3 = INT64_MAX;
        int pf_rate[3];
        int64_t pf_dist[3];

        mbmi->interp_filter = EIGHTTAP;
        vp9_build_inter_predictors_sby(xd, mi_row, mi_col, bsize);
        model_rd_for_sb_y(cpi, bsize, x, xd, &pf_rate[EIGHTTAP],
                          &pf_dist[EIGHTTAP]);
        tmp_rdcost1 = RDCOST(x->rdmult, x->rddiv,
332
                             vp9_get_switchable_rate(cpi) + pf_rate[EIGHTTAP],
333 334 335 336 337 338
                             pf_dist[EIGHTTAP]);

        mbmi->interp_filter = EIGHTTAP_SHARP;
        vp9_build_inter_predictors_sby(xd, mi_row, mi_col, bsize);
        model_rd_for_sb_y(cpi, bsize, x, xd, &pf_rate[EIGHTTAP_SHARP],
                          &pf_dist[EIGHTTAP_SHARP]);
339 340 341
        tmp_rdcost2 = RDCOST(x->rdmult, x->rddiv, vp9_get_switchable_rate(cpi) +
                                 pf_rate[EIGHTTAP_SHARP],
                             pf_dist[EIGHTTAP_SHARP]);
342 343 344 345 346

        mbmi->interp_filter = EIGHTTAP_SMOOTH;
        vp9_build_inter_predictors_sby(xd, mi_row, mi_col, bsize);
        model_rd_for_sb_y(cpi, bsize, x, xd, &pf_rate[EIGHTTAP_SMOOTH],
                          &pf_dist[EIGHTTAP_SMOOTH]);
347 348 349
        tmp_rdcost3 = RDCOST(x->rdmult, x->rddiv, vp9_get_switchable_rate(cpi) +
                                 pf_rate[EIGHTTAP_SMOOTH],
                             pf_dist[EIGHTTAP_SMOOTH]);
350 351 352 353 354 355 356 357 358 359 360 361 362 363 364 365 366 367 368 369 370

        if (tmp_rdcost2 < tmp_rdcost1) {
          if (tmp_rdcost2 < tmp_rdcost3)
            mbmi->interp_filter = EIGHTTAP_SHARP;
          else
            mbmi->interp_filter = EIGHTTAP_SMOOTH;
        } else {
          if (tmp_rdcost1 < tmp_rdcost3)
            mbmi->interp_filter = EIGHTTAP;
          else
            mbmi->interp_filter = EIGHTTAP_SMOOTH;
        }

        rate = pf_rate[mbmi->interp_filter];
        dist = pf_dist[mbmi->interp_filter];
      } else {
        mbmi->interp_filter = (filter_ref == SWITCHABLE) ? EIGHTTAP: filter_ref;
        vp9_build_inter_predictors_sby(xd, mi_row, mi_col, bsize);
        model_rd_for_sb_y(cpi, bsize, x, xd, &rate, &dist);
      }

371
      rate += rate_mv;
372
      rate += cpi->inter_mode_cost[mbmi->mode_context[ref_frame]]
373 374
                                [INTER_OFFSET(this_mode)];
      this_rd = RDCOST(x->rdmult, x->rddiv, rate, dist);
375 376 377

      if (this_rd < best_rd) {
        best_rd = this_rd;
378 379
        *returnrate = rate;
        *returndistortion = dist;
380
        best_mode = this_mode;
381
        best_pred_filter = mbmi->interp_filter;
382
        best_ref_frame = ref_frame;
383 384 385 386
      }
    }
  }

387
  mbmi->mode = best_mode;
388
  mbmi->interp_filter = best_pred_filter;
389 390
  mbmi->ref_frame[0] = best_ref_frame;
  mbmi->mv[0].as_int = frame_mv[best_mode][best_ref_frame].as_int;
391
  xd->mi[0]->bmi[0].as_mv[0].as_int = mbmi->mv[0].as_int;
392

393 394
  // Perform intra prediction search, if the best SAD is above a certain
  // threshold.
395
  if (best_rd > inter_mode_thresh && bsize < cpi->sf.max_intra_bsize) {
396
    for (this_mode = DC_PRED; this_mode <= DC_PRED; ++this_mode) {
397 398 399 400 401
      vp9_predict_intra_block(xd, 0, b_width_log2(bsize),
                              mbmi->tx_size, this_mode,
                              &p->src.buf[0], p->src.stride,
                              &pd->dst.buf[0], pd->dst.stride, 0, 0, 0);

402
      model_rd_for_sb_y(cpi, bsize, x, xd, &rate, &dist);
403
      rate += cpi->mbmode_cost[this_mode];
404
      rate += intra_cost_penalty;
405
      this_rd = RDCOST(x->rdmult, x->rddiv, rate, dist);
406 407 408

      if (this_rd + intra_mode_cost < best_rd) {
        best_rd = this_rd;
409 410
        *returnrate = rate;
        *returndistortion = dist;
411 412 413
        mbmi->mode = this_mode;
        mbmi->ref_frame[0] = INTRA_FRAME;
        mbmi->uv_mode = this_mode;
414
        mbmi->mv[0].as_int = INVALID_MV;
415 416 417
      }
    }
  }
418

419 420
  return INT64_MAX;
}