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
read_golomb(MACROBLOCKD * xd,aom_reader * r)22 static int read_golomb(MACROBLOCKD *xd, aom_reader *r) {
23 int x = 1;
24 int length = 0;
25 int i = 0;
26
27 while (!i) {
28 i = aom_read_bit(r, ACCT_STR);
29 ++length;
30 if (length > 20) {
31 aom_internal_error(xd->error_info, AOM_CODEC_CORRUPT_FRAME,
32 "Invalid length in read_golomb");
33 break;
34 }
35 }
36
37 for (i = 0; i < length - 1; ++i) {
38 x <<= 1;
39 x += aom_read_bit(r, ACCT_STR);
40 }
41
42 return x - 1;
43 }
44
rec_eob_pos(const int eob_token,const int extra)45 static INLINE int rec_eob_pos(const int eob_token, const int extra) {
46 int eob = k_eob_group_start[eob_token];
47 if (eob > 2) {
48 eob += extra;
49 }
50 return eob;
51 }
52
get_dqv(const int16_t * dequant,int coeff_idx,const qm_val_t * iqmatrix)53 static INLINE int get_dqv(const int16_t *dequant, int coeff_idx,
54 const qm_val_t *iqmatrix) {
55 int dqv = dequant[!!coeff_idx];
56 if (iqmatrix != NULL)
57 dqv =
58 ((iqmatrix[coeff_idx] * dqv) + (1 << (AOM_QM_BITS - 1))) >> AOM_QM_BITS;
59 return dqv;
60 }
61
read_coeffs_reverse_2d(aom_reader * r,TX_SIZE tx_size,int start_si,int end_si,const int16_t * scan,int bwl,uint8_t * levels,base_cdf_arr base_cdf,br_cdf_arr br_cdf)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 br_cdf_arr br_cdf) {
68 for (int c = end_si; c >= start_si; --c) {
69 const int pos = scan[c];
70 const int coeff_ctx = get_lower_levels_ctx_2d(levels, pos, bwl, tx_size);
71 const int nsymbs = 4;
72 int level = aom_read_symbol(r, base_cdf[coeff_ctx], nsymbs, ACCT_STR);
73 if (level > NUM_BASE_LEVELS) {
74 const int br_ctx = get_br_ctx_2d(levels, pos, bwl);
75 aom_cdf_prob *cdf = br_cdf[br_ctx];
76 for (int idx = 0; idx < COEFF_BASE_RANGE; idx += BR_CDF_SIZE - 1) {
77 const int k = aom_read_symbol(r, cdf, BR_CDF_SIZE, ACCT_STR);
78 level += k;
79 if (k < BR_CDF_SIZE - 1) break;
80 }
81 }
82 levels[get_padded_idx(pos, bwl)] = level;
83 }
84 }
85
read_coeffs_reverse(aom_reader * r,TX_SIZE tx_size,TX_CLASS tx_class,int start_si,int end_si,const int16_t * scan,int bwl,uint8_t * levels,base_cdf_arr base_cdf,br_cdf_arr br_cdf)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 br_cdf_arr br_cdf) {
91 for (int c = end_si; c >= start_si; --c) {
92 const int pos = scan[c];
93 const int coeff_ctx =
94 get_lower_levels_ctx(levels, pos, bwl, tx_size, tx_class);
95 const int nsymbs = 4;
96 int level = aom_read_symbol(r, base_cdf[coeff_ctx], nsymbs, ACCT_STR);
97 if (level > NUM_BASE_LEVELS) {
98 const int br_ctx = get_br_ctx(levels, pos, bwl, tx_class);
99 aom_cdf_prob *cdf = br_cdf[br_ctx];
100 for (int idx = 0; idx < COEFF_BASE_RANGE; idx += BR_CDF_SIZE - 1) {
101 const int k = aom_read_symbol(r, cdf, BR_CDF_SIZE, ACCT_STR);
102 level += k;
103 if (k < BR_CDF_SIZE - 1) break;
104 }
105 }
106 levels[get_padded_idx(pos, bwl)] = level;
107 }
108 }
109
av1_read_coeffs_txb(const AV1_COMMON * const cm,MACROBLOCKD * const xd,aom_reader * const r,const int blk_row,const int blk_col,const int plane,const TXB_CTX * const txb_ctx,const TX_SIZE tx_size)110 uint8_t av1_read_coeffs_txb(const AV1_COMMON *const cm, MACROBLOCKD *const xd,
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 const TX_SIZE tx_size) {
115 FRAME_CONTEXT *const ec_ctx = xd->tile_ctx;
116 const int32_t max_value = (1 << (7 + xd->bd)) - 1;
117 const int32_t min_value = -(1 << (7 + xd->bd));
118 const TX_SIZE txs_ctx = get_txsize_entropy_ctx(tx_size);
119 const PLANE_TYPE plane_type = get_plane_type(plane);
120 MB_MODE_INFO *const mbmi = xd->mi[0];
121 struct macroblockd_plane *const pd = &xd->plane[plane];
122 const int16_t *const dequant = pd->seg_dequant_QTX[mbmi->segment_id];
123 tran_low_t *const tcoeffs = pd->dqcoeff_block + xd->cb_offset[plane];
124 const int shift = av1_get_tx_scale(tx_size);
125 const int bwl = get_txb_bwl(tx_size);
126 const int width = get_txb_wide(tx_size);
127 const int height = get_txb_high(tx_size);
128 int cul_level = 0;
129 int dc_val = 0;
130 uint8_t levels_buf[TX_PAD_2D];
131 uint8_t *const levels = set_levels(levels_buf, width);
132 const int all_zero = aom_read_symbol(
133 r, ec_ctx->txb_skip_cdf[txs_ctx][txb_ctx->txb_skip_ctx], 2, ACCT_STR);
134 eob_info *eob_data = pd->eob_data + xd->txb_offset[plane];
135 uint16_t *const eob = &(eob_data->eob);
136 uint16_t *const max_scan_line = &(eob_data->max_scan_line);
137 *max_scan_line = 0;
138 *eob = 0;
139 if (all_zero) {
140 *max_scan_line = 0;
141 if (plane == 0) {
142 const int txk_type_idx =
143 av1_get_txk_type_index(mbmi->sb_type, blk_row, blk_col);
144 mbmi->txk_type[txk_type_idx] = DCT_DCT;
145 }
146 return 0;
147 }
148
149 memset(levels_buf, 0,
150 sizeof(*levels_buf) *
151 ((width + TX_PAD_HOR) * (height + TX_PAD_VER) + TX_PAD_END));
152 if (plane == AOM_PLANE_Y) {
153 // only y plane's tx_type is transmitted
154 av1_read_tx_type(cm, xd, blk_row, blk_col, tx_size, r);
155 }
156 const TX_TYPE tx_type = av1_get_tx_type(plane_type, xd, blk_row, blk_col,
157 tx_size, cm->reduced_tx_set_used);
158 const TX_CLASS tx_class = tx_type_to_class[tx_type];
159 const TX_SIZE qm_tx_size = av1_get_adjusted_tx_size(tx_size);
160 const qm_val_t *iqmatrix =
161 IS_2D_TRANSFORM(tx_type)
162 ? pd->seg_iqmatrix[mbmi->segment_id][qm_tx_size]
163 : cm->giqmatrix[NUM_QM_LEVELS - 1][0][qm_tx_size];
164 const SCAN_ORDER *const scan_order = get_scan(tx_size, tx_type);
165 const int16_t *const scan = scan_order->scan;
166 int eob_extra = 0;
167 int eob_pt = 1;
168
169 const int eob_multi_size = txsize_log2_minus4[tx_size];
170 const int eob_multi_ctx = (tx_class == TX_CLASS_2D) ? 0 : 1;
171 switch (eob_multi_size) {
172 case 0:
173 eob_pt =
174 aom_read_symbol(r, ec_ctx->eob_flag_cdf16[plane_type][eob_multi_ctx],
175 5, ACCT_STR) +
176 1;
177 break;
178 case 1:
179 eob_pt =
180 aom_read_symbol(r, ec_ctx->eob_flag_cdf32[plane_type][eob_multi_ctx],
181 6, ACCT_STR) +
182 1;
183 break;
184 case 2:
185 eob_pt =
186 aom_read_symbol(r, ec_ctx->eob_flag_cdf64[plane_type][eob_multi_ctx],
187 7, ACCT_STR) +
188 1;
189 break;
190 case 3:
191 eob_pt =
192 aom_read_symbol(r, ec_ctx->eob_flag_cdf128[plane_type][eob_multi_ctx],
193 8, ACCT_STR) +
194 1;
195 break;
196 case 4:
197 eob_pt =
198 aom_read_symbol(r, ec_ctx->eob_flag_cdf256[plane_type][eob_multi_ctx],
199 9, ACCT_STR) +
200 1;
201 break;
202 case 5:
203 eob_pt =
204 aom_read_symbol(r, ec_ctx->eob_flag_cdf512[plane_type][eob_multi_ctx],
205 10, ACCT_STR) +
206 1;
207 break;
208 case 6:
209 default:
210 eob_pt = aom_read_symbol(
211 r, ec_ctx->eob_flag_cdf1024[plane_type][eob_multi_ctx], 11,
212 ACCT_STR) +
213 1;
214 break;
215 }
216
217 if (k_eob_offset_bits[eob_pt] > 0) {
218 const int eob_ctx = eob_pt - 3;
219 int bit = aom_read_symbol(
220 r, ec_ctx->eob_extra_cdf[txs_ctx][plane_type][eob_ctx], 2, ACCT_STR);
221 if (bit) {
222 eob_extra += (1 << (k_eob_offset_bits[eob_pt] - 1));
223 }
224
225 for (int i = 1; i < k_eob_offset_bits[eob_pt]; i++) {
226 bit = aom_read_bit(r, ACCT_STR);
227 if (bit) {
228 eob_extra += (1 << (k_eob_offset_bits[eob_pt] - 1 - i));
229 }
230 }
231 }
232 *eob = rec_eob_pos(eob_pt, eob_extra);
233
234 {
235 // Read the non-zero coefficient with scan index eob-1
236 // TODO(angiebird): Put this into a function
237 const int c = *eob - 1;
238 const int pos = scan[c];
239 const int coeff_ctx = get_lower_levels_ctx_eob(bwl, height, c);
240 const int nsymbs = 3;
241 aom_cdf_prob *cdf =
242 ec_ctx->coeff_base_eob_cdf[txs_ctx][plane_type][coeff_ctx];
243 int level = aom_read_symbol(r, cdf, nsymbs, ACCT_STR) + 1;
244 if (level > NUM_BASE_LEVELS) {
245 const int br_ctx = get_br_ctx(levels, pos, bwl, tx_class);
246 for (int idx = 0; idx < COEFF_BASE_RANGE; idx += BR_CDF_SIZE - 1) {
247 const int k = aom_read_symbol(
248 r,
249 ec_ctx->coeff_br_cdf[AOMMIN(txs_ctx, TX_32X32)][plane_type][br_ctx],
250 BR_CDF_SIZE, ACCT_STR);
251 level += k;
252 if (k < BR_CDF_SIZE - 1) break;
253 }
254 }
255 levels[get_padded_idx(pos, bwl)] = level;
256 }
257 if (*eob > 1) {
258 base_cdf_arr base_cdf = ec_ctx->coeff_base_cdf[txs_ctx][plane_type];
259 br_cdf_arr br_cdf =
260 ec_ctx->coeff_br_cdf[AOMMIN(txs_ctx, TX_32X32)][plane_type];
261 if (tx_class == TX_CLASS_2D) {
262 read_coeffs_reverse_2d(r, tx_size, 1, *eob - 1 - 1, scan, bwl, levels,
263 base_cdf, br_cdf);
264 read_coeffs_reverse(r, tx_size, tx_class, 0, 0, scan, bwl, levels,
265 base_cdf, br_cdf);
266 } else {
267 read_coeffs_reverse(r, tx_size, tx_class, 0, *eob - 1 - 1, scan, bwl,
268 levels, base_cdf, br_cdf);
269 }
270 }
271
272 int16_t num_zero_coeffs = 0;
273 for (int c = 0; c < *eob; ++c) {
274 const int pos = scan[c];
275 num_zero_coeffs = AOMMAX(num_zero_coeffs, pos);
276 }
277 memset(tcoeffs, 0, (num_zero_coeffs + 1) * sizeof(tcoeffs[0]));
278
279 for (int c = 0; c < *eob; ++c) {
280 const int pos = scan[c];
281 uint8_t sign;
282 tran_low_t level = levels[get_padded_idx(pos, bwl)];
283 if (level) {
284 *max_scan_line = AOMMAX(*max_scan_line, pos);
285 if (c == 0) {
286 const int dc_sign_ctx = txb_ctx->dc_sign_ctx;
287 sign = aom_read_symbol(r, ec_ctx->dc_sign_cdf[plane_type][dc_sign_ctx],
288 2, ACCT_STR);
289 } else {
290 sign = aom_read_bit(r, ACCT_STR);
291 }
292 if (level >= MAX_BASE_BR_RANGE) {
293 level += read_golomb(xd, r);
294 }
295
296 if (c == 0) dc_val = sign ? -level : level;
297
298 // Bitmasking to clamp level to valid range:
299 // The valid range for 8/10/12 bit vdieo is at most 14/16/18 bit
300 level &= 0xfffff;
301 cul_level += level;
302 tran_low_t dq_coeff;
303 // Bitmasking to clamp dq_coeff to valid range:
304 // The valid range for 8/10/12 bit video is at most 17/19/21 bit
305 dq_coeff = (tran_low_t)(
306 (int64_t)level * get_dqv(dequant, scan[c], iqmatrix) & 0xffffff);
307 dq_coeff = dq_coeff >> shift;
308 if (sign) {
309 dq_coeff = -dq_coeff;
310 }
311 tcoeffs[pos] = clamp(dq_coeff, min_value, max_value);
312 }
313 }
314
315 cul_level = AOMMIN(COEFF_CONTEXT_MASK, cul_level);
316
317 // DC value
318 set_dc_sign(&cul_level, dc_val);
319
320 return cul_level;
321 }
322
av1_read_coeffs_txb_facade(const AV1_COMMON * const cm,MACROBLOCKD * const xd,aom_reader * const r,const int plane,const int row,const int col,const TX_SIZE tx_size)323 void av1_read_coeffs_txb_facade(const AV1_COMMON *const cm,
324 MACROBLOCKD *const xd, aom_reader *const r,
325 const int plane, const int row, const int col,
326 const TX_SIZE tx_size) {
327 #if TXCOEFF_TIMER
328 struct aom_usec_timer timer;
329 aom_usec_timer_start(&timer);
330 #endif
331 MB_MODE_INFO *const mbmi = xd->mi[0];
332 struct macroblockd_plane *const pd = &xd->plane[plane];
333
334 const BLOCK_SIZE bsize = mbmi->sb_type;
335 const BLOCK_SIZE plane_bsize =
336 get_plane_block_size(bsize, pd->subsampling_x, pd->subsampling_y);
337
338 TXB_CTX txb_ctx;
339 get_txb_ctx(plane_bsize, tx_size, plane, pd->above_context + col,
340 pd->left_context + row, &txb_ctx);
341 const uint8_t cul_level =
342 av1_read_coeffs_txb(cm, xd, r, row, col, plane, &txb_ctx, tx_size);
343 av1_set_contexts(xd, pd, plane, plane_bsize, tx_size, cul_level, col, row);
344
345 if (is_inter_block(mbmi)) {
346 PLANE_TYPE plane_type = get_plane_type(plane);
347 // tx_type will be read out in av1_read_coeffs_txb_facade
348 const TX_TYPE tx_type = av1_get_tx_type(plane_type, xd, row, col, tx_size,
349 cm->reduced_tx_set_used);
350
351 if (plane == 0)
352 update_txk_array(mbmi->txk_type, mbmi->sb_type, row, col, tx_size,
353 tx_type);
354 }
355
356 #if TXCOEFF_TIMER
357 aom_usec_timer_mark(&timer);
358 const int64_t elapsed_time = aom_usec_timer_elapsed(&timer);
359 cm->txcoeff_timer += elapsed_time;
360 ++cm->txb_count;
361 #endif
362 }
363