| /* |
| * Copyright (c) 2017, Alliance for Open Media. All rights reserved |
| * |
| * 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. |
| */ |
| |
| #include "aom_ports/mem.h" |
| #include "av1/common/scan.h" |
| #include "av1/common/idct.h" |
| #include "av1/common/txb_common.h" |
| #include "av1/decoder/decodemv.h" |
| #include "av1/decoder/decodetxb.h" |
| #include "av1/decoder/symbolrate.h" |
| |
| #define ACCT_STR __func__ |
| |
| static int read_golomb(MACROBLOCKD *xd, aom_reader *r, FRAME_COUNTS *counts) { |
| #if !CONFIG_SYMBOLRATE |
| (void)counts; |
| #endif |
| int x = 1; |
| int length = 0; |
| int i = 0; |
| |
| while (!i) { |
| i = av1_read_record_bit(counts, r, ACCT_STR); |
| ++length; |
| if (length >= 32) { |
| aom_internal_error(xd->error_info, AOM_CODEC_CORRUPT_FRAME, |
| "Invalid length in read_golomb"); |
| break; |
| } |
| } |
| |
| for (i = 0; i < length - 1; ++i) { |
| x <<= 1; |
| x += av1_read_record_bit(counts, r, ACCT_STR); |
| } |
| |
| return x - 1; |
| } |
| |
| static INLINE int rec_eob_pos(const int eob_token, const int extra) { |
| int eob = k_eob_group_start[eob_token]; |
| if (eob > 2) { |
| eob += extra; |
| } |
| return eob; |
| } |
| |
| uint8_t 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, |
| int16_t *const max_scan_line, int *const eob) { |
| FRAME_CONTEXT *const ec_ctx = xd->tile_ctx; |
| #if TXCOEFF_TIMER |
| FRAME_COUNTS *const counts = NULL; |
| #else |
| FRAME_COUNTS *const counts = xd->counts; |
| #endif |
| const TX_SIZE txs_ctx = get_txsize_entropy_ctx(tx_size); |
| const PLANE_TYPE plane_type = get_plane_type(plane); |
| MB_MODE_INFO *const mbmi = &xd->mi[0]->mbmi; |
| const int seg_eob = av1_get_max_eob(tx_size); |
| int c = 0; |
| int num_updates = 0; |
| struct macroblockd_plane *const pd = &xd->plane[plane]; |
| const int16_t *const dequant = pd->seg_dequant_QTX[mbmi->segment_id]; |
| tran_low_t *const tcoeffs = pd->dqcoeff; |
| #if !CONFIG_DAALA_TX |
| const int shift = av1_get_tx_scale(tx_size); |
| #endif |
| const int bwl = get_txb_bwl(tx_size); |
| const int width = get_txb_wide(tx_size); |
| const int height = get_txb_high(tx_size); |
| int cul_level = 0; |
| uint8_t levels_buf[TX_PAD_2D]; |
| uint8_t *const levels = set_levels(levels_buf, width); |
| DECLARE_ALIGNED(16, uint8_t, level_counts[MAX_TX_SQUARE]); |
| int8_t signs[MAX_TX_SQUARE]; |
| uint16_t update_pos[MAX_TX_SQUARE]; |
| |
| const int all_zero = av1_read_record_bin( |
| counts, r, ec_ctx->txb_skip_cdf[txs_ctx][txb_ctx->txb_skip_ctx], 2, |
| ACCT_STR); |
| // printf("txb_skip: %d %2d\n", txs_ctx, txb_ctx->txb_skip_ctx); |
| if (xd->counts) |
| ++xd->counts->txb_skip[txs_ctx][txb_ctx->txb_skip_ctx][all_zero]; |
| *eob = 0; |
| if (all_zero) { |
| *max_scan_line = 0; |
| #if CONFIG_TXK_SEL |
| if (plane == 0) |
| mbmi->txk_type[(blk_row << MAX_MIB_SIZE_LOG2) + blk_col] = DCT_DCT; |
| #endif |
| return 0; |
| } |
| |
| memset(levels_buf, 0, |
| sizeof(*levels_buf) * |
| ((width + TX_PAD_HOR) * (height + TX_PAD_VER) + TX_PAD_END)); |
| |
| (void)blk_row; |
| (void)blk_col; |
| #if CONFIG_TXK_SEL |
| av1_read_tx_type(cm, xd, blk_row, blk_col, plane, get_min_tx_size(tx_size), |
| r); |
| #endif |
| const TX_TYPE tx_type = |
| av1_get_tx_type(plane_type, xd, blk_row, blk_col, tx_size); |
| const SCAN_ORDER *const scan_order = get_scan(cm, tx_size, tx_type, mbmi); |
| const int16_t *const scan = scan_order->scan; |
| int dummy; |
| const int max_eob_pt = get_eob_pos_token(seg_eob, &dummy); |
| int eob_extra = 0; |
| int eob_pt = 1; |
| |
| for (eob_pt = 1; eob_pt < max_eob_pt; eob_pt++) { |
| const int eob_pos_ctx = av1_get_eob_pos_ctx(tx_type, eob_pt); |
| const int is_equal = av1_read_record_bin( |
| counts, r, ec_ctx->eob_flag_cdf[txs_ctx][plane_type][eob_pos_ctx], 2, |
| ACCT_STR); |
| // printf("eob_flag_cdf: %d %d %2d\n", txs_ctx, plane_type, eob_pos_ctx); |
| // aom_read_symbol(r, |
| // ec_ctx->eob_flag_cdf[AOMMIN(txs_ctx,3)][plane_type][eob_pos_ctx], 2, |
| // ACCT_STR); |
| if (counts) ++counts->eob_flag[txs_ctx][plane_type][eob_pos_ctx][is_equal]; |
| |
| if (is_equal) { |
| break; |
| } |
| } |
| |
| // printf("Dec: "); |
| if (k_eob_offset_bits[eob_pt] > 0) { |
| int bit = av1_read_record_bin( |
| counts, r, ec_ctx->eob_extra_cdf[txs_ctx][plane_type][eob_pt], 2, |
| ACCT_STR); |
| // printf("eob_extra_cdf: %d %d %2d\n", txs_ctx, plane_type, eob_pt); |
| if (counts) ++counts->eob_extra[txs_ctx][plane_type][eob_pt][bit]; |
| if (bit) { |
| eob_extra += (1 << (k_eob_offset_bits[eob_pt] - 1)); |
| } |
| |
| for (int i = 1; i < k_eob_offset_bits[eob_pt]; i++) { |
| bit = av1_read_record_bit(counts, r, ACCT_STR); |
| // printf("eob_bit:\n"); |
| if (bit) { |
| eob_extra += (1 << (k_eob_offset_bits[eob_pt] - 1 - i)); |
| } |
| // printf("%d ", bit); |
| } |
| } |
| *eob = rec_eob_pos(eob_pt, eob_extra); |
| // printf("=>[%d, %d], (%d, %d)\n", seg_eob, *eob, eob_pt, eob_extra); |
| |
| for (int i = 0; i < *eob; ++i) { |
| c = *eob - 1 - i; |
| const int pos = scan[c]; |
| #if CONFIG_LV_MAP_MULTI |
| const int coeff_ctx = get_nz_map_ctx(levels, pos, bwl, height, c, |
| c == *eob - 1, tx_size, tx_type); |
| #if USE_BASE_EOB_ALPHABET |
| aom_cdf_prob *cdf; |
| int nsymbs; |
| if (c == *eob - 1) { |
| cdf = ec_ctx->coeff_base_eob_cdf[txs_ctx][plane_type] |
| [coeff_ctx - SIG_COEF_CONTEXTS + |
| SIG_COEF_CONTEXTS_EOB]; |
| nsymbs = 3; |
| } else { |
| cdf = ec_ctx->coeff_base_cdf[txs_ctx][plane_type][coeff_ctx]; |
| nsymbs = 4; |
| } |
| const int level = av1_read_record_symbol(counts, r, cdf, nsymbs, ACCT_STR) + |
| (c == *eob - 1); |
| if (counts) { |
| if (c < *eob - 1) { |
| ++counts->nz_map[txs_ctx][plane_type][coeff_ctx][level != 0]; |
| } |
| if (level != 0) { |
| for (int k = 0; k < NUM_BASE_LEVELS; ++k) { |
| ++counts |
| ->coeff_base[txs_ctx][plane_type][k][coeff_ctx][level > k + 1]; |
| if (level == k + 1) break; |
| } |
| } |
| } |
| #else |
| const int level = av1_read_record_symbol( |
| counts, r, ec_ctx->coeff_base_cdf[txs_ctx][plane_type][coeff_ctx], 4, |
| ACCT_STR); |
| #endif |
| // printf("base_cdf: %d %d %2d\n", txs_ctx, plane_type, coeff_ctx); |
| // printf("base_cdf: %d %d %2d : %3d %3d %3d\n", txs_ctx, plane_type, |
| // coeff_ctx, |
| // ec_ctx->coeff_base_cdf[txs_ctx][plane_type][coeff_ctx][0]>>7, |
| // ec_ctx->coeff_base_cdf[txs_ctx][plane_type][coeff_ctx][1]>>7, |
| // ec_ctx->coeff_base_cdf[txs_ctx][plane_type][coeff_ctx][2]>>7); |
| if (level) { |
| levels[get_padded_idx(pos, bwl)] = level; |
| *max_scan_line = AOMMAX(*max_scan_line, pos); |
| if (level < 3) { |
| cul_level += level; |
| #if CONFIG_DAALA_TX |
| tcoeffs[pos] = level * dequant[!!c]; |
| #else |
| tcoeffs[pos] = (level * dequant[!!c]) >> shift; |
| #endif |
| } else { |
| update_pos[num_updates++] = pos; |
| } |
| } |
| #else |
| int is_nz; |
| const int coeff_ctx = get_nz_map_ctx(levels, pos, bwl, tx_size, tx_type); |
| |
| if (c < *eob - 1) { |
| is_nz = av1_read_record_bin( |
| counts, r, ec_ctx->nz_map_cdf[txs_ctx][plane_type][coeff_ctx], 2, |
| ACCT_STR); |
| } else { |
| is_nz = 1; |
| } |
| |
| if (is_nz) { |
| int k; |
| for (k = 0; k < NUM_BASE_LEVELS; ++k) { |
| const int ctx = coeff_ctx; |
| const int is_k = av1_read_record_bin( |
| counts, r, ec_ctx->coeff_base_cdf[txs_ctx][plane_type][k][ctx], 2, |
| ACCT_STR); |
| if (counts) ++counts->coeff_base[txs_ctx][plane_type][k][ctx][is_k]; |
| |
| // semantic: is_k = 1 if level > (k+1) |
| if (is_k == 0) { |
| cul_level += k + 1; |
| #if CONFIG_DAALA_TX |
| tcoeffs[pos] = (k + 1) * dequant[!!c]; |
| #else |
| tcoeffs[pos] = ((k + 1) * dequant[!!c]) >> shift; |
| #endif |
| break; |
| } |
| } |
| levels[get_padded_idx(pos, bwl)] = k + 1; |
| *max_scan_line = AOMMAX(*max_scan_line, pos); |
| if (k == NUM_BASE_LEVELS) { |
| update_pos[num_updates++] = pos; |
| } |
| } |
| #endif |
| } |
| |
| // Loop to decode all signs in the transform block, |
| // starting with the sign of the DC (if applicable) |
| for (c = 0; c < *eob; ++c) { |
| const int pos = scan[c]; |
| int8_t *const sign = &signs[pos]; |
| if (levels[get_padded_idx(pos, bwl)] == 0) continue; |
| if (c == 0) { |
| const int dc_sign_ctx = txb_ctx->dc_sign_ctx; |
| #if LV_MAP_PROB |
| *sign = av1_read_record_bin( |
| counts, r, ec_ctx->dc_sign_cdf[plane_type][dc_sign_ctx], 2, ACCT_STR); |
| // printf("dc_sign: %d %d\n", plane_type, dc_sign_ctx); |
| #else |
| *sign = aom_read(r, ec_ctx->dc_sign[plane_type][dc_sign_ctx], ACCT_STR); |
| #endif |
| if (counts) ++counts->dc_sign[plane_type][dc_sign_ctx][*sign]; |
| } else { |
| *sign = av1_read_record_bit(counts, r, ACCT_STR); |
| } |
| if (*sign) tcoeffs[pos] = -tcoeffs[pos]; |
| } |
| |
| if (num_updates) { |
| #if !CONFIG_LV_MAP_MULTI |
| av1_get_br_level_counts(levels, width, height, level_counts); |
| #endif |
| for (c = 0; c < num_updates; ++c) { |
| const int pos = update_pos[c]; |
| uint8_t *const level = &levels[get_padded_idx(pos, bwl)]; |
| int idx; |
| int ctx; |
| |
| assert(*level > NUM_BASE_LEVELS); |
| |
| #if !CONFIG_LV_MAP_MULTI |
| ctx = get_br_ctx(levels, pos, bwl, level_counts[pos]); |
| #endif |
| |
| #if CONFIG_LV_MAP_MULTI |
| #if USE_CAUSAL_BR_CTX |
| ctx = get_br_ctx(levels, pos, bwl, level_counts[pos], tx_type); |
| #else |
| ctx = get_br_ctx(levels, pos, bwl, level_counts[pos]); |
| #endif |
| for (idx = 0; idx < COEFF_BASE_RANGE / (BR_CDF_SIZE - 1); ++idx) { |
| int k = av1_read_record_symbol( |
| counts, r, |
| ec_ctx->coeff_br_cdf[AOMMIN(txs_ctx, TX_16X16)][plane_type][ctx], |
| BR_CDF_SIZE, ACCT_STR); |
| *level += k; |
| if (counts) { |
| for (int lps = 0; lps < BR_CDF_SIZE - 1; lps++) { |
| ++counts->coeff_lps[AOMMIN(txs_ctx, TX_16X16)][plane_type][lps][ctx] |
| [lps == k]; |
| if (lps == k) break; |
| } |
| } |
| if (k < BR_CDF_SIZE - 1) break; |
| } |
| if (*level <= NUM_BASE_LEVELS + COEFF_BASE_RANGE) { |
| cul_level += *level; |
| #if CONFIG_DAALA_TX |
| tran_low_t t = *level * dequant[!!pos]; |
| #else |
| tran_low_t t = (*level * dequant[!!pos]) >> shift; |
| #endif |
| if (signs[pos]) t = -t; |
| tcoeffs[pos] = t; |
| continue; |
| } |
| #else |
| for (idx = 0; idx < BASE_RANGE_SETS; ++idx) { |
| // printf("br: %d %d %d %d\n", txs_ctx, plane_type, idx, ctx); |
| if (av1_read_record_bin( |
| counts, r, ec_ctx->coeff_br_cdf[txs_ctx][plane_type][idx][ctx], |
| 2, ACCT_STR)) { |
| const int extra_bits = (1 << br_extra_bits[idx]) - 1; |
| // int br_offset = aom_read_literal(r, extra_bits, ACCT_STR); |
| int br_offset = 0; |
| int tok; |
| if (counts) ++counts->coeff_br[txs_ctx][plane_type][idx][ctx][1]; |
| for (tok = 0; tok < extra_bits; ++tok) { |
| if (av1_read_record_bin( |
| counts, r, ec_ctx->coeff_lps_cdf[txs_ctx][plane_type][ctx], |
| 2, ACCT_STR)) { |
| br_offset = tok; |
| if (counts) ++counts->coeff_lps[txs_ctx][plane_type][ctx][1]; |
| break; |
| } |
| if (counts) ++counts->coeff_lps[txs_ctx][plane_type][ctx][0]; |
| } |
| if (tok == extra_bits) br_offset = extra_bits; |
| |
| const int br_base = br_index_to_coeff[idx]; |
| |
| *level = NUM_BASE_LEVELS + 1 + br_base + br_offset; |
| cul_level += *level; |
| #if CONFIG_DAALA_TX |
| tran_low_t t = *level * dequant[!!pos]; |
| #else |
| tran_low_t t = (*level * dequant[!!pos]) >> shift; |
| #endif |
| if (signs[pos]) t = -t; |
| tcoeffs[pos] = t; |
| break; |
| } |
| if (counts) ++counts->coeff_br[txs_ctx][plane_type][idx][ctx][0]; |
| } |
| |
| if (idx < BASE_RANGE_SETS) continue; |
| #endif |
| // decode 0-th order Golomb code |
| *level = COEFF_BASE_RANGE + 1 + NUM_BASE_LEVELS; |
| // Save golomb in tcoeffs because adding it to level may incur overflow |
| tran_low_t t = *level + read_golomb(xd, r, counts); |
| cul_level += t; |
| #if CONFIG_DAALA_TX |
| t = t * dequant[!!pos]; |
| #else |
| t = (t * dequant[!!pos]) >> shift; |
| #endif |
| if (signs[pos]) t = -t; |
| tcoeffs[pos] = t; |
| } |
| } |
| |
| cul_level = AOMMIN(63, cul_level); |
| |
| // DC value |
| set_dc_sign(&cul_level, tcoeffs[0]); |
| |
| return cul_level; |
| } |
| |
| uint8_t av1_read_coeffs_txb_facade(const AV1_COMMON *const cm, |
| MACROBLOCKD *const xd, aom_reader *const r, |
| const int row, const int col, |
| const int plane, const TX_SIZE tx_size, |
| int16_t *const max_scan_line, |
| int *const eob) { |
| MB_MODE_INFO *const mbmi = &xd->mi[0]->mbmi; |
| struct macroblockd_plane *const pd = &xd->plane[plane]; |
| |
| const BLOCK_SIZE bsize = mbmi->sb_type; |
| const BLOCK_SIZE plane_bsize = |
| AOMMAX(BLOCK_4X4, get_plane_block_size(bsize, pd)); |
| |
| TXB_CTX txb_ctx; |
| get_txb_ctx(plane_bsize, tx_size, plane, pd->above_context + col, |
| pd->left_context + row, &txb_ctx); |
| uint8_t cul_level = av1_read_coeffs_txb(cm, xd, r, row, col, plane, &txb_ctx, |
| tx_size, max_scan_line, eob); |
| #if CONFIG_ADAPT_SCAN |
| PLANE_TYPE plane_type = get_plane_type(plane); |
| TX_TYPE tx_type = av1_get_tx_type(plane_type, xd, row, col, tx_size); |
| const int mi_row = -xd->mb_to_top_edge >> (3 + MI_SIZE_LOG2); |
| |
| if (xd->counts && *eob > 0) |
| av1_update_scan_count_facade(cm, xd, mi_row, tx_size, tx_type, pd->dqcoeff, |
| *eob); |
| #endif |
| av1_set_contexts(xd, pd, plane, tx_size, cul_level, col, row); |
| return cul_level; |
| } |
