decodeframe.c 192 KB
Newer Older
Jingning Han's avatar
Jingning Han committed
1
/*
Yaowu Xu's avatar
Yaowu Xu committed
2
 * Copyright (c) 2016, Alliance for Open Media. All rights reserved
Jingning Han's avatar
Jingning Han committed
3
 *
Yaowu Xu's avatar
Yaowu Xu committed
4 5 6 7 8 9
 * This source code is subject to the terms of the BSD 2 Clause License and
 * the Alliance for Open Media Patent License 1.0. If the BSD 2 Clause License
 * was not distributed with this source code in the LICENSE file, you can
 * obtain it at www.aomedia.org/license/software. If the Alliance for Open
 * Media Patent License 1.0 was not distributed with this source code in the
 * PATENTS file, you can obtain it at www.aomedia.org/license/patent.
Jingning Han's avatar
Jingning Han committed
10 11 12 13 14
 */

#include <assert.h>
#include <stdlib.h>  // qsort()

Yaowu Xu's avatar
Yaowu Xu committed
15 16 17
#include "./aom_config.h"
#include "./aom_dsp_rtcd.h"
#include "./aom_scale_rtcd.h"
Jingning Han's avatar
Jingning Han committed
18
#include "./av1_rtcd.h"
Jingning Han's avatar
Jingning Han committed
19

20
#include "aom/aom_codec.h"
Yaowu Xu's avatar
Yaowu Xu committed
21
#include "aom_dsp/aom_dsp_common.h"
Jingning Han's avatar
Jingning Han committed
22 23
#include "aom_dsp/bitreader.h"
#include "aom_dsp/bitreader_buffer.h"
24
#include "aom_dsp/binary_codes_reader.h"
Yaowu Xu's avatar
Yaowu Xu committed
25
#include "aom_mem/aom_mem.h"
26 27
#include "aom_ports/mem.h"
#include "aom_ports/mem_ops.h"
Yaowu Xu's avatar
Yaowu Xu committed
28 29
#include "aom_scale/aom_scale.h"
#include "aom_util/aom_thread.h"
Jingning Han's avatar
Jingning Han committed
30

31 32 33 34
#if CONFIG_BITSTREAM_DEBUG
#include "aom_util/debug_util.h"
#endif  // CONFIG_BITSTREAM_DEBUG

35
#include "av1/common/alloccommon.h"
36
#if CONFIG_CDEF
37
#include "av1/common/cdef.h"
38
#include "av1/common/clpf.h"
39
#endif
40 41 42
#if CONFIG_INSPECTION
#include "av1/decoder/inspection.h"
#endif
43 44 45
#include "av1/common/common.h"
#include "av1/common/entropy.h"
#include "av1/common/entropymode.h"
46
#include "av1/common/entropymv.h"
47 48 49 50
#include "av1/common/idct.h"
#include "av1/common/pred_common.h"
#include "av1/common/quant_common.h"
#include "av1/common/reconinter.h"
Jingning Han's avatar
Jingning Han committed
51
#include "av1/common/reconintra.h"
52
#include "av1/common/seg_common.h"
Jingning Han's avatar
Jingning Han committed
53
#include "av1/common/thread_common.h"
54
#include "av1/common/tile_common.h"
55

56 57 58
#include "av1/decoder/decodeframe.h"
#include "av1/decoder/decodemv.h"
#include "av1/decoder/decoder.h"
Angie Chiang's avatar
Angie Chiang committed
59 60 61
#if CONFIG_LV_MAP
#include "av1/decoder/decodetxb.h"
#endif
Jingning Han's avatar
Jingning Han committed
62
#include "av1/decoder/detokenize.h"
63
#include "av1/decoder/dsubexp.h"
Jingning Han's avatar
Jingning Han committed
64

65
#if CONFIG_WARPED_MOTION || CONFIG_GLOBAL_MOTION
Yue Chen's avatar
Yue Chen committed
66
#include "av1/common/warped_motion.h"
67
#endif  // CONFIG_WARPED_MOTION || CONFIG_GLOBAL_MOTION
Yue Chen's avatar
Yue Chen committed
68

Yaowu Xu's avatar
Yaowu Xu committed
69
#define MAX_AV1_HEADER_SIZE 80
Michael Bebenita's avatar
Michael Bebenita committed
70
#define ACCT_STR __func__
Jingning Han's avatar
Jingning Han committed
71

72 73
#if CONFIG_PVQ
#include "av1/common/partition.h"
hui su's avatar
hui su committed
74 75
#include "av1/common/pvq.h"
#include "av1/common/scan.h"
76
#include "av1/decoder/decint.h"
hui su's avatar
hui su committed
77 78
#include "av1/decoder/pvq_decoder.h"
#include "av1/encoder/encodemb.h"
79 80 81
#include "av1/encoder/hybrid_fwd_txfm.h"
#endif

82 83 84 85
#if CONFIG_CFL
#include "av1/common/cfl.h"
#endif

86 87 88 89 90 91 92 93
static struct aom_read_bit_buffer *init_read_bit_buffer(
    AV1Decoder *pbi, struct aom_read_bit_buffer *rb, const uint8_t *data,
    const uint8_t *data_end, uint8_t clear_data[MAX_AV1_HEADER_SIZE]);
static int read_compressed_header(AV1Decoder *pbi, const uint8_t *data,
                                  size_t partition_size);
static size_t read_uncompressed_header(AV1Decoder *pbi,
                                       struct aom_read_bit_buffer *rb);

Yaowu Xu's avatar
Yaowu Xu committed
94
static int is_compound_reference_allowed(const AV1_COMMON *cm) {
95 96 97
#if CONFIG_LOWDELAY_COMPOUND  // Normative in decoder
  return !frame_is_intra_only(cm);
#else
Jingning Han's avatar
Jingning Han committed
98
  int i;
99
  if (frame_is_intra_only(cm)) return 0;
100
  for (i = 1; i < INTER_REFS_PER_FRAME; ++i)
101
    if (cm->ref_frame_sign_bias[i + 1] != cm->ref_frame_sign_bias[1]) return 1;
Jingning Han's avatar
Jingning Han committed
102 103

  return 0;
104
#endif
Jingning Han's avatar
Jingning Han committed
105 106
}

Yaowu Xu's avatar
Yaowu Xu committed
107
static void setup_compound_reference_mode(AV1_COMMON *cm) {
108
#if CONFIG_EXT_REFS
109
  cm->comp_fwd_ref[0] = LAST_FRAME;
110 111 112 113
  cm->comp_fwd_ref[1] = LAST2_FRAME;
  cm->comp_fwd_ref[2] = LAST3_FRAME;
  cm->comp_fwd_ref[3] = GOLDEN_FRAME;

114 115
  cm->comp_bwd_ref[0] = BWDREF_FRAME;
  cm->comp_bwd_ref[1] = ALTREF_FRAME;
116
#else
Jingning Han's avatar
Jingning Han committed
117
  if (cm->ref_frame_sign_bias[LAST_FRAME] ==
118
      cm->ref_frame_sign_bias[GOLDEN_FRAME]) {
Jingning Han's avatar
Jingning Han committed
119 120 121 122
    cm->comp_fixed_ref = ALTREF_FRAME;
    cm->comp_var_ref[0] = LAST_FRAME;
    cm->comp_var_ref[1] = GOLDEN_FRAME;
  } else if (cm->ref_frame_sign_bias[LAST_FRAME] ==
123
             cm->ref_frame_sign_bias[ALTREF_FRAME]) {
Jingning Han's avatar
Jingning Han committed
124 125 126 127 128 129 130 131
    cm->comp_fixed_ref = GOLDEN_FRAME;
    cm->comp_var_ref[0] = LAST_FRAME;
    cm->comp_var_ref[1] = ALTREF_FRAME;
  } else {
    cm->comp_fixed_ref = LAST_FRAME;
    cm->comp_var_ref[0] = GOLDEN_FRAME;
    cm->comp_var_ref[1] = ALTREF_FRAME;
  }
132
#endif  // CONFIG_EXT_REFS
Jingning Han's avatar
Jingning Han committed
133 134 135 136 137 138
}

static int read_is_valid(const uint8_t *start, size_t len, const uint8_t *end) {
  return len != 0 && len <= (size_t)(end - start);
}

Yaowu Xu's avatar
Yaowu Xu committed
139 140
static int decode_unsigned_max(struct aom_read_bit_buffer *rb, int max) {
  const int data = aom_rb_read_literal(rb, get_unsigned_bits(max));
Jingning Han's avatar
Jingning Han committed
141 142 143
  return data > max ? max : data;
}

144 145 146
static TX_MODE read_tx_mode(AV1_COMMON *cm, MACROBLOCKD *xd,
                            struct aom_read_bit_buffer *rb) {
  int i, all_lossless = 1;
147
#if CONFIG_TX64X64
148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164
  TX_MODE tx_mode;
#endif

  if (cm->seg.enabled) {
    for (i = 0; i < MAX_SEGMENTS; ++i) {
      if (!xd->lossless[i]) {
        all_lossless = 0;
        break;
      }
    }
  } else {
    all_lossless = xd->lossless[0];
  }

  if (all_lossless) return ONLY_4X4;
#if CONFIG_TX64X64
  tx_mode = aom_rb_read_bit(rb) ? TX_MODE_SELECT : aom_rb_read_literal(rb, 2);
165 166 167
  if (tx_mode == ALLOW_32X32) tx_mode += aom_rb_read_bit(rb);
  return tx_mode;
#else
Yaowu Xu's avatar
Yaowu Xu committed
168
  return aom_rb_read_bit(rb) ? TX_MODE_SELECT : aom_rb_read_literal(rb, 2);
169
#endif  // CONFIG_TX64X64
170
}
Jingning Han's avatar
Jingning Han committed
171

172
#if !CONFIG_EC_ADAPT
Yaowu Xu's avatar
Yaowu Xu committed
173 174 175 176 177 178 179
static void read_tx_size_probs(FRAME_CONTEXT *fc, aom_reader *r) {
  int i, j, k;
  for (i = 0; i < MAX_TX_DEPTH; ++i)
    for (j = 0; j < TX_SIZE_CONTEXTS; ++j)
      for (k = 0; k < i + 1; ++k)
        av1_diff_update_prob(r, &fc->tx_size_probs[i][j][k], ACCT_STR);
}
180
#endif
Yaowu Xu's avatar
Yaowu Xu committed
181

182
#if !CONFIG_EC_ADAPT
Yaowu Xu's avatar
Yaowu Xu committed
183
static void read_switchable_interp_probs(FRAME_CONTEXT *fc, aom_reader *r) {
Jingning Han's avatar
Jingning Han committed
184
  int i, j;
185
  for (j = 0; j < SWITCHABLE_FILTER_CONTEXTS; ++j) {
Jingning Han's avatar
Jingning Han committed
186
    for (i = 0; i < SWITCHABLE_FILTERS - 1; ++i)
Michael Bebenita's avatar
Michael Bebenita committed
187
      av1_diff_update_prob(r, &fc->switchable_interp_prob[j][i], ACCT_STR);
188
  }
Jingning Han's avatar
Jingning Han committed
189
}
190
#endif
Jingning Han's avatar
Jingning Han committed
191

Yaowu Xu's avatar
Yaowu Xu committed
192
static void read_inter_mode_probs(FRAME_CONTEXT *fc, aom_reader *r) {
193
#if CONFIG_REF_MV
194
  int i;
195
  for (i = 0; i < NEWMV_MODE_CONTEXTS; ++i)
Michael Bebenita's avatar
Michael Bebenita committed
196
    av1_diff_update_prob(r, &fc->newmv_prob[i], ACCT_STR);
197
  for (i = 0; i < ZEROMV_MODE_CONTEXTS; ++i)
Michael Bebenita's avatar
Michael Bebenita committed
198
    av1_diff_update_prob(r, &fc->zeromv_prob[i], ACCT_STR);
199
  for (i = 0; i < REFMV_MODE_CONTEXTS; ++i)
Michael Bebenita's avatar
Michael Bebenita committed
200
    av1_diff_update_prob(r, &fc->refmv_prob[i], ACCT_STR);
201
  for (i = 0; i < DRL_MODE_CONTEXTS; ++i)
Michael Bebenita's avatar
Michael Bebenita committed
202
    av1_diff_update_prob(r, &fc->drl_prob[i], ACCT_STR);
203
#else
204
#if !CONFIG_EC_ADAPT
205
  int i, j;
206
  for (i = 0; i < INTER_MODE_CONTEXTS; ++i) {
Jingning Han's avatar
Jingning Han committed
207
    for (j = 0; j < INTER_MODES - 1; ++j)
Michael Bebenita's avatar
Michael Bebenita committed
208
      av1_diff_update_prob(r, &fc->inter_mode_probs[i][j], ACCT_STR);
209
  }
210 211 212
#else
  (void)fc;
  (void)r;
213
#endif
214
#endif
Jingning Han's avatar
Jingning Han committed
215 216
}

217
#if CONFIG_EXT_INTER
Yaowu Xu's avatar
Yaowu Xu committed
218
static void read_inter_compound_mode_probs(FRAME_CONTEXT *fc, aom_reader *r) {
219
  int i, j;
Michael Bebenita's avatar
Michael Bebenita committed
220
  if (aom_read(r, GROUP_DIFF_UPDATE_PROB, ACCT_STR)) {
221 222
    for (j = 0; j < INTER_MODE_CONTEXTS; ++j) {
      for (i = 0; i < INTER_COMPOUND_MODES - 1; ++i) {
Michael Bebenita's avatar
Michael Bebenita committed
223
        av1_diff_update_prob(r, &fc->inter_compound_mode_probs[j][i], ACCT_STR);
224 225 226 227 228
      }
    }
  }
}
#endif  // CONFIG_EXT_INTER
229
#if !CONFIG_EC_ADAPT
230 231 232 233 234 235 236 237 238 239 240 241 242 243 244 245 246 247 248 249
#if !CONFIG_EXT_TX
static void read_ext_tx_probs(FRAME_CONTEXT *fc, aom_reader *r) {
  int i, j, k;
  if (aom_read(r, GROUP_DIFF_UPDATE_PROB, ACCT_STR)) {
    for (i = TX_4X4; i < EXT_TX_SIZES; ++i) {
      for (j = 0; j < TX_TYPES; ++j) {
        for (k = 0; k < TX_TYPES - 1; ++k)
          av1_diff_update_prob(r, &fc->intra_ext_tx_prob[i][j][k], ACCT_STR);
      }
    }
  }
  if (aom_read(r, GROUP_DIFF_UPDATE_PROB, ACCT_STR)) {
    for (i = TX_4X4; i < EXT_TX_SIZES; ++i) {
      for (k = 0; k < TX_TYPES - 1; ++k)
        av1_diff_update_prob(r, &fc->inter_ext_tx_prob[i][k], ACCT_STR);
    }
  }
}
#endif
#endif
250

251
static REFERENCE_MODE read_frame_reference_mode(
Yaowu Xu's avatar
Yaowu Xu committed
252
    const AV1_COMMON *cm, struct aom_read_bit_buffer *rb) {
253
  if (is_compound_reference_allowed(cm)) {
Zoe Liu's avatar
Zoe Liu committed
254 255 256
#if CONFIG_REF_ADAPT
    return aom_rb_read_bit(rb) ? REFERENCE_MODE_SELECT : SINGLE_REFERENCE;
#else
Yaowu Xu's avatar
Yaowu Xu committed
257
    return aom_rb_read_bit(rb)
258
               ? REFERENCE_MODE_SELECT
Yaowu Xu's avatar
Yaowu Xu committed
259
               : (aom_rb_read_bit(rb) ? COMPOUND_REFERENCE : SINGLE_REFERENCE);
Zoe Liu's avatar
Zoe Liu committed
260
#endif  // CONFIG_REF_ADAPT
261 262 263 264
  } else {
    return SINGLE_REFERENCE;
  }
}
Jingning Han's avatar
Jingning Han committed
265

Yaowu Xu's avatar
Yaowu Xu committed
266
static void read_frame_reference_mode_probs(AV1_COMMON *cm, aom_reader *r) {
Jingning Han's avatar
Jingning Han committed
267
  FRAME_CONTEXT *const fc = cm->fc;
268
  int i, j;
Jingning Han's avatar
Jingning Han committed
269 270 271

  if (cm->reference_mode == REFERENCE_MODE_SELECT)
    for (i = 0; i < COMP_INTER_CONTEXTS; ++i)
Michael Bebenita's avatar
Michael Bebenita committed
272
      av1_diff_update_prob(r, &fc->comp_inter_prob[i], ACCT_STR);
Jingning Han's avatar
Jingning Han committed
273

274
  if (cm->reference_mode != COMPOUND_REFERENCE) {
Jingning Han's avatar
Jingning Han committed
275
    for (i = 0; i < REF_CONTEXTS; ++i) {
276
      for (j = 0; j < (SINGLE_REFS - 1); ++j) {
Michael Bebenita's avatar
Michael Bebenita committed
277
        av1_diff_update_prob(r, &fc->single_ref_prob[i][j], ACCT_STR);
278
      }
Jingning Han's avatar
Jingning Han committed
279
    }
280
  }
Jingning Han's avatar
Jingning Han committed
281

282 283
  if (cm->reference_mode != SINGLE_REFERENCE) {
    for (i = 0; i < REF_CONTEXTS; ++i) {
284
#if CONFIG_EXT_REFS
285
      for (j = 0; j < (FWD_REFS - 1); ++j)
Michael Bebenita's avatar
Michael Bebenita committed
286
        av1_diff_update_prob(r, &fc->comp_ref_prob[i][j], ACCT_STR);
287
      for (j = 0; j < (BWD_REFS - 1); ++j)
Michael Bebenita's avatar
Michael Bebenita committed
288
        av1_diff_update_prob(r, &fc->comp_bwdref_prob[i][j], ACCT_STR);
289 290
#else
      for (j = 0; j < (COMP_REFS - 1); ++j)
Michael Bebenita's avatar
Michael Bebenita committed
291
        av1_diff_update_prob(r, &fc->comp_ref_prob[i][j], ACCT_STR);
292
#endif  // CONFIG_EXT_REFS
293 294
    }
  }
Jingning Han's avatar
Jingning Han committed
295 296
}

Yaowu Xu's avatar
Yaowu Xu committed
297
static void update_mv_probs(aom_prob *p, int n, aom_reader *r) {
Jingning Han's avatar
Jingning Han committed
298
  int i;
Michael Bebenita's avatar
Michael Bebenita committed
299
  for (i = 0; i < n; ++i) av1_diff_update_prob(r, &p[i], ACCT_STR);
Jingning Han's avatar
Jingning Han committed
300 301
}

Yaowu Xu's avatar
Yaowu Xu committed
302
static void read_mv_probs(nmv_context *ctx, int allow_hp, aom_reader *r) {
303
  int i;
Jingning Han's avatar
Jingning Han committed
304

Yaowu Xu's avatar
Yaowu Xu committed
305
#if !CONFIG_EC_ADAPT
306
  int j;
Jingning Han's avatar
Jingning Han committed
307 308 309 310 311 312 313 314 315 316 317
  update_mv_probs(ctx->joints, MV_JOINTS - 1, r);

  for (i = 0; i < 2; ++i) {
    nmv_component *const comp_ctx = &ctx->comps[i];
    update_mv_probs(&comp_ctx->sign, 1, r);
    update_mv_probs(comp_ctx->classes, MV_CLASSES - 1, r);
    update_mv_probs(comp_ctx->class0, CLASS0_SIZE - 1, r);
    update_mv_probs(comp_ctx->bits, MV_OFFSET_BITS, r);
  }
  for (i = 0; i < 2; ++i) {
    nmv_component *const comp_ctx = &ctx->comps[i];
318
    for (j = 0; j < CLASS0_SIZE; ++j) {
Jingning Han's avatar
Jingning Han committed
319
      update_mv_probs(comp_ctx->class0_fp[j], MV_FP_SIZE - 1, r);
320
    }
321
    update_mv_probs(comp_ctx->fp, MV_FP_SIZE - 1, r);
Jingning Han's avatar
Jingning Han committed
322
  }
323
#endif  // !CONFIG_EC_ADAPT
Jingning Han's avatar
Jingning Han committed
324 325 326 327 328 329 330 331 332 333

  if (allow_hp) {
    for (i = 0; i < 2; ++i) {
      nmv_component *const comp_ctx = &ctx->comps[i];
      update_mv_probs(&comp_ctx->class0_hp, 1, r);
      update_mv_probs(&comp_ctx->hp, 1, r);
    }
  }
}

334
static void inverse_transform_block(MACROBLOCKD *xd, int plane,
335
                                    const TX_TYPE tx_type,
336
                                    const TX_SIZE tx_size, uint8_t *dst,
337
                                    int stride, int16_t scan_line, int eob) {
Jingning Han's avatar
Jingning Han committed
338
  struct macroblockd_plane *const pd = &xd->plane[plane];
339
  tran_low_t *const dqcoeff = pd->dqcoeff;
Angie Chiang's avatar
Angie Chiang committed
340
  av1_inverse_transform_block(xd, dqcoeff, tx_type, tx_size, dst, stride, eob);
341
  memset(dqcoeff, 0, (scan_line + 1) * sizeof(dqcoeff[0]));
Jingning Han's avatar
Jingning Han committed
342 343
}

344
#if CONFIG_PVQ
Thomas Daede's avatar
Thomas Daede committed
345
static int av1_pvq_decode_helper(MACROBLOCKD *xd, tran_low_t *ref_coeff,
346
                                 tran_low_t *dqcoeff, int16_t *quant, int pli,
347
                                 int bs, TX_TYPE tx_type, int xdec,
ltrudeau's avatar
ltrudeau committed
348
                                 PVQ_SKIP_TYPE ac_dc_coded) {
349 350 351 352
  unsigned int flags;  // used for daala's stream analyzer.
  int off;
  const int is_keyframe = 0;
  const int has_dc_skip = 1;
353
  int coeff_shift = 3 - av1_get_tx_scale(bs);
Thomas Daede's avatar
Thomas Daede committed
354
  int hbd_downshift = 0;
355
  int rounding_mask;
356 357 358 359 360 361
  // DC quantizer for PVQ
  int pvq_dc_quant;
  int lossless = (quant[0] == 0);
  const int blk_size = tx_size_wide[bs];
  int eob = 0;
  int i;
Thomas Daede's avatar
Thomas Daede committed
362
  od_dec_ctx *dec = &xd->daala_dec;
363
  int use_activity_masking = dec->use_activity_masking;
364 365
  DECLARE_ALIGNED(16, tran_low_t, dqcoeff_pvq[OD_TXSIZE_MAX * OD_TXSIZE_MAX]);
  DECLARE_ALIGNED(16, tran_low_t, ref_coeff_pvq[OD_TXSIZE_MAX * OD_TXSIZE_MAX]);
366

367 368
  od_coeff ref_int32[OD_TXSIZE_MAX * OD_TXSIZE_MAX];
  od_coeff out_int32[OD_TXSIZE_MAX * OD_TXSIZE_MAX];
369

Thomas Daede's avatar
Thomas Daede committed
370 371
  hbd_downshift = xd->bd - 8;

372 373 374
  od_raster_to_coding_order(ref_coeff_pvq, blk_size, tx_type, ref_coeff,
                            blk_size);

Thomas Daede's avatar
Thomas Daede committed
375
  assert(OD_COEFF_SHIFT >= 4);
376 377 378
  if (lossless)
    pvq_dc_quant = 1;
  else {
379 380
    if (use_activity_masking)
      pvq_dc_quant = OD_MAXI(
Thomas Daede's avatar
Thomas Daede committed
381
          1, (quant[0] << (OD_COEFF_SHIFT - 3) >> hbd_downshift) *
382 383 384
                     dec->state.pvq_qm_q4[pli][od_qm_get_index(bs, 0)] >>
                 4);
    else
Thomas Daede's avatar
Thomas Daede committed
385 386
      pvq_dc_quant =
          OD_MAXI(1, quant[0] << (OD_COEFF_SHIFT - 3) >> hbd_downshift);
387 388 389 390 391
  }

  off = od_qm_offset(bs, xdec);

  // copy int16 inputs to int32
392
  for (i = 0; i < blk_size * blk_size; i++) {
393
    ref_int32[i] =
Thomas Daede's avatar
Thomas Daede committed
394 395
        AOM_SIGNED_SHL(ref_coeff_pvq[i], OD_COEFF_SHIFT - coeff_shift) >>
        hbd_downshift;
396
  }
397

Thomas Daede's avatar
Thomas Daede committed
398 399
  od_pvq_decode(dec, ref_int32, out_int32,
                OD_MAXI(1, quant[1] << (OD_COEFF_SHIFT - 3) >> hbd_downshift),
400
                pli, bs, OD_PVQ_BETA[use_activity_masking][pli][bs],
401 402
                is_keyframe, &flags, ac_dc_coded, dec->state.qm + off,
                dec->state.qm_inv + off);
403

404 405
  if (!has_dc_skip || out_int32[0]) {
    out_int32[0] =
406
        has_dc_skip + generic_decode(dec->r, &dec->state.adapt->model_dc[pli],
407 408
                                     &dec->state.adapt->ex_dc[pli][bs][0], 2,
                                     "dc:mag");
409 410 411 412 413 414 415 416
    if (out_int32[0]) out_int32[0] *= aom_read_bit(dec->r, "dc:sign") ? -1 : 1;
  }
  out_int32[0] = out_int32[0] * pvq_dc_quant + ref_int32[0];

  // copy int32 result back to int16
  assert(OD_COEFF_SHIFT > coeff_shift);
  rounding_mask = (1 << (OD_COEFF_SHIFT - coeff_shift - 1)) - 1;
  for (i = 0; i < blk_size * blk_size; i++) {
Thomas Daede's avatar
Thomas Daede committed
417
    out_int32[i] = AOM_SIGNED_SHL(out_int32[i], hbd_downshift);
418 419
    dqcoeff_pvq[i] = (out_int32[i] + (out_int32[i] < 0) + rounding_mask) >>
                     (OD_COEFF_SHIFT - coeff_shift);
420 421 422 423 424 425 426 427 428
  }

  od_coding_order_to_raster(dqcoeff, blk_size, tx_type, dqcoeff_pvq, blk_size);

  eob = blk_size * blk_size;

  return eob;
}

ltrudeau's avatar
ltrudeau committed
429 430 431 432 433 434
static PVQ_SKIP_TYPE read_pvq_skip(AV1_COMMON *cm, MACROBLOCKD *const xd,
                                   int plane, TX_SIZE tx_size) {
  // decode ac/dc coded flag. bit0: DC coded, bit1 : AC coded
  // NOTE : we don't use 5 symbols for luma here in aom codebase,
  // since block partition is taken care of by aom.
  // So, only AC/DC skip info is coded
435
  const int ac_dc_coded = aom_read_symbol(
ltrudeau's avatar
ltrudeau committed
436
      xd->daala_dec.r,
437
      xd->daala_dec.state.adapt->skip_cdf[2 * tx_size + (plane != 0)], 4,
438
      "skip");
ltrudeau's avatar
ltrudeau committed
439 440 441 442 443 444 445 446
  if (ac_dc_coded < 0 || ac_dc_coded > 3) {
    aom_internal_error(&cm->error, AOM_CODEC_INVALID_PARAM,
                       "Invalid PVQ Skip Type");
  }
  return ac_dc_coded;
}

static int av1_pvq_decode_helper2(AV1_COMMON *cm, MACROBLOCKD *const xd,
447 448
                                  MB_MODE_INFO *const mbmi, int plane, int row,
                                  int col, TX_SIZE tx_size, TX_TYPE tx_type) {
449 450 451 452 453 454 455 456 457 458
  struct macroblockd_plane *const pd = &xd->plane[plane];
  // transform block size in pixels
  int tx_blk_size = tx_size_wide[tx_size];
  int i, j;
  tran_low_t *pvq_ref_coeff = pd->pvq_ref_coeff;
  const int diff_stride = tx_blk_size;
  int16_t *pred = pd->pred;
  tran_low_t *const dqcoeff = pd->dqcoeff;
  uint8_t *dst;
  int eob;
ltrudeau's avatar
ltrudeau committed
459
  const PVQ_SKIP_TYPE ac_dc_coded = read_pvq_skip(cm, xd, plane, tx_size);
460 461 462 463 464 465 466 467 468

  eob = 0;
  dst = &pd->dst.buf[4 * row * pd->dst.stride + 4 * col];

  if (ac_dc_coded) {
    int xdec = pd->subsampling_x;
    int seg_id = mbmi->segment_id;
    int16_t *quant;
    FWD_TXFM_PARAM fwd_txfm_param;
469 470
    // ToDo(yaowu): correct this with optimal number from decoding process.
    const int max_scan_line = tx_size_2d[tx_size];
471
#if CONFIG_HIGHBITDEPTH
Thomas Daede's avatar
Thomas Daede committed
472 473 474 475 476 477 478 479 480 481
    if (xd->cur_buf->flags & YV12_FLAG_HIGHBITDEPTH) {
      for (j = 0; j < tx_blk_size; j++)
        for (i = 0; i < tx_blk_size; i++)
          pred[diff_stride * j + i] =
              CONVERT_TO_SHORTPTR(dst)[pd->dst.stride * j + i];
    } else {
#endif
      for (j = 0; j < tx_blk_size; j++)
        for (i = 0; i < tx_blk_size; i++)
          pred[diff_stride * j + i] = dst[pd->dst.stride * j + i];
482
#if CONFIG_HIGHBITDEPTH
Thomas Daede's avatar
Thomas Daede committed
483 484
    }
#endif
485 486 487 488 489

    fwd_txfm_param.tx_type = tx_type;
    fwd_txfm_param.tx_size = tx_size;
    fwd_txfm_param.lossless = xd->lossless[seg_id];

490
#if CONFIG_HIGHBITDEPTH
Thomas Daede's avatar
Thomas Daede committed
491 492
    if (xd->cur_buf->flags & YV12_FLAG_HIGHBITDEPTH) {
      fwd_txfm_param.bd = xd->bd;
hui su's avatar
hui su committed
493
      av1_highbd_fwd_txfm(pred, pvq_ref_coeff, diff_stride, &fwd_txfm_param);
Thomas Daede's avatar
Thomas Daede committed
494
    } else {
495
#endif  // CONFIG_HIGHBITDEPTH
hui su's avatar
hui su committed
496
      av1_fwd_txfm(pred, pvq_ref_coeff, diff_stride, &fwd_txfm_param);
497
#if CONFIG_HIGHBITDEPTH
Thomas Daede's avatar
Thomas Daede committed
498
    }
499
#endif  // CONFIG_HIGHBITDEPTH
500 501 502

    quant = &pd->seg_dequant[seg_id][0];  // aom's quantizer

Thomas Daede's avatar
Thomas Daede committed
503 504
    eob = av1_pvq_decode_helper(xd, pvq_ref_coeff, dqcoeff, quant, plane,
                                tx_size, tx_type, xdec, ac_dc_coded);
505

506 507
    inverse_transform_block(xd, plane, tx_type, tx_size, dst, pd->dst.stride,
                            max_scan_line, eob);
508 509 510 511 512 513
  }

  return eob;
}
#endif

514 515 516 517 518 519 520 521 522 523 524 525 526 527 528 529 530 531 532 533
static int get_block_idx(const MACROBLOCKD *xd, int plane, int row, int col) {
  const int bsize = xd->mi[0]->mbmi.sb_type;
  const struct macroblockd_plane *pd = &xd->plane[plane];
#if CONFIG_CB4X4
#if CONFIG_CHROMA_2X2
  const BLOCK_SIZE plane_bsize = get_plane_block_size(bsize, pd);
#else
  const BLOCK_SIZE plane_bsize =
      AOMMAX(BLOCK_4X4, get_plane_block_size(bsize, pd));
#endif  // CONFIG_CHROMA_2X2
#else
  const BLOCK_SIZE plane_bsize =
      get_plane_block_size(AOMMAX(BLOCK_8X8, bsize), pd);
#endif
  const int max_blocks_wide = max_block_wide(xd, plane_bsize, plane);
  const TX_SIZE tx_size = get_tx_size(plane, xd);
  const uint8_t txh_unit = tx_size_high_unit[tx_size];
  return row * max_blocks_wide + col * txh_unit;
}

534 535 536
static void predict_and_reconstruct_intra_block(
    AV1_COMMON *cm, MACROBLOCKD *const xd, aom_reader *const r,
    MB_MODE_INFO *const mbmi, int plane, int row, int col, TX_SIZE tx_size) {
537
  PLANE_TYPE plane_type = get_plane_type(plane);
538
  const int block_idx = get_block_idx(xd, plane, row, col);
539 540 541
#if CONFIG_PVQ
  (void)r;
#endif
542
  av1_predict_intra_block_facade(xd, plane, block_idx, col, row, tx_size);
Jingning Han's avatar
Jingning Han committed
543 544

  if (!mbmi->skip) {
545
#if !CONFIG_PVQ
546
    struct macroblockd_plane *const pd = &xd->plane[plane];
Angie Chiang's avatar
Angie Chiang committed
547 548
#if CONFIG_LV_MAP
    int16_t max_scan_line = 0;
549 550 551
    int eob;
    av1_read_coeffs_txb_facade(cm, xd, r, row, col, block_idx, plane,
                               pd->dqcoeff, &max_scan_line, &eob);
552 553
    // tx_type will be read out in av1_read_coeffs_txb_facade
    TX_TYPE tx_type = get_tx_type(plane_type, xd, block_idx, tx_size);
Angie Chiang's avatar
Angie Chiang committed
554
#else   // CONFIG_LV_MAP
555
    TX_TYPE tx_type = get_tx_type(plane_type, xd, block_idx, tx_size);
Angie Chiang's avatar
Angie Chiang committed
556
    const SCAN_ORDER *scan_order = get_scan(cm, tx_size, tx_type, 0);
557 558
    int16_t max_scan_line = 0;
    const int eob =
559
        av1_decode_block_tokens(cm, xd, plane, scan_order, col, row, tx_size,
560
                                tx_type, &max_scan_line, r, mbmi->segment_id);
Angie Chiang's avatar
Angie Chiang committed
561
#endif  // CONFIG_LV_MAP
562 563 564
    if (eob) {
      uint8_t *dst =
          &pd->dst.buf[(row * pd->dst.stride + col) << tx_size_wide_log2[0]];
565 566
      inverse_transform_block(xd, plane, tx_type, tx_size, dst, pd->dst.stride,
                              max_scan_line, eob);
567
    }
568
#else
569
    TX_TYPE tx_type = get_tx_type(plane_type, xd, block_idx, tx_size);
ltrudeau's avatar
ltrudeau committed
570
    av1_pvq_decode_helper2(cm, xd, mbmi, plane, row, col, tx_size, tx_type);
571
#endif
Jingning Han's avatar
Jingning Han committed
572
  }
573 574 575 576 577 578 579 580
#if CONFIG_CFL
  if (plane == AOM_PLANE_Y) {
    struct macroblockd_plane *const pd = &xd->plane[plane];
    uint8_t *dst =
        &pd->dst.buf[(row * pd->dst.stride + col) << tx_size_wide_log2[0]];
    cfl_store(xd->cfl, dst, pd->dst.stride, row, col, tx_size);
  }
#endif
Jingning Han's avatar
Jingning Han committed
581 582
}

583
#if CONFIG_VAR_TX && !CONFIG_COEF_INTERLEAVE
Angie Chiang's avatar
Angie Chiang committed
584 585
static void decode_reconstruct_tx(AV1_COMMON *cm, MACROBLOCKD *const xd,
                                  aom_reader *r, MB_MODE_INFO *const mbmi,
586
                                  int plane, BLOCK_SIZE plane_bsize,
587 588
                                  int blk_row, int blk_col, TX_SIZE tx_size,
                                  int *eob_total) {
589
  const struct macroblockd_plane *const pd = &xd->plane[plane];
590
  const BLOCK_SIZE bsize = txsize_to_bsize[tx_size];
591 592
  const int tx_row = blk_row >> (1 - pd->subsampling_y);
  const int tx_col = blk_col >> (1 - pd->subsampling_x);
593
  const TX_SIZE plane_tx_size =
594
      plane ? uv_txsize_lookup[bsize][mbmi->inter_tx_size[tx_row][tx_col]][0][0]
595
            : mbmi->inter_tx_size[tx_row][tx_col];
596
  // Scale to match transform block unit.
597 598
  const int max_blocks_high = max_block_high(xd, plane_bsize, plane);
  const int max_blocks_wide = max_block_wide(xd, plane_bsize, plane);
599

600
  if (blk_row >= max_blocks_high || blk_col >= max_blocks_wide) return;
601

602
  if (tx_size == plane_tx_size) {
603
    PLANE_TYPE plane_type = get_plane_type(plane);
604
    int block_idx = get_block_idx(xd, plane, blk_row, blk_col);
Angie Chiang's avatar
Angie Chiang committed
605 606 607
#if CONFIG_LV_MAP
    (void)segment_id;
    int16_t max_scan_line = 0;
608 609 610
    int eob;
    av1_read_coeffs_txb_facade(cm, xd, r, row, col, block_idx, plane,
                               pd->dqcoeff, &max_scan_line, &eob);
611 612
    // tx_type will be read out in av1_read_coeffs_txb_facade
    TX_TYPE tx_type = get_tx_type(plane_type, xd, block_idx, plane_tx_size);
Angie Chiang's avatar
Angie Chiang committed
613
#else   // CONFIG_LV_MAP
614
    TX_TYPE tx_type = get_tx_type(plane_type, xd, block_idx, plane_tx_size);
Angie Chiang's avatar
Angie Chiang committed
615
    const SCAN_ORDER *sc = get_scan(cm, plane_tx_size, tx_type, 1);
616
    int16_t max_scan_line = 0;
617 618 619
    const int eob = av1_decode_block_tokens(
        cm, xd, plane, sc, blk_col, blk_row, plane_tx_size, tx_type,
        &max_scan_line, r, mbmi->segment_id);
Angie Chiang's avatar
Angie Chiang committed
620
#endif  // CONFIG_LV_MAP
621 622 623 624
    inverse_transform_block(xd, plane, tx_type, plane_tx_size,
                            &pd->dst.buf[(blk_row * pd->dst.stride + blk_col)
                                         << tx_size_wide_log2[0]],
                            pd->dst.stride, max_scan_line, eob);
625 626
    *eob_total += eob;
  } else {
627 628
    const TX_SIZE sub_txs = sub_tx_size_map[tx_size];
    const int bsl = tx_size_wide_unit[sub_txs];
629 630 631 632 633
    int i;

    assert(bsl > 0);

    for (i = 0; i < 4; ++i) {
634 635
      const int offsetr = blk_row + (i >> 1) * bsl;
      const int offsetc = blk_col + (i & 0x01) * bsl;
636

637
      if (offsetr >= max_blocks_high || offsetc >= max_blocks_wide) continue;
638

639
      decode_reconstruct_tx(cm, xd, r, mbmi, plane, plane_bsize, offsetr,
640
                            offsetc, sub_txs, eob_total);
641 642 643
    }
  }
}
644
#endif  // CONFIG_VAR_TX
645

646
#if !CONFIG_VAR_TX || CONFIG_SUPERTX || CONFIG_COEF_INTERLEAVE || \
647
    (!CONFIG_VAR_TX && CONFIG_EXT_TX && CONFIG_RECT_TX)
Angie Chiang's avatar
Angie Chiang committed
648
static int reconstruct_inter_block(AV1_COMMON *cm, MACROBLOCKD *const xd,
649 650
                                   aom_reader *const r, int segment_id,
                                   int plane, int row, int col,
651
                                   TX_SIZE tx_size) {
652
  PLANE_TYPE plane_type = get_plane_type(plane);
653
  int block_idx = get_block_idx(xd, plane, row, col);
654 655 656
#if CONFIG_PVQ
  int eob;
  (void)r;
657 658 659
  (void)segment_id;
#else
  struct macroblockd_plane *const pd = &xd->plane[plane];
660 661 662
#endif

#if !CONFIG_PVQ
Angie Chiang's avatar
Angie Chiang committed
663 664
#if CONFIG_LV_MAP
  (void)segment_id;
665
  int16_t max_scan_line = 0;
666 667 668
  int eob;
  av1_read_coeffs_txb_facade(cm, xd, r, row, col, block_idx, plane, pd->dqcoeff,
                             &max_scan_line, &eob);
669 670
  // tx_type will be read out in av1_read_coeffs_txb_facade
  TX_TYPE tx_type = get_tx_type(plane_type, xd, block_idx, tx_size);
Angie Chiang's avatar
Angie Chiang committed
671 672
#else   // CONFIG_LV_MAP
  int16_t max_scan_line = 0;
673
  TX_TYPE tx_type = get_tx_type(plane_type, xd, block_idx, tx_size);
Angie Chiang's avatar
Angie Chiang committed
674
  const SCAN_ORDER *scan_order = get_scan(cm, tx_size, tx_type, 1);
675
  const int eob =
676 677
      av1_decode_block_tokens(cm, xd, plane, scan_order, col, row, tx_size,
                              tx_type, &max_scan_line, r, segment_id);
Angie Chiang's avatar
Angie Chiang committed
678
#endif  // CONFIG_LV_MAP
679 680
  uint8_t *dst =
      &pd->dst.buf[(row * pd->dst.stride + col) << tx_size_wide_log2[0]];
681
  if (eob)
682 683
    inverse_transform_block(xd, plane, tx_type, tx_size, dst, pd->dst.stride,
                            max_scan_line, eob);
684
#else
685
  TX_TYPE tx_type = get_tx_type(plane_type, xd, block_idx, tx_size);
ltrudeau's avatar
ltrudeau committed
686 687
  eob = av1_pvq_decode_helper2(cm, xd, &xd->mi[0]->mbmi, plane, row, col,
                               tx_size, tx_type);
688
#endif
Jingning Han's avatar
Jingning Han committed
689 690
  return eob;
}
691
#endif  // !CONFIG_VAR_TX || CONFIG_SUPER_TX
Jingning Han's avatar
Jingning Han committed
692

693 694 695
static void set_offsets(AV1_COMMON *const cm, MACROBLOCKD *const xd,
                        BLOCK_SIZE bsize, int mi_row, int mi_col, int bw,
                        int bh, int x_mis, int y_mis) {
Jingning Han's avatar
Jingning Han committed
696 697 698 699 700 701 702 703 704
  const int offset = mi_row * cm->mi_stride + mi_col;
  int x, y;
  const TileInfo *const tile = &xd->tile;

  xd->mi = cm->mi_grid_visible + offset;
  xd->mi[0] = &cm->mi[offset];
  // TODO(slavarnway): Generate sb_type based on bwl and bhl, instead of
  // passing bsize from decode_partition().
  xd->mi[0]->mbmi.sb_type = bsize;
705 706 707 708
#if CONFIG_RD_DEBUG
  xd->mi[0]->mbmi.mi_row = mi_row;
  xd->mi[0]->mbmi.mi_col = mi_col;
#endif
Jingning Han's avatar
Jingning Han committed
709
  for (y = 0; y < y_mis; ++y)
Jingning Han's avatar
Jingning Han committed
710
    for (x = !y; x < x_mis; ++x) xd->mi[y * cm->mi_stride + x] = xd->mi[0];
Jingning Han's avatar
Jingning Han committed
711

Jingning Han's avatar
Jingning Han committed
712
  set_plane_n4(xd, bw, bh);
Jingning Han's avatar
Jingning Han committed
713 714
  set_skip_context(xd, mi_row, mi_col);

715 716 717 718
#if CONFIG_VAR_TX
  xd->max_tx_size = max_txsize_lookup[bsize];
#endif

719 720 721
  // Distance of Mb to the various image edges. These are specified to 8th pel
  // as they are always compared to values that are in 1/8th pel units
  set_mi_row_col(xd, tile, mi_row, bh, mi_col, bw,
722
#if CONFIG_DEPENDENT_HORZTILES
723 724 725
                 cm->dependent_horz_tiles,
#endif  // CONFIG_DEPENDENT_HORZTILES
                 cm->mi_rows, cm->mi_cols);
Jingning Han's avatar
Jingning Han committed
726

727 728
  av1_setup_dst_planes(xd->plane, bsize, get_frame_new_buffer(cm), mi_row,
                       mi_col);
Jingning Han's avatar
Jingning Han committed
729 730
}

731
#if CONFIG_SUPERTX
Yaowu Xu's avatar
Yaowu Xu committed
732
static MB_MODE_INFO *set_offsets_extend(AV1_COMMON *const cm,
733 734
                                        MACROBLOCKD *const xd,
                                        const TileInfo *const tile,
735 736 737
                                        BLOCK_SIZE bsize_pred, int mi_row_pred,
                                        int mi_col_pred, int mi_row_ori,
                                        int mi_col_ori) {
738 739 740
  // Used in supertx
  // (mi_row_ori, mi_col_ori): location for mv
  // (mi_row_pred, mi_col_pred, bsize_pred): region to predict
Jingning Han's avatar
Jingning Han committed
741 742
  const int bw = mi_size_wide[bsize_pred];
  const int bh = mi_size_high[bsize_pred];
743 744 745
  const int offset = mi_row_ori * cm->mi_stride + mi_col_ori;
  xd->mi = cm->mi_grid_visible + offset;
  xd->mi[0] = cm->mi + offset;
746
  set_mi_row_col(xd, tile, mi_row_pred, bh, mi_col_pred, bw,
747
#if CONFIG_DEPENDENT_HORZTILES
748 749 750
                 cm->dependent_horz_tiles,
#endif  // CONFIG_DEPENDENT_HORZTILES
                 cm->mi_rows, cm->mi_cols);
751

752 753
  xd->up_available = (mi_row_ori > tile->mi_row_start);
  xd->left_available = (mi_col_ori > tile->mi_col_start);
754

Jingning Han's avatar
Jingning Han committed
755
  set_plane_n4(xd, bw, bh);
756 757 758 759

  return &xd->mi[0]->mbmi;
}

Angie Chiang's avatar
Angie Chiang committed
760
#if CONFIG_SUPERTX
Yaowu Xu's avatar
Yaowu Xu committed
761 762 763
static MB_MODE_INFO *set_mb_offsets(AV1_COMMON *const cm, MACROBLOCKD *const xd,
                                    BLOCK_SIZE bsize, int mi_row, int mi_col,
                                    int bw, int bh, int x_mis, int y_mis) {
764 765 766 767 768 769 770 771
  const int offset = mi_row * cm->mi_stride + mi_col;
  const TileInfo *const tile = &xd->tile;
  int x, y;

  xd->mi = cm->mi_grid_visible + offset;
  xd->mi[0] = cm->mi + offset;
  xd->mi[0]->mbmi.sb_type = bsize;
  for (y = 0; y < y_mis; ++y)
772
    for (x = !y; x < x_mis; ++x) xd->mi[y * cm->mi_stride + x] = xd->mi[0];
773

774
  set_mi_row_col(xd, tile, mi_row, bh, mi_col, bw,
775
#if CONFIG_DEPENDENT_HORZTILES
776 777 778
                 cm->dependent_horz_tiles,
#endif  // CONFIG_DEPENDENT_HORZTILES
                 cm->mi_rows, cm->mi_cols);
779 780
  return &xd->mi[0]->mbmi;
}
Angie Chiang's avatar
Angie Chiang committed
781
#endif
782

Yaowu Xu's avatar
Yaowu Xu committed
783
static void set_offsets_topblock(AV1_COMMON *const cm, MACROBLOCKD *const xd,
784 785
                                 const TileInfo *const tile, BLOCK_SIZE bsize,
                                 int mi_row, int mi_col) {
Jingning Han's avatar
Jingning Han committed
786 787
  const int bw = mi_size_wide[bsize];
  const int bh = mi_size_high[bsize];
788 789 790 791 792
  const int offset = mi_row * cm->mi_stride + mi_col;

  xd->mi = cm->mi_grid_visible + offset;
  xd->mi[0] = cm->mi + offset;

Jingning Han's avatar
Jingning Han committed
793
  set_plane_n4(xd, bw, bh);
794

795
  set_mi_row_col(xd, tile, mi_row, bh, mi_col, bw,
796
#if CONFIG_DEPENDENT_HORZTILES
797 798 799
                 cm->dependent_horz_tiles,
#endif  // CONFIG_DEPENDENT_HORZTILES
                 cm->mi_rows, cm->mi_cols);
Debargha Mukherjee's avatar