media/libvpx/vp9/common/vp9_entropy.c

changeset 0
6474c204b198
     1.1 --- /dev/null	Thu Jan 01 00:00:00 1970 +0000
     1.2 +++ b/media/libvpx/vp9/common/vp9_entropy.c	Wed Dec 31 06:09:35 2014 +0100
     1.3 @@ -0,0 +1,424 @@
     1.4 +/*
     1.5 + *  Copyright (c) 2010 The WebM project authors. All Rights Reserved.
     1.6 + *
     1.7 + *  Use of this source code is governed by a BSD-style license
     1.8 + *  that can be found in the LICENSE file in the root of the source
     1.9 + *  tree. An additional intellectual property rights grant can be found
    1.10 + *  in the file PATENTS.  All contributing project authors may
    1.11 + *  be found in the AUTHORS file in the root of the source tree.
    1.12 + */
    1.13 +
    1.14 +#include "vp9/common/vp9_entropy.h"
    1.15 +#include "vp9/common/vp9_blockd.h"
    1.16 +#include "vp9/common/vp9_onyxc_int.h"
    1.17 +#include "vp9/common/vp9_entropymode.h"
    1.18 +#include "vpx_mem/vpx_mem.h"
    1.19 +#include "vpx/vpx_integer.h"
    1.20 +
    1.21 +#define MODEL_NODES (ENTROPY_NODES - UNCONSTRAINED_NODES)
    1.22 +
    1.23 +DECLARE_ALIGNED(16, const uint8_t, vp9_norm[256]) = {
    1.24 +  0, 7, 6, 6, 5, 5, 5, 5, 4, 4, 4, 4, 4, 4, 4, 4,
    1.25 +  3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3,
    1.26 +  2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2,
    1.27 +  2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2,
    1.28 +  1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
    1.29 +  1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
    1.30 +  1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
    1.31 +  1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
    1.32 +  0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
    1.33 +  0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
    1.34 +  0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
    1.35 +  0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
    1.36 +  0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
    1.37 +  0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
    1.38 +  0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
    1.39 +  0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0
    1.40 +};
    1.41 +
    1.42 +DECLARE_ALIGNED(16, const uint8_t,
    1.43 +                vp9_coefband_trans_8x8plus[1024]) = {
    1.44 +  0, 1, 1, 2, 2, 2, 3, 3, 3, 3, 4, 4, 4, 4, 4, 4,
    1.45 +  4, 4, 4, 4, 4, 5,
    1.46 +  // beyond MAXBAND_INDEX+1 all values are filled as 5
    1.47 +                    5, 5, 5, 5, 5, 5, 5, 5, 5, 5,
    1.48 +  5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5,
    1.49 +  5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5,
    1.50 +  5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5,
    1.51 +  5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5,
    1.52 +  5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5,
    1.53 +  5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5,
    1.54 +  5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5,
    1.55 +  5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5,
    1.56 +  5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5,
    1.57 +  5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5,
    1.58 +  5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5,
    1.59 +  5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5,
    1.60 +  5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5,
    1.61 +  5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5,
    1.62 +  5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5,
    1.63 +  5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5,
    1.64 +  5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5,
    1.65 +  5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5,
    1.66 +  5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5,
    1.67 +  5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5,
    1.68 +  5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5,
    1.69 +  5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5,
    1.70 +  5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5,
    1.71 +  5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5,
    1.72 +  5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5,
    1.73 +  5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5,
    1.74 +  5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5,
    1.75 +  5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5,
    1.76 +  5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5,
    1.77 +  5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5,
    1.78 +  5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5,
    1.79 +  5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5,
    1.80 +  5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5,
    1.81 +  5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5,
    1.82 +  5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5,
    1.83 +  5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5,
    1.84 +  5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5,
    1.85 +  5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5,
    1.86 +  5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5,
    1.87 +  5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5,
    1.88 +  5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5,
    1.89 +  5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5,
    1.90 +  5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5,
    1.91 +  5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5,
    1.92 +  5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5,
    1.93 +  5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5,
    1.94 +  5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5,
    1.95 +  5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5,
    1.96 +  5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5,
    1.97 +  5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5,
    1.98 +  5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5,
    1.99 +  5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5,
   1.100 +  5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5,
   1.101 +  5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5,
   1.102 +  5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5,
   1.103 +  5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5,
   1.104 +  5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5,
   1.105 +  5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5,
   1.106 +  5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5,
   1.107 +  5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5,
   1.108 +  5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5,
   1.109 +  5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5,
   1.110 +};
   1.111 +
   1.112 +DECLARE_ALIGNED(16, const uint8_t,
   1.113 +                vp9_coefband_trans_4x4[16]) = {
   1.114 +  0, 1, 1, 2, 2, 2, 3, 3, 3, 3, 4, 4, 4, 5, 5, 5,
   1.115 +};
   1.116 +
   1.117 +DECLARE_ALIGNED(16, const uint8_t, vp9_pt_energy_class[MAX_ENTROPY_TOKENS]) = {
   1.118 +  0, 1, 2, 3, 3, 4, 4, 5, 5, 5, 5, 5
   1.119 +};
   1.120 +
   1.121 +
   1.122 +
   1.123 +/* Array indices are identical to previously-existing CONTEXT_NODE indices */
   1.124 +
   1.125 +const vp9_tree_index vp9_coef_tree[TREE_SIZE(MAX_ENTROPY_TOKENS)] = {
   1.126 +  -DCT_EOB_TOKEN, 2,                          /* 0 = EOB */
   1.127 +  -ZERO_TOKEN, 4,                             /* 1 = ZERO */
   1.128 +  -ONE_TOKEN, 6,                              /* 2 = ONE */
   1.129 +  8, 12,                                      /* 3 = LOW_VAL */
   1.130 +  -TWO_TOKEN, 10,                            /* 4 = TWO */
   1.131 +  -THREE_TOKEN, -FOUR_TOKEN,                /* 5 = THREE */
   1.132 +  14, 16,                                   /* 6 = HIGH_LOW */
   1.133 +  -DCT_VAL_CATEGORY1, -DCT_VAL_CATEGORY2,   /* 7 = CAT_ONE */
   1.134 +  18, 20,                                   /* 8 = CAT_THREEFOUR */
   1.135 +  -DCT_VAL_CATEGORY3, -DCT_VAL_CATEGORY4,   /* 9 = CAT_THREE */
   1.136 +  -DCT_VAL_CATEGORY5, -DCT_VAL_CATEGORY6    /* 10 = CAT_FIVE */
   1.137 +};
   1.138 +
   1.139 +struct vp9_token vp9_coef_encodings[MAX_ENTROPY_TOKENS];
   1.140 +
   1.141 +/* Trees for extra bits.  Probabilities are constant and
   1.142 +   do not depend on previously encoded bits */
   1.143 +
   1.144 +static const vp9_prob Pcat1[] = { 159};
   1.145 +static const vp9_prob Pcat2[] = { 165, 145};
   1.146 +static const vp9_prob Pcat3[] = { 173, 148, 140};
   1.147 +static const vp9_prob Pcat4[] = { 176, 155, 140, 135};
   1.148 +static const vp9_prob Pcat5[] = { 180, 157, 141, 134, 130};
   1.149 +static const vp9_prob Pcat6[] = {
   1.150 +  254, 254, 254, 252, 249, 243, 230, 196, 177, 153, 140, 133, 130, 129
   1.151 +};
   1.152 +
   1.153 +const vp9_tree_index vp9_coefmodel_tree[6] = {
   1.154 +  -DCT_EOB_MODEL_TOKEN, 2,                      /* 0 = EOB */
   1.155 +  -ZERO_TOKEN, 4,                               /* 1 = ZERO */
   1.156 +  -ONE_TOKEN, -TWO_TOKEN,
   1.157 +};
   1.158 +
   1.159 +// Model obtained from a 2-sided zero-centerd distribuition derived
   1.160 +// from a Pareto distribution. The cdf of the distribution is:
   1.161 +// cdf(x) = 0.5 + 0.5 * sgn(x) * [1 - {alpha/(alpha + |x|)} ^ beta]
   1.162 +//
   1.163 +// For a given beta and a given probablity of the 1-node, the alpha
   1.164 +// is first solved, and then the {alpha, beta} pair is used to generate
   1.165 +// the probabilities for the rest of the nodes.
   1.166 +
   1.167 +// beta = 8
   1.168 +static const vp9_prob modelcoefprobs_pareto8[COEFPROB_MODELS][MODEL_NODES] = {
   1.169 +  {  3,  86, 128,   6,  86,  23,  88,  29},
   1.170 +  {  9,  86, 129,  17,  88,  61,  94,  76},
   1.171 +  { 15,  87, 129,  28,  89,  93, 100, 110},
   1.172 +  { 20,  88, 130,  38,  91, 118, 106, 136},
   1.173 +  { 26,  89, 131,  48,  92, 139, 111, 156},
   1.174 +  { 31,  90, 131,  58,  94, 156, 117, 171},
   1.175 +  { 37,  90, 132,  66,  95, 171, 122, 184},
   1.176 +  { 42,  91, 132,  75,  97, 183, 127, 194},
   1.177 +  { 47,  92, 133,  83,  98, 193, 132, 202},
   1.178 +  { 52,  93, 133,  90, 100, 201, 137, 208},
   1.179 +  { 57,  94, 134,  98, 101, 208, 142, 214},
   1.180 +  { 62,  94, 135, 105, 103, 214, 146, 218},
   1.181 +  { 66,  95, 135, 111, 104, 219, 151, 222},
   1.182 +  { 71,  96, 136, 117, 106, 224, 155, 225},
   1.183 +  { 76,  97, 136, 123, 107, 227, 159, 228},
   1.184 +  { 80,  98, 137, 129, 109, 231, 162, 231},
   1.185 +  { 84,  98, 138, 134, 110, 234, 166, 233},
   1.186 +  { 89,  99, 138, 140, 112, 236, 170, 235},
   1.187 +  { 93, 100, 139, 145, 113, 238, 173, 236},
   1.188 +  { 97, 101, 140, 149, 115, 240, 176, 238},
   1.189 +  {101, 102, 140, 154, 116, 242, 179, 239},
   1.190 +  {105, 103, 141, 158, 118, 243, 182, 240},
   1.191 +  {109, 104, 141, 162, 119, 244, 185, 241},
   1.192 +  {113, 104, 142, 166, 120, 245, 187, 242},
   1.193 +  {116, 105, 143, 170, 122, 246, 190, 243},
   1.194 +  {120, 106, 143, 173, 123, 247, 192, 244},
   1.195 +  {123, 107, 144, 177, 125, 248, 195, 244},
   1.196 +  {127, 108, 145, 180, 126, 249, 197, 245},
   1.197 +  {130, 109, 145, 183, 128, 249, 199, 245},
   1.198 +  {134, 110, 146, 186, 129, 250, 201, 246},
   1.199 +  {137, 111, 147, 189, 131, 251, 203, 246},
   1.200 +  {140, 112, 147, 192, 132, 251, 205, 247},
   1.201 +  {143, 113, 148, 194, 133, 251, 207, 247},
   1.202 +  {146, 114, 149, 197, 135, 252, 208, 248},
   1.203 +  {149, 115, 149, 199, 136, 252, 210, 248},
   1.204 +  {152, 115, 150, 201, 138, 252, 211, 248},
   1.205 +  {155, 116, 151, 204, 139, 253, 213, 249},
   1.206 +  {158, 117, 151, 206, 140, 253, 214, 249},
   1.207 +  {161, 118, 152, 208, 142, 253, 216, 249},
   1.208 +  {163, 119, 153, 210, 143, 253, 217, 249},
   1.209 +  {166, 120, 153, 212, 144, 254, 218, 250},
   1.210 +  {168, 121, 154, 213, 146, 254, 220, 250},
   1.211 +  {171, 122, 155, 215, 147, 254, 221, 250},
   1.212 +  {173, 123, 155, 217, 148, 254, 222, 250},
   1.213 +  {176, 124, 156, 218, 150, 254, 223, 250},
   1.214 +  {178, 125, 157, 220, 151, 254, 224, 251},
   1.215 +  {180, 126, 157, 221, 152, 254, 225, 251},
   1.216 +  {183, 127, 158, 222, 153, 254, 226, 251},
   1.217 +  {185, 128, 159, 224, 155, 255, 227, 251},
   1.218 +  {187, 129, 160, 225, 156, 255, 228, 251},
   1.219 +  {189, 131, 160, 226, 157, 255, 228, 251},
   1.220 +  {191, 132, 161, 227, 159, 255, 229, 251},
   1.221 +  {193, 133, 162, 228, 160, 255, 230, 252},
   1.222 +  {195, 134, 163, 230, 161, 255, 231, 252},
   1.223 +  {197, 135, 163, 231, 162, 255, 231, 252},
   1.224 +  {199, 136, 164, 232, 163, 255, 232, 252},
   1.225 +  {201, 137, 165, 233, 165, 255, 233, 252},
   1.226 +  {202, 138, 166, 233, 166, 255, 233, 252},
   1.227 +  {204, 139, 166, 234, 167, 255, 234, 252},
   1.228 +  {206, 140, 167, 235, 168, 255, 235, 252},
   1.229 +  {207, 141, 168, 236, 169, 255, 235, 252},
   1.230 +  {209, 142, 169, 237, 171, 255, 236, 252},
   1.231 +  {210, 144, 169, 237, 172, 255, 236, 252},
   1.232 +  {212, 145, 170, 238, 173, 255, 237, 252},
   1.233 +  {214, 146, 171, 239, 174, 255, 237, 253},
   1.234 +  {215, 147, 172, 240, 175, 255, 238, 253},
   1.235 +  {216, 148, 173, 240, 176, 255, 238, 253},
   1.236 +  {218, 149, 173, 241, 177, 255, 239, 253},
   1.237 +  {219, 150, 174, 241, 179, 255, 239, 253},
   1.238 +  {220, 152, 175, 242, 180, 255, 240, 253},
   1.239 +  {222, 153, 176, 242, 181, 255, 240, 253},
   1.240 +  {223, 154, 177, 243, 182, 255, 240, 253},
   1.241 +  {224, 155, 178, 244, 183, 255, 241, 253},
   1.242 +  {225, 156, 178, 244, 184, 255, 241, 253},
   1.243 +  {226, 158, 179, 244, 185, 255, 242, 253},
   1.244 +  {228, 159, 180, 245, 186, 255, 242, 253},
   1.245 +  {229, 160, 181, 245, 187, 255, 242, 253},
   1.246 +  {230, 161, 182, 246, 188, 255, 243, 253},
   1.247 +  {231, 163, 183, 246, 189, 255, 243, 253},
   1.248 +  {232, 164, 184, 247, 190, 255, 243, 253},
   1.249 +  {233, 165, 185, 247, 191, 255, 244, 253},
   1.250 +  {234, 166, 185, 247, 192, 255, 244, 253},
   1.251 +  {235, 168, 186, 248, 193, 255, 244, 253},
   1.252 +  {236, 169, 187, 248, 194, 255, 244, 253},
   1.253 +  {236, 170, 188, 248, 195, 255, 245, 253},
   1.254 +  {237, 171, 189, 249, 196, 255, 245, 254},
   1.255 +  {238, 173, 190, 249, 197, 255, 245, 254},
   1.256 +  {239, 174, 191, 249, 198, 255, 245, 254},
   1.257 +  {240, 175, 192, 249, 199, 255, 246, 254},
   1.258 +  {240, 177, 193, 250, 200, 255, 246, 254},
   1.259 +  {241, 178, 194, 250, 201, 255, 246, 254},
   1.260 +  {242, 179, 195, 250, 202, 255, 246, 254},
   1.261 +  {242, 181, 196, 250, 203, 255, 247, 254},
   1.262 +  {243, 182, 197, 251, 204, 255, 247, 254},
   1.263 +  {244, 184, 198, 251, 205, 255, 247, 254},
   1.264 +  {244, 185, 199, 251, 206, 255, 247, 254},
   1.265 +  {245, 186, 200, 251, 207, 255, 247, 254},
   1.266 +  {246, 188, 201, 252, 207, 255, 248, 254},
   1.267 +  {246, 189, 202, 252, 208, 255, 248, 254},
   1.268 +  {247, 191, 203, 252, 209, 255, 248, 254},
   1.269 +  {247, 192, 204, 252, 210, 255, 248, 254},
   1.270 +  {248, 194, 205, 252, 211, 255, 248, 254},
   1.271 +  {248, 195, 206, 252, 212, 255, 249, 254},
   1.272 +  {249, 197, 207, 253, 213, 255, 249, 254},
   1.273 +  {249, 198, 208, 253, 214, 255, 249, 254},
   1.274 +  {250, 200, 210, 253, 215, 255, 249, 254},
   1.275 +  {250, 201, 211, 253, 215, 255, 249, 254},
   1.276 +  {250, 203, 212, 253, 216, 255, 249, 254},
   1.277 +  {251, 204, 213, 253, 217, 255, 250, 254},
   1.278 +  {251, 206, 214, 254, 218, 255, 250, 254},
   1.279 +  {252, 207, 216, 254, 219, 255, 250, 254},
   1.280 +  {252, 209, 217, 254, 220, 255, 250, 254},
   1.281 +  {252, 211, 218, 254, 221, 255, 250, 254},
   1.282 +  {253, 213, 219, 254, 222, 255, 250, 254},
   1.283 +  {253, 214, 221, 254, 223, 255, 250, 254},
   1.284 +  {253, 216, 222, 254, 224, 255, 251, 254},
   1.285 +  {253, 218, 224, 254, 225, 255, 251, 254},
   1.286 +  {254, 220, 225, 254, 225, 255, 251, 254},
   1.287 +  {254, 222, 227, 255, 226, 255, 251, 254},
   1.288 +  {254, 224, 228, 255, 227, 255, 251, 254},
   1.289 +  {254, 226, 230, 255, 228, 255, 251, 254},
   1.290 +  {255, 228, 231, 255, 230, 255, 251, 254},
   1.291 +  {255, 230, 233, 255, 231, 255, 252, 254},
   1.292 +  {255, 232, 235, 255, 232, 255, 252, 254},
   1.293 +  {255, 235, 237, 255, 233, 255, 252, 254},
   1.294 +  {255, 238, 240, 255, 235, 255, 252, 255},
   1.295 +  {255, 241, 243, 255, 236, 255, 252, 254},
   1.296 +  {255, 246, 247, 255, 239, 255, 253, 255}
   1.297 +};
   1.298 +
   1.299 +static void extend_model_to_full_distribution(vp9_prob p,
   1.300 +                                              vp9_prob *tree_probs) {
   1.301 +  const int l = (p - 1) / 2;
   1.302 +  const vp9_prob (*model)[MODEL_NODES] = modelcoefprobs_pareto8;
   1.303 +  if (p & 1) {
   1.304 +    vpx_memcpy(tree_probs + UNCONSTRAINED_NODES,
   1.305 +               model[l], MODEL_NODES * sizeof(vp9_prob));
   1.306 +  } else {
   1.307 +    // interpolate
   1.308 +    int i;
   1.309 +    for (i = UNCONSTRAINED_NODES; i < ENTROPY_NODES; ++i)
   1.310 +      tree_probs[i] = (model[l][i - UNCONSTRAINED_NODES] +
   1.311 +                       model[l + 1][i - UNCONSTRAINED_NODES]) >> 1;
   1.312 +  }
   1.313 +}
   1.314 +
   1.315 +void vp9_model_to_full_probs(const vp9_prob *model, vp9_prob *full) {
   1.316 +  if (full != model)
   1.317 +    vpx_memcpy(full, model, sizeof(vp9_prob) * UNCONSTRAINED_NODES);
   1.318 +  extend_model_to_full_distribution(model[PIVOT_NODE], full);
   1.319 +}
   1.320 +
   1.321 +static vp9_tree_index cat1[2], cat2[4], cat3[6], cat4[8], cat5[10], cat6[28];
   1.322 +
   1.323 +static void init_bit_tree(vp9_tree_index *p, int n) {
   1.324 +  int i = 0;
   1.325 +
   1.326 +  while (++i < n) {
   1.327 +    p[0] = p[1] = i << 1;
   1.328 +    p += 2;
   1.329 +  }
   1.330 +
   1.331 +  p[0] = p[1] = 0;
   1.332 +}
   1.333 +
   1.334 +static void init_bit_trees() {
   1.335 +  init_bit_tree(cat1, 1);
   1.336 +  init_bit_tree(cat2, 2);
   1.337 +  init_bit_tree(cat3, 3);
   1.338 +  init_bit_tree(cat4, 4);
   1.339 +  init_bit_tree(cat5, 5);
   1.340 +  init_bit_tree(cat6, 14);
   1.341 +}
   1.342 +
   1.343 +const vp9_extra_bit vp9_extra_bits[MAX_ENTROPY_TOKENS] = {
   1.344 +  { 0, 0, 0, 0},
   1.345 +  { 0, 0, 0, 1},
   1.346 +  { 0, 0, 0, 2},
   1.347 +  { 0, 0, 0, 3},
   1.348 +  { 0, 0, 0, 4},
   1.349 +  { cat1, Pcat1, 1, 5},
   1.350 +  { cat2, Pcat2, 2, 7},
   1.351 +  { cat3, Pcat3, 3, 11},
   1.352 +  { cat4, Pcat4, 4, 19},
   1.353 +  { cat5, Pcat5, 5, 35},
   1.354 +  { cat6, Pcat6, 14, 67},
   1.355 +  { 0, 0, 0, 0}
   1.356 +};
   1.357 +
   1.358 +#include "vp9/common/vp9_default_coef_probs.h"
   1.359 +
   1.360 +void vp9_default_coef_probs(VP9_COMMON *cm) {
   1.361 +  vp9_copy(cm->fc.coef_probs[TX_4X4], default_coef_probs_4x4);
   1.362 +  vp9_copy(cm->fc.coef_probs[TX_8X8], default_coef_probs_8x8);
   1.363 +  vp9_copy(cm->fc.coef_probs[TX_16X16], default_coef_probs_16x16);
   1.364 +  vp9_copy(cm->fc.coef_probs[TX_32X32], default_coef_probs_32x32);
   1.365 +}
   1.366 +
   1.367 +void vp9_coef_tree_initialize() {
   1.368 +  init_bit_trees();
   1.369 +  vp9_tokens_from_tree(vp9_coef_encodings, vp9_coef_tree);
   1.370 +}
   1.371 +
   1.372 +// #define COEF_COUNT_TESTING
   1.373 +
   1.374 +#define COEF_COUNT_SAT 24
   1.375 +#define COEF_MAX_UPDATE_FACTOR 112
   1.376 +#define COEF_COUNT_SAT_KEY 24
   1.377 +#define COEF_MAX_UPDATE_FACTOR_KEY 112
   1.378 +#define COEF_COUNT_SAT_AFTER_KEY 24
   1.379 +#define COEF_MAX_UPDATE_FACTOR_AFTER_KEY 128
   1.380 +
   1.381 +static void adapt_coef_probs(VP9_COMMON *cm, TX_SIZE tx_size,
   1.382 +                             unsigned int count_sat,
   1.383 +                             unsigned int update_factor) {
   1.384 +  const FRAME_CONTEXT *pre_fc = &cm->frame_contexts[cm->frame_context_idx];
   1.385 +
   1.386 +  vp9_coeff_probs_model *dst_coef_probs = cm->fc.coef_probs[tx_size];
   1.387 +  const vp9_coeff_probs_model *pre_coef_probs = pre_fc->coef_probs[tx_size];
   1.388 +  vp9_coeff_count_model *coef_counts = cm->counts.coef[tx_size];
   1.389 +  unsigned int (*eob_branch_count)[REF_TYPES][COEF_BANDS][PREV_COEF_CONTEXTS] =
   1.390 +      cm->counts.eob_branch[tx_size];
   1.391 +  int i, j, k, l, m;
   1.392 +  unsigned int branch_ct[UNCONSTRAINED_NODES][2];
   1.393 +
   1.394 +  for (i = 0; i < BLOCK_TYPES; ++i)
   1.395 +    for (j = 0; j < REF_TYPES; ++j)
   1.396 +      for (k = 0; k < COEF_BANDS; ++k)
   1.397 +        for (l = 0; l < PREV_COEF_CONTEXTS; ++l) {
   1.398 +          if (l >= 3 && k == 0)
   1.399 +            continue;
   1.400 +          vp9_tree_probs_from_distribution(vp9_coefmodel_tree, branch_ct,
   1.401 +                                           coef_counts[i][j][k][l]);
   1.402 +          branch_ct[0][1] = eob_branch_count[i][j][k][l] - branch_ct[0][0];
   1.403 +          for (m = 0; m < UNCONSTRAINED_NODES; ++m)
   1.404 +            dst_coef_probs[i][j][k][l][m] = merge_probs(
   1.405 +                                                pre_coef_probs[i][j][k][l][m],
   1.406 +                                                branch_ct[m],
   1.407 +                                                count_sat, update_factor);
   1.408 +        }
   1.409 +}
   1.410 +
   1.411 +void vp9_adapt_coef_probs(VP9_COMMON *cm) {
   1.412 +  TX_SIZE t;
   1.413 +  unsigned int count_sat, update_factor;
   1.414 +
   1.415 +  if (frame_is_intra_only(cm)) {
   1.416 +    update_factor = COEF_MAX_UPDATE_FACTOR_KEY;
   1.417 +    count_sat = COEF_COUNT_SAT_KEY;
   1.418 +  } else if (cm->last_frame_type == KEY_FRAME) {
   1.419 +    update_factor = COEF_MAX_UPDATE_FACTOR_AFTER_KEY;  /* adapt quickly */
   1.420 +    count_sat = COEF_COUNT_SAT_AFTER_KEY;
   1.421 +  } else {
   1.422 +    update_factor = COEF_MAX_UPDATE_FACTOR;
   1.423 +    count_sat = COEF_COUNT_SAT;
   1.424 +  }
   1.425 +  for (t = TX_4X4; t <= TX_32X32; t++)
   1.426 +    adapt_coef_probs(cm, t, count_sat, update_factor);
   1.427 +}

mercurial