michael@0: /*
michael@0:  *  Copyright (c) 2010 The WebM project authors. All Rights Reserved.
michael@0:  *
michael@0:  *  Use of this source code is governed by a BSD-style license
michael@0:  *  that can be found in the LICENSE file in the root of the source
michael@0:  *  tree. An additional intellectual property rights grant can be found
michael@0:  *  in the file PATENTS.  All contributing project authors may
michael@0:  *  be found in the AUTHORS file in the root of the source tree.
michael@0:  */
michael@0: 
michael@0: #include "vpx_mem/vpx_mem.h"
michael@0: #include "vpx_ports/mem.h"
michael@0: 
michael@0: #include "vp9/common/vp9_blockd.h"
michael@0: #include "vp9/common/vp9_common.h"
michael@0: #include "vp9/common/vp9_seg_common.h"
michael@0: 
michael@0: #include "vp9/decoder/vp9_dboolhuff.h"
michael@0: #include "vp9/decoder/vp9_detokenize.h"
michael@0: #include "vp9/decoder/vp9_onyxd_int.h"
michael@0: #include "vp9/decoder/vp9_treereader.h"
michael@0: 
michael@0: #define EOB_CONTEXT_NODE            0
michael@0: #define ZERO_CONTEXT_NODE           1
michael@0: #define ONE_CONTEXT_NODE            2
michael@0: #define LOW_VAL_CONTEXT_NODE        3
michael@0: #define TWO_CONTEXT_NODE            4
michael@0: #define THREE_CONTEXT_NODE          5
michael@0: #define HIGH_LOW_CONTEXT_NODE       6
michael@0: #define CAT_ONE_CONTEXT_NODE        7
michael@0: #define CAT_THREEFOUR_CONTEXT_NODE  8
michael@0: #define CAT_THREE_CONTEXT_NODE      9
michael@0: #define CAT_FIVE_CONTEXT_NODE       10
michael@0: 
michael@0: #define CAT1_MIN_VAL    5
michael@0: #define CAT2_MIN_VAL    7
michael@0: #define CAT3_MIN_VAL   11
michael@0: #define CAT4_MIN_VAL   19
michael@0: #define CAT5_MIN_VAL   35
michael@0: #define CAT6_MIN_VAL   67
michael@0: #define CAT1_PROB0    159
michael@0: #define CAT2_PROB0    145
michael@0: #define CAT2_PROB1    165
michael@0: 
michael@0: #define CAT3_PROB0 140
michael@0: #define CAT3_PROB1 148
michael@0: #define CAT3_PROB2 173
michael@0: 
michael@0: #define CAT4_PROB0 135
michael@0: #define CAT4_PROB1 140
michael@0: #define CAT4_PROB2 155
michael@0: #define CAT4_PROB3 176
michael@0: 
michael@0: #define CAT5_PROB0 130
michael@0: #define CAT5_PROB1 134
michael@0: #define CAT5_PROB2 141
michael@0: #define CAT5_PROB3 157
michael@0: #define CAT5_PROB4 180
michael@0: 
michael@0: static const vp9_prob cat6_prob[15] = {
michael@0:   254, 254, 254, 252, 249, 243, 230, 196, 177, 153, 140, 133, 130, 129, 0
michael@0: };
michael@0: 
michael@0: static const int token_to_counttoken[MAX_ENTROPY_TOKENS] = {
michael@0:   ZERO_TOKEN, ONE_TOKEN, TWO_TOKEN, TWO_TOKEN,
michael@0:   TWO_TOKEN, TWO_TOKEN, TWO_TOKEN, TWO_TOKEN,
michael@0:   TWO_TOKEN, TWO_TOKEN, TWO_TOKEN, DCT_EOB_MODEL_TOKEN
michael@0: };
michael@0: 
michael@0: #define INCREMENT_COUNT(token)                              \
michael@0:   do {                                                      \
michael@0:      if (!cm->frame_parallel_decoding_mode)                 \
michael@0:        ++coef_counts[band][pt][token_to_counttoken[token]]; \
michael@0:   } while (0)
michael@0: 
michael@0: 
michael@0: #define WRITE_COEF_CONTINUE(val, token)                  \
michael@0:   {                                                      \
michael@0:     dqcoeff_ptr[scan[c]] = (vp9_read_bit(r) ? -val : val) * \
michael@0:                             dq[c > 0] / (1 + (tx_size == TX_32X32)); \
michael@0:     INCREMENT_COUNT(token);                              \
michael@0:     token_cache[scan[c]] = vp9_pt_energy_class[token];   \
michael@0:     ++c;                                                 \
michael@0:     continue;                                            \
michael@0:   }
michael@0: 
michael@0: #define ADJUST_COEF(prob, bits_count)                   \
michael@0:   do {                                                  \
michael@0:     val += (vp9_read(r, prob) << bits_count);           \
michael@0:   } while (0)
michael@0: 
michael@0: static int decode_coefs(VP9_COMMON *cm, const MACROBLOCKD *xd,
michael@0:                         vp9_reader *r, int block_idx,
michael@0:                         PLANE_TYPE type, int seg_eob, int16_t *dqcoeff_ptr,
michael@0:                         TX_SIZE tx_size, const int16_t *dq, int pt,
michael@0:                         uint8_t *token_cache) {
michael@0:   const FRAME_CONTEXT *const fc = &cm->fc;
michael@0:   FRAME_COUNTS *const counts = &cm->counts;
michael@0:   const int ref = is_inter_block(&xd->mi_8x8[0]->mbmi);
michael@0:   int band, c = 0;
michael@0:   const vp9_prob (*coef_probs)[PREV_COEF_CONTEXTS][UNCONSTRAINED_NODES] =
michael@0:       fc->coef_probs[tx_size][type][ref];
michael@0:   vp9_prob coef_probs_full[COEF_BANDS][PREV_COEF_CONTEXTS][ENTROPY_NODES];
michael@0:   uint8_t load_map[COEF_BANDS][PREV_COEF_CONTEXTS] = { { 0 } };
michael@0:   const vp9_prob *prob;
michael@0:   unsigned int (*coef_counts)[PREV_COEF_CONTEXTS][UNCONSTRAINED_NODES + 1] =
michael@0:       counts->coef[tx_size][type][ref];
michael@0:   unsigned int (*eob_branch_count)[PREV_COEF_CONTEXTS] =
michael@0:       counts->eob_branch[tx_size][type][ref];
michael@0:   const int16_t *scan, *nb;
michael@0:   const uint8_t *cat6;
michael@0:   const uint8_t *band_translate = get_band_translate(tx_size);
michael@0:   get_scan(xd, tx_size, type, block_idx, &scan, &nb);
michael@0: 
michael@0:   while (c < seg_eob) {
michael@0:     int val;
michael@0:     if (c)
michael@0:       pt = get_coef_context(nb, token_cache, c);
michael@0:     band = *band_translate++;
michael@0:     prob = coef_probs[band][pt];
michael@0:     if (!cm->frame_parallel_decoding_mode)
michael@0:       ++eob_branch_count[band][pt];
michael@0:     if (!vp9_read(r, prob[EOB_CONTEXT_NODE]))
michael@0:       break;
michael@0:     goto DECODE_ZERO;
michael@0: 
michael@0:   SKIP_START:
michael@0:     if (c >= seg_eob)
michael@0:       break;
michael@0:     if (c)
michael@0:       pt = get_coef_context(nb, token_cache, c);
michael@0:     band = *band_translate++;
michael@0:     prob = coef_probs[band][pt];
michael@0: 
michael@0:   DECODE_ZERO:
michael@0:     if (!vp9_read(r, prob[ZERO_CONTEXT_NODE])) {
michael@0:       INCREMENT_COUNT(ZERO_TOKEN);
michael@0:       token_cache[scan[c]] = vp9_pt_energy_class[ZERO_TOKEN];
michael@0:       ++c;
michael@0:       goto SKIP_START;
michael@0:     }
michael@0: 
michael@0:     // ONE_CONTEXT_NODE_0_
michael@0:     if (!vp9_read(r, prob[ONE_CONTEXT_NODE])) {
michael@0:       WRITE_COEF_CONTINUE(1, ONE_TOKEN);
michael@0:     }
michael@0:     // Load full probabilities if not already loaded
michael@0:     if (!load_map[band][pt]) {
michael@0:       vp9_model_to_full_probs(coef_probs[band][pt],
michael@0:                               coef_probs_full[band][pt]);
michael@0:       load_map[band][pt] = 1;
michael@0:     }
michael@0:     prob = coef_probs_full[band][pt];
michael@0:     // LOW_VAL_CONTEXT_NODE_0_
michael@0:     if (!vp9_read(r, prob[LOW_VAL_CONTEXT_NODE])) {
michael@0:       if (!vp9_read(r, prob[TWO_CONTEXT_NODE])) {
michael@0:         WRITE_COEF_CONTINUE(2, TWO_TOKEN);
michael@0:       }
michael@0:       if (!vp9_read(r, prob[THREE_CONTEXT_NODE])) {
michael@0:         WRITE_COEF_CONTINUE(3, THREE_TOKEN);
michael@0:       }
michael@0:       WRITE_COEF_CONTINUE(4, FOUR_TOKEN);
michael@0:     }
michael@0:     // HIGH_LOW_CONTEXT_NODE_0_
michael@0:     if (!vp9_read(r, prob[HIGH_LOW_CONTEXT_NODE])) {
michael@0:       if (!vp9_read(r, prob[CAT_ONE_CONTEXT_NODE])) {
michael@0:         val = CAT1_MIN_VAL;
michael@0:         ADJUST_COEF(CAT1_PROB0, 0);
michael@0:         WRITE_COEF_CONTINUE(val, DCT_VAL_CATEGORY1);
michael@0:       }
michael@0:       val = CAT2_MIN_VAL;
michael@0:       ADJUST_COEF(CAT2_PROB1, 1);
michael@0:       ADJUST_COEF(CAT2_PROB0, 0);
michael@0:       WRITE_COEF_CONTINUE(val, DCT_VAL_CATEGORY2);
michael@0:     }
michael@0:     // CAT_THREEFOUR_CONTEXT_NODE_0_
michael@0:     if (!vp9_read(r, prob[CAT_THREEFOUR_CONTEXT_NODE])) {
michael@0:       if (!vp9_read(r, prob[CAT_THREE_CONTEXT_NODE])) {
michael@0:         val = CAT3_MIN_VAL;
michael@0:         ADJUST_COEF(CAT3_PROB2, 2);
michael@0:         ADJUST_COEF(CAT3_PROB1, 1);
michael@0:         ADJUST_COEF(CAT3_PROB0, 0);
michael@0:         WRITE_COEF_CONTINUE(val, DCT_VAL_CATEGORY3);
michael@0:       }
michael@0:       val = CAT4_MIN_VAL;
michael@0:       ADJUST_COEF(CAT4_PROB3, 3);
michael@0:       ADJUST_COEF(CAT4_PROB2, 2);
michael@0:       ADJUST_COEF(CAT4_PROB1, 1);
michael@0:       ADJUST_COEF(CAT4_PROB0, 0);
michael@0:       WRITE_COEF_CONTINUE(val, DCT_VAL_CATEGORY4);
michael@0:     }
michael@0:     // CAT_FIVE_CONTEXT_NODE_0_:
michael@0:     if (!vp9_read(r, prob[CAT_FIVE_CONTEXT_NODE])) {
michael@0:       val = CAT5_MIN_VAL;
michael@0:       ADJUST_COEF(CAT5_PROB4, 4);
michael@0:       ADJUST_COEF(CAT5_PROB3, 3);
michael@0:       ADJUST_COEF(CAT5_PROB2, 2);
michael@0:       ADJUST_COEF(CAT5_PROB1, 1);
michael@0:       ADJUST_COEF(CAT5_PROB0, 0);
michael@0:       WRITE_COEF_CONTINUE(val, DCT_VAL_CATEGORY5);
michael@0:     }
michael@0:     val = 0;
michael@0:     cat6 = cat6_prob;
michael@0:     while (*cat6)
michael@0:       val = (val << 1) | vp9_read(r, *cat6++);
michael@0:     val += CAT6_MIN_VAL;
michael@0: 
michael@0:     WRITE_COEF_CONTINUE(val, DCT_VAL_CATEGORY6);
michael@0:   }
michael@0: 
michael@0:   if (c < seg_eob) {
michael@0:     if (!cm->frame_parallel_decoding_mode)
michael@0:       ++coef_counts[band][pt][DCT_EOB_MODEL_TOKEN];
michael@0:   }
michael@0: 
michael@0:   return c;
michael@0: }
michael@0: 
michael@0: int vp9_decode_block_tokens(VP9_COMMON *cm, MACROBLOCKD *xd,
michael@0:                             int plane, int block, BLOCK_SIZE plane_bsize,
michael@0:                             int x, int y, TX_SIZE tx_size, vp9_reader *r,
michael@0:                             uint8_t *token_cache) {
michael@0:   struct macroblockd_plane *const pd = &xd->plane[plane];
michael@0:   const int seg_eob = get_tx_eob(&cm->seg, xd->mi_8x8[0]->mbmi.segment_id,
michael@0:                                  tx_size);
michael@0:   const int pt = get_entropy_context(tx_size, pd->above_context + x,
michael@0:                                               pd->left_context + y);
michael@0:   const int eob = decode_coefs(cm, xd, r, block, pd->plane_type, seg_eob,
michael@0:                                BLOCK_OFFSET(pd->dqcoeff, block), tx_size,
michael@0:                                pd->dequant, pt, token_cache);
michael@0:   set_contexts(xd, pd, plane_bsize, tx_size, eob > 0, x, y);
michael@0:   pd->eobs[block] = eob;
michael@0:   return eob;
michael@0: }
michael@0: 
michael@0: