1*77c1e3ccSAndroid Build Coastguard Worker /*
2*77c1e3ccSAndroid Build Coastguard Worker * Copyright (c) 2017, Alliance for Open Media. All rights reserved.
3*77c1e3ccSAndroid Build Coastguard Worker *
4*77c1e3ccSAndroid Build Coastguard Worker * This source code is subject to the terms of the BSD 2 Clause License and
5*77c1e3ccSAndroid Build Coastguard Worker * the Alliance for Open Media Patent License 1.0. If the BSD 2 Clause License
6*77c1e3ccSAndroid Build Coastguard Worker * was not distributed with this source code in the LICENSE file, you can
7*77c1e3ccSAndroid Build Coastguard Worker * obtain it at www.aomedia.org/license/software. If the Alliance for Open
8*77c1e3ccSAndroid Build Coastguard Worker * Media Patent License 1.0 was not distributed with this source code in the
9*77c1e3ccSAndroid Build Coastguard Worker * PATENTS file, you can obtain it at www.aomedia.org/license/patent.
10*77c1e3ccSAndroid Build Coastguard Worker */
11*77c1e3ccSAndroid Build Coastguard Worker
12*77c1e3ccSAndroid Build Coastguard Worker #include "aom_dsp/binary_codes_reader.h"
13*77c1e3ccSAndroid Build Coastguard Worker #include "aom_dsp/recenter.h"
14*77c1e3ccSAndroid Build Coastguard Worker
15*77c1e3ccSAndroid Build Coastguard Worker #define read_primitive_quniform(r, n, ACCT_STR_NAME) \
16*77c1e3ccSAndroid Build Coastguard Worker read_primitive_quniform_(r, n ACCT_STR_ARG(ACCT_STR_NAME))
17*77c1e3ccSAndroid Build Coastguard Worker #define read_primitive_subexpfin(r, n, k, ACCT_STR_NAME) \
18*77c1e3ccSAndroid Build Coastguard Worker read_primitive_subexpfin_(r, n, k ACCT_STR_ARG(ACCT_STR_NAME))
19*77c1e3ccSAndroid Build Coastguard Worker
read_primitive_quniform_(aom_reader * r,uint16_t n ACCT_STR_PARAM)20*77c1e3ccSAndroid Build Coastguard Worker static uint16_t read_primitive_quniform_(aom_reader *r,
21*77c1e3ccSAndroid Build Coastguard Worker uint16_t n ACCT_STR_PARAM) {
22*77c1e3ccSAndroid Build Coastguard Worker if (n <= 1) return 0;
23*77c1e3ccSAndroid Build Coastguard Worker const int l = get_msb(n) + 1;
24*77c1e3ccSAndroid Build Coastguard Worker const int m = (1 << l) - n;
25*77c1e3ccSAndroid Build Coastguard Worker const int v = aom_read_literal(r, l - 1, ACCT_STR_NAME);
26*77c1e3ccSAndroid Build Coastguard Worker return v < m ? v : (v << 1) - m + aom_read_bit(r, ACCT_STR_NAME);
27*77c1e3ccSAndroid Build Coastguard Worker }
28*77c1e3ccSAndroid Build Coastguard Worker
29*77c1e3ccSAndroid Build Coastguard Worker // Decode finite subexponential code that for a symbol v in [0, n-1] with
30*77c1e3ccSAndroid Build Coastguard Worker // parameter k
read_primitive_subexpfin_(aom_reader * r,uint16_t n,uint16_t k ACCT_STR_PARAM)31*77c1e3ccSAndroid Build Coastguard Worker static uint16_t read_primitive_subexpfin_(aom_reader *r, uint16_t n,
32*77c1e3ccSAndroid Build Coastguard Worker uint16_t k ACCT_STR_PARAM) {
33*77c1e3ccSAndroid Build Coastguard Worker int i = 0;
34*77c1e3ccSAndroid Build Coastguard Worker int mk = 0;
35*77c1e3ccSAndroid Build Coastguard Worker
36*77c1e3ccSAndroid Build Coastguard Worker while (1) {
37*77c1e3ccSAndroid Build Coastguard Worker int b = (i ? k + i - 1 : k);
38*77c1e3ccSAndroid Build Coastguard Worker int a = (1 << b);
39*77c1e3ccSAndroid Build Coastguard Worker
40*77c1e3ccSAndroid Build Coastguard Worker if (n <= mk + 3 * a) {
41*77c1e3ccSAndroid Build Coastguard Worker return read_primitive_quniform(r, n - mk, ACCT_STR_NAME) + mk;
42*77c1e3ccSAndroid Build Coastguard Worker }
43*77c1e3ccSAndroid Build Coastguard Worker
44*77c1e3ccSAndroid Build Coastguard Worker if (!aom_read_bit(r, ACCT_STR_NAME)) {
45*77c1e3ccSAndroid Build Coastguard Worker return aom_read_literal(r, b, ACCT_STR_NAME) + mk;
46*77c1e3ccSAndroid Build Coastguard Worker }
47*77c1e3ccSAndroid Build Coastguard Worker
48*77c1e3ccSAndroid Build Coastguard Worker i = i + 1;
49*77c1e3ccSAndroid Build Coastguard Worker mk += a;
50*77c1e3ccSAndroid Build Coastguard Worker }
51*77c1e3ccSAndroid Build Coastguard Worker
52*77c1e3ccSAndroid Build Coastguard Worker assert(0);
53*77c1e3ccSAndroid Build Coastguard Worker return 0;
54*77c1e3ccSAndroid Build Coastguard Worker }
55*77c1e3ccSAndroid Build Coastguard Worker
aom_read_primitive_refsubexpfin_(aom_reader * r,uint16_t n,uint16_t k,uint16_t ref ACCT_STR_PARAM)56*77c1e3ccSAndroid Build Coastguard Worker uint16_t aom_read_primitive_refsubexpfin_(aom_reader *r, uint16_t n, uint16_t k,
57*77c1e3ccSAndroid Build Coastguard Worker uint16_t ref ACCT_STR_PARAM) {
58*77c1e3ccSAndroid Build Coastguard Worker return inv_recenter_finite_nonneg(
59*77c1e3ccSAndroid Build Coastguard Worker n, ref, read_primitive_subexpfin(r, n, k, ACCT_STR_NAME));
60*77c1e3ccSAndroid Build Coastguard Worker }
61