vp9_blockd.h 28.8 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 {
73 74 75 76 77 78 79 80 81 82 83
  DC_PRED,         // Average of above and left pixels
  V_PRED,          // Vertical
  H_PRED,          // Horizontal
  D45_PRED,        // Directional 45  deg = round(arctan(1/1) * 180/pi)
  D135_PRED,       // Directional 135 deg = 180 - 45
  D117_PRED,       // Directional 117 deg = 180 - 63
  D153_PRED,       // Directional 153 deg = 180 - 27
  D27_PRED,        // Directional 27  deg = round(arctan(1/2) * 180/pi)
  D63_PRED,        // Directional 63  deg = round(arctan(2/1) * 180/pi)
  TM_PRED,         // True-motion
  I4X4_PRED,       // Each 4x4 subblock 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
static INLINE int mi_width_log2(BLOCK_SIZE_TYPE sb_type) {
  int a = b_width_log2(sb_type) - 1;
225 226 227 228 229
#if CONFIG_AB4X4
  // align 4x4 block to mode_info
  if (a < 0)
    a = 0;
#endif
230 231
  assert(a >= 0);
  return a;
232 233
}

234 235
static INLINE int mi_height_log2(BLOCK_SIZE_TYPE sb_type) {
  int a = b_height_log2(sb_type) - 1;
236 237 238 239
#if CONFIG_AB4X4
  if (a < 0)
    a = 0;
#endif
240 241
  assert(a >= 0);
  return a;
242 243
}

John Koleszar's avatar
John Koleszar committed
244 245 246 247
typedef struct {
  MB_PREDICTION_MODE mode, uv_mode;
  MV_REFERENCE_FRAME ref_frame, second_ref_frame;
  TX_SIZE txfm_size;
248
  int_mv mv[2]; // for each reference frame used
Paul Wilkins's avatar
Paul Wilkins committed
249
  int_mv ref_mvs[MAX_REF_FRAMES][MAX_MV_REF_CANDIDATES];
Deb Mukherjee's avatar
Deb Mukherjee committed
250
  int_mv best_mv, best_second_mv;
Paul Wilkins's avatar
Paul Wilkins committed
251

Paul Wilkins's avatar
Paul Wilkins committed
252 253
  int mb_mode_context[MAX_REF_FRAMES];

John Koleszar's avatar
John Koleszar committed
254 255 256
  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;
257
  unsigned char segment_id;           // Segment id for current frame
John Koleszar's avatar
John Koleszar committed
258 259 260 261 262 263 264 265 266

  // 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;
267

268
  INTERPOLATIONFILTERTYPE interp_filter;
269

Ronald S. Bultje's avatar
Ronald S. Bultje committed
270
  BLOCK_SIZE_TYPE sb_type;
John Koleszar's avatar
John Koleszar committed
271 272
} MB_MODE_INFO;

John Koleszar's avatar
John Koleszar committed
273 274
typedef struct {
  MB_MODE_INFO mbmi;
275
  union b_mode_info bmi[16];
John Koleszar's avatar
John Koleszar committed
276 277
} MODE_INFO;

278 279 280 281
struct scale_factors {
  int x_num;
  int x_den;
  int x_offset_q4;
282
  int x_step_q4;
283 284 285
  int y_num;
  int y_den;
  int y_offset_q4;
286
  int y_step_q4;
287 288 289 290 291 292 293 294 295 296 297

  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);

298
  convolve_fn_t predict[2][2][2];  // horiz, vert, avg
299 300
};

301 302
enum { MAX_MB_PLANE = 3 };

303 304 305 306 307
struct buf_2d {
  uint8_t *buf;
  int stride;
};

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

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

Jim Bankoski's avatar
Jim Bankoski committed
325
typedef struct macroblockd {
326
  struct macroblockd_plane plane[MAX_MB_PLANE];
327

328 329
  struct scale_factors scale_factor[2];
  struct scale_factors scale_factor_uv[2];
John Koleszar's avatar
John Koleszar committed
330

John Koleszar's avatar
John Koleszar committed
331 332 333
  MODE_INFO *prev_mode_info_context;
  MODE_INFO *mode_info_context;
  int mode_info_stride;
John Koleszar's avatar
John Koleszar committed
334

John Koleszar's avatar
John Koleszar committed
335
  FRAME_TYPE frame_type;
John Koleszar's avatar
John Koleszar committed
336

John Koleszar's avatar
John Koleszar committed
337 338
  int up_available;
  int left_available;
339
  int right_available;
John Koleszar's avatar
John Koleszar committed
340

341 342 343 344
  // partition contexts
  PARTITION_CONTEXT *above_seg_context;
  PARTITION_CONTEXT *left_seg_context;

345
  /* 0 (disable) 1 (enable) segmentation */
John Koleszar's avatar
John Koleszar committed
346
  unsigned char segmentation_enabled;
John Koleszar's avatar
John Koleszar committed
347

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

351 352 353 354
#if CONFIG_IMPLICIT_SEGMENTATION
  unsigned char allow_implicit_segment_update;
#endif

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

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

John Koleszar's avatar
John Koleszar committed
361 362
  /* 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 */
363

John Koleszar's avatar
John Koleszar committed
364
  // Probability Tree used to code Segment number
365
  vp9_prob mb_segment_tree_probs[MB_SEG_TREE_PROBS];
366

John Koleszar's avatar
John Koleszar committed
367 368 369
  // 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
370

John Koleszar's avatar
John Koleszar committed
371 372 373 374 375
  /* 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
376 377 378 379 380 381 382 383
  /* 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
384 385 386 387 388 389 390 391 392

  /* 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;
393

Yaowu Xu's avatar
Yaowu Xu committed
394
  int lossless;
395
  /* Inverse transform function pointers. */
Yaowu Xu's avatar
Yaowu Xu committed
396 397
  void (*inv_txm4x4_1)(int16_t *input, int16_t *output, int pitch);
  void (*inv_txm4x4)(int16_t *input, int16_t *output, int pitch);
398 399 400 401 402
  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);
403

404
  struct subpix_fn_table  subpix;
405

John Koleszar's avatar
John Koleszar committed
406
  int allow_high_precision_mv;
John Koleszar's avatar
John Koleszar committed
407

John Koleszar's avatar
John Koleszar committed
408
  int corrupted;
409

410 411 412
  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
413 414 415
#if CONFIG_AB4X4
  int ab_index;   // index of 4x4 block inside the 8x8 block
#endif
Jingning Han's avatar
Jingning Han committed
416 417
  int q_index;

John Koleszar's avatar
John Koleszar committed
418 419
} MACROBLOCKD;

420 421 422 423 424 425 426 427 428 429 430 431 432 433 434 435 436 437 438 439 440 441 442 443 444 445
static int *get_sb_index(MACROBLOCKD *xd, BLOCK_SIZE_TYPE subsize) {
  switch (subsize) {
    case BLOCK_SIZE_SB64X32:
    case BLOCK_SIZE_SB32X64:
    case BLOCK_SIZE_SB32X32:
      return &xd->sb_index;
    case BLOCK_SIZE_SB32X16:
    case BLOCK_SIZE_SB16X32:
    case BLOCK_SIZE_MB16X16:
      return &xd->mb_index;
    case BLOCK_SIZE_SB16X8:
    case BLOCK_SIZE_SB8X16:
    case BLOCK_SIZE_SB8X8:
      return &xd->b_index;
#if CONFIG_AB4X4
    case BLOCK_SIZE_SB8X4:
    case BLOCK_SIZE_SB4X8:
    case BLOCK_SIZE_AB4X4:
      return &xd->ab_index;
#endif
    default:
      assert(0);
      return NULL;
  }
}

446 447 448
static INLINE void update_partition_context(MACROBLOCKD *xd,
                                            BLOCK_SIZE_TYPE sb_type,
                                            BLOCK_SIZE_TYPE sb_size) {
449
  int bsl = mi_width_log2(sb_size), bs = 1 << bsl;
450 451 452
  int bwl = mi_width_log2(sb_type);
  int bhl = mi_height_log2(sb_type);
  int boffset = mi_width_log2(BLOCK_SIZE_SB64X64) - bsl;
453
  int i;
454 455 456

#if !CONFIG_AB4X4
  // skip 8x8 block partition
457 458
  if (bsl == 0)
    return;
459
#endif
460 461 462 463 464 465 466 467 468 469 470 471 472 473 474 475 476 477 478 479 480 481 482 483 484 485 486 487 488 489 490

  // 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) {
491
  int bsl = mi_width_log2(sb_type), bs = 1 << bsl;
492
  int above = 0, left = 0, i;
493
  int boffset = mi_width_log2(BLOCK_SIZE_SB64X64) - bsl;
494

495
  assert(mi_width_log2(sb_type) == mi_height_log2(sb_type));
496 497 498 499 500 501 502 503 504 505 506
  assert(bsl >= 0);
  assert(boffset >= 0);

  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);

507 508 509
#if CONFIG_AB4X4
  return (left * 2 + above) + bsl * PARTITION_PLOFFSET;
#else
510
  return (left * 2 + above) + (bsl - 1) * PARTITION_PLOFFSET;
511
#endif
512 513
}

514 515 516 517 518 519 520 521 522 523 524 525 526 527
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;
528 529 530 531
#if CONFIG_AB4X4
      else if (bsize == BLOCK_SIZE_SB8X8)
        subsize = BLOCK_SIZE_SB8X4;
#endif
532 533 534 535 536 537 538 539 540 541
      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;
542 543 544 545
#if CONFIG_AB4X4
      else if (bsize == BLOCK_SIZE_SB8X8)
        subsize = BLOCK_SIZE_SB4X8;
#endif
546 547 548 549 550 551 552 553 554 555
      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;
556 557 558 559
#if CONFIG_AB4X4
      else if (bsize == BLOCK_SIZE_SB8X8)
        subsize = BLOCK_SIZE_AB4X4;
#endif
560 561 562 563 564 565 566 567 568
      else
        assert(0);
      break;
    default:
      assert(0);
  }
  return subsize;
}

569
#define ACTIVE_HT   110                // quantization stepsize threshold
Deb Mukherjee's avatar
Deb Mukherjee committed
570

571
#define ACTIVE_HT8  300
Deb Mukherjee's avatar
Deb Mukherjee committed
572

573
#define ACTIVE_HT16 300
Deb Mukherjee's avatar
Deb Mukherjee committed
574 575 576 577

// 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
578
    case DC_PRED: return B_DC_PRED;
Yaowu Xu's avatar
Yaowu Xu committed
579 580
    case V_PRED: return B_V_PRED;
    case H_PRED: return B_H_PRED;
Dmitry Kovalev's avatar
Dmitry Kovalev committed
581
    case TM_PRED: return B_TM_PRED;
Yaowu Xu's avatar
Yaowu Xu committed
582 583 584 585 586 587
    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
588 589 590
    default:
       assert(0);
       return B_MODE_COUNT;  // Dummy value
Deb Mukherjee's avatar
Deb Mukherjee committed
591 592 593
  }
}

594
// transform mapping
Deb Mukherjee's avatar
Deb Mukherjee committed
595
static TX_TYPE txfm_map(B_PREDICTION_MODE bmode) {
596 597
  switch (bmode) {
    case B_TM_PRED :
Yaowu Xu's avatar
Yaowu Xu committed
598
    case B_D135_PRED :
Dmitry Kovalev's avatar
Dmitry Kovalev committed
599
      return ADST_ADST;
600

Yaowu Xu's avatar
Yaowu Xu committed
601 602
    case B_V_PRED :
    case B_D117_PRED :
Yaowu Xu's avatar
Yaowu Xu committed
603
    case B_D63_PRED:
Dmitry Kovalev's avatar
Dmitry Kovalev committed
604
      return ADST_DCT;
605

Yaowu Xu's avatar
Yaowu Xu committed
606 607 608
    case B_H_PRED :
    case B_D153_PRED :
    case B_D27_PRED :
Dmitry Kovalev's avatar
Dmitry Kovalev committed
609
      return DCT_ADST;
610

Dmitry Kovalev's avatar
Dmitry Kovalev committed
611 612
    default:
      return DCT_DCT;
613
  }
Deb Mukherjee's avatar
Deb Mukherjee committed
614 615
}

616 617
#define USE_ADST_FOR_I16X16_8X8   1
#define USE_ADST_FOR_I16X16_4X4   1
618 619
#define USE_ADST_FOR_I8X8_4X4     1
#define USE_ADST_PERIPHERY_ONLY   1
620 621
#define USE_ADST_FOR_SB           1
#define USE_ADST_FOR_REMOTE_EDGE  0
622

623
static TX_TYPE get_tx_type_4x4(const MACROBLOCKD *xd, int ib) {
624 625
  // TODO(debargha): explore different patterns for ADST usage when blocksize
  // is smaller than the prediction size
Deb Mukherjee's avatar
Deb Mukherjee committed
626
  TX_TYPE tx_type = DCT_DCT;
627
  const BLOCK_SIZE_TYPE sb_type = xd->mode_info_context->mbmi.sb_type;
628
  const int wb = b_width_log2(sb_type);
629
#if !USE_ADST_FOR_SB
630
  if (sb_type > BLOCK_SIZE_MB16X16)
631
    return tx_type;
632
#endif
Yaowu Xu's avatar
Yaowu Xu committed
633
  if (xd->lossless)
Yaowu Xu's avatar
Yaowu Xu committed
634
    return DCT_DCT;
Yaowu Xu's avatar
Yaowu Xu committed
635
  if (xd->mode_info_context->mbmi.mode == I4X4_PRED &&
Deb Mukherjee's avatar
Deb Mukherjee committed
636
      xd->q_index < ACTIVE_HT) {
637
    tx_type = txfm_map(
Scott LaVarnway's avatar
Scott LaVarnway committed
638
        xd->mode_info_context->bmi[ib].as_mode.first);
639
  } else if (xd->mode_info_context->mbmi.mode <= TM_PRED &&
640 641 642
             xd->q_index < ACTIVE_HT) {
#if USE_ADST_FOR_I16X16_4X4
#if USE_ADST_PERIPHERY_ONLY
643
    const int hmax = 1 << wb;
644 645 646 647 648 649 650 651 652 653 654 655 656 657 658
    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
659 660 661 662 663 664 665 666
#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
667 668 669 670
  }
  return tx_type;
}

671
static TX_TYPE get_tx_type_8x8(const MACROBLOCKD *xd, int ib) {
672 673
  // TODO(debargha): explore different patterns for ADST usage when blocksize
  // is smaller than the prediction size
Deb Mukherjee's avatar
Deb Mukherjee committed
674
  TX_TYPE tx_type = DCT_DCT;
675
  const BLOCK_SIZE_TYPE sb_type = xd->mode_info_context->mbmi.sb_type;
676
  const int wb = b_width_log2(sb_type), hb = b_height_log2(sb_type);
677
#if !USE_ADST_FOR_SB
678
  if (sb_type > BLOCK_SIZE_MB16X16)
679
    return tx_type;
680
#endif
681
  if (ib >= (1 << (wb + hb)))  // no chroma adst
682
    return tx_type;
683 684
  if (xd->mode_info_context->mbmi.mode <= TM_PRED &&
      xd->q_index < ACTIVE_HT8) {
685
#if USE_ADST_FOR_I16X16_8X8
686
#if USE_ADST_PERIPHERY_ONLY
687
    const int hmax = 1 << wb;
688 689 690 691 692 693 694 695 696 697 698 699 700 701 702
    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
703 704 705 706 707 708 709 710
#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
711 712
  }
  return tx_type;
713
}
714

715
static TX_TYPE get_tx_type_16x16(const MACROBLOCKD *xd, int ib) {
Deb Mukherjee's avatar
Deb Mukherjee committed
716
  TX_TYPE tx_type = DCT_DCT;
717
  const BLOCK_SIZE_TYPE sb_type = xd->mode_info_context->mbmi.sb_type;
718
  const int wb = b_width_log2(sb_type), hb = b_height_log2(sb_type);
719
#if !USE_ADST_FOR_SB
720
  if (sb_type > BLOCK_SIZE_MB16X16)
721 722
    return tx_type;
#endif
723
  if (ib >= (1 << (wb + hb)))
724
    return tx_type;
725
  if (xd->mode_info_context->mbmi.mode <= TM_PRED &&
Deb Mukherjee's avatar
Deb Mukherjee committed
726 727
      xd->q_index < ACTIVE_HT16) {
    tx_type = txfm_map(pred_mode_conv(xd->mode_info_context->mbmi.mode));
728
#if USE_ADST_PERIPHERY_ONLY
729
    if (sb_type > BLOCK_SIZE_MB16X16) {
730
      const int hmax = 1 << wb;
731 732 733 734 735 736 737 738 739 740 741 742 743 744 745 746
#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
747 748 749 750
  }
  return tx_type;
}

751 752
void vp9_setup_block_dptrs(MACROBLOCKD *xd,
                           int subsampling_x, int subsampling_y);
John Koleszar's avatar
John Koleszar committed
753

754
static TX_SIZE get_uv_tx_size(const MACROBLOCKD *xd) {
Dmitry Kovalev's avatar
Dmitry Kovalev committed
755 756 757 758 759 760
  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;
761 762
    case BLOCK_SIZE_SB64X32:
    case BLOCK_SIZE_SB32X64:
Dmitry Kovalev's avatar
Dmitry Kovalev committed
763 764 765 766 767
    case BLOCK_SIZE_SB32X32:
      if (size == TX_32X32)
        return TX_16X16;
      else
        return size;
768 769 770 771 772 773 774 775 776
    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;
777
  }
Dmitry Kovalev's avatar
Dmitry Kovalev committed
778 779

  return size;
780
}
781

782 783 784 785 786 787 788
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
789
static INLINE struct plane_block_idx plane_block_idx(int y_blocks,
790
                                                     int b_idx) {
John Koleszar's avatar
John Koleszar committed
791
  const int v_offset = y_blocks * 5 / 4;
792 793
  struct plane_block_idx res;

John Koleszar's avatar
John Koleszar committed
794
  if (b_idx < y_blocks) {
795 796 797 798
    res.plane = 0;
    res.block = b_idx;
  } else if (b_idx < v_offset) {
    res.plane = 1;
John Koleszar's avatar
John Koleszar committed
799
    res.block = b_idx - y_blocks;
800
  } else {
John Koleszar's avatar
John Koleszar committed
801
    assert(b_idx < y_blocks * 3 / 2);
802 803 804 805 806 807
    res.plane = 2;
    res.block = b_idx - v_offset;
  }
  return res;
}

John Koleszar's avatar
John Koleszar committed
808
typedef void (*foreach_transformed_block_visitor)(int plane, int block,
809
                                                  BLOCK_SIZE_TYPE bsize,
John Koleszar's avatar
John Koleszar committed
810 811 812
                                                  int ss_txfrm_size,
                                                  void *arg);
static INLINE void foreach_transformed_block_in_plane(
813
    const MACROBLOCKD* const xd, BLOCK_SIZE_TYPE bsize, int plane,
814
    foreach_transformed_block_visitor visit, void *arg) {
815 816
  const int bw = b_width_log2(bsize), bh = b_height_log2(bsize);

John Koleszar's avatar
John Koleszar committed
817 818
  // block and transform sizes, in number of 4x4 blocks log 2 ("*_b")
  // 4x4=0, 8x8=2, 16x16=4, 32x32=6, 64x64=8
John Koleszar's avatar
John Koleszar committed
819 820 821
  // transform size varies per plane, look it up in a common way.
  const TX_SIZE tx_size = plane ? get_uv_tx_size(xd)
                                : xd->mode_info_context->mbmi.txfm_size;
822
  const int block_size_b = bw + bh;
John Koleszar's avatar
John Koleszar committed
823 824 825 826 827 828 829
  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;

John Koleszar's avatar
John Koleszar committed
830
  const int step = 1 << txfrm_size_b;
John Koleszar's avatar
John Koleszar committed
831 832 833 834

  int i;

  assert(txfrm_size_b <= block_size_b);
John Koleszar's avatar
John Koleszar committed
835
  assert(txfrm_size_b <= ss_block_size);
836
  for (i = 0; i < (1 << ss_block_size); i += step) {
John Koleszar's avatar
John Koleszar committed
837
    visit(plane, i, bsize, txfrm_size_b, arg);
John Koleszar's avatar
John Koleszar committed
838 839 840 841
  }
}

static INLINE void foreach_transformed_block(
842
    const MACROBLOCKD* const xd, BLOCK_SIZE_TYPE bsize,
John Koleszar's avatar
John Koleszar committed
843 844 845 846
    foreach_transformed_block_visitor visit, void *arg) {
  int plane;

  for (plane = 0; plane < MAX_MB_PLANE; plane++) {
847
    foreach_transformed_block_in_plane(xd, bsize, plane,
John Koleszar's avatar
John Koleszar committed
848 849 850
                                       visit, arg);
  }
}
851

852
static INLINE void foreach_transformed_block_uv(
853
    const MACROBLOCKD* const xd, BLOCK_SIZE_TYPE bsize,
854 855 856 857
    foreach_transformed_block_visitor visit, void *arg) {
  int plane;

  for (plane = 1; plane < MAX_MB_PLANE; plane++) {
858
    foreach_transformed_block_in_plane(xd, bsize, plane,
859 860 861 862
                                       visit, arg);
  }
}

863 864 865 866 867 868 869 870 871 872
// 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) {
873 874
  int i, x, y;
  const MB_PREDICTION_MODE mode = xd->mode_info_context->mbmi.mode;
875 876 877 878

  // 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
879 880
  const int bw = b_width_log2(bsize) - xd->plane[plane].subsampling_x;
  const int bh = b_height_log2(bsize) - xd->plane[plane].subsampling_y;
881 882

  // size of the predictor to use.
883 884 885
  int pred_w, pred_h;

  if (mode == SPLITMV) {
886 887
    pred_w = 0;
    pred_h = 0;
888 889 890 891 892 893 894 895 896 897 898 899 900 901 902 903
  } 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);
904 905 906 907 908 909 910 911 912 913 914 915 916 917 918 919 920 921 922 923
  }
}
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);
  }
}
924
static int raster_block_offset(MACROBLOCKD *xd, BLOCK_SIZE_TYPE bsize,
John Koleszar's avatar
John Koleszar committed
925
                               int plane, int block, int stride) {
926 927 928 929 930 931 932
  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
933 934 935 936 937 938 939 940 941
  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);
942
}
943

John Koleszar's avatar
John Koleszar committed
944 945 946 947 948 949 950 951 952 953 954 955 956 957 958 959 960 961 962 963 964 965 966 967 968 969 970 971
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);
}

972
#endif  // VP9_COMMON_VP9_BLOCKD_H_