vp9_tokenize.c 35 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 13 14 15
 */


#include <math.h>
#include <stdio.h>
#include <string.h>
#include <assert.h>
16 17
#include "vp9/encoder/vp9_onyx_int.h"
#include "vp9/encoder/vp9_tokenize.h"
John Koleszar's avatar
John Koleszar committed
18 19
#include "vpx_mem/vpx_mem.h"

20 21 22
#include "vp9/common/vp9_pred_common.h"
#include "vp9/common/vp9_seg_common.h"
#include "vp9/common/vp9_entropy.h"
23

John Koleszar's avatar
John Koleszar committed
24
/* Global event counters used for accumulating statistics across several
25
   compressions, then generating vp9_context.c = initial stats. */
John Koleszar's avatar
John Koleszar committed
26 27

#ifdef ENTROPY_STATS
28 29 30
vp9_coeff_accum context_counters_4x4[BLOCK_TYPES];
vp9_coeff_accum context_counters_8x8[BLOCK_TYPES];
vp9_coeff_accum context_counters_16x16[BLOCK_TYPES];
31
vp9_coeff_accum context_counters_32x32[BLOCK_TYPES];
32

33 34 35
extern vp9_coeff_stats tree_update_hist_4x4[BLOCK_TYPES];
extern vp9_coeff_stats tree_update_hist_8x8[BLOCK_TYPES];
extern vp9_coeff_stats tree_update_hist_16x16[BLOCK_TYPES];
36
extern vp9_coeff_stats tree_update_hist_32x32[BLOCK_TYPES];
37 38
#endif  /* ENTROPY_STATS */

39 40 41 42 43 44 45 46 47 48 49 50 51 52 53
#if CONFIG_CODE_NONZEROCOUNT
#ifdef NZC_STATS
unsigned int nzc_counts_4x4[MAX_NZC_CONTEXTS][REF_TYPES][BLOCK_TYPES]
                           [NZC4X4_TOKENS];
unsigned int nzc_counts_8x8[MAX_NZC_CONTEXTS][REF_TYPES][BLOCK_TYPES]
                           [NZC8X8_TOKENS];
unsigned int nzc_counts_16x16[MAX_NZC_CONTEXTS][REF_TYPES][BLOCK_TYPES]
                             [NZC16X16_TOKENS];
unsigned int nzc_counts_32x32[MAX_NZC_CONTEXTS][REF_TYPES][BLOCK_TYPES]
                             [NZC32X32_TOKENS];
unsigned int nzc_pcat_counts[MAX_NZC_CONTEXTS][NZC_TOKENS_EXTRA]
                            [NZC_BITS_EXTRA][2];
#endif
#endif

John Koleszar's avatar
John Koleszar committed
54
static TOKENVALUE dct_value_tokens[DCT_MAX_VALUE * 2];
55
const TOKENVALUE *vp9_dct_value_tokens_ptr;
John Koleszar's avatar
John Koleszar committed
56
static int dct_value_cost[DCT_MAX_VALUE * 2];
57
const int *vp9_dct_value_cost_ptr;
58

John Koleszar's avatar
John Koleszar committed
59
static void fill_value_tokens() {
John Koleszar's avatar
John Koleszar committed
60

John Koleszar's avatar
John Koleszar committed
61
  TOKENVALUE *const t = dct_value_tokens + DCT_MAX_VALUE;
62
  vp9_extra_bit_struct *const e = vp9_extra_bits;
John Koleszar's avatar
John Koleszar committed
63

John Koleszar's avatar
John Koleszar committed
64 65
  int i = -DCT_MAX_VALUE;
  int sign = 1;
John Koleszar's avatar
John Koleszar committed
66

John Koleszar's avatar
John Koleszar committed
67 68 69
  do {
    if (!i)
      sign = 0;
John Koleszar's avatar
John Koleszar committed
70

John Koleszar's avatar
John Koleszar committed
71 72 73
    {
      const int a = sign ? -i : i;
      int eb = sign;
John Koleszar's avatar
John Koleszar committed
74

John Koleszar's avatar
John Koleszar committed
75 76
      if (a > 4) {
        int j = 4;
John Koleszar's avatar
John Koleszar committed
77

John Koleszar's avatar
John Koleszar committed
78
        while (++j < 11  &&  e[j].base_val <= a) {}
John Koleszar's avatar
John Koleszar committed
79

John Koleszar's avatar
John Koleszar committed
80 81 82 83
        t[i].Token = --j;
        eb |= (a - e[j].base_val) << 1;
      } else
        t[i].Token = a;
John Koleszar's avatar
John Koleszar committed
84

John Koleszar's avatar
John Koleszar committed
85 86
      t[i].Extra = eb;
    }
John Koleszar's avatar
John Koleszar committed
87

John Koleszar's avatar
John Koleszar committed
88 89 90
    // initialize the cost for extra bits for all possible coefficient value.
    {
      int cost = 0;
91
      vp9_extra_bit_struct *p = vp9_extra_bits + t[i].Token;
John Koleszar's avatar
John Koleszar committed
92

John Koleszar's avatar
John Koleszar committed
93 94 95
      if (p->base_val) {
        const int extra = t[i].Extra;
        const int Length = p->Len;
John Koleszar's avatar
John Koleszar committed
96

John Koleszar's avatar
John Koleszar committed
97
        if (Length)
98
          cost += treed_cost(p->tree, p->prob, extra >> 1, Length);
John Koleszar's avatar
John Koleszar committed
99

100
        cost += vp9_cost_bit(vp9_prob_half, extra & 1); /* sign */
John Koleszar's avatar
John Koleszar committed
101 102
        dct_value_cost[i + DCT_MAX_VALUE] = cost;
      }
John Koleszar's avatar
John Koleszar committed
103 104 105

    }

John Koleszar's avatar
John Koleszar committed
106 107
  } while (++i < DCT_MAX_VALUE);

108 109
  vp9_dct_value_tokens_ptr = dct_value_tokens + DCT_MAX_VALUE;
  vp9_dct_value_cost_ptr   = dct_value_cost + DCT_MAX_VALUE;
John Koleszar's avatar
John Koleszar committed
110
}
111

112
static void tokenize_b(VP9_COMP *cpi,
113
                       MACROBLOCKD *xd,
Yaowu Xu's avatar
Yaowu Xu committed
114
                       const int ib,
115 116 117 118
                       TOKENEXTRA **tp,
                       PLANE_TYPE type,
                       TX_SIZE tx_size,
                       int dry_run) {
119
  MB_MODE_INFO *mbmi = &xd->mode_info_context->mbmi;
Daniel Kang's avatar
Daniel Kang committed
120
  int pt; /* near block/prev token context index */
121
  int c = 0;
122
  int recent_energy = 0;
123
  const int eob = xd->eobs[ib];     /* one beyond last nonzero coeff */
Daniel Kang's avatar
Daniel Kang committed
124
  TOKENEXTRA *t = *tp;        /* store tokens starting here */
125
  int16_t *qcoeff_ptr = xd->qcoeff + 16 * ib;
126
  int seg_eob;
127 128
  const int segment_id = mbmi->segment_id;
  const BLOCK_SIZE_TYPE sb_type = mbmi->sb_type;
129
  const int *scan;
130 131
  vp9_coeff_count *counts;
  vp9_coeff_probs *probs;
132 133
  const int ref = mbmi->ref_frame != INTRA_FRAME;
  ENTROPY_CONTEXT *a, *l, *a1, *l1, *a2, *l2, *a3, *l3, a_ec, l_ec;
134 135 136 137 138
#if CONFIG_CODE_NONZEROCOUNT
  int zerosleft, nzc = 0;
  if (eob == 0)
    assert(xd->nzcs[ib] == 0);
#endif
139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160

  if (sb_type == BLOCK_SIZE_SB64X64) {
    a = (ENTROPY_CONTEXT *)xd->above_context +
                                             vp9_block2above_sb64[tx_size][ib];
    l = (ENTROPY_CONTEXT *)xd->left_context + vp9_block2left_sb64[tx_size][ib];
    a1 = a + sizeof(ENTROPY_CONTEXT_PLANES) / sizeof(ENTROPY_CONTEXT);
    l1 = l + sizeof(ENTROPY_CONTEXT_PLANES) / sizeof(ENTROPY_CONTEXT);
    a2 = a1 + sizeof(ENTROPY_CONTEXT_PLANES) / sizeof(ENTROPY_CONTEXT);
    l2 = l1 + sizeof(ENTROPY_CONTEXT_PLANES) / sizeof(ENTROPY_CONTEXT);
    a3 = a2 + sizeof(ENTROPY_CONTEXT_PLANES) / sizeof(ENTROPY_CONTEXT);
    l3 = l2 + sizeof(ENTROPY_CONTEXT_PLANES) / sizeof(ENTROPY_CONTEXT);
  } else if (sb_type == BLOCK_SIZE_SB32X32) {
    a = (ENTROPY_CONTEXT *)xd->above_context + vp9_block2above_sb[tx_size][ib];
    l = (ENTROPY_CONTEXT *)xd->left_context + vp9_block2left_sb[tx_size][ib];
    a1 = a + sizeof(ENTROPY_CONTEXT_PLANES) / sizeof(ENTROPY_CONTEXT);
    l1 = l + sizeof(ENTROPY_CONTEXT_PLANES) / sizeof(ENTROPY_CONTEXT);
    a2 = a3 = l2 = l3 = NULL;
  } else {
    a = (ENTROPY_CONTEXT *)xd->above_context + vp9_block2above[tx_size][ib];
    l = (ENTROPY_CONTEXT *)xd->left_context + vp9_block2left[tx_size][ib];
    a1 = l1 = a2 = l2 = a3 = l3 = NULL;
  }
161

162 163
  switch (tx_size) {
    default:
164 165 166
    case TX_4X4: {
      const TX_TYPE tx_type = (type == PLANE_TYPE_Y_WITH_DC) ?
                              get_tx_type_4x4(xd, ib) : DCT_DCT;
167 168
      a_ec = *a;
      l_ec = *l;
169
      seg_eob = 16;
170
      scan = vp9_default_zig_zag1d_4x4;
171 172
      if (tx_type != DCT_DCT) {
        if (tx_type == ADST_DCT) {
173
          scan = vp9_row_scan_4x4;
174
        } else if (tx_type == DCT_ADST) {
175
          scan = vp9_col_scan_4x4;
176 177
        }
      }
178 179
      counts = cpi->coef_counts_4x4;
      probs = cpi->common.fc.coef_probs_4x4;
180
      break;
181
    }
182
    case TX_8X8:
183 184 185 186
      a_ec = (a[0] + a[1]) != 0;
      l_ec = (l[0] + l[1]) != 0;
      seg_eob = 64;
      scan = vp9_default_zig_zag1d_8x8;
187 188
      counts = cpi->coef_counts_8x8;
      probs = cpi->common.fc.coef_probs_8x8;
189 190
      break;
    case TX_16X16:
Yaowu Xu's avatar
Yaowu Xu committed
191 192 193
      if (type != PLANE_TYPE_UV) {
        a_ec = (a[0] + a[1] + a[2] + a[3]) != 0;
        l_ec = (l[0] + l[1] + l[2] + l[3]) != 0;
194 195 196
      } else {
        a_ec = (a[0] + a[1] + a1[0] + a1[1]) != 0;
        l_ec = (l[0] + l[1] + l1[0] + l1[1]) != 0;
Yaowu Xu's avatar
Yaowu Xu committed
197
      }
198
      seg_eob = 256;
199
      scan = vp9_default_zig_zag1d_16x16;
200 201
      counts = cpi->coef_counts_16x16;
      probs = cpi->common.fc.coef_probs_16x16;
202 203
      break;
    case TX_32X32:
204 205 206 207 208 209 210 211 212 213 214
      if (type != PLANE_TYPE_UV) {
        a_ec = (a[0] + a[1] + a[2] + a[3] +
                a1[0] + a1[1] + a1[2] + a1[3]) != 0;
        l_ec = (l[0] + l[1] + l[2] + l[3] +
                l1[0] + l1[1] + l1[2] + l1[3]) != 0;
      } else {
        a_ec = (a[0] + a[1] + a1[0] + a1[1] +
                a2[0] + a2[1] + a3[0] + a3[1]) != 0;
        l_ec = (l[0] + l[1] + l1[0] + l1[1] +
                l2[0] + l2[1] + l3[0] + l3[1]) != 0;
      }
215 216 217 218
      seg_eob = 1024;
      scan = vp9_default_zig_zag1d_32x32;
      counts = cpi->coef_counts_32x32;
      probs = cpi->common.fc.coef_probs_32x32;
219
      break;
John Koleszar's avatar
John Koleszar committed
220 221
  }

Yaowu Xu's avatar
Yaowu Xu committed
222 223
  VP9_COMBINEENTROPYCONTEXTS(pt, a_ec, l_ec);

Paul Wilkins's avatar
Paul Wilkins committed
224 225
  if (vp9_segfeature_active(xd, segment_id, SEG_LVL_SKIP))
    seg_eob = 0;
John Koleszar's avatar
John Koleszar committed
226

227
  do {
228
    const int band = get_coef_band(tx_size, c);
229
    int token;
230 231 232 233
    int v = 0;
#if CONFIG_CODE_NONZEROCOUNT
    zerosleft = seg_eob - xd->nzcs[ib] - c + nzc;
#endif
234
    if (c < eob) {
235
      const int rc = scan[c];
236
      v = qcoeff_ptr[rc];
237 238
      assert(-DCT_MAX_VALUE <= v  &&  v < DCT_MAX_VALUE);

239 240
      t->Extra = vp9_dct_value_tokens_ptr[v].Extra;
      token    = vp9_dct_value_tokens_ptr[v].Token;
241
    } else {
242 243 244
#if CONFIG_CODE_NONZEROCOUNT
      break;
#else
245
      token = DCT_EOB_TOKEN;
246
#endif
247
    }
John Koleszar's avatar
John Koleszar committed
248 249

    t->Token = token;
250
    t->context_tree = probs[type][ref][band][pt];
251 252 253 254
#if CONFIG_CODE_NONZEROCOUNT
    // Skip zero node if there are no zeros left
    t->skip_eob_node = 1 + (zerosleft == 0);
#else
255
    t->skip_eob_node = (pt == 0) && (band > 0);
256
#endif
257
    assert(vp9_coef_encodings[t->Token].Len - t->skip_eob_node > 0);
258
    if (!dry_run) {
259
      ++counts[type][ref][band][pt][token];
260
    }
261 262 263
#if CONFIG_CODE_NONZEROCOUNT
    nzc += (v != 0);
#endif
264 265

    pt = vp9_get_coef_context(&recent_energy, token);
266
    ++t;
267
  } while (c < eob && ++c < seg_eob);
268 269 270
#if CONFIG_CODE_NONZEROCOUNT
  assert(nzc == xd->nzcs[ib]);
#endif
271 272

  *tp = t;
273
  a_ec = l_ec = (c > 0); /* 0 <-> all coeff data is zero */
Yaowu Xu's avatar
Yaowu Xu committed
274 275 276
  a[0] = a_ec;
  l[0] = l_ec;

277
  if (tx_size == TX_8X8) {
Yaowu Xu's avatar
Yaowu Xu committed
278 279 280
    a[1] = a_ec;
    l[1] = l_ec;
  } else if (tx_size == TX_16X16) {
Yaowu Xu's avatar
Yaowu Xu committed
281 282 283 284
    if (type != PLANE_TYPE_UV) {
      a[1] = a[2] = a[3] = a_ec;
      l[1] = l[2] = l[3] = l_ec;
    } else {
285 286
      a1[0] = a1[1] = a[1] = a_ec;
      l1[0] = l1[1] = l[1] = l_ec;
Yaowu Xu's avatar
Yaowu Xu committed
287
    }
288
  } else if (tx_size == TX_32X32) {
289 290 291 292 293 294 295 296 297 298 299
    if (type != PLANE_TYPE_UV) {
      a[1] = a[2] = a[3] = a_ec;
      l[1] = l[2] = l[3] = l_ec;
      a1[0] = a1[1] = a1[2] = a1[3] = a_ec;
      l1[0] = l1[1] = l1[2] = l1[3] = l_ec;
    } else {
      a[1] = a1[0] = a1[1] = a_ec;
      l[1] = l1[0] = l1[1] = l_ec;
      a2[0] = a2[1] = a3[0] = a3[1] = a_ec;
      l2[0] = l2[1] = l3[0] = l3[1] = l_ec;
    }
Yaowu Xu's avatar
Yaowu Xu committed
300
  }
John Koleszar's avatar
John Koleszar committed
301 302
}

303
int vp9_mby_is_skippable_4x4(MACROBLOCKD *xd) {
John Koleszar's avatar
John Koleszar committed
304 305
  int skip = 1;
  int i = 0;
306

307
  for (i = 0; i < 16; i++)
308
    skip &= (!xd->eobs[i]);
309

John Koleszar's avatar
John Koleszar committed
310
  return skip;
311
}
312

313
int vp9_mbuv_is_skippable_4x4(MACROBLOCKD *xd) {
John Koleszar's avatar
John Koleszar committed
314 315
  int skip = 1;
  int i;
316

John Koleszar's avatar
John Koleszar committed
317
  for (i = 16; i < 24; i++)
318
    skip &= (!xd->eobs[i]);
John Koleszar's avatar
John Koleszar committed
319
  return skip;
320 321
}

322 323
static int mb_is_skippable_4x4(MACROBLOCKD *xd) {
  return (vp9_mby_is_skippable_4x4(xd) &
324
          vp9_mbuv_is_skippable_4x4(xd));
325
}
Yaowu Xu's avatar
Yaowu Xu committed
326

327
int vp9_mby_is_skippable_8x8(MACROBLOCKD *xd) {
John Koleszar's avatar
John Koleszar committed
328 329
  int skip = 1;
  int i = 0;
330

331
  for (i = 0; i < 16; i += 4)
332
    skip &= (!xd->eobs[i]);
333

John Koleszar's avatar
John Koleszar committed
334
  return skip;
335
}
336

337
int vp9_mbuv_is_skippable_8x8(MACROBLOCKD *xd) {
338
  return (!xd->eobs[16]) & (!xd->eobs[20]);
339
}
340

341 342
static int mb_is_skippable_8x8(MACROBLOCKD *xd) {
  return (vp9_mby_is_skippable_8x8(xd) &
343
          vp9_mbuv_is_skippable_8x8(xd));
Deb Mukherjee's avatar
Deb Mukherjee committed
344 345
}

346 347
static int mb_is_skippable_8x8_4x4uv(MACROBLOCKD *xd) {
  return (vp9_mby_is_skippable_8x8(xd) &
348
          vp9_mbuv_is_skippable_4x4(xd));
349
}
Yaowu Xu's avatar
Yaowu Xu committed
350

351
int vp9_mby_is_skippable_16x16(MACROBLOCKD *xd) {
352
  return (!xd->eobs[0]);
Daniel Kang's avatar
Daniel Kang committed
353 354
}

355 356
static int mb_is_skippable_16x16(MACROBLOCKD *xd) {
  return (vp9_mby_is_skippable_16x16(xd) & vp9_mbuv_is_skippable_8x8(xd));
Daniel Kang's avatar
Daniel Kang committed
357
}
358

359
int vp9_sby_is_skippable_32x32(MACROBLOCKD *xd) {
360
  return (!xd->eobs[0]);
361 362 363
}

int vp9_sbuv_is_skippable_16x16(MACROBLOCKD *xd) {
364
  return (!xd->eobs[64]) & (!xd->eobs[80]);
365 366 367 368 369 370 371
}

static int sb_is_skippable_32x32(MACROBLOCKD *xd) {
  return vp9_sby_is_skippable_32x32(xd) &&
         vp9_sbuv_is_skippable_16x16(xd);
}

372
int vp9_sby_is_skippable_16x16(MACROBLOCKD *xd) {
373 374 375 376 377 378 379 380 381 382
  int skip = 1;
  int i = 0;

  for (i = 0; i < 64; i += 16)
    skip &= (!xd->eobs[i]);

  return skip;
}

static int sb_is_skippable_16x16(MACROBLOCKD *xd) {
383
  return vp9_sby_is_skippable_16x16(xd) & vp9_sbuv_is_skippable_16x16(xd);
384 385
}

386
int vp9_sby_is_skippable_8x8(MACROBLOCKD *xd) {
387 388 389 390 391 392 393 394 395
  int skip = 1;
  int i = 0;

  for (i = 0; i < 64; i += 4)
    skip &= (!xd->eobs[i]);

  return skip;
}

396
int vp9_sbuv_is_skippable_8x8(MACROBLOCKD *xd) {
397 398 399 400 401 402 403 404 405 406
  int skip = 1;
  int i = 0;

  for (i = 64; i < 96; i += 4)
    skip &= (!xd->eobs[i]);

  return skip;
}

static int sb_is_skippable_8x8(MACROBLOCKD *xd) {
407
  return vp9_sby_is_skippable_8x8(xd) & vp9_sbuv_is_skippable_8x8(xd);
408 409
}

410
int vp9_sby_is_skippable_4x4(MACROBLOCKD *xd) {
411 412 413 414 415 416 417 418 419
  int skip = 1;
  int i = 0;

  for (i = 0; i < 64; i++)
    skip &= (!xd->eobs[i]);

  return skip;
}

420
int vp9_sbuv_is_skippable_4x4(MACROBLOCKD *xd) {
421 422 423 424 425 426 427 428 429 430
  int skip = 1;
  int i = 0;

  for (i = 64; i < 96; i++)
    skip &= (!xd->eobs[i]);

  return skip;
}

static int sb_is_skippable_4x4(MACROBLOCKD *xd) {
431
  return vp9_sby_is_skippable_4x4(xd) & vp9_sbuv_is_skippable_4x4(xd);
432 433
}

434 435 436 437 438 439 440 441 442
void vp9_tokenize_sb(VP9_COMP *cpi,
                     MACROBLOCKD *xd,
                     TOKENEXTRA **t,
                     int dry_run) {
  VP9_COMMON * const cm = &cpi->common;
  MB_MODE_INFO * const mbmi = &xd->mode_info_context->mbmi;
  TOKENEXTRA *t_backup = *t;
  const int mb_skip_context = vp9_get_pred_context(cm, xd, PRED_MBSKIP);
  const int segment_id = mbmi->segment_id;
Paul Wilkins's avatar
Paul Wilkins committed
443
  const int skip_inc = !vp9_segfeature_active(xd, segment_id, SEG_LVL_SKIP);
444 445
  int b;

446 447 448 449 450 451 452 453 454 455 456 457 458 459 460
  switch (mbmi->txfm_size) {
    case TX_32X32:
      mbmi->mb_skip_coeff = sb_is_skippable_32x32(xd);
      break;
    case TX_16X16:
      mbmi->mb_skip_coeff = sb_is_skippable_16x16(xd);
      break;
    case TX_8X8:
      mbmi->mb_skip_coeff = sb_is_skippable_8x8(xd);
      break;
    case TX_4X4:
      mbmi->mb_skip_coeff = sb_is_skippable_4x4(xd);
      break;
    default: assert(0);
  }
461 462 463 464 465 466 467

  if (mbmi->mb_skip_coeff) {
    if (!dry_run)
      cpi->skip_true_count[mb_skip_context] += skip_inc;
    if (!cm->mb_no_coeff_skip) {
      vp9_stuff_sb(cpi, xd, t, dry_run);
    } else {
468
      vp9_reset_sb_tokens_context(xd);
469 470 471 472 473 474 475 476 477
    }
    if (dry_run)
      *t = t_backup;
    return;
  }

  if (!dry_run)
    cpi->skip_false_count[mb_skip_context] += skip_inc;

478 479 480 481 482 483 484 485 486 487 488 489 490 491 492 493 494 495 496 497 498 499 500 501 502 503 504 505 506 507 508 509 510 511
  switch (mbmi->txfm_size) {
    case TX_32X32:
      tokenize_b(cpi, xd, 0, t, PLANE_TYPE_Y_WITH_DC,
                 TX_32X32, dry_run);
      for (b = 64; b < 96; b += 16)
        tokenize_b(cpi, xd, b, t, PLANE_TYPE_UV,
                   TX_16X16, dry_run);
      break;
    case TX_16X16:
      for (b = 0; b < 64; b += 16)
        tokenize_b(cpi, xd, b, t, PLANE_TYPE_Y_WITH_DC,
                   TX_16X16, dry_run);
      for (b = 64; b < 96; b += 16)
        tokenize_b(cpi, xd, b, t, PLANE_TYPE_UV,
                   TX_16X16, dry_run);
      break;
    case TX_8X8:
      for (b = 0; b < 64; b += 4)
        tokenize_b(cpi, xd, b, t, PLANE_TYPE_Y_WITH_DC,
                   TX_8X8, dry_run);
      for (b = 64; b < 96; b += 4)
        tokenize_b(cpi, xd, b, t, PLANE_TYPE_UV,
                   TX_8X8, dry_run);
      break;
    case TX_4X4:
      for (b = 0; b < 64; b++)
        tokenize_b(cpi, xd, b, t, PLANE_TYPE_Y_WITH_DC,
                   TX_4X4, dry_run);
      for (b = 64; b < 96; b++)
        tokenize_b(cpi, xd, b, t, PLANE_TYPE_UV,
                   TX_4X4, dry_run);
      break;
    default: assert(0);
  }
512

513 514 515 516
  if (dry_run)
    *t = t_backup;
}

517
int vp9_sb64y_is_skippable_32x32(MACROBLOCKD *xd) {
518 519 520 521 522 523 524 525 526 527 528 529 530 531
  int skip = 1;
  int i = 0;

  for (i = 0; i < 256; i += 64)
    skip &= (!xd->eobs[i]);

  return skip;
}

int vp9_sb64uv_is_skippable_32x32(MACROBLOCKD *xd) {
  return (!xd->eobs[256]) & (!xd->eobs[320]);
}

static int sb64_is_skippable_32x32(MACROBLOCKD *xd) {
532
  return vp9_sb64y_is_skippable_32x32(xd) & vp9_sb64uv_is_skippable_32x32(xd);
533 534
}

535
int vp9_sb64y_is_skippable_16x16(MACROBLOCKD *xd) {
536 537 538 539 540 541 542 543 544
  int skip = 1;
  int i = 0;

  for (i = 0; i < 256; i += 16)
    skip &= (!xd->eobs[i]);

  return skip;
}

545
int vp9_sb64uv_is_skippable_16x16(MACROBLOCKD *xd) {
546 547 548 549 550 551 552 553 554 555
  int skip = 1;
  int i = 0;

  for (i = 256; i < 384; i += 16)
    skip &= (!xd->eobs[i]);

  return skip;
}

static int sb64_is_skippable_16x16(MACROBLOCKD *xd) {
556
  return vp9_sb64y_is_skippable_16x16(xd) & vp9_sb64uv_is_skippable_16x16(xd);
557 558
}

559
int vp9_sb64y_is_skippable_8x8(MACROBLOCKD *xd) {
560 561 562 563 564 565 566 567 568
  int skip = 1;
  int i = 0;

  for (i = 0; i < 256; i += 4)
    skip &= (!xd->eobs[i]);

  return skip;
}

569
int vp9_sb64uv_is_skippable_8x8(MACROBLOCKD *xd) {
570 571 572 573 574 575 576 577 578 579
  int skip = 1;
  int i = 0;

  for (i = 256; i < 384; i += 4)
    skip &= (!xd->eobs[i]);

  return skip;
}

static int sb64_is_skippable_8x8(MACROBLOCKD *xd) {
580
  return vp9_sb64y_is_skippable_8x8(xd) & vp9_sb64uv_is_skippable_8x8(xd);
581 582
}

583
int vp9_sb64y_is_skippable_4x4(MACROBLOCKD *xd) {
584 585 586 587 588 589 590 591 592
  int skip = 1;
  int i = 0;

  for (i = 0; i < 256; i++)
    skip &= (!xd->eobs[i]);

  return skip;
}

593
int vp9_sb64uv_is_skippable_4x4(MACROBLOCKD *xd) {
594 595 596 597 598 599 600 601 602 603
  int skip = 1;
  int i = 0;

  for (i = 256; i < 384; i++)
    skip &= (!xd->eobs[i]);

  return skip;
}

static int sb64_is_skippable_4x4(MACROBLOCKD *xd) {
604
  return vp9_sb64y_is_skippable_4x4(xd) & vp9_sb64uv_is_skippable_4x4(xd);
605 606 607 608 609 610 611 612 613 614 615 616 617 618 619 620 621 622 623 624 625 626 627 628 629 630 631 632 633 634 635 636 637 638 639 640 641 642 643 644 645 646 647 648 649 650 651 652 653 654 655 656 657 658 659 660 661 662 663 664 665 666 667 668 669 670 671 672 673 674 675 676 677 678 679 680 681 682 683 684
}

void vp9_tokenize_sb64(VP9_COMP *cpi,
                       MACROBLOCKD *xd,
                       TOKENEXTRA **t,
                       int dry_run) {
  VP9_COMMON * const cm = &cpi->common;
  MB_MODE_INFO * const mbmi = &xd->mode_info_context->mbmi;
  TOKENEXTRA *t_backup = *t;
  const int mb_skip_context = vp9_get_pred_context(cm, xd, PRED_MBSKIP);
  const int segment_id = mbmi->segment_id;
  const int skip_inc = !vp9_segfeature_active(xd, segment_id, SEG_LVL_SKIP);
  int b;

  switch (mbmi->txfm_size) {
    case TX_32X32:
      mbmi->mb_skip_coeff = sb64_is_skippable_32x32(xd);
      break;
    case TX_16X16:
      mbmi->mb_skip_coeff = sb64_is_skippable_16x16(xd);
      break;
    case TX_8X8:
      mbmi->mb_skip_coeff = sb64_is_skippable_8x8(xd);
      break;
    case TX_4X4:
      mbmi->mb_skip_coeff = sb64_is_skippable_4x4(xd);
      break;
    default: assert(0);
  }

  if (mbmi->mb_skip_coeff) {
    if (!dry_run)
      cpi->skip_true_count[mb_skip_context] += skip_inc;
    if (!cm->mb_no_coeff_skip) {
      vp9_stuff_sb64(cpi, xd, t, dry_run);
    } else {
      vp9_reset_sb64_tokens_context(xd);
    }
    if (dry_run)
      *t = t_backup;
    return;
  }

  if (!dry_run)
    cpi->skip_false_count[mb_skip_context] += skip_inc;

  switch (mbmi->txfm_size) {
    case TX_32X32:
      for (b = 0; b < 256; b += 64)
        tokenize_b(cpi, xd, b, t, PLANE_TYPE_Y_WITH_DC,
                   TX_32X32, dry_run);
      for (b = 256; b < 384; b += 64)
        tokenize_b(cpi, xd, b, t, PLANE_TYPE_UV,
                   TX_32X32, dry_run);
      break;
    case TX_16X16:
      for (b = 0; b < 256; b += 16)
        tokenize_b(cpi, xd, b, t, PLANE_TYPE_Y_WITH_DC,
                   TX_16X16, dry_run);
      for (b = 256; b < 384; b += 16)
        tokenize_b(cpi, xd, b, t, PLANE_TYPE_UV,
                   TX_16X16, dry_run);
      break;
    case TX_8X8:
      for (b = 0; b < 256; b += 4)
        tokenize_b(cpi, xd, b, t, PLANE_TYPE_Y_WITH_DC,
                   TX_8X8, dry_run);
      for (b = 256; b < 384; b += 4)
        tokenize_b(cpi, xd, b, t, PLANE_TYPE_UV,
                   TX_8X8, dry_run);
      break;
    case TX_4X4:
      for (b = 0; b < 256; b++)
        tokenize_b(cpi, xd, b, t, PLANE_TYPE_Y_WITH_DC,
                   TX_4X4, dry_run);
      for (b = 256; b < 384; b++)
        tokenize_b(cpi, xd, b, t, PLANE_TYPE_UV,
                   TX_4X4, dry_run);
      break;
    default: assert(0);
685
  }
686

687 688 689 690
  if (dry_run)
    *t = t_backup;
}

691
void vp9_tokenize_mb(VP9_COMP *cpi,
Paul Wilkins's avatar
Paul Wilkins committed
692
                     MACROBLOCKD *xd,
693 694
                     TOKENEXTRA **t,
                     int dry_run) {
John Koleszar's avatar
John Koleszar committed
695
  int b;
696
  int tx_size = xd->mode_info_context->mbmi.txfm_size;
Paul Wilkins's avatar
Paul Wilkins committed
697
  int mb_skip_context = vp9_get_pred_context(&cpi->common, xd, PRED_MBSKIP);
698
  TOKENEXTRA *t_backup = *t;
Paul Wilkins's avatar
Paul Wilkins committed
699

John Koleszar's avatar
John Koleszar committed
700 701 702 703
  // If the MB is going to be skipped because of a segment level flag
  // exclude this from the skip count stats used to calculate the
  // transmitted skip probability;
  int skip_inc;
Paul Wilkins's avatar
Paul Wilkins committed
704
  int segment_id = xd->mode_info_context->mbmi.segment_id;
705

Paul Wilkins's avatar
Paul Wilkins committed
706
  if (!vp9_segfeature_active(xd, segment_id, SEG_LVL_SKIP)) {
John Koleszar's avatar
John Koleszar committed
707 708 709
    skip_inc = 1;
  } else
    skip_inc = 0;
710

711
  switch (tx_size) {
Daniel Kang's avatar
Daniel Kang committed
712
    case TX_16X16:
713

Daniel Kang's avatar
Daniel Kang committed
714
      xd->mode_info_context->mbmi.mb_skip_coeff = mb_is_skippable_16x16(xd);
Daniel Kang's avatar
Daniel Kang committed
715 716
      break;
    case TX_8X8:
717 718
      if (xd->mode_info_context->mbmi.mode == I8X8_PRED ||
          xd->mode_info_context->mbmi.mode == SPLITMV)
719
        xd->mode_info_context->mbmi.mb_skip_coeff =
720
            mb_is_skippable_8x8_4x4uv(xd);
Deb Mukherjee's avatar
Deb Mukherjee committed
721
      else
722
        xd->mode_info_context->mbmi.mb_skip_coeff =
723
            mb_is_skippable_8x8(xd);
Daniel Kang's avatar
Daniel Kang committed
724
      break;
Deb Mukherjee's avatar
Deb Mukherjee committed
725

Daniel Kang's avatar
Daniel Kang committed
726
    default:
727
      xd->mode_info_context->mbmi.mb_skip_coeff =
728
          mb_is_skippable_4x4(xd);
Daniel Kang's avatar
Daniel Kang committed
729 730
      break;
  }
731

Paul Wilkins's avatar
Paul Wilkins committed
732
  if (xd->mode_info_context->mbmi.mb_skip_coeff) {
733 734
    if (!dry_run)
      cpi->skip_true_count[mb_skip_context] += skip_inc;
John Koleszar's avatar
John Koleszar committed
735
    if (!cpi->common.mb_no_coeff_skip) {
736
      vp9_stuff_mb(cpi, xd, t, dry_run);
John Koleszar's avatar
John Koleszar committed
737
    } else {
Yaowu Xu's avatar
Yaowu Xu committed
738
      vp9_reset_mb_tokens_context(xd);
John Koleszar's avatar
John Koleszar committed
739
    }
Yaowu Xu's avatar
Yaowu Xu committed
740

741 742
    if (dry_run)
      *t = t_backup;
John Koleszar's avatar
John Koleszar committed
743 744 745
    return;
  }

746 747
  if (!dry_run)
    cpi->skip_false_count[mb_skip_context] += skip_inc;
748

749
  if (tx_size == TX_16X16) {
Yaowu Xu's avatar
Yaowu Xu committed
750
    tokenize_b(cpi, xd, 0, t, PLANE_TYPE_Y_WITH_DC, TX_16X16, dry_run);
Daniel Kang's avatar
Daniel Kang committed
751
    for (b = 16; b < 24; b += 4) {
Yaowu Xu's avatar
Yaowu Xu committed
752
      tokenize_b(cpi, xd, b, t, PLANE_TYPE_UV, TX_8X8, dry_run);
Daniel Kang's avatar
Daniel Kang committed
753
    }
754
  } else if (tx_size == TX_8X8) {
John Koleszar's avatar
John Koleszar committed
755
    for (b = 0; b < 16; b += 4) {
756
      tokenize_b(cpi, xd, b, t, PLANE_TYPE_Y_WITH_DC, TX_8X8, dry_run);
John Koleszar's avatar
John Koleszar committed
757
    }
758 759
    if (xd->mode_info_context->mbmi.mode == I8X8_PRED ||
        xd->mode_info_context->mbmi.mode == SPLITMV) {
760
      for (b = 16; b < 24; b++) {
Yaowu Xu's avatar
Yaowu Xu committed
761
        tokenize_b(cpi, xd, b, t, PLANE_TYPE_UV, TX_4X4, dry_run);
762
      }
763
    } else {
Deb Mukherjee's avatar
Deb Mukherjee committed
764
      for (b = 16; b < 24; b += 4) {
Yaowu Xu's avatar
Yaowu Xu committed
765
        tokenize_b(cpi, xd, b, t, PLANE_TYPE_UV, TX_8X8, dry_run);
Deb Mukherjee's avatar
Deb Mukherjee committed
766
      }
767
    }
Jingning Han's avatar
Jingning Han committed
768
  } else {
769 770 771 772
    for (b = 0; b < 16; b++)
      tokenize_b(cpi, xd, b, t, PLANE_TYPE_Y_WITH_DC, TX_4X4, dry_run);
    for (b = 16; b < 24; b++)
      tokenize_b(cpi, xd, b, t, PLANE_TYPE_UV, TX_4X4, dry_run);
Jingning Han's avatar
Jingning Han committed
773
  }
774 775
  if (dry_run)
    *t = t_backup;
John Koleszar's avatar
John Koleszar committed
776
}
777

John Koleszar's avatar
John Koleszar committed
778
#ifdef ENTROPY_STATS
John Koleszar's avatar
John Koleszar committed
779 780 781
void init_context_counters(void) {
  FILE *f = fopen("context.bin", "rb");
  if (!f) {
782
    vpx_memset(context_counters_4x4, 0, sizeof(context_counters_4x4));
John Koleszar's avatar
John Koleszar committed
783
    vpx_memset(context_counters_8x8, 0, sizeof(context_counters_8x8));
Daniel Kang's avatar
Daniel Kang committed
784
    vpx_memset(context_counters_16x16, 0, sizeof(context_counters_16x16));
785
    vpx_memset(context_counters_32x32, 0, sizeof(context_counters_32x32));
John Koleszar's avatar
John Koleszar committed
786
  } else {
787
    fread(context_counters_4x4, sizeof(context_counters_4x4), 1, f);
John Koleszar's avatar
John Koleszar committed
788
    fread(context_counters_8x8, sizeof(context_counters_8x8), 1, f);
Daniel Kang's avatar
Daniel Kang committed
789
    fread(context_counters_16x16, sizeof(context_counters_16x16), 1, f);
790
    fread(context_counters_32x32, sizeof(context_counters_32x32), 1, f);
John Koleszar's avatar
John Koleszar committed
791 792 793 794 795
    fclose(f);
  }

  f = fopen("treeupdate.bin", "rb");
  if (!f) {
796
    vpx_memset(tree_update_hist_4x4, 0, sizeof(tree_update_hist_4x4));
John Koleszar's avatar
John Koleszar committed
797
    vpx_memset(tree_update_hist_8x8, 0, sizeof(tree_update_hist_8x8));
Daniel Kang's avatar
Daniel Kang committed
798
    vpx_memset(tree_update_hist_16x16, 0, sizeof(tree_update_hist_16x16));
799
    vpx_memset(tree_update_hist_32x32, 0, sizeof(tree_update_hist_32x32));
John Koleszar's avatar
John Koleszar committed
800
  } else {
801
    fread(tree_update_hist_4x4, sizeof(tree_update_hist_4x4), 1, f);
John Koleszar's avatar
John Koleszar committed
802
    fread(tree_update_hist_8x8, sizeof(tree_update_hist_8x8), 1, f);
Daniel Kang's avatar
Daniel Kang committed
803
    fread(tree_update_hist_16x16, sizeof(tree_update_hist_16x16), 1, f);
804
    fread(tree_update_hist_32x32, sizeof(tree_update_hist_32x32), 1, f);
John Koleszar's avatar
John Koleszar committed
805 806
    fclose(f);
  }
John Koleszar's avatar
John Koleszar committed
807 808
}

809 810
static void print_counter(FILE *f, vp9_coeff_accum *context_counters,
                          int block_types, const char *header) {
811
  int type, ref, band, pt, t;
John Koleszar's avatar
John Koleszar committed
812

813
  fprintf(f, "static const vp9_coeff_count %s = {\n", header);
John Koleszar's avatar
John Koleszar committed
814

815
#define Comma(X) (X ? "," : "")
John Koleszar's avatar
John Koleszar committed
816 817
  type = 0;
  do {
818
    ref = 0;
John Koleszar's avatar
John Koleszar committed
819 820
    fprintf(f, "%s\n  { /* block Type %d */", Comma(type), type);
    do {
821 822
      fprintf(f, "%s\n    { /* %s */", Comma(type), ref ? "Inter" : "Intra");
      band = 0;
John Koleszar's avatar
John Koleszar committed
823
      do {
824 825
        fprintf(f, "%s\n      { /* Coeff Band %d */", Comma(band), band);
        pt = 0;
John Koleszar's avatar
John Koleszar committed
826
        do {
827 828 829 830 831 832 833 834 835 836 837 838 839 840
          fprintf(f, "%s\n        {", Comma(pt));

          t = 0;
          do {
            const int64_t x = context_counters[type][ref][band][pt][t];
            const int y = (int) x;

            assert(x == (int64_t) y);  /* no overflow handling yet */
            fprintf(f, "%s %d", Comma(t), y);
          } while (++t < MAX_ENTROPY_TOKENS);
          fprintf(f, "}");
        } while (++pt < PREV_COEF_CONTEXTS);
        fprintf(f, "\n      }");
      } while (++band < COEF_BANDS);
John Koleszar's avatar
John Koleszar committed
841
      fprintf(f, "\n    }");
842
    } while (++ref < REF_TYPES);
John Koleszar's avatar
John Koleszar committed
843
    fprintf(f, "\n  }");
844
  } while (++type < block_types);
Daniel Kang's avatar
Daniel Kang committed
845
  fprintf(f, "\n};\n");
846
}
847

848 849
static void print_probs(FILE *f, vp9_coeff_accum *context_counters,
                        int block_types, const char *header) {
850
  int type, ref, band, pt, t;
Daniel Kang's avatar
Daniel Kang committed
851

852
  fprintf(f, "static const vp9_coeff_probs %s = {", header);
853

John Koleszar's avatar
John Koleszar committed
854
  type = 0;
855
#define Newline(x, spaces) (x ? " " : "\n" spaces)
John Koleszar's avatar
John Koleszar committed
856
  do {
857 858
    fprintf(f, "%s%s{ /* block Type %d */",
            Comma(type), Newline(type, "  "), type);
859
    ref = 0;
John Koleszar's avatar
John Koleszar committed
860
    do {
861 862 863
      fprintf(f, "%s%s{ /* %s */",
              Comma(band), Newline(band, "    "), ref ? "Inter" : "Intra");
      band = 0;
John Koleszar's avatar
John Koleszar committed
864
      do {
865 866 867
        fprintf(f, "%s%s{ /* Coeff Band %d */",
                Comma(band), Newline(band, "      "), band);
        pt = 0;