vp9_blockd.h 29.7 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
#ifndef VP9_COMMON_VP9_BLOCKD_H_
#define VP9_COMMON_VP9_BLOCKD_H_
John Koleszar's avatar
John Koleszar committed
14

15
#include "./vpx_config.h"
John Koleszar's avatar
John Koleszar committed
16
#include "vpx_scale/yv12config.h"
17
#include "vp9/common/vp9_convolve.h"
18 19
#include "vp9/common/vp9_mv.h"
#include "vp9/common/vp9_treecoder.h"
John Koleszar's avatar
John Koleszar committed
20
#include "vpx_ports/mem.h"
21
#include "vp9/common/vp9_common.h"
22
#include "vp9/common/vp9_enums.h"
John Koleszar's avatar
John Koleszar committed
23

John Koleszar's avatar
John Koleszar committed
24
// #define MODE_STATS
25

26 27
#define MAX_MB_SEGMENTS     8
#define MB_SEG_TREE_PROBS   (MAX_MB_SEGMENTS-1)
28
#define PREDICTION_PROBS 3
29

30 31 32 33
#define DEFAULT_PRED_PROB_0 120
#define DEFAULT_PRED_PROB_1 80
#define DEFAULT_PRED_PROB_2 40

34 35
#define MBSKIP_CONTEXTS 3

John Koleszar's avatar
John Koleszar committed
36 37 38
#define MAX_REF_LF_DELTAS       4
#define MAX_MODE_LF_DELTAS      4

39
/* Segment Feature Masks */
John Koleszar's avatar
John Koleszar committed
40 41
#define SEGMENT_DELTADATA   0
#define SEGMENT_ABSDATA     1
Paul Wilkins's avatar
Paul Wilkins committed
42
#define MAX_MV_REFS 9
Paul Wilkins's avatar
Paul Wilkins committed
43
#define MAX_MV_REF_CANDIDATES 2
John Koleszar's avatar
John Koleszar committed
44

Dmitry Kovalev's avatar
Dmitry Kovalev committed
45
typedef enum {
46
  PLANE_TYPE_Y_WITH_DC,
47
  PLANE_TYPE_UV,
48
} PLANE_TYPE;
Scott LaVarnway's avatar
Scott LaVarnway committed
49

50
typedef char ENTROPY_CONTEXT;
John Koleszar's avatar
John Koleszar committed
51

52 53
typedef char PARTITION_CONTEXT;

54 55 56 57
static INLINE int combine_entropy_contexts(ENTROPY_CONTEXT a,
                                           ENTROPY_CONTEXT b) {
  return (a != 0) + (b != 0);
}
John Koleszar's avatar
John Koleszar committed
58

John Koleszar's avatar
John Koleszar committed
59 60 61
typedef enum {
  KEY_FRAME = 0,
  INTER_FRAME = 1
John Koleszar's avatar
John Koleszar committed
62 63
} FRAME_TYPE;

Dmitry Kovalev's avatar
Dmitry Kovalev committed
64
typedef enum {
Adrian Grange's avatar
Adrian Grange committed
65 66 67 68
  EIGHTTAP_SMOOTH,
  EIGHTTAP,
  EIGHTTAP_SHARP,
  BILINEAR,
69 70 71
  SWITCHABLE  /* should be the last one */
} INTERPOLATIONFILTERTYPE;

Dmitry Kovalev's avatar
Dmitry Kovalev committed
72
typedef enum {
John Koleszar's avatar
John Koleszar committed
73 74 75 76 77 78 79 80 81 82
  DC_PRED,            /* average of above and left pixels */
  V_PRED,             /* vertical prediction */
  H_PRED,             /* horizontal prediction */
  D45_PRED,           /* Directional 45 deg prediction  [anti-clockwise from 0 deg hor] */
  D135_PRED,          /* Directional 135 deg prediction [anti-clockwise from 0 deg hor] */
  D117_PRED,          /* Directional 112 deg prediction [anti-clockwise from 0 deg hor] */
  D153_PRED,          /* Directional 157 deg prediction [anti-clockwise from 0 deg hor] */
  D27_PRED,           /* Directional 22 deg prediction  [anti-clockwise from 0 deg hor] */
  D63_PRED,           /* Directional 67 deg prediction  [anti-clockwise from 0 deg hor] */
  TM_PRED,            /* Truemotion prediction */
Yaowu Xu's avatar
Yaowu Xu committed
83
  I4X4_PRED,          /* 4x4 based prediction, each 4x4 has its own mode */
John Koleszar's avatar
John Koleszar committed
84 85 86 87 88 89
  NEARESTMV,
  NEARMV,
  ZEROMV,
  NEWMV,
  SPLITMV,
  MB_MODE_COUNT
John Koleszar's avatar
John Koleszar committed
90 91
} MB_PREDICTION_MODE;

92 93 94 95 96
static INLINE int is_inter_mode(MB_PREDICTION_MODE mode) {
  return mode >= NEARESTMV && mode <= SPLITMV;
}


Paul Wilkins's avatar
Paul Wilkins committed
97
// Segment level features.
John Koleszar's avatar
John Koleszar committed
98 99 100 101
typedef enum {
  SEG_LVL_ALT_Q = 0,               // Use alternate Quantizer ....
  SEG_LVL_ALT_LF = 1,              // Use alternate loop filter value...
  SEG_LVL_REF_FRAME = 2,           // Optional Segment reference frame
Paul Wilkins's avatar
Paul Wilkins committed
102
  SEG_LVL_SKIP = 3,                // Optional Segment (0,0) + skip mode
Paul Wilkins's avatar
Paul Wilkins committed
103
  SEG_LVL_MAX = 4                  // Number of MB level features supported
Paul Wilkins's avatar
Paul Wilkins committed
104
} SEG_LVL_FEATURES;
John Koleszar's avatar
John Koleszar committed
105

Paul Wilkins's avatar
Paul Wilkins committed
106
// Segment level features.
John Koleszar's avatar
John Koleszar committed
107
typedef enum {
Yaowu Xu's avatar
Yaowu Xu committed
108 109 110 111 112 113
  TX_4X4 = 0,                      // 4x4 dct transform
  TX_8X8 = 1,                      // 8x8 dct transform
  TX_16X16 = 2,                    // 16x16 dct transform
  TX_SIZE_MAX_MB = 3,              // Number of different transforms available
  TX_32X32 = TX_SIZE_MAX_MB,       // 32x32 dct transform
  TX_SIZE_MAX_SB,                  // Number of transforms available to SBs
Paul Wilkins's avatar
Paul Wilkins committed
114 115
} TX_SIZE;

Jingning Han's avatar
Jingning Han committed
116 117
typedef enum {
  DCT_DCT   = 0,                      // DCT  in both horizontal and vertical
118 119
  ADST_DCT  = 1,                      // ADST in vertical, DCT in horizontal
  DCT_ADST  = 2,                      // DCT  in vertical, ADST in horizontal
Jingning Han's avatar
Jingning Han committed
120 121 122
  ADST_ADST = 3                       // ADST in both directions
} TX_TYPE;

Yaowu Xu's avatar
Yaowu Xu committed
123
#define VP9_YMODES  (I4X4_PRED + 1)
124 125
#define VP9_UV_MODES (TM_PRED + 1)
#define VP9_I32X32_MODES (TM_PRED + 1)
John Koleszar's avatar
John Koleszar committed
126

127
#define VP9_MVREFS (1 + SPLITMV - NEARESTMV)
John Koleszar's avatar
John Koleszar committed
128

Yaowu Xu's avatar
Yaowu Xu committed
129
#define WHT_UPSCALE_FACTOR 2
130

John Koleszar's avatar
John Koleszar committed
131 132
typedef enum {
  B_DC_PRED,          /* average of above and left pixels */
Yaowu Xu's avatar
Yaowu Xu committed
133 134 135 136 137 138 139 140
  B_V_PRED,          /* vertical prediction */
  B_H_PRED,          /* horizontal prediction */
  B_D45_PRED,
  B_D135_PRED,
  B_D117_PRED,
  B_D153_PRED,
  B_D27_PRED,
  B_D63_PRED,
John Koleszar's avatar
John Koleszar committed
141
  B_TM_PRED,
John Koleszar's avatar
John Koleszar committed
142

John Koleszar's avatar
John Koleszar committed
143 144 145 146
  LEFT4X4,
  ABOVE4X4,
  ZERO4X4,
  NEW4X4,
John Koleszar's avatar
John Koleszar committed
147

John Koleszar's avatar
John Koleszar committed
148
  B_MODE_COUNT
John Koleszar's avatar
John Koleszar committed
149 150
} B_PREDICTION_MODE;

151
#define VP9_BINTRAMODES (LEFT4X4)
152
#define VP9_SUBMVREFS (1 + NEW4X4 - LEFT4X4)
John Koleszar's avatar
John Koleszar committed
153

154 155 156
#define VP9_KF_BINTRAMODES (VP9_BINTRAMODES)   /* 10 */
#define VP9_NKF_BINTRAMODES (VP9_BINTRAMODES)  /* 10 */

John Koleszar's avatar
John Koleszar committed
157 158 159 160
/* For keyframes, intra block modes are predicted by the (already decoded)
   modes for the Y blocks to the left and above us; for interframes, there
   is a single probability table. */

John Koleszar's avatar
John Koleszar committed
161 162 163 164
union b_mode_info {
  struct {
    B_PREDICTION_MODE first;
  } as_mode;
165
  int_mv as_mv[2];  // first, second inter predictor motion vectors
166
};
John Koleszar's avatar
John Koleszar committed
167

John Koleszar's avatar
John Koleszar committed
168
typedef enum {
169
  NONE = -1,
John Koleszar's avatar
John Koleszar committed
170 171 172 173 174
  INTRA_FRAME = 0,
  LAST_FRAME = 1,
  GOLDEN_FRAME = 2,
  ALTREF_FRAME = 3,
  MAX_REF_FRAMES = 4
John Koleszar's avatar
John Koleszar committed
175 176
} MV_REFERENCE_FRAME;

177
static INLINE int b_width_log2(BLOCK_SIZE_TYPE sb_type) {
178
  switch (sb_type) {
179 180 181
#if CONFIG_AB4X4
    case BLOCK_SIZE_SB4X8:
#endif
182
    case BLOCK_SIZE_AB4X4: return 0;
183 184 185
#if CONFIG_AB4X4
    case BLOCK_SIZE_SB8X4:
#endif
186 187
    case BLOCK_SIZE_SB8X8:
    case BLOCK_SIZE_SB8X16: return 1;
188
    case BLOCK_SIZE_SB16X8:
189
    case BLOCK_SIZE_MB16X16:
190
    case BLOCK_SIZE_SB16X32: return 2;
191
    case BLOCK_SIZE_SB32X16:
192 193
    case BLOCK_SIZE_SB32X32:
    case BLOCK_SIZE_SB32X64: return 3;
194
    case BLOCK_SIZE_SB64X32:
195
    case BLOCK_SIZE_SB64X64: return 4;
196 197 198 199
    default: assert(0);
  }
}

200
static INLINE int b_height_log2(BLOCK_SIZE_TYPE sb_type) {
201
  switch (sb_type) {
202 203 204
#if CONFIG_AB4X4
    case BLOCK_SIZE_SB8X4:
#endif
205
    case BLOCK_SIZE_AB4X4: return 0;
206 207 208
#if CONFIG_AB4X4
    case BLOCK_SIZE_SB4X8:
#endif
209 210
    case BLOCK_SIZE_SB8X8:
    case BLOCK_SIZE_SB16X8: return 1;
211
    case BLOCK_SIZE_SB8X16:
212
    case BLOCK_SIZE_MB16X16:
213
    case BLOCK_SIZE_SB32X16: return 2;
214
    case BLOCK_SIZE_SB16X32:
215 216
    case BLOCK_SIZE_SB32X32:
    case BLOCK_SIZE_SB64X32: return 3;
217
    case BLOCK_SIZE_SB32X64:
218
    case BLOCK_SIZE_SB64X64: return 4;
219 220 221
    default: assert(0);
  }
}
Ronald S. Bultje's avatar
Ronald S. Bultje committed
222

223 224 225 226
static INLINE int mi_width_log2(BLOCK_SIZE_TYPE sb_type) {
  int a = b_width_log2(sb_type) - 1;
  assert(a >= 0);
  return a;
227 228
}

229 230 231 232
static INLINE int mi_height_log2(BLOCK_SIZE_TYPE sb_type) {
  int a = b_height_log2(sb_type) - 1;
  assert(a >= 0);
  return a;
233 234
}

John Koleszar's avatar
John Koleszar committed
235 236 237 238
typedef struct {
  MB_PREDICTION_MODE mode, uv_mode;
  MV_REFERENCE_FRAME ref_frame, second_ref_frame;
  TX_SIZE txfm_size;
239
  int_mv mv[2]; // for each reference frame used
Paul Wilkins's avatar
Paul Wilkins committed
240
  int_mv ref_mvs[MAX_REF_FRAMES][MAX_MV_REF_CANDIDATES];
Deb Mukherjee's avatar
Deb Mukherjee committed
241
  int_mv best_mv, best_second_mv;
Paul Wilkins's avatar
Paul Wilkins committed
242

Paul Wilkins's avatar
Paul Wilkins committed
243 244
  int mb_mode_context[MAX_REF_FRAMES];

John Koleszar's avatar
John Koleszar committed
245 246 247
  unsigned char mb_skip_coeff;                                /* does this mb has coefficients at all, 1=no coefficients, 0=need decode tokens */
  unsigned char need_to_clamp_mvs;
  unsigned char need_to_clamp_secondmv;
248
  unsigned char segment_id;           // Segment id for current frame
John Koleszar's avatar
John Koleszar committed
249 250 251 252 253 254 255 256 257

  // Flags used for prediction status of various bistream signals
  unsigned char seg_id_predicted;
  unsigned char ref_predicted;

  // Indicates if the mb is part of the image (1) vs border (0)
  // This can be useful in determining whether the MB provides
  // a valid predictor
  unsigned char mb_in_image;
258

259
  INTERPOLATIONFILTERTYPE interp_filter;
260

Ronald S. Bultje's avatar
Ronald S. Bultje committed
261
  BLOCK_SIZE_TYPE sb_type;
John Koleszar's avatar
John Koleszar committed
262 263
} MB_MODE_INFO;

John Koleszar's avatar
John Koleszar committed
264 265
typedef struct {
  MB_MODE_INFO mbmi;
Jingning Han's avatar
Jingning Han committed
266
  union b_mode_info bmi[4];
John Koleszar's avatar
John Koleszar committed
267 268
} MODE_INFO;

269 270 271 272
struct scale_factors {
  int x_num;
  int x_den;
  int x_offset_q4;
273
  int x_step_q4;
274 275 276
  int y_num;
  int y_den;
  int y_offset_q4;
277
  int y_step_q4;
278 279 280 281 282 283 284 285 286 287 288

  int (*scale_value_x)(int val, const struct scale_factors *scale);
  int (*scale_value_y)(int val, const struct scale_factors *scale);
  void (*set_scaled_offsets)(struct scale_factors *scale, int row, int col);
  int_mv32 (*scale_motion_vector_q3_to_q4)(const int_mv *src_mv,
                                           const struct scale_factors *scale);
  int32_t (*scale_motion_vector_component_q4)(int mv_q4,
                                              int num,
                                              int den,
                                              int offset_q4);

289
  convolve_fn_t predict[2][2][2];  // horiz, vert, avg
290 291
};

292 293
enum { MAX_MB_PLANE = 3 };

294 295 296 297 298
struct buf_2d {
  uint8_t *buf;
  int stride;
};

299
struct macroblockd_plane {
300 301
  DECLARE_ALIGNED(16, int16_t,  qcoeff[64 * 64]);
  DECLARE_ALIGNED(16, int16_t,  dqcoeff[64 * 64]);
John Koleszar's avatar
John Koleszar committed
302
  DECLARE_ALIGNED(16, uint16_t, eobs[256]);
303
  DECLARE_ALIGNED(16, int16_t,  diff[64 * 64]);
304 305 306
  PLANE_TYPE plane_type;
  int subsampling_x;
  int subsampling_y;
307 308
  struct buf_2d dst;
  struct buf_2d pre[2];
309
  int16_t *dequant;
310 311
  ENTROPY_CONTEXT *above_context;
  ENTROPY_CONTEXT *left_context;
312 313 314 315 316 317 318 319 320
};

#define BLOCK_OFFSET(x, i, n) ((x) + (i) * (n))

#define MB_SUBBLOCK_FIELD(x, field, i) (\
  ((i) < 16) ? BLOCK_OFFSET((x)->plane[0].field, (i), 16) : \
  ((i) < 20) ? BLOCK_OFFSET((x)->plane[1].field, ((i) - 16), 16) : \
  BLOCK_OFFSET((x)->plane[2].field, ((i) - 20), 16))

Jim Bankoski's avatar
Jim Bankoski committed
321
typedef struct macroblockd {
322
  struct macroblockd_plane plane[MAX_MB_PLANE];
323

324 325
  struct scale_factors scale_factor[2];
  struct scale_factors scale_factor_uv[2];
John Koleszar's avatar
John Koleszar committed
326

John Koleszar's avatar
John Koleszar committed
327 328 329
  MODE_INFO *prev_mode_info_context;
  MODE_INFO *mode_info_context;
  int mode_info_stride;
John Koleszar's avatar
John Koleszar committed
330

John Koleszar's avatar
John Koleszar committed
331
  FRAME_TYPE frame_type;
John Koleszar's avatar
John Koleszar committed
332

John Koleszar's avatar
John Koleszar committed
333 334
  int up_available;
  int left_available;
335
  int right_available;
John Koleszar's avatar
John Koleszar committed
336

337 338 339 340
  // partition contexts
  PARTITION_CONTEXT *above_seg_context;
  PARTITION_CONTEXT *left_seg_context;

341
  /* 0 (disable) 1 (enable) segmentation */
John Koleszar's avatar
John Koleszar committed
342
  unsigned char segmentation_enabled;
John Koleszar's avatar
John Koleszar committed
343

John Koleszar's avatar
John Koleszar committed
344 345
  /* 0 (do not update) 1 (update) the macroblock segmentation map. */
  unsigned char update_mb_segmentation_map;
John Koleszar's avatar
John Koleszar committed
346

347 348 349 350
#if CONFIG_IMPLICIT_SEGMENTATION
  unsigned char allow_implicit_segment_update;
#endif

John Koleszar's avatar
John Koleszar committed
351 352
  /* 0 (do not update) 1 (update) the macroblock segmentation feature data. */
  unsigned char update_mb_segmentation_data;
John Koleszar's avatar
John Koleszar committed
353

John Koleszar's avatar
John Koleszar committed
354 355
  /* 0 (do not update) 1 (update) the macroblock segmentation feature data. */
  unsigned char mb_segment_abs_delta;
John Koleszar's avatar
John Koleszar committed
356

John Koleszar's avatar
John Koleszar committed
357 358
  /* Per frame flags that define which MB level features (such as quantizer or loop filter level) */
  /* are enabled and when enabled the proabilities used to decode the per MB flags in MB_MODE_INFO */
359

John Koleszar's avatar
John Koleszar committed
360
  // Probability Tree used to code Segment number
361
  vp9_prob mb_segment_tree_probs[MB_SEG_TREE_PROBS];
362

John Koleszar's avatar
John Koleszar committed
363 364 365
  // Segment features
  signed char segment_feature_data[MAX_MB_SEGMENTS][SEG_LVL_MAX];
  unsigned int segment_feature_mask[MAX_MB_SEGMENTS];
John Koleszar's avatar
John Koleszar committed
366

John Koleszar's avatar
John Koleszar committed
367 368 369 370 371
  /* mode_based Loop filter adjustment */
  unsigned char mode_ref_lf_delta_enabled;
  unsigned char mode_ref_lf_delta_update;

  /* Delta values have the range +/- MAX_LOOP_FILTER */
Yaowu Xu's avatar
Yaowu Xu committed
372 373 374 375 376 377 378 379
  /* 0 = Intra, Last, GF, ARF */
  signed char last_ref_lf_deltas[MAX_REF_LF_DELTAS];
  /* 0 = Intra, Last, GF, ARF */
  signed char ref_lf_deltas[MAX_REF_LF_DELTAS];
  /* 0 = I4X4_PRED, ZERO_MV, MV, SPLIT */
  signed char last_mode_lf_deltas[MAX_MODE_LF_DELTAS];
  /* 0 = I4X4_PRED, ZERO_MV, MV, SPLIT */
  signed char mode_lf_deltas[MAX_MODE_LF_DELTAS];
John Koleszar's avatar
John Koleszar committed
380 381 382 383 384 385 386 387 388

  /* Distance of MB away from frame edges */
  int mb_to_left_edge;
  int mb_to_right_edge;
  int mb_to_top_edge;
  int mb_to_bottom_edge;

  unsigned int frames_since_golden;
  unsigned int frames_till_alt_ref_frame;
389

Yaowu Xu's avatar
Yaowu Xu committed
390
  int lossless;
391
  /* Inverse transform function pointers. */
Yaowu Xu's avatar
Yaowu Xu committed
392 393
  void (*inv_txm4x4_1)(int16_t *input, int16_t *output, int pitch);
  void (*inv_txm4x4)(int16_t *input, int16_t *output, int pitch);
394 395 396 397 398
  void (*itxm_add)(int16_t *input, uint8_t *dest, int stride, int eob);
  void (*itxm_add_y_block)(int16_t *q, uint8_t *dst, int stride,
    struct macroblockd *xd);
  void (*itxm_add_uv_block)(int16_t *q, uint8_t *dst, int stride,
    uint16_t *eobs);
399

400
  struct subpix_fn_table  subpix;
401

John Koleszar's avatar
John Koleszar committed
402
  int allow_high_precision_mv;
John Koleszar's avatar
John Koleszar committed
403

John Koleszar's avatar
John Koleszar committed
404
  int corrupted;
405

406 407 408
  int sb_index;   // index of 32x32 block inside the 64x64 block
  int mb_index;   // index of 16x16 block inside the 32x32 block
  int b_index;    // index of 8x8 block inside the 16x16 block
Jingning Han's avatar
Jingning Han committed
409 410
  int q_index;

John Koleszar's avatar
John Koleszar committed
411 412
} MACROBLOCKD;

413 414 415
static INLINE void update_partition_context(MACROBLOCKD *xd,
                                            BLOCK_SIZE_TYPE sb_type,
                                            BLOCK_SIZE_TYPE sb_size) {
416
  int bsl = mi_width_log2(sb_size), bs;
417 418 419
  int bwl = mi_width_log2(sb_type);
  int bhl = mi_height_log2(sb_type);
  int boffset = mi_width_log2(BLOCK_SIZE_SB64X64) - bsl;
420 421 422 423 424
  int i;
  // skip macroblock partition
  if (bsl == 0)
    return;

425 426
  bs = 1 << (bsl - 1);

427 428 429 430 431 432 433 434 435 436 437 438 439 440 441 442 443 444 445 446 447 448 449 450 451 452 453 454 455 456
  // update the partition context at the end notes. set partition bits
  // of block sizes larger than the current one to be one, and partition
  // bits of smaller block sizes to be zero.
  if ((bwl == bsl) && (bhl == bsl)) {
    for (i = 0; i < bs; i++)
      xd->left_seg_context[i] = ~(0xf << boffset);
    for (i = 0; i < bs; i++)
      xd->above_seg_context[i] = ~(0xf << boffset);
  } else if ((bwl == bsl) && (bhl < bsl)) {
    for (i = 0; i < bs; i++)
      xd->left_seg_context[i] = ~(0xe << boffset);
    for (i = 0; i < bs; i++)
      xd->above_seg_context[i] = ~(0xf << boffset);
  }  else if ((bwl < bsl) && (bhl == bsl)) {
    for (i = 0; i < bs; i++)
      xd->left_seg_context[i] = ~(0xf << boffset);
    for (i = 0; i < bs; i++)
      xd->above_seg_context[i] = ~(0xe << boffset);
  } else if ((bwl < bsl) && (bhl < bsl)) {
    for (i = 0; i < bs; i++)
      xd->left_seg_context[i] = ~(0xe << boffset);
    for (i = 0; i < bs; i++)
      xd->above_seg_context[i] = ~(0xe << boffset);
  } else {
    assert(0);
  }
}

static INLINE int partition_plane_context(MACROBLOCKD *xd,
                                          BLOCK_SIZE_TYPE sb_type) {
457
  int bsl = mi_width_log2(sb_type), bs;
458
  int above = 0, left = 0, i;
459
  int boffset = mi_width_log2(BLOCK_SIZE_SB64X64) - bsl;
460

461 462
  bs = 1 << (bsl - 1);

463
  assert(mi_width_log2(sb_type) == mi_height_log2(sb_type));
464 465 466
  assert(bsl >= 0);
  assert(boffset >= 0);

467 468
  bs = 1 << (bsl - 1);

469 470 471 472 473 474 475 476 477 478 479
  for (i = 0; i < bs; i++)
    above |= (xd->above_seg_context[i] & (1 << boffset));
  for (i = 0; i < bs; i++)
    left |= (xd->left_seg_context[i] & (1 << boffset));

  above = (above > 0);
  left  = (left > 0);

  return (left * 2 + above) + (bsl - 1) * PARTITION_PLOFFSET;
}

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 512 513 514 515 516 517 518 519 520 521 522
static BLOCK_SIZE_TYPE get_subsize(BLOCK_SIZE_TYPE bsize,
                                   PARTITION_TYPE partition) {
  BLOCK_SIZE_TYPE subsize;
  switch (partition) {
    case PARTITION_NONE:
      subsize = bsize;
      break;
    case PARTITION_HORZ:
      if (bsize == BLOCK_SIZE_SB64X64)
        subsize = BLOCK_SIZE_SB64X32;
      else if (bsize == BLOCK_SIZE_SB32X32)
        subsize = BLOCK_SIZE_SB32X16;
      else if (bsize == BLOCK_SIZE_MB16X16)
        subsize = BLOCK_SIZE_SB16X8;
      else
        assert(0);
      break;
    case PARTITION_VERT:
      if (bsize == BLOCK_SIZE_SB64X64)
        subsize = BLOCK_SIZE_SB32X64;
      else if (bsize == BLOCK_SIZE_SB32X32)
        subsize = BLOCK_SIZE_SB16X32;
      else if (bsize == BLOCK_SIZE_MB16X16)
        subsize = BLOCK_SIZE_SB8X16;
      else
        assert(0);
      break;
    case PARTITION_SPLIT:
      if (bsize == BLOCK_SIZE_SB64X64)
        subsize = BLOCK_SIZE_SB32X32;
      else if (bsize == BLOCK_SIZE_SB32X32)
        subsize = BLOCK_SIZE_MB16X16;
      else if (bsize == BLOCK_SIZE_MB16X16)
        subsize = BLOCK_SIZE_SB8X8;
      else
        assert(0);
      break;
    default:
      assert(0);
  }
  return subsize;
}

523
#define ACTIVE_HT   110                // quantization stepsize threshold
Deb Mukherjee's avatar
Deb Mukherjee committed
524

525
#define ACTIVE_HT8  300
Deb Mukherjee's avatar
Deb Mukherjee committed
526

527
#define ACTIVE_HT16 300
Deb Mukherjee's avatar
Deb Mukherjee committed
528 529 530 531

// convert MB_PREDICTION_MODE to B_PREDICTION_MODE
static B_PREDICTION_MODE pred_mode_conv(MB_PREDICTION_MODE mode) {
  switch (mode) {
Dmitry Kovalev's avatar
Dmitry Kovalev committed
532
    case DC_PRED: return B_DC_PRED;
Yaowu Xu's avatar
Yaowu Xu committed
533 534
    case V_PRED: return B_V_PRED;
    case H_PRED: return B_H_PRED;
Dmitry Kovalev's avatar
Dmitry Kovalev committed
535
    case TM_PRED: return B_TM_PRED;
Yaowu Xu's avatar
Yaowu Xu committed
536 537 538 539 540 541
    case D45_PRED: return B_D45_PRED;
    case D135_PRED: return B_D135_PRED;
    case D117_PRED: return B_D117_PRED;
    case D153_PRED: return B_D153_PRED;
    case D27_PRED: return B_D27_PRED;
    case D63_PRED: return B_D63_PRED;
Dmitry Kovalev's avatar
Dmitry Kovalev committed
542 543 544
    default:
       assert(0);
       return B_MODE_COUNT;  // Dummy value
Deb Mukherjee's avatar
Deb Mukherjee committed
545 546 547
  }
}

548
// transform mapping
Deb Mukherjee's avatar
Deb Mukherjee committed
549
static TX_TYPE txfm_map(B_PREDICTION_MODE bmode) {
550 551
  switch (bmode) {
    case B_TM_PRED :
Yaowu Xu's avatar
Yaowu Xu committed
552
    case B_D135_PRED :
Dmitry Kovalev's avatar
Dmitry Kovalev committed
553
      return ADST_ADST;
554

Yaowu Xu's avatar
Yaowu Xu committed
555 556
    case B_V_PRED :
    case B_D117_PRED :
Dmitry Kovalev's avatar
Dmitry Kovalev committed
557
      return ADST_DCT;
558

Yaowu Xu's avatar
Yaowu Xu committed
559 560 561
    case B_H_PRED :
    case B_D153_PRED :
    case B_D27_PRED :
Dmitry Kovalev's avatar
Dmitry Kovalev committed
562
      return DCT_ADST;
563

Dmitry Kovalev's avatar
Dmitry Kovalev committed
564 565
    default:
      return DCT_DCT;
566
  }
Deb Mukherjee's avatar
Deb Mukherjee committed
567 568
}

569 570
#define USE_ADST_FOR_I16X16_8X8   1
#define USE_ADST_FOR_I16X16_4X4   1
571 572
#define USE_ADST_FOR_I8X8_4X4     1
#define USE_ADST_PERIPHERY_ONLY   1
573 574
#define USE_ADST_FOR_SB           1
#define USE_ADST_FOR_REMOTE_EDGE  0
575

576
static TX_TYPE get_tx_type_4x4(const MACROBLOCKD *xd, int ib) {
577 578
  // TODO(debargha): explore different patterns for ADST usage when blocksize
  // is smaller than the prediction size
Deb Mukherjee's avatar
Deb Mukherjee committed
579
  TX_TYPE tx_type = DCT_DCT;
580
  const BLOCK_SIZE_TYPE sb_type = xd->mode_info_context->mbmi.sb_type;
581
  const int wb = b_width_log2(sb_type), hb = b_height_log2(sb_type);
582
#if !USE_ADST_FOR_SB
583
  if (sb_type > BLOCK_SIZE_MB16X16)
584
    return tx_type;
585
#endif
586
  if (ib >= (1 << (wb + hb)))  // no chroma adst
587
    return tx_type;
Yaowu Xu's avatar
Yaowu Xu committed
588
  if (xd->lossless)
Yaowu Xu's avatar
Yaowu Xu committed
589
    return DCT_DCT;
Yaowu Xu's avatar
Yaowu Xu committed
590
  if (xd->mode_info_context->mbmi.mode == I4X4_PRED &&
Deb Mukherjee's avatar
Deb Mukherjee committed
591
      xd->q_index < ACTIVE_HT) {
592
    tx_type = txfm_map(
Scott LaVarnway's avatar
Scott LaVarnway committed
593
        xd->mode_info_context->bmi[ib].as_mode.first);
594
  } else if (xd->mode_info_context->mbmi.mode <= TM_PRED &&
595 596 597
             xd->q_index < ACTIVE_HT) {
#if USE_ADST_FOR_I16X16_4X4
#if USE_ADST_PERIPHERY_ONLY
598
    const int hmax = 1 << wb;
599 600 601 602 603 604 605 606 607 608 609 610 611 612 613
    tx_type = txfm_map(pred_mode_conv(xd->mode_info_context->mbmi.mode));
#if USE_ADST_FOR_REMOTE_EDGE
    if ((ib & (hmax - 1)) != 0 && ib >= hmax)
      tx_type = DCT_DCT;
#else
    if (ib >= 1 && ib < hmax) {
      if (tx_type == ADST_ADST) tx_type = ADST_DCT;
      else if (tx_type == DCT_ADST) tx_type = DCT_DCT;
    } else if (ib >= 1 && (ib & (hmax - 1)) == 0) {
      if (tx_type == ADST_ADST) tx_type = DCT_ADST;
      else if (tx_type == ADST_DCT) tx_type = DCT_DCT;
    } else if (ib != 0) {
      tx_type = DCT_DCT;
    }
#endif
614 615 616 617 618 619 620 621
#else
    // Use ADST
    tx_type = txfm_map(pred_mode_conv(xd->mode_info_context->mbmi.mode));
#endif
#else
    // Use 2D DCT
    tx_type = DCT_DCT;
#endif
Deb Mukherjee's avatar
Deb Mukherjee committed
622 623 624 625
  }
  return tx_type;
}

626
static TX_TYPE get_tx_type_8x8(const MACROBLOCKD *xd, int ib) {
627 628
  // TODO(debargha): explore different patterns for ADST usage when blocksize
  // is smaller than the prediction size
Deb Mukherjee's avatar
Deb Mukherjee committed
629
  TX_TYPE tx_type = DCT_DCT;
630
  const BLOCK_SIZE_TYPE sb_type = xd->mode_info_context->mbmi.sb_type;
631
  const int wb = b_width_log2(sb_type), hb = b_height_log2(sb_type);
632
#if !USE_ADST_FOR_SB
633
  if (sb_type > BLOCK_SIZE_MB16X16)
634
    return tx_type;
635
#endif
636
  if (ib >= (1 << (wb + hb)))  // no chroma adst
637
    return tx_type;
638 639
  if (xd->mode_info_context->mbmi.mode <= TM_PRED &&
      xd->q_index < ACTIVE_HT8) {
640
#if USE_ADST_FOR_I16X16_8X8
641
#if USE_ADST_PERIPHERY_ONLY
642
    const int hmax = 1 << wb;
643 644 645 646 647 648 649 650 651 652 653 654 655 656 657
    tx_type = txfm_map(pred_mode_conv(xd->mode_info_context->mbmi.mode));
#if USE_ADST_FOR_REMOTE_EDGE
    if ((ib & (hmax - 1)) != 0 && ib >= hmax)
      tx_type = DCT_DCT;
#else
    if (ib >= 1 && ib < hmax) {
      if (tx_type == ADST_ADST) tx_type = ADST_DCT;
      else if (tx_type == DCT_ADST) tx_type = DCT_DCT;
    } else if (ib >= 1 && (ib & (hmax - 1)) == 0) {
      if (tx_type == ADST_ADST) tx_type = DCT_ADST;
      else if (tx_type == ADST_DCT) tx_type = DCT_DCT;
    } else if (ib != 0) {
      tx_type = DCT_DCT;
    }
#endif
658 659 660 661 662 663 664 665
#else
    // Use ADST
    tx_type = txfm_map(pred_mode_conv(xd->mode_info_context->mbmi.mode));
#endif
#else
    // Use 2D DCT
    tx_type = DCT_DCT;
#endif
Deb Mukherjee's avatar
Deb Mukherjee committed
666 667
  }
  return tx_type;
668
}
669

670
static TX_TYPE get_tx_type_16x16(const MACROBLOCKD *xd, int ib) {
Deb Mukherjee's avatar
Deb Mukherjee committed
671
  TX_TYPE tx_type = DCT_DCT;
672
  const BLOCK_SIZE_TYPE sb_type = xd->mode_info_context->mbmi.sb_type;
673
  const int wb = b_width_log2(sb_type), hb = b_height_log2(sb_type);
674
#if !USE_ADST_FOR_SB
675
  if (sb_type > BLOCK_SIZE_MB16X16)
676 677
    return tx_type;
#endif
678
  if (ib >= (1 << (wb + hb)))
679
    return tx_type;
680
  if (xd->mode_info_context->mbmi.mode <= TM_PRED &&
Deb Mukherjee's avatar
Deb Mukherjee committed
681 682
      xd->q_index < ACTIVE_HT16) {
    tx_type = txfm_map(pred_mode_conv(xd->mode_info_context->mbmi.mode));
683
#if USE_ADST_PERIPHERY_ONLY
684
    if (sb_type > BLOCK_SIZE_MB16X16) {
685
      const int hmax = 1 << wb;
686 687 688 689 690 691 692 693 694 695 696 697 698 699 700 701
#if USE_ADST_FOR_REMOTE_EDGE
      if ((ib & (hmax - 1)) != 0 && ib >= hmax)
        tx_type = DCT_DCT;
#else
      if (ib >= 1 && ib < hmax) {
        if (tx_type == ADST_ADST) tx_type = ADST_DCT;
        else if (tx_type == DCT_ADST) tx_type = DCT_DCT;
      } else if (ib >= 1 && (ib & (hmax - 1)) == 0) {
        if (tx_type == ADST_ADST) tx_type = DCT_ADST;
        else if (tx_type == ADST_DCT) tx_type = DCT_DCT;
      } else if (ib != 0) {
        tx_type = DCT_DCT;
      }
#endif
    }
#endif
Deb Mukherjee's avatar
Deb Mukherjee committed
702 703 704 705
  }
  return tx_type;
}

Dmitry Kovalev's avatar
Dmitry Kovalev committed
706
void vp9_setup_block_dptrs(MACROBLOCKD *xd);
John Koleszar's avatar
John Koleszar committed
707

708
static TX_SIZE get_uv_tx_size(const MACROBLOCKD *xd) {
Dmitry Kovalev's avatar
Dmitry Kovalev committed
709 710 711 712 713 714
  MB_MODE_INFO *mbmi = &xd->mode_info_context->mbmi;
  const TX_SIZE size = mbmi->txfm_size;

  switch (mbmi->sb_type) {
    case BLOCK_SIZE_SB64X64:
      return size;
715 716
    case BLOCK_SIZE_SB64X32:
    case BLOCK_SIZE_SB32X64:
Dmitry Kovalev's avatar
Dmitry Kovalev committed
717 718 719 720 721
    case BLOCK_SIZE_SB32X32:
      if (size == TX_32X32)
        return TX_16X16;
      else
        return size;
722 723 724 725 726 727 728 729 730
    case BLOCK_SIZE_SB32X16:
    case BLOCK_SIZE_SB16X32:
    case BLOCK_SIZE_MB16X16:
      if (size == TX_16X16)
        return TX_8X8;
      else
        return size;
    default:
      return TX_4X4;
731
  }
Dmitry Kovalev's avatar
Dmitry Kovalev committed
732 733

  return size;
734
}
735

736 737 738 739 740 741 742
struct plane_block_idx {
  int plane;
  int block;
};

// TODO(jkoleszar): returning a struct so it can be used in a const context,
// expect to refactor this further later.
John Koleszar's avatar
John Koleszar committed
743
static INLINE struct plane_block_idx plane_block_idx(int y_blocks,
744
                                                     int b_idx) {
John Koleszar's avatar
John Koleszar committed
745
  const int v_offset = y_blocks * 5 / 4;
746 747
  struct plane_block_idx res;

John Koleszar's avatar
John Koleszar committed
748
  if (b_idx < y_blocks) {
749 750 751 752
    res.plane = 0;
    res.block = b_idx;
  } else if (b_idx < v_offset) {
    res.plane = 1;
John Koleszar's avatar
John Koleszar committed
753
    res.block = b_idx - y_blocks;
754
  } else {
John Koleszar's avatar
John Koleszar committed
755
    assert(b_idx < y_blocks * 3 / 2);
756 757 758 759 760 761
    res.plane = 2;
    res.block = b_idx - v_offset;
  }
  return res;
}

John Koleszar's avatar
John Koleszar committed
762
typedef void (*foreach_transformed_block_visitor)(int plane, int block,
763
                                                  BLOCK_SIZE_TYPE bsize,
John Koleszar's avatar
John Koleszar committed
764 765 766
                                                  int ss_txfrm_size,
                                                  void *arg);
static INLINE void foreach_transformed_block_in_plane(
767
    const MACROBLOCKD* const xd, BLOCK_SIZE_TYPE bsize, int plane,
768
    foreach_transformed_block_visitor visit, void *arg) {
769 770
  const int bw = b_width_log2(bsize), bh = b_height_log2(bsize);

John Koleszar's avatar
John Koleszar committed
771 772 773
  // block and transform sizes, in number of 4x4 blocks log 2 ("*_b")
  // 4x4=0, 8x8=2, 16x16=4, 32x32=6, 64x64=8
  const TX_SIZE tx_size = xd->mode_info_context->mbmi.txfm_size;
774
  const int block_size_b = bw + bh;
John Koleszar's avatar
John Koleszar committed
775 776 777 778 779 780 781 782 783 784 785
  const int txfrm_size_b = tx_size * 2;

  // subsampled size of the block
  const int ss_sum = xd->plane[plane].subsampling_x +
                     xd->plane[plane].subsampling_y;
  const int ss_block_size = block_size_b - ss_sum;

  // size of the transform to use. scale the transform down if it's larger
  // than the size of the subsampled data, or forced externally by the mb mode.
  const int ss_max = MAX(xd->plane[plane].subsampling_x,
                         xd->plane[plane].subsampling_y);
786
  const int ss_txfrm_size = txfrm_size_b > ss_block_size
John Koleszar's avatar
John Koleszar committed
787 788
                                ? txfrm_size_b - ss_max * 2
                                : txfrm_size_b;
789
  const int step = 1 << ss_txfrm_size;
John Koleszar's avatar
John Koleszar committed
790 791 792 793 794

  int i;

  assert(txfrm_size_b <= block_size_b);
  assert(ss_txfrm_size <= ss_block_size);
795 796
  for (i = 0; i < (1 << ss_block_size); i += step) {
    visit(plane, i, bsize, ss_txfrm_size, arg);
John Koleszar's avatar
John Koleszar committed
797 798 799 800
  }
}

static INLINE void foreach_transformed_block(
801
    const MACROBLOCKD* const xd, BLOCK_SIZE_TYPE bsize,
John Koleszar's avatar
John Koleszar committed
802 803 804 805
    foreach_transformed_block_visitor visit, void *arg) {
  int plane;

  for (plane = 0; plane < MAX_MB_PLANE; plane++) {
806
    foreach_transformed_block_in_plane(xd, bsize, plane,
John Koleszar's avatar
John Koleszar committed
807 808 809
                                       visit, arg);
  }
}
810

811
static INLINE void foreach_transformed_block_uv(
812
    const MACROBLOCKD* const xd, BLOCK_SIZE_TYPE bsize,
813 814 815 816
    foreach_transformed_block_visitor visit, void *arg) {
  int plane;

  for (plane = 1; plane < MAX_MB_PLANE; plane++) {
817
    foreach_transformed_block_in_plane(xd, bsize, plane,
818 819 820 821
                                       visit, arg);
  }
}

822 823 824 825 826 827 828 829 830 831
// TODO(jkoleszar): In principle, pred_w, pred_h are unnecessary, as we could
// calculate the subsampled BLOCK_SIZE_TYPE, but that type isn't defined for
// sizes smaller than 16x16 yet.
typedef void (*foreach_predicted_block_visitor)(int plane, int block,
                                                BLOCK_SIZE_TYPE bsize,
                                                int pred_w, int pred_h,
                                                void *arg);
static INLINE void foreach_predicted_block_in_plane(
    const MACROBLOCKD* const xd, BLOCK_SIZE_TYPE bsize, int plane,
    foreach_predicted_block_visitor visit, void *arg) {
832 833
  int i, x, y;
  const MB_PREDICTION_MODE mode = xd->mode_info_context->mbmi.mode;
834 835 836 837

  // block sizes in number of 4x4 blocks log 2 ("*_b")
  // 4x4=0, 8x8=2, 16x16=4, 32x32=6, 64x64=8
  // subsampled size of the block
838 839
  const int bw = b_width_log2(bsize) - xd->plane[plane].subsampling_x;
  const int bh = b_height_log2(bsize) - xd->plane[plane].subsampling_y;
840 841

  // size of the predictor to use.
842 843 844
  int pred_w, pred_h;

  if (mode == SPLITMV) {
845 846
    pred_w = 0;
    pred_h = 0;
847 848 849 850 851 852 853 854 855 856 857 858 859 860 861 862
  } else {
    pred_w = bw;
    pred_h = bh;
  }
  assert(pred_w <= bw);
  assert(pred_h <= bh);

  // visit each subblock in raster order
  i = 0;
  for (y = 0; y < 1 << bh; y += 1 << pred_h) {
    for (x = 0; x < 1 << bw; x += 1 << pred_w) {
      visit(plane, i, bsize, pred_w, pred_h, arg);
      i += 1 << pred_w;
    }
    i -= 1 << bw;
    i += 1 << (bw + pred_h);
863 864 865 866 867 868 869 870 871 872 873 874 875 876 877 878 879 880 881 882
  }
}
static INLINE void foreach_predicted_block(
    const MACROBLOCKD* const xd, BLOCK_SIZE_TYPE bsize,
    foreach_predicted_block_visitor visit, void *arg) {
  int plane;

  for (plane = 0; plane < MAX_MB_PLANE; plane++) {
    foreach_predicted_block_in_plane(xd, bsize, plane, visit, arg);
  }
}
static INLINE void foreach_predicted_block_uv(
    const MACROBLOCKD* const xd, BLOCK_SIZE_TYPE bsize,
    foreach_predicted_block_visitor visit, void *arg) {
  int plane;

  for (plane = 1; plane < MAX_MB_PLANE; plane++) {
    foreach_predicted_block_in_plane(xd, bsize, plane, visit, arg);
  }
}
883
static int raster_block_offset(MACROBLOCKD *xd, BLOCK_SIZE_TYPE bsize,
John Koleszar's avatar
John Koleszar committed
884
                               int plane, int block, int stride) {
885 886 887 888 889 890 891
  const int bw = b_width_log2(bsize) - xd->plane[plane].subsampling_x;
  const int y = 4 * (block >> bw), x = 4 * (block & ((1 << bw) - 1));
  return y * stride + x;
}
static int16_t* raster_block_offset_int16(MACROBLOCKD *xd,
                                         BLOCK_SIZE_TYPE bsize,
                                         int plane, int block, int16_t *base) {
John Koleszar's avatar
John Koleszar committed
892 893 894 895 896 897 898 899 900
  const int bw = b_width_log2(bsize) - xd->plane[plane].subsampling_x;
  const int stride = 4 << bw;
  return base + raster_block_offset(xd, bsize, plane, block, stride);
}
static uint8_t* raster_block_offset_uint8(MACROBLOCKD *xd,
                                         BLOCK_SIZE_TYPE bsize,
                                         int plane, int block,
                                         uint8_t *base, int stride) {
  return base + raster_block_offset(xd, bsize, plane, block, stride);
901
}
902

John Koleszar's avatar
John Koleszar committed
903 904 905 906 907 908 909 910 911 912 913 914 915 916 917 918 919 920 921 922 923 924 925 926 927 928 929 930 931 932 933 934 935 936 937 938 939 940 941 942 943 944 945 946 947 948 949 950 951 952 953 954 955 956 957 958 959 960
static int txfrm_block_to_raster_block(MACROBLOCKD *xd,
                                       BLOCK_SIZE_TYPE bsize,
                                       int plane, int block,
                                       int ss_txfrm_size) {
  const int bwl = b_width_log2(bsize) - xd->plane[plane].subsampling_x;
  const int txwl = ss_txfrm_size / 2;
  const int tx_cols_lg2 = bwl - txwl;
  const int tx_cols = 1 << tx_cols_lg2;
  const int raster_mb = block >> ss_txfrm_size;
  const int x = (raster_mb & (tx_cols - 1)) << (txwl);
  const int y = raster_mb >> tx_cols_lg2 << (txwl);
  return x + (y << bwl);
}

static void txfrm_block_to_raster_xy(MACROBLOCKD *xd,
                                     BLOCK_SIZE_TYPE bsize,
                                     int plane, int block,
                                     int ss_txfrm_size,
                                     int *x, int *y) {
  const int bwl = b_width_log2(bsize) - xd->plane[plane].subsampling_x;
  const int txwl = ss_txfrm_size / 2;
  const int tx_cols_lg2 = bwl - txwl;
  const int tx_cols = 1 << tx_cols_lg2;
  const int raster_mb = block >> ss_txfrm_size;
  *x = (raster_mb & (tx_cols - 1)) << (txwl);
  *y = raster_mb >> tx_cols_lg2 << (txwl);
}

static TX_SIZE tx_size_for_plane(MACROBLOCKD *xd, BLOCK_SIZE_TYPE bsize,
                                 int plane) {
  // TODO(jkoleszar): This duplicates a ton of code, but we're going to be
  // moving this to a per-plane lookup shortly, and this will go away then.
  if (!plane) {
    return xd->mode_info_context->mbmi.txfm_size;
  } else {
    const int bw = b_width_log2(bsize), bh = b_height_log2(bsize);
    // block and transform sizes, in number of 4x4 blocks log 2 ("*_b")
    // 4x4=0, 8x8=2, 16x16=4, 32x32=6, 64x64=8
    const TX_SIZE tx_size = xd->mode_info_context->mbmi.txfm_size;
    const int block_size_b = bw + bh;
    const int txfrm_size_b = tx_size * 2;

    // subsampled size of the block
    const int ss_sum = xd->plane[plane].subsampling_x +
                       xd->plane[plane].subsampling_y;
    const int ss_block_size = block_size_b - ss_sum;

    // size of the transform to use. scale the transform down if it's larger
    // than the size of the subsampled data, or forced externally by the mb mode
    const int ss_max = MAX(xd->plane[plane].subsampling_x,
                           xd->plane[plane].subsampling_y);
    const int ss_txfrm_size = txfrm_size_b > ss_block_size
                                  ? txfrm_size_b - ss_max * 2
                                  : txfrm_size_b;
    return (TX_SIZE)(ss_txfrm_size / 2);
  }
}

961
#endif  // VP9_COMMON_VP9_BLOCKD_H_