Coverage Report

Created: 2022-08-24 06:17

/src/aom/av1/decoder/decodetxb.c
Line
Count
Source (jump to first uncovered line)
1
/*
2
 * Copyright (c) 2017, Alliance for Open Media. All rights reserved
3
 *
4
 * This source code is subject to the terms of the BSD 2 Clause License and
5
 * the Alliance for Open Media Patent License 1.0. If the BSD 2 Clause License
6
 * was not distributed with this source code in the LICENSE file, you can
7
 * obtain it at www.aomedia.org/license/software. If the Alliance for Open
8
 * Media Patent License 1.0 was not distributed with this source code in the
9
 * PATENTS file, you can obtain it at www.aomedia.org/license/patent.
10
 */
11
12
#include "av1/decoder/decodetxb.h"
13
14
#include "aom_ports/mem.h"
15
#include "av1/common/idct.h"
16
#include "av1/common/scan.h"
17
#include "av1/common/txb_common.h"
18
#include "av1/decoder/decodemv.h"
19
20
#define ACCT_STR __func__
21
22
1.25M
static int read_golomb(MACROBLOCKD *xd, aom_reader *r) {
23
1.25M
  int x = 1;
24
1.25M
  int length = 0;
25
1.25M
  int i = 0;
26
27
3.75M
  while (!i) {
28
2.49M
    i = aom_read_bit(r, ACCT_STR);
29
2.49M
    ++length;
30
2.49M
    if (length > 20) {
31
20
      aom_internal_error(xd->error_info, AOM_CODEC_CORRUPT_FRAME,
32
20
                         "Invalid length in read_golomb");
33
20
      break;
34
20
    }
35
2.49M
  }
36
37
2.49M
  for (i = 0; i < length - 1; ++i) {
38
1.23M
    x <<= 1;
39
1.23M
    x += aom_read_bit(r, ACCT_STR);
40
1.23M
  }
41
42
1.25M
  return x - 1;
43
1.25M
}
44
45
771k
static INLINE int rec_eob_pos(const int eob_token, const int extra) {
46
771k
  int eob = av1_eob_group_start[eob_token];
47
771k
  if (eob > 2) {
48
503k
    eob += extra;
49
503k
  }
50
771k
  return eob;
51
771k
}
52
53
static INLINE int get_dqv(const int16_t *dequant, int coeff_idx,
54
7.16M
                          const qm_val_t *iqmatrix) {
55
7.16M
  int dqv = dequant[!!coeff_idx];
56
7.16M
  if (iqmatrix != NULL)
57
868k
    dqv =
58
868k
        ((iqmatrix[coeff_idx] * dqv) + (1 << (AOM_QM_BITS - 1))) >> AOM_QM_BITS;
59
7.16M
  return dqv;
60
7.16M
}
61
62
static INLINE void read_coeffs_reverse_2d(aom_reader *r, TX_SIZE tx_size,
63
                                          int start_si, int end_si,
64
                                          const int16_t *scan, int bwl,
65
                                          uint8_t *levels,
66
                                          base_cdf_arr base_cdf,
67
495k
                                          br_cdf_arr br_cdf) {
68
15.2M
  for (int c = end_si; c >= start_si; --c) {
69
14.7M
    const int pos = scan[c];
70
14.7M
    const int coeff_ctx = get_lower_levels_ctx_2d(levels, pos, bwl, tx_size);
71
14.7M
    const int nsymbs = 4;
72
14.7M
    int level = aom_read_symbol(r, base_cdf[coeff_ctx], nsymbs, ACCT_STR);
73
14.7M
    if (level > NUM_BASE_LEVELS) {
74
2.32M
      const int br_ctx = get_br_ctx_2d(levels, pos, bwl);
75
2.32M
      aom_cdf_prob *cdf = br_cdf[br_ctx];
76
8.01M
      for (int idx = 0; idx < COEFF_BASE_RANGE; idx += BR_CDF_SIZE - 1) {
77
6.76M
        const int k = aom_read_symbol(r, cdf, BR_CDF_SIZE, ACCT_STR);
78
6.76M
        level += k;
79
6.76M
        if (k < BR_CDF_SIZE - 1) break;
80
6.76M
      }
81
2.32M
    }
82
14.7M
    levels[get_padded_idx(pos, bwl)] = level;
83
14.7M
  }
84
495k
}
85
86
static INLINE void read_coeffs_reverse(aom_reader *r, TX_SIZE tx_size,
87
                                       TX_CLASS tx_class, int start_si,
88
                                       int end_si, const int16_t *scan, int bwl,
89
                                       uint8_t *levels, base_cdf_arr base_cdf,
90
540k
                                       br_cdf_arr br_cdf) {
91
1.95M
  for (int c = end_si; c >= start_si; --c) {
92
1.41M
    const int pos = scan[c];
93
1.41M
    const int coeff_ctx =
94
1.41M
        get_lower_levels_ctx(levels, pos, bwl, tx_size, tx_class);
95
1.41M
    const int nsymbs = 4;
96
1.41M
    int level = aom_read_symbol(r, base_cdf[coeff_ctx], nsymbs, ACCT_STR);
97
1.41M
    if (level > NUM_BASE_LEVELS) {
98
196k
      const int br_ctx = get_br_ctx(levels, pos, bwl, tx_class);
99
196k
      aom_cdf_prob *cdf = br_cdf[br_ctx];
100
300k
      for (int idx = 0; idx < COEFF_BASE_RANGE; idx += BR_CDF_SIZE - 1) {
101
291k
        const int k = aom_read_symbol(r, cdf, BR_CDF_SIZE, ACCT_STR);
102
291k
        level += k;
103
291k
        if (k < BR_CDF_SIZE - 1) break;
104
291k
      }
105
196k
    }
106
1.41M
    levels[get_padded_idx(pos, bwl)] = level;
107
1.41M
  }
108
540k
}
109
110
uint8_t av1_read_coeffs_txb(const AV1_COMMON *const cm, DecoderCodingBlock *dcb,
111
                            aom_reader *const r, const int blk_row,
112
                            const int blk_col, const int plane,
113
                            const TXB_CTX *const txb_ctx,
114
1.70M
                            const TX_SIZE tx_size) {
115
1.70M
  MACROBLOCKD *const xd = &dcb->xd;
116
1.70M
  FRAME_CONTEXT *const ec_ctx = xd->tile_ctx;
117
1.70M
  const int32_t max_value = (1 << (7 + xd->bd)) - 1;
118
1.70M
  const int32_t min_value = -(1 << (7 + xd->bd));
119
1.70M
  const TX_SIZE txs_ctx = get_txsize_entropy_ctx(tx_size);
120
1.70M
  const PLANE_TYPE plane_type = get_plane_type(plane);
121
1.70M
  MB_MODE_INFO *const mbmi = xd->mi[0];
122
1.70M
  struct macroblockd_plane *const pd = &xd->plane[plane];
123
1.70M
  const int16_t *const dequant = pd->seg_dequant_QTX[mbmi->segment_id];
124
1.70M
  tran_low_t *const tcoeffs = dcb->dqcoeff_block[plane] + dcb->cb_offset[plane];
125
1.70M
  const int shift = av1_get_tx_scale(tx_size);
126
1.70M
  const int bwl = get_txb_bwl(tx_size);
127
1.70M
  const int width = get_txb_wide(tx_size);
128
1.70M
  const int height = get_txb_high(tx_size);
129
1.70M
  int cul_level = 0;
130
1.70M
  int dc_val = 0;
131
1.70M
  uint8_t levels_buf[TX_PAD_2D];
132
1.70M
  uint8_t *const levels = set_levels(levels_buf, width);
133
1.70M
  const int all_zero = aom_read_symbol(
134
1.70M
      r, ec_ctx->txb_skip_cdf[txs_ctx][txb_ctx->txb_skip_ctx], 2, ACCT_STR);
135
1.70M
  eob_info *eob_data = dcb->eob_data[plane] + dcb->txb_offset[plane];
136
1.70M
  uint16_t *const eob = &(eob_data->eob);
137
1.70M
  uint16_t *const max_scan_line = &(eob_data->max_scan_line);
138
1.70M
  *max_scan_line = 0;
139
1.70M
  *eob = 0;
140
141
#if CONFIG_INSPECTION
142
  if (plane == 0) {
143
    const int txk_type_idx =
144
        av1_get_txk_type_index(mbmi->bsize, blk_row, blk_col);
145
    mbmi->tx_skip[txk_type_idx] = all_zero;
146
  }
147
#endif
148
149
1.70M
  if (all_zero) {
150
932k
    *max_scan_line = 0;
151
932k
    if (plane == 0) {
152
178k
      xd->tx_type_map[blk_row * xd->tx_type_map_stride + blk_col] = DCT_DCT;
153
178k
    }
154
932k
    return 0;
155
932k
  }
156
157
771k
  if (plane == AOM_PLANE_Y) {
158
    // only y plane's tx_type is transmitted
159
585k
    av1_read_tx_type(cm, xd, blk_row, blk_col, tx_size, r);
160
585k
  }
161
771k
  const TX_TYPE tx_type =
162
771k
      av1_get_tx_type(xd, plane_type, blk_row, blk_col, tx_size,
163
771k
                      cm->features.reduced_tx_set_used);
164
771k
  const TX_CLASS tx_class = tx_type_to_class[tx_type];
165
771k
  const qm_val_t *iqmatrix =
166
771k
      av1_get_iqmatrix(&cm->quant_params, xd, plane, tx_size, tx_type);
167
771k
  const SCAN_ORDER *const scan_order = get_scan(tx_size, tx_type);
168
771k
  const int16_t *const scan = scan_order->scan;
169
771k
  int eob_extra = 0;
170
771k
  int eob_pt = 1;
171
172
771k
  const int eob_multi_size = txsize_log2_minus4[tx_size];
173
771k
  const int eob_multi_ctx = (tx_class == TX_CLASS_2D) ? 0 : 1;
174
771k
  switch (eob_multi_size) {
175
146k
    case 0:
176
146k
      eob_pt =
177
146k
          aom_read_symbol(r, ec_ctx->eob_flag_cdf16[plane_type][eob_multi_ctx],
178
146k
                          5, ACCT_STR) +
179
146k
          1;
180
146k
      break;
181
77.1k
    case 1:
182
77.1k
      eob_pt =
183
77.1k
          aom_read_symbol(r, ec_ctx->eob_flag_cdf32[plane_type][eob_multi_ctx],
184
77.1k
                          6, ACCT_STR) +
185
77.1k
          1;
186
77.1k
      break;
187
156k
    case 2:
188
156k
      eob_pt =
189
156k
          aom_read_symbol(r, ec_ctx->eob_flag_cdf64[plane_type][eob_multi_ctx],
190
156k
                          7, ACCT_STR) +
191
156k
          1;
192
156k
      break;
193
79.3k
    case 3:
194
79.3k
      eob_pt =
195
79.3k
          aom_read_symbol(r, ec_ctx->eob_flag_cdf128[plane_type][eob_multi_ctx],
196
79.3k
                          8, ACCT_STR) +
197
79.3k
          1;
198
79.3k
      break;
199
134k
    case 4:
200
134k
      eob_pt =
201
134k
          aom_read_symbol(r, ec_ctx->eob_flag_cdf256[plane_type][eob_multi_ctx],
202
134k
                          9, ACCT_STR) +
203
134k
          1;
204
134k
      break;
205
58.0k
    case 5:
206
58.0k
      eob_pt =
207
58.0k
          aom_read_symbol(r, ec_ctx->eob_flag_cdf512[plane_type][eob_multi_ctx],
208
58.0k
                          10, ACCT_STR) +
209
58.0k
          1;
210
58.0k
      break;
211
120k
    case 6:
212
120k
    default:
213
120k
      eob_pt = aom_read_symbol(
214
120k
                   r, ec_ctx->eob_flag_cdf1024[plane_type][eob_multi_ctx], 11,
215
120k
                   ACCT_STR) +
216
120k
               1;
217
120k
      break;
218
771k
  }
219
220
771k
  const int eob_offset_bits = av1_eob_offset_bits[eob_pt];
221
771k
  if (eob_offset_bits > 0) {
222
503k
    const int eob_ctx = eob_pt - 3;
223
503k
    int bit = aom_read_symbol(
224
503k
        r, ec_ctx->eob_extra_cdf[txs_ctx][plane_type][eob_ctx], 2, ACCT_STR);
225
503k
    if (bit) {
226
183k
      eob_extra += (1 << (eob_offset_bits - 1));
227
183k
    }
228
229
1.61M
    for (int i = 1; i < eob_offset_bits; i++) {
230
1.11M
      bit = aom_read_bit(r, ACCT_STR);
231
1.11M
      if (bit) {
232
555k
        eob_extra += (1 << (eob_offset_bits - 1 - i));
233
555k
      }
234
1.11M
    }
235
503k
  }
236
771k
  *eob = rec_eob_pos(eob_pt, eob_extra);
237
238
771k
  if (*eob > 1) {
239
540k
    memset(levels_buf, 0,
240
540k
           sizeof(*levels_buf) *
241
540k
               ((width + TX_PAD_HOR) * (height + TX_PAD_VER) + TX_PAD_END));
242
540k
  }
243
244
771k
  {
245
    // Read the non-zero coefficient with scan index eob-1
246
    // TODO(angiebird): Put this into a function
247
771k
    const int c = *eob - 1;
248
771k
    const int pos = scan[c];
249
771k
    const int coeff_ctx = get_lower_levels_ctx_eob(bwl, height, c);
250
771k
    const int nsymbs = 3;
251
771k
    aom_cdf_prob *cdf =
252
771k
        ec_ctx->coeff_base_eob_cdf[txs_ctx][plane_type][coeff_ctx];
253
771k
    int level = aom_read_symbol(r, cdf, nsymbs, ACCT_STR) + 1;
254
771k
    if (level > NUM_BASE_LEVELS) {
255
17.2k
      const int br_ctx = get_br_ctx_eob(pos, bwl, tx_class);
256
17.2k
      cdf = ec_ctx->coeff_br_cdf[AOMMIN(txs_ctx, TX_32X32)][plane_type][br_ctx];
257
22.7k
      for (int idx = 0; idx < COEFF_BASE_RANGE; idx += BR_CDF_SIZE - 1) {
258
22.5k
        const int k = aom_read_symbol(r, cdf, BR_CDF_SIZE, ACCT_STR);
259
22.5k
        level += k;
260
22.5k
        if (k < BR_CDF_SIZE - 1) break;
261
22.5k
      }
262
17.2k
    }
263
771k
    levels[get_padded_idx(pos, bwl)] = level;
264
771k
  }
265
771k
  if (*eob > 1) {
266
540k
    base_cdf_arr base_cdf = ec_ctx->coeff_base_cdf[txs_ctx][plane_type];
267
540k
    br_cdf_arr br_cdf =
268
540k
        ec_ctx->coeff_br_cdf[AOMMIN(txs_ctx, TX_32X32)][plane_type];
269
540k
    if (tx_class == TX_CLASS_2D) {
270
495k
      read_coeffs_reverse_2d(r, tx_size, 1, *eob - 1 - 1, scan, bwl, levels,
271
495k
                             base_cdf, br_cdf);
272
495k
      read_coeffs_reverse(r, tx_size, tx_class, 0, 0, scan, bwl, levels,
273
495k
                          base_cdf, br_cdf);
274
495k
    } else {
275
44.4k
      read_coeffs_reverse(r, tx_size, tx_class, 0, *eob - 1 - 1, scan, bwl,
276
44.4k
                          levels, base_cdf, br_cdf);
277
44.4k
    }
278
540k
  }
279
280
17.7M
  for (int c = 0; c < *eob; ++c) {
281
16.9M
    const int pos = scan[c];
282
16.9M
    uint8_t sign;
283
16.9M
    tran_low_t level = levels[get_padded_idx(pos, bwl)];
284
16.9M
    if (level) {
285
7.16M
      *max_scan_line = AOMMAX(*max_scan_line, pos);
286
7.16M
      if (c == 0) {
287
629k
        const int dc_sign_ctx = txb_ctx->dc_sign_ctx;
288
629k
        sign = aom_read_symbol(r, ec_ctx->dc_sign_cdf[plane_type][dc_sign_ctx],
289
629k
                               2, ACCT_STR);
290
6.53M
      } else {
291
6.53M
        sign = aom_read_bit(r, ACCT_STR);
292
6.53M
      }
293
7.16M
      if (level >= MAX_BASE_BR_RANGE) {
294
1.25M
        level += read_golomb(xd, r);
295
1.25M
      }
296
297
7.16M
      if (c == 0) dc_val = sign ? -level : level;
298
299
      // Bitmasking to clamp level to valid range:
300
      //   The valid range for 8/10/12 bit vdieo is at most 14/16/18 bit
301
7.16M
      level &= 0xfffff;
302
7.16M
      cul_level += level;
303
7.16M
      tran_low_t dq_coeff;
304
      // Bitmasking to clamp dq_coeff to valid range:
305
      //   The valid range for 8/10/12 bit video is at most 17/19/21 bit
306
7.16M
      dq_coeff = (tran_low_t)(
307
7.16M
          (int64_t)level * get_dqv(dequant, scan[c], iqmatrix) & 0xffffff);
308
7.16M
      dq_coeff = dq_coeff >> shift;
309
7.16M
      if (sign) {
310
3.59M
        dq_coeff = -dq_coeff;
311
3.59M
      }
312
7.16M
      tcoeffs[pos] = clamp(dq_coeff, min_value, max_value);
313
7.16M
    }
314
16.9M
  }
315
316
771k
  cul_level = AOMMIN(COEFF_CONTEXT_MASK, cul_level);
317
318
  // DC value
319
771k
  set_dc_sign(&cul_level, dc_val);
320
321
771k
  return cul_level;
322
771k
}
323
324
void av1_read_coeffs_txb_facade(const AV1_COMMON *const cm,
325
                                DecoderCodingBlock *dcb, aom_reader *const r,
326
                                const int plane, const int row, const int col,
327
1.70M
                                const TX_SIZE tx_size) {
328
#if TXCOEFF_TIMER
329
  struct aom_usec_timer timer;
330
  aom_usec_timer_start(&timer);
331
#endif
332
1.70M
  MACROBLOCKD *const xd = &dcb->xd;
333
1.70M
  MB_MODE_INFO *const mbmi = xd->mi[0];
334
1.70M
  struct macroblockd_plane *const pd = &xd->plane[plane];
335
336
1.70M
  const BLOCK_SIZE bsize = mbmi->bsize;
337
1.70M
  assert(bsize < BLOCK_SIZES_ALL);
338
1.70M
  const BLOCK_SIZE plane_bsize =
339
1.70M
      get_plane_block_size(bsize, pd->subsampling_x, pd->subsampling_y);
340
341
1.70M
  TXB_CTX txb_ctx;
342
1.70M
  get_txb_ctx(plane_bsize, tx_size, plane, pd->above_entropy_context + col,
343
1.70M
              pd->left_entropy_context + row, &txb_ctx);
344
1.70M
  const uint8_t cul_level =
345
1.70M
      av1_read_coeffs_txb(cm, dcb, r, row, col, plane, &txb_ctx, tx_size);
346
1.70M
  av1_set_entropy_contexts(xd, pd, plane, plane_bsize, tx_size, cul_level, col,
347
1.70M
                           row);
348
349
1.70M
  if (is_inter_block(mbmi)) {
350
29.8k
    const PLANE_TYPE plane_type = get_plane_type(plane);
351
    // tx_type will be read out in av1_read_coeffs_txb_facade
352
29.8k
    const TX_TYPE tx_type = av1_get_tx_type(xd, plane_type, row, col, tx_size,
353
29.8k
                                            cm->features.reduced_tx_set_used);
354
355
29.8k
    if (plane == 0) {
356
16.6k
      const int txw = tx_size_wide_unit[tx_size];
357
16.6k
      const int txh = tx_size_high_unit[tx_size];
358
      // The 16x16 unit is due to the constraint from tx_64x64 which sets the
359
      // maximum tx size for chroma as 32x32. Coupled with 4x1 transform block
360
      // size, the constraint takes effect in 32x16 / 16x32 size too. To solve
361
      // the intricacy, cover all the 16x16 units inside a 64 level transform.
362
16.6k
      if (txw == tx_size_wide_unit[TX_64X64] ||
363
16.6k
          txh == tx_size_high_unit[TX_64X64]) {
364
498
        const int tx_unit = tx_size_wide_unit[TX_16X16];
365
498
        const int stride = xd->tx_type_map_stride;
366
2.13k
        for (int idy = 0; idy < txh; idy += tx_unit) {
367
7.62k
          for (int idx = 0; idx < txw; idx += tx_unit) {
368
5.99k
            xd->tx_type_map[(row + idy) * stride + col + idx] = tx_type;
369
5.99k
          }
370
1.63k
        }
371
498
      }
372
16.6k
    }
373
29.8k
  }
374
375
#if TXCOEFF_TIMER
376
  aom_usec_timer_mark(&timer);
377
  const int64_t elapsed_time = aom_usec_timer_elapsed(&timer);
378
  cm->txcoeff_timer += elapsed_time;
379
  ++cm->txb_count;
380
#endif
381
1.70M
}