vp9_encodeframe.c 74.9 KB
Newer Older
John Koleszar's avatar
John Koleszar committed
1
/*
2
 *  Copyright (c) 2010 The WebM project authors. All Rights Reserved.
John Koleszar's avatar
John Koleszar committed
3
 *
4
 *  Use of this source code is governed by a BSD-style license
5 6
 *  that can be found in the LICENSE file in the root of the source
 *  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.
John Koleszar's avatar
John Koleszar committed
9 10 11
 */


12
#include "./vpx_config.h"
13
#include "vp9/encoder/vp9_encodeframe.h"
14 15
#include "vp9/encoder/vp9_encodemb.h"
#include "vp9/encoder/vp9_encodemv.h"
16
#include "vp9/common/vp9_common.h"
17
#include "vp9/encoder/vp9_onyx_int.h"
18
#include "vp9/common/vp9_extend.h"
Yaowu Xu's avatar
Yaowu Xu committed
19
#include "vp9/common/vp9_entropy.h"
20 21
#include "vp9/common/vp9_entropymode.h"
#include "vp9/common/vp9_quant_common.h"
22
#include "vp9/encoder/vp9_segmentation.h"
23
#include "vp9/common/vp9_setupintrarecon.h"
24
#include "vp9/encoder/vp9_encodeintra.h"
25 26
#include "vp9/common/vp9_reconinter.h"
#include "vp9/common/vp9_invtrans.h"
27
#include "vp9/encoder/vp9_rdopt.h"
28 29 30
#include "vp9/common/vp9_findnearmv.h"
#include "vp9/common/vp9_reconintra.h"
#include "vp9/common/vp9_seg_common.h"
31
#include "vp9/common/vp9_tile_common.h"
32
#include "vp9/encoder/vp9_tokenize.h"
Yaowu Xu's avatar
Yaowu Xu committed
33
#include "./vp9_rtcd.h"
John Koleszar's avatar
John Koleszar committed
34
#include <stdio.h>
35
#include <math.h>
John Koleszar's avatar
John Koleszar committed
36 37
#include <limits.h>
#include "vpx_ports/vpx_timer.h"
38 39
#include "vp9/common/vp9_pred_common.h"
#include "vp9/common/vp9_mvref_common.h"
Paul Wilkins's avatar
Paul Wilkins committed
40

Yaowu Xu's avatar
Yaowu Xu committed
41
#define DBG_PRNT_SEGMAP 0
42

43
// #define ENC_DEBUG
44
#ifdef ENC_DEBUG
John Koleszar's avatar
John Koleszar committed
45
int enc_debug = 0;
46 47
#endif

Dmitry Kovalev's avatar
Dmitry Kovalev committed
48
void vp9_select_interp_filter_type(VP9_COMP *cpi);
49

Ronald S. Bultje's avatar
Ronald S. Bultje committed
50 51
static void encode_macroblock(VP9_COMP *cpi, TOKENEXTRA **t,
                              int output_enabled, int mb_row, int mb_col);
52

53 54 55
static void encode_superblock(VP9_COMP *cpi, TOKENEXTRA **t,
                              int output_enabled, int mb_row, int mb_col,
                              BLOCK_SIZE_TYPE bsize);
56

57
static void adjust_act_zbin(VP9_COMP *cpi, MACROBLOCK *x);
58

John Koleszar's avatar
John Koleszar committed
59
#ifdef MODE_STATS
60
unsigned int inter_y_modes[MB_MODE_COUNT];
61
unsigned int inter_uv_modes[VP9_UV_MODES];
62
unsigned int inter_b_modes[B_MODE_COUNT];
63 64 65 66
unsigned int y_modes[VP9_YMODES];
unsigned int i8x8_modes[VP9_I8X8_MODES];
unsigned int uv_modes[VP9_UV_MODES];
unsigned int uv_modes_y[VP9_YMODES][VP9_UV_MODES];
67
unsigned int b_modes[B_MODE_COUNT];
John Koleszar's avatar
John Koleszar committed
68 69 70
#endif


71 72 73
/* activity_avg must be positive, or flat regions could get a zero weight
 *  (infinite lambda), which confounds analysis.
 * This also avoids the need for divide by zero checks in
74
 *  vp9_activity_masking().
75
 */
76
#define VP9_ACTIVITY_AVG_MIN (64)
77 78 79 80 81 82

/* This is used as a reference when computing the source variance for the
 *  purposes of activity masking.
 * Eventually this should be replaced by custom no-reference routines,
 *  which will be faster.
 */
83
static const uint8_t VP9_VAR_OFFS[16] = {
John Koleszar's avatar
John Koleszar committed
84
  128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128
85 86
};

87 88

// Original activity measure from Tim T's code.
89
static unsigned int tt_activity_measure(VP9_COMP *cpi, MACROBLOCK *x) {
John Koleszar's avatar
John Koleszar committed
90 91 92 93 94 95 96 97 98
  unsigned int act;
  unsigned int sse;
  /* TODO: This could also be done over smaller areas (8x8), but that would
   *  require extensive changes elsewhere, as lambda is assumed to be fixed
   *  over an entire MB in most of the code.
   * Another option is to compute four 8x8 variances, and pick a single
   *  lambda using a non-linear combination (e.g., the smallest, or second
   *  smallest, etc.).
   */
99
  act = vp9_variance16x16(x->src.y_buffer, x->src.y_stride, VP9_VAR_OFFS, 0,
100
                          &sse);
Dmitry Kovalev's avatar
Dmitry Kovalev committed
101
  act <<= 4;
John Koleszar's avatar
John Koleszar committed
102 103 104 105 106 107

  /* If the region is flat, lower the activity some more. */
  if (act < 8 << 12)
    act = act < 5 << 12 ? act : 5 << 12;

  return act;
108 109
}

110
// Stub for alternative experimental activity measures.
111
static unsigned int alt_activity_measure(VP9_COMP *cpi,
John Koleszar's avatar
John Koleszar committed
112
                                         MACROBLOCK *x, int use_dc_pred) {
113
  return vp9_encode_intra(cpi, x, use_dc_pred);
114 115 116 117 118
}


// Measure the activity of the current macroblock
// What we measure here is TBD so abstracted to this function
119
#define ALT_ACT_MEASURE 1
120
static unsigned int mb_activity_measure(VP9_COMP *cpi, MACROBLOCK *x,
John Koleszar's avatar
John Koleszar committed
121 122
                                        int mb_row, int mb_col) {
  unsigned int mb_activity;
123

John Koleszar's avatar
John Koleszar committed
124 125
  if (ALT_ACT_MEASURE) {
    int use_dc_pred = (mb_col || mb_row) && (!mb_col || !mb_row);
126

John Koleszar's avatar
John Koleszar committed
127 128 129 130 131 132
    // Or use and alternative.
    mb_activity = alt_activity_measure(cpi, x, use_dc_pred);
  } else {
    // Original activity measure from Tim T's code.
    mb_activity = tt_activity_measure(cpi, x);
  }
133

134 135
  if (mb_activity < VP9_ACTIVITY_AVG_MIN)
    mb_activity = VP9_ACTIVITY_AVG_MIN;
136

John Koleszar's avatar
John Koleszar committed
137
  return mb_activity;
138 139 140
}

// Calculate an "average" mb activity value for the frame
141
#define ACT_MEDIAN 0
142
static void calc_av_activity(VP9_COMP *cpi, int64_t activity_sum) {
143
#if ACT_MEDIAN
John Koleszar's avatar
John Koleszar committed
144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172
  // Find median: Simple n^2 algorithm for experimentation
  {
    unsigned int median;
    unsigned int i, j;
    unsigned int *sortlist;
    unsigned int tmp;

    // Create a list to sort to
    CHECK_MEM_ERROR(sortlist,
    vpx_calloc(sizeof(unsigned int),
    cpi->common.MBs));

    // Copy map to sort list
    vpx_memcpy(sortlist, cpi->mb_activity_map,
    sizeof(unsigned int) * cpi->common.MBs);


    // Ripple each value down to its correct position
    for (i = 1; i < cpi->common.MBs; i ++) {
      for (j = i; j > 0; j --) {
        if (sortlist[j] < sortlist[j - 1]) {
          // Swap values
          tmp = sortlist[j - 1];
          sortlist[j - 1] = sortlist[j];
          sortlist[j] = tmp;
        } else
          break;
      }
    }
173

John Koleszar's avatar
John Koleszar committed
174 175 176
    // Even number MBs so estimate median as mean of two either side.
    median = (1 + sortlist[cpi->common.MBs >> 1] +
              sortlist[(cpi->common.MBs >> 1) + 1]) >> 1;
177

John Koleszar's avatar
John Koleszar committed
178
    cpi->activity_avg = median;
179

John Koleszar's avatar
John Koleszar committed
180 181
    vpx_free(sortlist);
  }
182
#else
John Koleszar's avatar
John Koleszar committed
183 184
  // Simple mean for now
  cpi->activity_avg = (unsigned int)(activity_sum / cpi->common.MBs);
185 186
#endif

187 188
  if (cpi->activity_avg < VP9_ACTIVITY_AVG_MIN)
    cpi->activity_avg = VP9_ACTIVITY_AVG_MIN;
189

John Koleszar's avatar
John Koleszar committed
190 191 192
  // Experimental code: return fixed value normalized for several clips
  if (ALT_ACT_MEASURE)
    cpi->activity_avg = 100000;
193 194
}

195
#define USE_ACT_INDEX   0
196
#define OUTPUT_NORM_ACT_STATS   0
197 198

#if USE_ACT_INDEX
199
// Calculate an activity index for each mb
200 201
static void calc_activity_index(VP9_COMP *cpi, MACROBLOCK *x) {
  VP9_COMMON *const cm = &cpi->common;
John Koleszar's avatar
John Koleszar committed
202
  int mb_row, mb_col;
203

John Koleszar's avatar
John Koleszar committed
204 205 206
  int64_t act;
  int64_t a;
  int64_t b;
207 208

#if OUTPUT_NORM_ACT_STATS
John Koleszar's avatar
John Koleszar committed
209 210
  FILE *f = fopen("norm_act.stt", "a");
  fprintf(f, "\n%12d\n", cpi->activity_avg);
211 212
#endif

John Koleszar's avatar
John Koleszar committed
213 214
  // Reset pointers to start of activity map
  x->mb_activity_ptr = cpi->mb_activity_map;
215

John Koleszar's avatar
John Koleszar committed
216 217 218 219 220 221
  // Calculate normalized mb activity number.
  for (mb_row = 0; mb_row < cm->mb_rows; mb_row++) {
    // for each macroblock col in image
    for (mb_col = 0; mb_col < cm->mb_cols; mb_col++) {
      // Read activity from the map
      act = *(x->mb_activity_ptr);
222

John Koleszar's avatar
John Koleszar committed
223 224 225
      // Calculate a normalized activity number
      a = act + 4 * cpi->activity_avg;
      b = 4 * act + cpi->activity_avg;
226

John Koleszar's avatar
John Koleszar committed
227 228 229 230
      if (b >= a)
        *(x->activity_ptr) = (int)((b + (a >> 1)) / a) - 1;
      else
        *(x->activity_ptr) = 1 - (int)((a + (b >> 1)) / b);
231 232

#if OUTPUT_NORM_ACT_STATS
John Koleszar's avatar
John Koleszar committed
233
      fprintf(f, " %6d", *(x->mb_activity_ptr));
234
#endif
John Koleszar's avatar
John Koleszar committed
235 236 237
      // Increment activity map pointers
      x->mb_activity_ptr++;
    }
238 239

#if OUTPUT_NORM_ACT_STATS
John Koleszar's avatar
John Koleszar committed
240
    fprintf(f, "\n");
241 242
#endif

John Koleszar's avatar
John Koleszar committed
243
  }
244 245

#if OUTPUT_NORM_ACT_STATS
John Koleszar's avatar
John Koleszar committed
246
  fclose(f);
247 248 249
#endif

}
250
#endif
251 252 253

// Loop through all MBs. Note activity of each, average activity and
// calculate a normalized activity for each
254
static void build_activity_map(VP9_COMP *cpi) {
John Koleszar's avatar
John Koleszar committed
255
  MACROBLOCK *const x = &cpi->mb;
John Koleszar's avatar
John Koleszar committed
256
  MACROBLOCKD *xd = &x->e_mbd;
257
  VP9_COMMON *const cm = &cpi->common;
258

259
#if ALT_ACT_MEASURE
John Koleszar's avatar
John Koleszar committed
260 261 262
  YV12_BUFFER_CONFIG *new_yv12 = &cm->yv12_fb[cm->new_fb_idx];
  int recon_yoffset;
  int recon_y_stride = new_yv12->y_stride;
263 264
#endif

John Koleszar's avatar
John Koleszar committed
265 266 267
  int mb_row, mb_col;
  unsigned int mb_activity;
  int64_t activity_sum = 0;
268

269 270
  x->mb_activity_ptr = cpi->mb_activity_map;

John Koleszar's avatar
John Koleszar committed
271 272
  // for each macroblock row in image
  for (mb_row = 0; mb_row < cm->mb_rows; mb_row++) {
273
#if ALT_ACT_MEASURE
John Koleszar's avatar
John Koleszar committed
274 275 276
    // reset above block coeffs
    xd->up_available = (mb_row != 0);
    recon_yoffset = (mb_row * recon_y_stride * 16);
277
#endif
John Koleszar's avatar
John Koleszar committed
278 279
    // for each macroblock col in image
    for (mb_col = 0; mb_col < cm->mb_cols; mb_col++) {
280
#if ALT_ACT_MEASURE
John Koleszar's avatar
John Koleszar committed
281 282 283
      xd->dst.y_buffer = new_yv12->y_buffer + recon_yoffset;
      xd->left_available = (mb_col != 0);
      recon_yoffset += 16;
284
#endif
285

John Koleszar's avatar
John Koleszar committed
286 287
      // measure activity
      mb_activity = mb_activity_measure(cpi, x, mb_row, mb_col);
288

John Koleszar's avatar
John Koleszar committed
289 290
      // Keep frame sum
      activity_sum += mb_activity;
291

John Koleszar's avatar
John Koleszar committed
292 293
      // Store MB level activity details.
      *x->mb_activity_ptr = mb_activity;
294

John Koleszar's avatar
John Koleszar committed
295 296
      // Increment activity map pointer
      x->mb_activity_ptr++;
297

John Koleszar's avatar
John Koleszar committed
298 299 300
      // adjust to the next column of source macroblocks
      x->src.y_buffer += 16;
    }
301

302

John Koleszar's avatar
John Koleszar committed
303 304
    // adjust to the next row of mbs
    x->src.y_buffer += 16 * x->src.y_stride - 16 * cm->mb_cols;
305 306

#if ALT_ACT_MEASURE
John Koleszar's avatar
John Koleszar committed
307
    // extend the recon for intra prediction
308
    vp9_extend_mb_row(new_yv12, xd->dst.y_buffer + 16,
John Koleszar's avatar
John Koleszar committed
309
                      xd->dst.u_buffer + 8, xd->dst.v_buffer + 8);
310 311
#endif

John Koleszar's avatar
John Koleszar committed
312
  }
313

John Koleszar's avatar
John Koleszar committed
314 315
  // Calculate an "average" MB activity
  calc_av_activity(cpi, activity_sum);
316

317
#if USE_ACT_INDEX
John Koleszar's avatar
John Koleszar committed
318 319
  // Calculate an activity index number of each mb
  calc_activity_index(cpi, x);
320 321
#endif

322 323
}

324
// Macroblock activity masking
325
void vp9_activity_masking(VP9_COMP *cpi, MACROBLOCK *x) {
326
#if USE_ACT_INDEX
John Koleszar's avatar
John Koleszar committed
327 328 329
  x->rdmult += *(x->mb_activity_ptr) * (x->rdmult >> 2);
  x->errorperbit = x->rdmult * 100 / (110 * x->rddiv);
  x->errorperbit += (x->errorperbit == 0);
330
#else
John Koleszar's avatar
John Koleszar committed
331 332 333
  int64_t a;
  int64_t b;
  int64_t act = *(x->mb_activity_ptr);
334

John Koleszar's avatar
John Koleszar committed
335 336 337
  // Apply the masking to the RD multiplier.
  a = act + (2 * cpi->activity_avg);
  b = (2 * act) + cpi->activity_avg;
338

John Koleszar's avatar
John Koleszar committed
339 340 341
  x->rdmult = (unsigned int)(((int64_t)x->rdmult * b + (a >> 1)) / a);
  x->errorperbit = x->rdmult * 100 / (110 * x->rddiv);
  x->errorperbit += (x->errorperbit == 0);
342
#endif
343

John Koleszar's avatar
John Koleszar committed
344 345
  // Activity based Zbin adjustment
  adjust_act_zbin(cpi, x);
346
}
John Koleszar's avatar
John Koleszar committed
347

Deb Mukherjee's avatar
Deb Mukherjee committed
348 349 350 351 352 353 354 355 356 357 358 359 360 361 362 363 364 365 366 367 368 369 370 371 372 373 374 375 376 377 378 379 380 381 382 383 384 385 386 387 388
#if CONFIG_NEW_MVREF
static int vp9_cost_mv_ref_id(vp9_prob * ref_id_probs, int mv_ref_id) {
  int cost;

  // Encode the index for the MV reference.
  switch (mv_ref_id) {
    case 0:
      cost = vp9_cost_zero(ref_id_probs[0]);
      break;
    case 1:
      cost = vp9_cost_one(ref_id_probs[0]);
      cost += vp9_cost_zero(ref_id_probs[1]);
      break;
    case 2:
      cost = vp9_cost_one(ref_id_probs[0]);
      cost += vp9_cost_one(ref_id_probs[1]);
      cost += vp9_cost_zero(ref_id_probs[2]);
      break;
    case 3:
      cost = vp9_cost_one(ref_id_probs[0]);
      cost += vp9_cost_one(ref_id_probs[1]);
      cost += vp9_cost_one(ref_id_probs[2]);
      break;

      // TRAP.. This should not happen
    default:
      assert(0);
      break;
  }
  return cost;
}

// Estimate the cost of each coding the vector using each reference candidate
static unsigned int pick_best_mv_ref(MACROBLOCK *x,
                                     MV_REFERENCE_FRAME ref_frame,
                                     int_mv target_mv,
                                     int_mv * mv_ref_list,
                                     int_mv * best_ref) {
  int i;
  int best_index = 0;
  int cost, cost2;
389
  int zero_seen = (mv_ref_list[0].as_int) ? 0 : 1;
Deb Mukherjee's avatar
Deb Mukherjee committed
390 391 392
  MACROBLOCKD *xd = &x->e_mbd;
  int max_mv = MV_MAX;

393
  cost = vp9_cost_mv_ref_id(xd->mb_mv_ref_probs[ref_frame], 0) +
Deb Mukherjee's avatar
Deb Mukherjee committed
394 395 396
         vp9_mv_bit_cost(&target_mv, &mv_ref_list[0], x->nmvjointcost,
                         x->mvcost, 96, xd->allow_high_precision_mv);

Paul Wilkins's avatar
Paul Wilkins committed
397
  for (i = 1; i < MAX_MV_REF_CANDIDATES; ++i) {
Deb Mukherjee's avatar
Deb Mukherjee committed
398 399
    // If we see a 0,0 reference vector for a second time we have reached
    // the end of the list of valid candidate vectors.
Jim Bankoski's avatar
Jim Bankoski committed
400
    if (!mv_ref_list[i].as_int) {
Deb Mukherjee's avatar
Deb Mukherjee committed
401 402 403
      if (zero_seen)
        break;
      else
404
        zero_seen = 1;
Jim Bankoski's avatar
Jim Bankoski committed
405
    }
Deb Mukherjee's avatar
Deb Mukherjee committed
406 407 408 409 410 411 412 413

    // Check for cases where the reference choice would give rise to an
    // uncodable/out of range residual for row or col.
    if ((abs(target_mv.as_mv.row - mv_ref_list[i].as_mv.row) > max_mv) ||
        (abs(target_mv.as_mv.col - mv_ref_list[i].as_mv.col) > max_mv)) {
      continue;
    }

414
    cost2 = vp9_cost_mv_ref_id(xd->mb_mv_ref_probs[ref_frame], i) +
Deb Mukherjee's avatar
Deb Mukherjee committed
415 416 417 418 419 420 421 422 423 424 425 426 427 428
            vp9_mv_bit_cost(&target_mv, &mv_ref_list[i], x->nmvjointcost,
                            x->mvcost, 96, xd->allow_high_precision_mv);

    if (cost2 < cost) {
      cost = cost2;
      best_index = i;
    }
  }
  best_ref->as_int = mv_ref_list[best_index].as_int;

  return best_index;
}
#endif

429
static void update_state(VP9_COMP *cpi,
430 431
                         PICK_MODE_CONTEXT *ctx,
                         BLOCK_SIZE_TYPE bsize,
Ronald S. Bultje's avatar
Ronald S. Bultje committed
432 433
                         int output_enabled) {
  int i, x_idx, y;
434
  VP9_COMMON *const cm = &cpi->common;
435 436
  MACROBLOCK *const x = &cpi->mb;
  MACROBLOCKD *const xd = &x->e_mbd;
John Koleszar's avatar
John Koleszar committed
437
  MODE_INFO *mi = &ctx->mic;
Ronald S. Bultje's avatar
Ronald S. Bultje committed
438
  MB_MODE_INFO *const mbmi = &xd->mode_info_context->mbmi;
John Koleszar's avatar
John Koleszar committed
439 440
  int mb_mode = mi->mbmi.mode;
  int mb_mode_index = ctx->best_mode_index;
Ronald S. Bultje's avatar
Ronald S. Bultje committed
441
  const int mis = cpi->common.mode_info_stride;
442
  const int bh = 1 << mb_height_log2(bsize), bw = 1 << mb_width_log2(bsize);
Adrian Grange's avatar
Adrian Grange committed
443 444

#if CONFIG_DEBUG
John Koleszar's avatar
John Koleszar committed
445 446 447
  assert(mb_mode < MB_MODE_COUNT);
  assert(mb_mode_index < MAX_MODES);
  assert(mi->mbmi.ref_frame < MAX_REF_FRAMES);
Adrian Grange's avatar
Adrian Grange committed
448
#endif
449
  assert(mi->mbmi.sb_type == bsize);
Adrian Grange's avatar
Adrian Grange committed
450

John Koleszar's avatar
John Koleszar committed
451 452
  // Restore the coding context of the MB to that that was in place
  // when the mode was picked for it
453 454 455 456
  for (y = 0; y < bh; y++) {
    for (x_idx = 0; x_idx < bw; x_idx++) {
      if ((xd->mb_to_right_edge >> 7) + bw > x_idx &&
          (xd->mb_to_bottom_edge >> 7) + bh > y) {
Ronald S. Bultje's avatar
Ronald S. Bultje committed
457 458 459 460
        MODE_INFO *mi_addr = xd->mode_info_context + x_idx + y * mis;

        vpx_memcpy(mi_addr, mi, sizeof(MODE_INFO));
      }
Ronald S. Bultje's avatar
Ronald S. Bultje committed
461
    }
Ronald S. Bultje's avatar
Ronald S. Bultje committed
462
  }
463
  if (bsize < BLOCK_SIZE_SB32X32) {
464
    ctx->txfm_rd_diff[ALLOW_32X32] = ctx->txfm_rd_diff[ALLOW_16X16];
Ronald S. Bultje's avatar
Ronald S. Bultje committed
465
  }
Adrian Grange's avatar
Adrian Grange committed
466

Yaowu Xu's avatar
Yaowu Xu committed
467
  if (mb_mode == I4X4_PRED) {
John Koleszar's avatar
John Koleszar committed
468 469
    for (i = 0; i < 16; i++) {
      xd->block[i].bmi.as_mode = xd->mode_info_context->bmi[i].as_mode;
Frank Galligan's avatar
Frank Galligan committed
470
      assert(xd->block[i].bmi.as_mode.first < B_MODE_COUNT);
Adrian Grange's avatar
Adrian Grange committed
471
    }
John Koleszar's avatar
John Koleszar committed
472 473 474
  } else if (mb_mode == I8X8_PRED) {
    for (i = 0; i < 16; i++) {
      xd->block[i].bmi = xd->mode_info_context->bmi[i];
Adrian Grange's avatar
Adrian Grange committed
475
    }
John Koleszar's avatar
John Koleszar committed
476 477 478 479
  } else if (mb_mode == SPLITMV) {
    vpx_memcpy(x->partition_info, &ctx->partition_info,
               sizeof(PARTITION_INFO));

480 481
    mbmi->mv[0].as_int = x->partition_info->bmi[15].mv.as_int;
    mbmi->mv[1].as_int = x->partition_info->bmi[15].second_mv.as_int;
John Koleszar's avatar
John Koleszar committed
482 483
  }

484
  x->skip = ctx->skip;
Ronald S. Bultje's avatar
Ronald S. Bultje committed
485 486 487
  if (!output_enabled)
    return;

488
  {
489
    int segment_id = mbmi->segment_id, ref_pred_flag;
Paul Wilkins's avatar
Paul Wilkins committed
490
    if (!vp9_segfeature_active(xd, segment_id, SEG_LVL_SKIP)) {
491 492 493 494
      for (i = 0; i < NB_TXFM_MODES; i++) {
        cpi->rd_tx_select_diff[i] += ctx->txfm_rd_diff[i];
      }
    }
495 496 497 498 499 500 501 502 503 504 505 506 507 508 509 510 511

    // Did the chosen reference frame match its predicted value.
    ref_pred_flag = ((xd->mode_info_context->mbmi.ref_frame ==
                      vp9_get_pred_ref(cm, xd)));
    vp9_set_pred_flag(xd, PRED_REF, ref_pred_flag);
    if (!xd->segmentation_enabled ||
        !vp9_segfeature_active(xd, segment_id, SEG_LVL_REF_FRAME) ||
        vp9_check_segref(xd, segment_id, INTRA_FRAME)  +
        vp9_check_segref(xd, segment_id, LAST_FRAME)   +
        vp9_check_segref(xd, segment_id, GOLDEN_FRAME) +
        vp9_check_segref(xd, segment_id, ALTREF_FRAME) > 1) {
      // Get the prediction context and status
      int pred_context = vp9_get_pred_context(cm, xd, PRED_REF);

      // Count prediction success
      cpi->ref_pred_count[pred_context][ref_pred_flag]++;
    }
512 513
  }

John Koleszar's avatar
John Koleszar committed
514 515
  if (cpi->common.frame_type == KEY_FRAME) {
    // Restore the coding modes to that held in the coding context
Yaowu Xu's avatar
Yaowu Xu committed
516
    // if (mb_mode == I4X4_PRED)
John Koleszar's avatar
John Koleszar committed
517 518 519 520 521 522
    //    for (i = 0; i < 16; i++)
    //    {
    //        xd->block[i].bmi.as_mode =
    //                          xd->mode_info_context->bmi[i].as_mode;
    //        assert(xd->mode_info_context->bmi[i].as_mode < MB_MODE_COUNT);
    //    }
523
#if CONFIG_INTERNAL_STATS
John Koleszar's avatar
John Koleszar committed
524 525 526 527 528 529 530 531 532 533 534 535
    static const int kf_mode_index[] = {
      THR_DC /*DC_PRED*/,
      THR_V_PRED /*V_PRED*/,
      THR_H_PRED /*H_PRED*/,
      THR_D45_PRED /*D45_PRED*/,
      THR_D135_PRED /*D135_PRED*/,
      THR_D117_PRED /*D117_PRED*/,
      THR_D153_PRED /*D153_PRED*/,
      THR_D27_PRED /*D27_PRED*/,
      THR_D63_PRED /*D63_PRED*/,
      THR_TM /*TM_PRED*/,
      THR_I8X8_PRED /*I8X8_PRED*/,
Yaowu Xu's avatar
Yaowu Xu committed
536
      THR_B_PRED /*I4X4_PRED*/,
John Koleszar's avatar
John Koleszar committed
537 538
    };
    cpi->mode_chosen_counts[kf_mode_index[mb_mode]]++;
539
#endif
John Koleszar's avatar
John Koleszar committed
540 541 542 543 544 545 546
  } else {
    /*
            // Reduce the activation RD thresholds for the best choice mode
            if ((cpi->rd_baseline_thresh[mb_mode_index] > 0) &&
                (cpi->rd_baseline_thresh[mb_mode_index] < (INT_MAX >> 2)))
            {
                int best_adjustment = (cpi->rd_thresh_mult[mb_mode_index] >> 2);
Adrian Grange's avatar
Adrian Grange committed
547

John Koleszar's avatar
John Koleszar committed
548 549 550 551 552 553 554 555
                cpi->rd_thresh_mult[mb_mode_index] =
                        (cpi->rd_thresh_mult[mb_mode_index]
                         >= (MIN_THRESHMULT + best_adjustment)) ?
                                cpi->rd_thresh_mult[mb_mode_index] - best_adjustment :
                                MIN_THRESHMULT;
                cpi->rd_threshes[mb_mode_index] =
                        (cpi->rd_baseline_thresh[mb_mode_index] >> 7)
                        * cpi->rd_thresh_mult[mb_mode_index];
Adrian Grange's avatar
Adrian Grange committed
556

John Koleszar's avatar
John Koleszar committed
557 558 559 560
            }
    */
    // Note how often each mode chosen as best
    cpi->mode_chosen_counts[mb_mode_index]++;
Deb Mukherjee's avatar
Deb Mukherjee committed
561 562 563
    if (mbmi->mode == SPLITMV || mbmi->mode == NEWMV) {
      int_mv best_mv, best_second_mv;
      MV_REFERENCE_FRAME rf = mbmi->ref_frame;
Jim Bankoski's avatar
Jim Bankoski committed
564 565
#if CONFIG_NEW_MVREF
      unsigned int best_index;
Deb Mukherjee's avatar
Deb Mukherjee committed
566
      MV_REFERENCE_FRAME sec_ref_frame = mbmi->second_ref_frame;
Jim Bankoski's avatar
Jim Bankoski committed
567
#endif
Deb Mukherjee's avatar
Deb Mukherjee committed
568 569 570 571 572 573 574 575 576
      best_mv.as_int = ctx->best_ref_mv.as_int;
      best_second_mv.as_int = ctx->second_best_ref_mv.as_int;
      if (mbmi->mode == NEWMV) {
        best_mv.as_int = mbmi->ref_mvs[rf][0].as_int;
        best_second_mv.as_int = mbmi->ref_mvs[mbmi->second_ref_frame][0].as_int;
#if CONFIG_NEW_MVREF
        best_index = pick_best_mv_ref(x, rf, mbmi->mv[0],
                                      mbmi->ref_mvs[rf], &best_mv);
        mbmi->best_index = best_index;
577
        ++cpi->mb_mv_ref_count[rf][best_index];
Deb Mukherjee's avatar
Deb Mukherjee committed
578

579
        if (mbmi->second_ref_frame > 0) {
Deb Mukherjee's avatar
Deb Mukherjee committed
580 581 582 583 584 585
          unsigned int best_index;
          best_index =
              pick_best_mv_ref(x, sec_ref_frame, mbmi->mv[1],
                               mbmi->ref_mvs[sec_ref_frame],
                               &best_second_mv);
          mbmi->best_second_index = best_index;
586
          ++cpi->mb_mv_ref_count[sec_ref_frame][best_index];
Deb Mukherjee's avatar
Deb Mukherjee committed
587 588 589 590 591 592 593
        }
#endif
      }
      mbmi->best_mv.as_int = best_mv.as_int;
      mbmi->best_second_mv.as_int = best_second_mv.as_int;
      vp9_update_nmv_count(cpi, x, &best_mv, &best_second_mv);
    }
594 595 596 597 598 599 600 601 602 603 604 605 606
#if CONFIG_COMP_INTERINTRA_PRED
    if (mbmi->mode >= NEARESTMV && mbmi->mode < SPLITMV &&
        mbmi->second_ref_frame <= INTRA_FRAME) {
      if (mbmi->second_ref_frame == INTRA_FRAME) {
        ++cpi->interintra_count[1];
        ++cpi->ymode_count[mbmi->interintra_mode];
#if SEPARATE_INTERINTRA_UV
        ++cpi->y_uv_mode_count[mbmi->interintra_mode][mbmi->interintra_uv_mode];
#endif
      } else {
        ++cpi->interintra_count[0];
      }
    }
607
#endif
608 609 610 611 612 613 614
    if (cpi->common.mcomp_filter_type == SWITCHABLE &&
        mbmi->mode >= NEARESTMV &&
        mbmi->mode <= SPLITMV) {
      ++cpi->switchable_interp_count
          [vp9_get_pred_context(&cpi->common, xd, PRED_SWITCHABLE_INTERP)]
          [vp9_switchable_interp_map[mbmi->interp_filter]];
    }
Adrian Grange's avatar
Adrian Grange committed
615

616 617 618
    cpi->rd_comp_pred_diff[SINGLE_PREDICTION_ONLY] += ctx->single_pred_diff;
    cpi->rd_comp_pred_diff[COMP_PREDICTION_ONLY]   += ctx->comp_pred_diff;
    cpi->rd_comp_pred_diff[HYBRID_PREDICTION]      += ctx->hybrid_pred_diff;
John Koleszar's avatar
John Koleszar committed
619
  }
Adrian Grange's avatar
Adrian Grange committed
620 621
}

622
static unsigned find_seg_id(uint8_t *buf, BLOCK_SIZE_TYPE bsize,
Ronald S. Bultje's avatar
Ronald S. Bultje committed
623
                            int start_y, int height, int start_x, int width) {
624 625 626
  const int bw = 1 << mb_width_log2(bsize), bh = 1 << mb_height_log2(bsize);
  const int end_x = MIN(start_x + bw, width);
  const int end_y = MIN(start_y + bh, height);
Ronald S. Bultje's avatar
Ronald S. Bultje committed
627 628 629 630 631 632 633 634 635 636 637 638 639 640
  int x, y;
  unsigned seg_id = -1;

  buf += width * start_y;
  for (y = start_y; y < end_y; y++, buf += width) {
    for (x = start_x; x < end_x; x++) {
      seg_id = MIN(seg_id, buf[x]);
    }
  }

  return seg_id;
}

static void set_offsets(VP9_COMP *cpi,
641
                        int mb_row, int mb_col, BLOCK_SIZE_TYPE bsize) {
Ronald S. Bultje's avatar
Ronald S. Bultje committed
642 643 644 645 646 647 648
  MACROBLOCK *const x = &cpi->mb;
  VP9_COMMON *const cm = &cpi->common;
  MACROBLOCKD *const xd = &x->e_mbd;
  MB_MODE_INFO *mbmi;
  const int dst_fb_idx = cm->new_fb_idx;
  const int idx_map = mb_row * cm->mb_cols + mb_col;
  const int idx_str = xd->mode_info_stride * mb_row + mb_col;
649
  const int bw = 1 << mb_width_log2(bsize), bh = 1 << mb_height_log2(bsize);
Ronald S. Bultje's avatar
Ronald S. Bultje committed
650 651 652 653 654 655 656 657 658 659 660 661 662 663 664 665 666 667 668

  // entropy context structures
  xd->above_context = cm->above_context + mb_col;
  xd->left_context  = cm->left_context + (mb_row & 3);

  // GF active flags data structure
  x->gf_active_ptr = (signed char *)&cpi->gf_active_flags[idx_map];

  // Activity map pointer
  x->mb_activity_ptr = &cpi->mb_activity_map[idx_map];
  x->active_ptr = cpi->active_map + idx_map;

  /* pointers to mode info contexts */
  x->partition_info          = x->pi + idx_str;
  xd->mode_info_context      = cm->mi + idx_str;
  mbmi = &xd->mode_info_context->mbmi;
  xd->prev_mode_info_context = cm->prev_mi + idx_str;

  // Set up destination pointers
669 670
  setup_pred_block(&xd->dst,
                   &cm->yv12_fb[dst_fb_idx],
671
                   mb_row, mb_col, NULL, NULL);
Ronald S. Bultje's avatar
Ronald S. Bultje committed
672 673 674 675 676 677

  /* Set up limit values for MV components to prevent them from
   * extending beyond the UMV borders assuming 16x16 block size */
  x->mv_row_min = -((mb_row * 16) + VP9BORDERINPIXELS - VP9_INTERP_EXTEND);
  x->mv_col_min = -((mb_col * 16) + VP9BORDERINPIXELS - VP9_INTERP_EXTEND);
  x->mv_row_max = ((cm->mb_rows - mb_row) * 16 +
678
                   (VP9BORDERINPIXELS - 16 * bh - VP9_INTERP_EXTEND));
Ronald S. Bultje's avatar
Ronald S. Bultje committed
679
  x->mv_col_max = ((cm->mb_cols - mb_col) * 16 +
680
                   (VP9BORDERINPIXELS - 16 * bw - VP9_INTERP_EXTEND));
Ronald S. Bultje's avatar
Ronald S. Bultje committed
681 682

  // Set up distance of MB to edge of frame in 1/8th pel units
683 684 685
  assert(!(mb_col & (bw - 1)) && !(mb_row & (bh - 1)));
  set_mb_row(cm, xd, mb_row, bh);
  set_mb_col(cm, xd, mb_col, bw);
Ronald S. Bultje's avatar
Ronald S. Bultje committed
686 687

  /* set up source buffers */
688
  setup_pred_block(&x->src, cpi->Source, mb_row, mb_col, NULL, NULL);
Ronald S. Bultje's avatar
Ronald S. Bultje committed
689 690 691 692 693 694 695 696

  /* R/D setup */
  x->rddiv = cpi->RDDIV;
  x->rdmult = cpi->RDMULT;

  /* segment ID */
  if (xd->segmentation_enabled) {
    if (xd->update_mb_segmentation_map) {
697
      mbmi->segment_id = find_seg_id(cpi->segmentation_map, bsize,
Ronald S. Bultje's avatar
Ronald S. Bultje committed
698 699
                                     mb_row, cm->mb_rows, mb_col, cm->mb_cols);
    } else {
700
      mbmi->segment_id = find_seg_id(cm->last_frame_seg_map, bsize,
Ronald S. Bultje's avatar
Ronald S. Bultje committed
701 702 703 704 705 706 707 708 709 710 711 712 713 714 715 716 717 718
                                     mb_row, cm->mb_rows, mb_col, cm->mb_cols);
    }
    assert(mbmi->segment_id <= 3);
    vp9_mb_init_quantizer(cpi, x);

    if (xd->segmentation_enabled && cpi->seg0_cnt > 0 &&
        !vp9_segfeature_active(xd, 0, SEG_LVL_REF_FRAME) &&
        vp9_segfeature_active(xd, 1, SEG_LVL_REF_FRAME) &&
        vp9_check_segref(xd, 1, INTRA_FRAME)  +
        vp9_check_segref(xd, 1, LAST_FRAME)   +
        vp9_check_segref(xd, 1, GOLDEN_FRAME) +
        vp9_check_segref(xd, 1, ALTREF_FRAME) == 1) {
      cpi->seg0_progress = (cpi->seg0_idx << 16) / cpi->seg0_cnt;
    } else {
      const int y = mb_row & ~3;
      const int x = mb_col & ~3;
      const int p16 = ((mb_row & 1) << 1) +  (mb_col & 1);
      const int p32 = ((mb_row & 2) << 2) + ((mb_col & 2) << 1);
719 720
      const int tile_progress = cm->cur_tile_mb_col_start * cm->mb_rows;
      const int mb_cols = cm->cur_tile_mb_col_end - cm->cur_tile_mb_col_start;
Ronald S. Bultje's avatar
Ronald S. Bultje committed
721 722

      cpi->seg0_progress =
723
          ((y * mb_cols + x * 4 + p32 + p16 + tile_progress) << 16) / cm->MBs;
Ronald S. Bultje's avatar
Ronald S. Bultje committed
724 725 726 727 728 729
    }
  } else {
    mbmi->segment_id = 0;
  }
}

730 731 732 733 734 735
static int pick_mb_mode(VP9_COMP *cpi,
                        int mb_row,
                        int mb_col,
                        TOKENEXTRA **tp,
                        int *totalrate,
                        int *totaldist) {
736 737 738
  VP9_COMMON *const cm = &cpi->common;
  MACROBLOCK *const x = &cpi->mb;
  MACROBLOCKD *const xd = &x->e_mbd;
739
  int splitmodes_used = 0;
740
  MB_MODE_INFO *mbmi;
741

742
  set_offsets(cpi, mb_row, mb_col, BLOCK_SIZE_MB16X16);
Ronald S. Bultje's avatar
Ronald S. Bultje committed
743

744 745
  if (cpi->oxcf.tuning == VP8_TUNE_SSIM)
    vp9_activity_masking(cpi, x);
John Koleszar's avatar
John Koleszar committed
746

747 748
  mbmi = &xd->mode_info_context->mbmi;
  mbmi->sb_type = BLOCK_SIZE_MB16X16;
John Koleszar's avatar
John Koleszar committed
749

750 751 752 753
  // Find best coding mode & reconstruct the MB so it is available
  // as a predictor for MBs that follow in the SB
  if (cm->frame_type == KEY_FRAME) {
    vp9_rd_pick_intra_mode(cpi, x, totalrate, totaldist);
754

755 756 757 758 759 760 761
    // Save the coding context
    vpx_memcpy(&x->mb_context[xd->sb_index][xd->mb_index].mic,
               xd->mode_info_context, sizeof(MODE_INFO));
  } else {
    vp9_pick_mode_inter_macroblock(cpi, x, mb_row, mb_col,
                                   totalrate, totaldist);
    splitmodes_used += (mbmi->mode == SPLITMV);
John Koleszar's avatar
John Koleszar committed
762

763
    if (cpi->mb.e_mbd.segmentation_enabled && mbmi->segment_id == 0) {
764
      cpi->seg0_idx++;
John Koleszar's avatar
John Koleszar committed
765 766 767
    }
  }

768
  return splitmodes_used;
Adrian Grange's avatar
Adrian Grange committed
769 770
}

771 772 773
static void pick_sb_modes(VP9_COMP *cpi, int mb_row, int mb_col,
                          TOKENEXTRA **tp, int *totalrate, int *totaldist,
                          BLOCK_SIZE_TYPE bsize, PICK_MODE_CONTEXT *ctx) {
774 775 776
  VP9_COMMON *const cm = &cpi->common;
  MACROBLOCK *const x = &cpi->mb;
  MACROBLOCKD *const xd = &x->e_mbd;
Ronald S. Bultje's avatar
Ronald S. Bultje committed
777

778 779
  set_offsets(cpi, mb_row, mb_col, bsize);
  xd->mode_info_context->mbmi.sb_type = bsize;
Ronald S. Bultje's avatar
Ronald S. Bultje committed
780
  if (cpi->oxcf.tuning == VP8_TUNE_SSIM)
781
    vp9_activity_masking(cpi, x);
Ronald S. Bultje's avatar
Ronald S. Bultje committed
782 783 784

  /* Find best coding mode & reconstruct the MB so it is available
   * as a predictor for MBs that follow in the SB */
Ronald S. Bultje's avatar
Ronald S. Bultje committed
785
  if (cm->frame_type == KEY_FRAME) {
786
    vp9_rd_pick_intra_mode_sb(cpi, x, totalrate, totaldist, bsize, ctx);
Ronald S. Bultje's avatar
Ronald S. Bultje committed
787
  } else {
788
    vp9_rd_pick_inter_mode_sb(cpi, x, mb_row, mb_col, totalrate, totaldist,
789
                              bsize, ctx);
Ronald S. Bultje's avatar
Ronald S. Bultje committed
790 791
  }
}
Adrian Grange's avatar
Adrian Grange committed
792

793
static void update_stats(VP9_COMP *cpi, int mb_row, int mb_col) {
Ronald S. Bultje's avatar
Ronald S. Bultje committed
794 795 796 797 798
  VP9_COMMON *const cm = &cpi->common;
  MACROBLOCK *const x = &cpi->mb;
  MACROBLOCKD *const xd = &x->e_mbd;
  MODE_INFO *mi = xd->mode_info_context;
  MB_MODE_INFO *const mbmi = &mi->mbmi;
Adrian Grange's avatar
Adrian Grange committed
799

Ronald S. Bultje's avatar
Ronald S. Bultje committed
800 801 802
  if (cm->frame_type == KEY_FRAME) {
#ifdef MODE_STATS
    y_modes[mbmi->mode]++;
Adrian Grange's avatar
Adrian Grange committed
803
#endif
Ronald S. Bultje's avatar
Ronald S. Bultje committed
804 805
  } else {
    int segment_id, seg_ref_active;
John Koleszar's avatar
John Koleszar committed
806

Ronald S. Bultje's avatar
Ronald S. Bultje committed
807 808
    if (mbmi->ref_frame) {
      int pred_context = vp9_get_pred_context(cm, xd, PRED_COMP);
809

Ronald S. Bultje's avatar
Ronald S. Bultje committed
810 811 812 813
      if (mbmi->second_ref_frame <= INTRA_FRAME)
        cpi->single_pred_count[pred_context]++;
      else
        cpi->comp_pred_count[pred_context]++;
Ronald S. Bultje's avatar
Ronald S. Bultje committed
814
    }
Adrian Grange's avatar
Adrian Grange committed
815

Ronald S. Bultje's avatar
Ronald S. Bultje committed
816 817
#ifdef MODE_STATS
    inter_y_modes[mbmi->mode]++;
Adrian Grange's avatar
Adrian Grange committed
818

Ronald S. Bultje's avatar
Ronald S. Bultje committed
819 820
    if (mbmi->mode == SPLITMV) {
      int b;
John Koleszar's avatar
John Koleszar committed
821

Ronald S. Bultje's avatar
Ronald S. Bultje committed
822 823 824 825
      for (b = 0; b < x->partition_info->count; b++) {
        inter_b_modes[x->partition_info->bmi[b].mode]++;
      }
    }
826
#endif
827

Ronald S. Bultje's avatar
Ronald S. Bultje committed
828 829 830 831 832 833 834 835 836 837 838 839 840 841
    // If we have just a single reference frame coded for a segment then
    // exclude from the reference frame counts used to work out
    // probabilities. NOTE: At the moment we dont support custom trees
    // for the reference frame coding for each segment but this is a
    // possible future action.
    segment_id = mbmi->segment_id;
    seg_ref_active = vp9_segfeature_active(xd, segment_id,
                                           SEG_LVL_REF_FRAME);
    if (!seg_ref_active ||
        ((vp9_check_segref(xd, segment_id, INTRA_FRAME) +
          vp9_check_segref(xd, segment_id, LAST_FRAME) +
          vp9_check_segref(xd, segment_id, GOLDEN_FRAME) +
          vp9_check_segref(xd, segment_id, ALTREF_FRAME)) > 1)) {
      cpi->count_mb_ref_frame_usage[mbmi->ref_frame]++;
Ronald S. Bultje's avatar
Ronald S. Bultje committed
842
    }
Ronald S. Bultje's avatar
Ronald S. Bultje committed
843 844 845 846
    // Count of last ref frame 0,0 usage
    if ((mbmi->mode == ZEROMV) && (mbmi->ref_frame == LAST_FRAME))
      cpi->inter_zz_count++;
  }
847 848 849
#if CONFIG_CODE_NONZEROCOUNT
  vp9_update_nzc_counts(&cpi->common, xd, mb_row, mb_col);
#endif
Ronald S. Bultje's avatar
Ronald S. Bultje committed
850
}
John Koleszar's avatar
John Koleszar committed
851

Ronald S. Bultje's avatar
Ronald S. Bultje committed
852 853 854 855 856
static void encode_sb(VP9_COMP *cpi,
                      int mb_row,
                      int mb_col,
                      int output_enabled,
                      TOKENEXTRA **tp, int is_sb) {
857 858 859
  VP9_COMMON *const cm = &cpi->common;
  MACROBLOCK *const x = &cpi->mb;
  MACROBLOCKD *const xd = &x->e_mbd;
860
  BLOCK_SIZE_TYPE bsize = BLOCK_SIZE_SB32X32;
John Koleszar's avatar
John Koleszar committed
861

Ronald S. Bultje's avatar
Ronald S. Bultje committed
862
  if (is_sb) {
863
    set_offsets(cpi, mb_row, mb_col, bsize);
864
    update_state(cpi, &x->sb32_context[xd->sb_index],
865
                 bsize, output_enabled);
866

867
    encode_superblock(cpi, tp,
868
                      output_enabled, mb_row, mb_col, bsize);
869 870
    if (output_enabled) {
      update_stats(cpi, mb_row, mb_col);
871
      cpi->partition_count[partition_plane(bsize)][PARTITION_NONE]++;
872
    }
873

Ronald S. Bultje's avatar
Ronald S. Bultje committed
874 875 876 877
    if (output_enabled) {
      (*tp<