/* * Copyright (c) 2010 The WebM project authors. All Rights Reserved. * * Use of this source code is governed by a BSD-style license * that can be found in the LICENSE file in the root of the source * tree. An additional intellectual property rights grant can be found * in the file PATENTS. All contributing project authors may * be found in the AUTHORS file in the root of the source tree. */ #include #include "vp9/common/vp9_entropy.h" #include "string.h" #include "vp9/common/vp9_blockd.h" #include "vp9/common/vp9_onyxc_int.h" #include "vp9/common/vp9_entropymode.h" #include "vpx_mem/vpx_mem.h" #include "vpx/vpx_integer.h" #include "vp9/common/vp9_coefupdateprobs.h" const int vp9_i8x8_block[4] = {0, 2, 8, 10}; DECLARE_ALIGNED(16, const uint8_t, vp9_norm[256]) = { 0, 7, 6, 6, 5, 5, 5, 5, 4, 4, 4, 4, 4, 4, 4, 4, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 }; // Unified coefficient band structure used by all block sizes DECLARE_ALIGNED(16, const int, vp9_coef_bands[32]) = { 0, 1, 2, 3, 3, 3, 4, 4, 4, 4, 4, 4, 4, 4, 4, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5 }; DECLARE_ALIGNED(16, const int, vp9_coef_bands4x4[16]) = { 0, 1, 2, 3, 3, 3, 4, 4, 4, 4, 5, 5, 5, 5, 5, 5 }; DECLARE_ALIGNED(16, const uint8_t, vp9_pt_energy_class[MAX_ENTROPY_TOKENS]) = { 0, 1, 2, 3, 3, 4, 4, 5, 5, 5, 5, 5 }; DECLARE_ALIGNED(16, const int, vp9_default_zig_zag1d_4x4[16]) = { 0, 1, 4, 8, 5, 2, 3, 6, 9, 12, 13, 10, 7, 11, 14, 15, }; DECLARE_ALIGNED(16, const int, vp9_col_scan_4x4[16]) = { 0, 4, 8, 12, 1, 5, 9, 13, 2, 6, 10, 14, 3, 7, 11, 15 }; DECLARE_ALIGNED(16, const int, vp9_row_scan_4x4[16]) = { 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15 }; DECLARE_ALIGNED(64, const int, vp9_default_zig_zag1d_8x8[64]) = { 0, 1, 8, 16, 9, 2, 3, 10, 17, 24, 32, 25, 18, 11, 4, 5, 12, 19, 26, 33, 40, 48, 41, 34, 27, 20, 13, 6, 7, 14, 21, 28, 35, 42, 49, 56, 57, 50, 43, 36, 29, 22, 15, 23, 30, 37, 44, 51, 58, 59, 52, 45, 38, 31, 39, 46, 53, 60, 61, 54, 47, 55, 62, 63, }; DECLARE_ALIGNED(16, const int, vp9_default_zig_zag1d_16x16[256]) = { 0, 1, 16, 32, 17, 2, 3, 18, 33, 48, 64, 49, 34, 19, 4, 5, 20, 35, 50, 65, 80, 96, 81, 66, 51, 36, 21, 6, 7, 22, 37, 52, 67, 82, 97, 112, 128, 113, 98, 83, 68, 53, 38, 23, 8, 9, 24, 39, 54, 69, 84, 99, 114, 129, 144, 160, 145, 130, 115, 100, 85, 70, 55, 40, 25, 10, 11, 26, 41, 56, 71, 86, 101, 116, 131, 146, 161, 176, 192, 177, 162, 147, 132, 117, 102, 87, 72, 57, 42, 27, 12, 13, 28, 43, 58, 73, 88, 103, 118, 133, 148, 163, 178, 193, 208, 224, 209, 194, 179, 164, 149, 134, 119, 104, 89, 74, 59, 44, 29, 14, 15, 30, 45, 60, 75, 90, 105, 120, 135, 150, 165, 180, 195, 210, 225, 240, 241, 226, 211, 196, 181, 166, 151, 136, 121, 106, 91, 76, 61, 46, 31, 47, 62, 77, 92, 107, 122, 137, 152, 167, 182, 197, 212, 227, 242, 243, 228, 213, 198, 183, 168, 153, 138, 123, 108, 93, 78, 63, 79, 94, 109, 124, 139, 154, 169, 184, 199, 214, 229, 244, 245, 230, 215, 200, 185, 170, 155, 140, 125, 110, 95, 111, 126, 141, 156, 171, 186, 201, 216, 231, 246, 247, 232, 217, 202, 187, 172, 157, 142, 127, 143, 158, 173, 188, 203, 218, 233, 248, 249, 234, 219, 204, 189, 174, 159, 175, 190, 205, 220, 235, 250, 251, 236, 221, 206, 191, 207, 222, 237, 252, 253, 238, 223, 239, 254, 255, }; DECLARE_ALIGNED(16, const int, vp9_default_zig_zag1d_32x32[1024]) = { 0, 1, 32, 64, 33, 2, 3, 34, 65, 96, 128, 97, 66, 35, 4, 5, 36, 67, 98, 129, 160, 192, 161, 130, 99, 68, 37, 6, 7, 38, 69, 100, 131, 162, 193, 224, 256, 225, 194, 163, 132, 101, 70, 39, 8, 9, 40, 71, 102, 133, 164, 195, 226, 257, 288, 320, 289, 258, 227, 196, 165, 134, 103, 72, 41, 10, 11, 42, 73, 104, 135, 166, 197, 228, 259, 290, 321, 352, 384, 353, 322, 291, 260, 229, 198, 167, 136, 105, 74, 43, 12, 13, 44, 75, 106, 137, 168, 199, 230, 261, 292, 323, 354, 385, 416, 448, 417, 386, 355, 324, 293, 262, 231, 200, 169, 138, 107, 76, 45, 14, 15, 46, 77, 108, 139, 170, 201, 232, 263, 294, 325, 356, 387, 418, 449, 480, 512, 481, 450, 419, 388, 357, 326, 295, 264, 233, 202, 171, 140, 109, 78, 47, 16, 17, 48, 79, 110, 141, 172, 203, 234, 265, 296, 327, 358, 389, 420, 451, 482, 513, 544, 576, 545, 514, 483, 452, 421, 390, 359, 328, 297, 266, 235, 204, 173, 142, 111, 80, 49, 18, 19, 50, 81, 112, 143, 174, 205, 236, 267, 298, 329, 360, 391, 422, 453, 484, 515, 546, 577, 608, 640, 609, 578, 547, 516, 485, 454, 423, 392, 361, 330, 299, 268, 237, 206, 175, 144, 113, 82, 51, 20, 21, 52, 83, 114, 145, 176, 207, 238, 269, 300, 331, 362, 393, 424, 455, 486, 517, 548, 579, 610, 641, 672, 704, 673, 642, 611, 580, 549, 518, 487, 456, 425, 394, 363, 332, 301, 270, 239, 208, 177, 146, 115, 84, 53, 22, 23, 54, 85, 116, 147, 178, 209, 240, 271, 302, 333, 364, 395, 426, 457, 488, 519, 550, 581, 612, 643, 674, 705, 736, 768, 737, 706, 675, 644, 613, 582, 551, 520, 489, 458, 427, 396, 365, 334, 303, 272, 241, 210, 179, 148, 117, 86, 55, 24, 25, 56, 87, 118, 149, 180, 211, 242, 273, 304, 335, 366, 397, 428, 459, 490, 521, 552, 583, 614, 645, 676, 707, 738, 769, 800, 832, 801, 770, 739, 708, 677, 646, 615, 584, 553, 522, 491, 460, 429, 398, 367, 336, 305, 274, 243, 212, 181, 150, 119, 88, 57, 26, 27, 58, 89, 120, 151, 182, 213, 244, 275, 306, 337, 368, 399, 430, 461, 492, 523, 554, 585, 616, 647, 678, 709, 740, 771, 802, 833, 864, 896, 865, 834, 803, 772, 741, 710, 679, 648, 617, 586, 555, 524, 493, 462, 431, 400, 369, 338, 307, 276, 245, 214, 183, 152, 121, 90, 59, 28, 29, 60, 91, 122, 153, 184, 215, 246, 277, 308, 339, 370, 401, 432, 463, 494, 525, 556, 587, 618, 649, 680, 711, 742, 773, 804, 835, 866, 897, 928, 960, 929, 898, 867, 836, 805, 774, 743, 712, 681, 650, 619, 588, 557, 526, 495, 464, 433, 402, 371, 340, 309, 278, 247, 216, 185, 154, 123, 92, 61, 30, 31, 62, 93, 124, 155, 186, 217, 248, 279, 310, 341, 372, 403, 434, 465, 496, 527, 558, 589, 620, 651, 682, 713, 744, 775, 806, 837, 868, 899, 930, 961, 992, 993, 962, 931, 900, 869, 838, 807, 776, 745, 714, 683, 652, 621, 590, 559, 528, 497, 466, 435, 404, 373, 342, 311, 280, 249, 218, 187, 156, 125, 94, 63, 95, 126, 157, 188, 219, 250, 281, 312, 343, 374, 405, 436, 467, 498, 529, 560, 591, 622, 653, 684, 715, 746, 777, 808, 839, 870, 901, 932, 963, 994, 995, 964, 933, 902, 871, 840, 809, 778, 747, 716, 685, 654, 623, 592, 561, 530, 499, 468, 437, 406, 375, 344, 313, 282, 251, 220, 189, 158, 127, 159, 190, 221, 252, 283, 314, 345, 376, 407, 438, 469, 500, 531, 562, 593, 624, 655, 686, 717, 748, 779, 810, 841, 872, 903, 934, 965, 996, 997, 966, 935, 904, 873, 842, 811, 780, 749, 718, 687, 656, 625, 594, 563, 532, 501, 470, 439, 408, 377, 346, 315, 284, 253, 222, 191, 223, 254, 285, 316, 347, 378, 409, 440, 471, 502, 533, 564, 595, 626, 657, 688, 719, 750, 781, 812, 843, 874, 905, 936, 967, 998, 999, 968, 937, 906, 875, 844, 813, 782, 751, 720, 689, 658, 627, 596, 565, 534, 503, 472, 441, 410, 379, 348, 317, 286, 255, 287, 318, 349, 380, 411, 442, 473, 504, 535, 566, 597, 628, 659, 690, 721, 752, 783, 814, 845, 876, 907, 938, 969, 1000, 1001, 970, 939, 908, 877, 846, 815, 784, 753, 722, 691, 660, 629, 598, 567, 536, 505, 474, 443, 412, 381, 350, 319, 351, 382, 413, 444, 475, 506, 537, 568, 599, 630, 661, 692, 723, 754, 785, 816, 847, 878, 909, 940, 971, 1002, 1003, 972, 941, 910, 879, 848, 817, 786, 755, 724, 693, 662, 631, 600, 569, 538, 507, 476, 445, 414, 383, 415, 446, 477, 508, 539, 570, 601, 632, 663, 694, 725, 756, 787, 818, 849, 880, 911, 942, 973, 1004, 1005, 974, 943, 912, 881, 850, 819, 788, 757, 726, 695, 664, 633, 602, 571, 540, 509, 478, 447, 479, 510, 541, 572, 603, 634, 665, 696, 727, 758, 789, 820, 851, 882, 913, 944, 975, 1006, 1007, 976, 945, 914, 883, 852, 821, 790, 759, 728, 697, 666, 635, 604, 573, 542, 511, 543, 574, 605, 636, 667, 698, 729, 760, 791, 822, 853, 884, 915, 946, 977, 1008, 1009, 978, 947, 916, 885, 854, 823, 792, 761, 730, 699, 668, 637, 606, 575, 607, 638, 669, 700, 731, 762, 793, 824, 855, 886, 917, 948, 979, 1010, 1011, 980, 949, 918, 887, 856, 825, 794, 763, 732, 701, 670, 639, 671, 702, 733, 764, 795, 826, 857, 888, 919, 950, 981, 1012, 1013, 982, 951, 920, 889, 858, 827, 796, 765, 734, 703, 735, 766, 797, 828, 859, 890, 921, 952, 983, 1014, 1015, 984, 953, 922, 891, 860, 829, 798, 767, 799, 830, 861, 892, 923, 954, 985, 1016, 1017, 986, 955, 924, 893, 862, 831, 863, 894, 925, 956, 987, 1018, 1019, 988, 957, 926, 895, 927, 958, 989, 1020, 1021, 990, 959, 991, 1022, 1023, }; /* Array indices are identical to previously-existing CONTEXT_NODE indices */ const vp9_tree_index vp9_coef_tree[ 22] = /* corresponding _CONTEXT_NODEs */ { -DCT_EOB_TOKEN, 2, /* 0 = EOB */ -ZERO_TOKEN, 4, /* 1 = ZERO */ -ONE_TOKEN, 6, /* 2 = ONE */ 8, 12, /* 3 = LOW_VAL */ -TWO_TOKEN, 10, /* 4 = TWO */ -THREE_TOKEN, -FOUR_TOKEN, /* 5 = THREE */ 14, 16, /* 6 = HIGH_LOW */ -DCT_VAL_CATEGORY1, -DCT_VAL_CATEGORY2, /* 7 = CAT_ONE */ 18, 20, /* 8 = CAT_THREEFOUR */ -DCT_VAL_CATEGORY3, -DCT_VAL_CATEGORY4, /* 9 = CAT_THREE */ -DCT_VAL_CATEGORY5, -DCT_VAL_CATEGORY6 /* 10 = CAT_FIVE */ }; struct vp9_token_struct vp9_coef_encodings[MAX_ENTROPY_TOKENS]; /* Trees for extra bits. Probabilities are constant and do not depend on previously encoded bits */ static const vp9_prob Pcat1[] = { 159}; static const vp9_prob Pcat2[] = { 165, 145}; static const vp9_prob Pcat3[] = { 173, 148, 140}; static const vp9_prob Pcat4[] = { 176, 155, 140, 135}; static const vp9_prob Pcat5[] = { 180, 157, 141, 134, 130}; static const vp9_prob Pcat6[] = { 254, 254, 254, 252, 249, 243, 230, 196, 177, 153, 140, 133, 130, 129 }; #if CONFIG_CODE_NONZEROCOUNT const vp9_tree_index vp9_nzc4x4_tree[2 * NZC4X4_NODES] = { -NZC_0, 2, 4, 6, -NZC_1, -NZC_2, -NZC_3TO4, 8, -NZC_5TO8, -NZC_9TO16, }; struct vp9_token_struct vp9_nzc4x4_encodings[NZC4X4_TOKENS]; const vp9_tree_index vp9_nzc8x8_tree[2 * NZC8X8_NODES] = { -NZC_0, 2, 4, 6, -NZC_1, -NZC_2, 8, 10, -NZC_3TO4, -NZC_5TO8, -NZC_9TO16, 12, -NZC_17TO32, -NZC_33TO64, }; struct vp9_token_struct vp9_nzc8x8_encodings[NZC8X8_TOKENS]; const vp9_tree_index vp9_nzc16x16_tree[2 * NZC16X16_NODES] = { -NZC_0, 2, 4, 6, -NZC_1, -NZC_2, 8, 10, -NZC_3TO4, -NZC_5TO8, 12, 14, -NZC_9TO16, -NZC_17TO32, -NZC_33TO64, 16, -NZC_65TO128, -NZC_129TO256, }; struct vp9_token_struct vp9_nzc16x16_encodings[NZC16X16_TOKENS]; const vp9_tree_index vp9_nzc32x32_tree[2 * NZC32X32_NODES] = { -NZC_0, 2, 4, 6, -NZC_1, -NZC_2, 8, 10, -NZC_3TO4, -NZC_5TO8, 12, 14, -NZC_9TO16, -NZC_17TO32, 16, 18, -NZC_33TO64, -NZC_65TO128, -NZC_129TO256, 20, -NZC_257TO512, -NZC_513TO1024, }; struct vp9_token_struct vp9_nzc32x32_encodings[NZC32X32_TOKENS]; const int vp9_extranzcbits[NZC32X32_TOKENS] = { 0, 0, 0, 1, 2, 3, 4, 5, 6, 7, 8, 9 }; const int vp9_basenzcvalue[NZC32X32_TOKENS] = { 0, 1, 2, 3, 5, 9, 17, 33, 65, 129, 257, 513 }; #endif // CONFIG_CODE_NONZEROCOUNT static vp9_tree_index cat1[2], cat2[4], cat3[6], cat4[8], cat5[10], cat6[28]; static void init_bit_tree(vp9_tree_index *p, int n) { int i = 0; while (++i < n) { p[0] = p[1] = i << 1; p += 2; } p[0] = p[1] = 0; } static void init_bit_trees() { init_bit_tree(cat1, 1); init_bit_tree(cat2, 2); init_bit_tree(cat3, 3); init_bit_tree(cat4, 4); init_bit_tree(cat5, 5); init_bit_tree(cat6, 14); } vp9_extra_bit_struct vp9_extra_bits[12] = { { 0, 0, 0, 0}, { 0, 0, 0, 1}, { 0, 0, 0, 2}, { 0, 0, 0, 3}, { 0, 0, 0, 4}, { cat1, Pcat1, 1, 5}, { cat2, Pcat2, 2, 7}, { cat3, Pcat3, 3, 11}, { cat4, Pcat4, 4, 19}, { cat5, Pcat5, 5, 35}, { cat6, Pcat6, 14, 67}, { 0, 0, 0, 0} }; #include "vp9/common/vp9_default_coef_probs.h" // This function updates and then returns n AC coefficient context // This is currently a placeholder function to allow experimentation // using various context models based on the energy earlier tokens // within the current block. // // For now it just returns the previously used context. int vp9_get_coef_context(int * recent_energy, int token) { // int token_energy; // int av_energy; /*token_energy = ((token != DCT_EOB_TOKEN) ? token : 0); if (!token_energy) { if (!(*recent_energy)) { av_energy = 0; } else { av_energy = 1; } } else { av_energy = ((token_energy + *recent_energy + 1) >> 1) + 1; if (av_energy > DCT_VAL_CATEGORY6) av_energy = DCT_VAL_CATEGORY6; } *recent_energy = token_energy;*/ return vp9_pt_energy_class[token]; }; void vp9_default_coef_probs(VP9_COMMON *pc) { #if CONFIG_CODE_NONZEROCOUNT #ifdef NZC_DEFAULT_COUNTS int h, g; for (h = 0; h < MAX_NZC_CONTEXTS; ++h) { for (g = 0; g < REF_TYPES; ++g) { int i; unsigned int branch_ct4x4[NZC4X4_NODES][2]; unsigned int branch_ct8x8[NZC8X8_NODES][2]; unsigned int branch_ct16x16[NZC16X16_NODES][2]; unsigned int branch_ct32x32[NZC32X32_NODES][2]; for (i = 0; i < BLOCK_TYPES; ++i) { vp9_tree_probs_from_distribution( vp9_nzc4x4_tree, pc->fc.nzc_probs_4x4[h][g][i], branch_ct4x4, default_nzc_counts_4x4[h][g][i], 0); } for (i = 0; i < BLOCK_TYPES; ++i) { vp9_tree_probs_from_distribution( vp9_nzc8x8_tree, pc->fc.nzc_probs_8x8[h][g][i], branch_ct8x8, default_nzc_counts_8x8[h][g][i], 0); } for (i = 0; i < BLOCK_TYPES; ++i) { vp9_tree_probs_from_distribution( vp9_nzc16x16_tree, pc->fc.nzc_probs_16x16[h][g][i], branch_ct16x16, default_nzc_counts_16x16[h][g][i], 0); } for (i = 0; i < BLOCK_TYPES; ++i) { vp9_tree_probs_from_distribution( vp9_nzc32x32_tree, pc->fc.nzc_probs_32x32[h][g][i], branch_ct32x32, default_nzc_counts_32x32[h][g][i], 0); } } } #else vpx_memcpy(pc->fc.nzc_probs_4x4, default_nzc_probs_4x4, sizeof(pc->fc.nzc_probs_4x4)); vpx_memcpy(pc->fc.nzc_probs_8x8, default_nzc_probs_8x8, sizeof(pc->fc.nzc_probs_8x8)); vpx_memcpy(pc->fc.nzc_probs_16x16, default_nzc_probs_16x16, sizeof(pc->fc.nzc_probs_16x16)); vpx_memcpy(pc->fc.nzc_probs_32x32, default_nzc_probs_32x32, sizeof(pc->fc.nzc_probs_32x32)); #endif vpx_memcpy(pc->fc.nzc_pcat_probs, default_nzc_pcat_probs, sizeof(pc->fc.nzc_pcat_probs)); #endif // CONFIG_CODE_NONZEROCOUNTyy vpx_memcpy(pc->fc.coef_probs_4x4, default_coef_probs_4x4, sizeof(pc->fc.coef_probs_4x4)); vpx_memcpy(pc->fc.coef_probs_8x8, default_coef_probs_8x8, sizeof(pc->fc.coef_probs_8x8)); vpx_memcpy(pc->fc.coef_probs_16x16, default_coef_probs_16x16, sizeof(pc->fc.coef_probs_16x16)); vpx_memcpy(pc->fc.coef_probs_32x32, default_coef_probs_32x32, sizeof(pc->fc.coef_probs_32x32)); } void vp9_coef_tree_initialize() { init_bit_trees(); vp9_tokens_from_tree(vp9_coef_encodings, vp9_coef_tree); #if CONFIG_CODE_NONZEROCOUNT vp9_tokens_from_tree(vp9_nzc4x4_encodings, vp9_nzc4x4_tree); vp9_tokens_from_tree(vp9_nzc8x8_encodings, vp9_nzc8x8_tree); vp9_tokens_from_tree(vp9_nzc16x16_encodings, vp9_nzc16x16_tree); vp9_tokens_from_tree(vp9_nzc32x32_encodings, vp9_nzc32x32_tree); #endif } #if CONFIG_CODE_NONZEROCOUNT #define mb_in_cur_tile(cm, mb_row, mb_col) \ ((mb_col) >= (cm)->cur_tile_mb_col_start && \ (mb_col) <= (cm)->cur_tile_mb_col_end && \ (mb_row) >= 0) #define choose_nzc_context(nzc_exp, t2, t1) \ ((nzc_exp) >= (t2) ? 2 : (nzc_exp) >= (t1) ? 1 : 0) #define NZC_T2_32X32 (16 << 6) #define NZC_T1_32X32 (4 << 6) #define NZC_T2_16X16 (12 << 6) #define NZC_T1_16X16 (3 << 6) #define NZC_T2_8X8 (8 << 6) #define NZC_T1_8X8 (2 << 6) #define NZC_T2_4X4 (4 << 6) #define NZC_T1_4X4 (1 << 6) // Transforms a mb16 block index to a sb64 block index static inline int mb16_to_sb64_index(int mb_row, int mb_col, int block) { int r = (mb_row & 3); int c = (mb_col & 3); int b; if (block < 16) { // Y int ib = block >> 2; int jb = block & 3; ib += r * 4; jb += c * 4; b = ib * 16 + jb; assert(b < 256); return b; } else { // UV int base = block - (block & 3); int ib = (block - base) >> 1; int jb = (block - base) & 1; ib += r * 2; jb += c * 2; b = base * 16 + ib * 8 + jb; assert(b >= 256 && b < 384); return b; } } // Transforms a mb16 block index to a sb32 block index static inline int mb16_to_sb32_index(int mb_row, int mb_col, int block) { int r = (mb_row & 1); int c = (mb_col & 1); int b; if (block < 16) { // Y int ib = block >> 2; int jb = block & 3; ib += r * 4; jb += c * 4; b = ib * 8 + jb; assert(b < 64); return b; } else { // UV int base = block - (block & 3); int ib = (block - base) >> 1; int jb = (block - base) & 1; ib += r * 2; jb += c * 2; b = base * 4 + ib * 4 + jb; assert(b >= 64 && b < 96); return b; } } static inline int block_to_txfm_index(int block, TX_SIZE tx_size, int s) { // s is the log of the number of 4x4 blocks in each row/col of larger block int b, ib, jb, nb; ib = block >> s; jb = block - (ib << s); ib >>= tx_size; jb >>= tx_size; nb = 1 << (s - tx_size); b = (ib * nb + jb) << (2 * tx_size); return b; } /* BEGIN - Helper functions to get the y nzcs */ static unsigned int get_nzc_4x4_y_sb64(MB_MODE_INFO *mi, int block) { int b; assert(block < 256); b = block_to_txfm_index(block, mi->txfm_size, 4); assert(b < 256); return mi->nzcs[b] << (6 - 2 * mi->txfm_size); } static unsigned int get_nzc_4x4_y_sb32(MB_MODE_INFO *mi, int block) { int b; assert(block < 64); b = block_to_txfm_index(block, mi->txfm_size, 3); assert(b < 64); return mi->nzcs[b] << (6 - 2 * mi->txfm_size); } static unsigned int get_nzc_4x4_y_mb16(MB_MODE_INFO *mi, int block) { int b; assert(block < 16); b = block_to_txfm_index(block, mi->txfm_size, 2); assert(b < 16); return mi->nzcs[b] << (6 - 2 * mi->txfm_size); } /* END - Helper functions to get the y nzcs */ /* Function to get y nzc where block index is in mb16 terms */ static unsigned int get_nzc_4x4_y(VP9_COMMON *cm, MODE_INFO *m, int mb_row, int mb_col, int block) { // NOTE: All values returned are at 64 times the true value at 4x4 scale MB_MODE_INFO *const mi = &m->mbmi; const int mis = cm->mode_info_stride; if (mi->mb_skip_coeff || !mb_in_cur_tile(cm, mb_row, mb_col)) return 0; if (mi->sb_type == BLOCK_SIZE_SB64X64) { int r = mb_row & 3; int c = mb_col & 3; m -= c + r * mis; if (m->mbmi.mb_skip_coeff || !mb_in_cur_tile(cm, mb_row - r, mb_col - c)) return 0; else return get_nzc_4x4_y_sb64( &m->mbmi, mb16_to_sb64_index(mb_row, mb_col, block)); } else if (mi->sb_type == BLOCK_SIZE_SB32X32) { int r = mb_row & 1; int c = mb_col & 1; m -= c + r * mis; if (m->mbmi.mb_skip_coeff || !mb_in_cur_tile(cm, mb_row - r, mb_col - c)) return 0; else return get_nzc_4x4_y_sb32( &m->mbmi, mb16_to_sb32_index(mb_row, mb_col, block)); } else { if (m->mbmi.mb_skip_coeff || !mb_in_cur_tile(cm, mb_row, mb_col)) return 0; return get_nzc_4x4_y_mb16(mi, block); } } /* BEGIN - Helper functions to get the uv nzcs */ static unsigned int get_nzc_4x4_uv_sb64(MB_MODE_INFO *mi, int block) { int b; int base, uvtxfm_size; assert(block >= 256 && block < 384); uvtxfm_size = mi->txfm_size; base = 256 + (block & 64); block -= base; b = base + block_to_txfm_index(block, uvtxfm_size, 3); assert(b >= 256 && b < 384); return mi->nzcs[b] << (6 - 2 * uvtxfm_size); } static unsigned int get_nzc_4x4_uv_sb32(MB_MODE_INFO *mi, int block) { int b; int base, uvtxfm_size; assert(block >= 64 && block < 96); if (mi->txfm_size == TX_32X32) uvtxfm_size = TX_16X16; else uvtxfm_size = mi->txfm_size; base = 64 + (block & 16); block -= base; b = base + block_to_txfm_index(block, uvtxfm_size, 2); assert(b >= 64 && b < 96); return mi->nzcs[b] << (6 - 2 * uvtxfm_size); } static unsigned int get_nzc_4x4_uv_mb16(MB_MODE_INFO *mi, int block) { int b; int base, uvtxfm_size; assert(block >= 16 && block < 24); if (mi->txfm_size == TX_8X8 && (mi->mode == SPLITMV || mi->mode == I8X8_PRED)) uvtxfm_size = TX_4X4; else if (mi->txfm_size == TX_16X16) uvtxfm_size = TX_8X8; else uvtxfm_size = mi->txfm_size; base = 16 + (block & 4); block -= base; b = base + block_to_txfm_index(block, uvtxfm_size, 1); assert(b >= 16 && b < 24); return mi->nzcs[b] << (6 - 2 * uvtxfm_size); } /* END - Helper functions to get the uv nzcs */ /* Function to get uv nzc where block index is in mb16 terms */ static unsigned int get_nzc_4x4_uv(VP9_COMMON *cm, MODE_INFO *m, int mb_row, int mb_col, int block) { // NOTE: All values returned are at 64 times the true value at 4x4 scale MB_MODE_INFO *const mi = &m->mbmi; const int mis = cm->mode_info_stride; if (mi->mb_skip_coeff || !mb_in_cur_tile(cm, mb_row, mb_col)) return 0; if (mi->sb_type == BLOCK_SIZE_SB64X64) { int r = mb_row & 3; int c = mb_col & 3; m -= c + r * mis; if (m->mbmi.mb_skip_coeff || !mb_in_cur_tile(cm, mb_row - r, mb_col - c)) return 0; else return get_nzc_4x4_uv_sb64( &m->mbmi, mb16_to_sb64_index(mb_row, mb_col, block)); } else if (mi->sb_type == BLOCK_SIZE_SB32X32) { int r = mb_row & 1; int c = mb_col & 1; m -= c + r * mis; if (m->mbmi.mb_skip_coeff || !mb_in_cur_tile(cm, mb_row - r, mb_col - c)) return 0; else return get_nzc_4x4_uv_sb32( &m->mbmi, mb16_to_sb32_index(mb_row, mb_col, block)); } else { return get_nzc_4x4_uv_mb16(mi, block); } } int vp9_get_nzc_context_y_sb64(VP9_COMMON *cm, MODE_INFO *cur, int mb_row, int mb_col, int block) { // returns an index in [0, MAX_NZC_CONTEXTS - 1] to reflect how busy // neighboring blocks are int mis = cm->mode_info_stride; int nzc_exp = 0; TX_SIZE txfm_size = cur->mbmi.txfm_size; assert(block < 256); switch (txfm_size) { case TX_32X32: assert((block & 63) == 0); if (block < 128) { int o = (block >> 6) * 2; nzc_exp = get_nzc_4x4_y(cm, cur - mis + o, mb_row - 1, mb_col + o, 12) + get_nzc_4x4_y(cm, cur - mis + o, mb_row - 1, mb_col + o, 13) + get_nzc_4x4_y(cm, cur - mis + o, mb_row - 1, mb_col + o, 14) + get_nzc_4x4_y(cm, cur - mis + o, mb_row - 1, mb_col + o, 15) + get_nzc_4x4_y(cm, cur - mis + o + 1, mb_row - 1, mb_col + o + 1, 12) + get_nzc_4x4_y(cm, cur - mis + o + 1, mb_row - 1, mb_col + o + 1, 13) + get_nzc_4x4_y(cm, cur - mis + o + 1, mb_row - 1, mb_col + o + 1, 14) + get_nzc_4x4_y(cm, cur - mis + o + 1, mb_row - 1, mb_col + o + 1, 15); } else { nzc_exp = cur->mbmi.nzcs[block - 128] << 3; } if ((block & 127) == 0) { int o = (block >> 7) * 2; nzc_exp += get_nzc_4x4_y(cm, cur - 1 + o * mis, mb_row + o, mb_col - 1, 3) + get_nzc_4x4_y(cm, cur - 1 + o * mis, mb_row + o, mb_col - 1, 7) + get_nzc_4x4_y(cm, cur - 1 + o * mis, mb_row + o, mb_col - 1, 11) + get_nzc_4x4_y(cm, cur - 1 + o * mis, mb_row + o, mb_col - 1, 15) + get_nzc_4x4_y(cm, cur - 1 + o * mis + mis, mb_row + o + 1, mb_col - 1, 3) + get_nzc_4x4_y(cm, cur - 1 + o * mis + mis, mb_row + o + 1, mb_col - 1, 7) + get_nzc_4x4_y(cm, cur - 1 + o * mis + mis, mb_row + o + 1, mb_col - 1, 11) + get_nzc_4x4_y(cm, cur - 1 + o * mis + mis, mb_row + o + 1, mb_col - 1, 15); } else { nzc_exp += cur->mbmi.nzcs[block - 64] << 3; } nzc_exp <<= 2; // Note nzc_exp is 64 times the average value expected at 32x32 scale return choose_nzc_context(nzc_exp, NZC_T2_32X32, NZC_T1_32X32); break; case TX_16X16: assert((block & 15) == 0); if (block < 64) { int o = block >> 4; nzc_exp = get_nzc_4x4_y(cm, cur - mis + o, mb_row - 1, mb_col + o, 12) + get_nzc_4x4_y(cm, cur - mis + o, mb_row - 1, mb_col + o, 13) + get_nzc_4x4_y(cm, cur - mis + o, mb_row - 1, mb_col + o, 14) + get_nzc_4x4_y(cm, cur - mis + o, mb_row - 1, mb_col + o, 15); } else { nzc_exp = cur->mbmi.nzcs[block - 64] << 4; } if ((block & 63) == 0) { int o = block >> 6; nzc_exp += get_nzc_4x4_y(cm, cur - 1 + o * mis, mb_row + o, mb_col - 1, 3) + get_nzc_4x4_y(cm, cur - 1 + o * mis, mb_row + o, mb_col - 1, 7) + get_nzc_4x4_y(cm, cur - 1 + o * mis, mb_row + o, mb_col - 1, 11) + get_nzc_4x4_y(cm, cur - 1 + o * mis, mb_row + o, mb_col - 1, 15); } else { nzc_exp += cur->mbmi.nzcs[block - 16] << 4; } nzc_exp <<= 1; // Note nzc_exp is 64 times the average value expected at 16x16 scale return choose_nzc_context(nzc_exp, NZC_T2_16X16, NZC_T1_16X16); break; case TX_8X8: assert((block & 3) == 0); if (block < 32) { int o = block >> 3; int p = ((block >> 2) & 1) ? 14 : 12; nzc_exp = get_nzc_4x4_y(cm, cur - mis + o, mb_row - 1, mb_col + o, p) + get_nzc_4x4_y(cm, cur - mis + o, mb_row - 1, mb_col + o, p + 1); } else { nzc_exp = cur->mbmi.nzcs[block - 32] << 5; } if ((block & 31) == 0) { int o = block >> 6; int p = ((block >> 5) & 1) ? 11 : 3; nzc_exp += get_nzc_4x4_y(cm, cur - 1 + o * mis, mb_row + o, mb_col - 1, p) + get_nzc_4x4_y(cm, cur - 1 + o * mis, mb_row + o, mb_col - 1, p + 4); } else { nzc_exp += cur->mbmi.nzcs[block - 4] << 5; } // Note nzc_exp is 64 times the average value expected at 8x8 scale return choose_nzc_context(nzc_exp, NZC_T2_8X8, NZC_T1_8X8); break; case TX_4X4: if (block < 16) { int o = block >> 2; int p = block & 3; nzc_exp = get_nzc_4x4_y(cm, cur - mis + o, mb_row - 1, mb_col + o, 12 + p); } else { nzc_exp = (cur->mbmi.nzcs[block - 16] << 6); } if ((block & 15) == 0) { int o = block >> 6; int p = (block >> 4) & 3; nzc_exp += get_nzc_4x4_y(cm, cur - 1 + o * mis, mb_row + o, mb_col - 1, 3 + 4 * p); } else { nzc_exp += (cur->mbmi.nzcs[block - 1] << 6); } nzc_exp >>= 1; // Note nzc_exp is 64 times the average value expected at 4x4 scale return choose_nzc_context(nzc_exp, NZC_T2_4X4, NZC_T1_4X4); break; default: return 0; } } int vp9_get_nzc_context_y_sb32(VP9_COMMON *cm, MODE_INFO *cur, int mb_row, int mb_col, int block) { // returns an index in [0, MAX_NZC_CONTEXTS - 1] to reflect how busy // neighboring blocks are int mis = cm->mode_info_stride; int nzc_exp = 0; TX_SIZE txfm_size = cur->mbmi.txfm_size; assert(block < 64); switch (txfm_size) { case TX_32X32: assert(block == 0); nzc_exp = (get_nzc_4x4_y(cm, cur - mis, mb_row - 1, mb_col, 12) + get_nzc_4x4_y(cm, cur - mis, mb_row - 1, mb_col, 13) + get_nzc_4x4_y(cm, cur - mis, mb_row - 1, mb_col, 14) + get_nzc_4x4_y(cm, cur - mis, mb_row - 1, mb_col, 15) + get_nzc_4x4_y(cm, cur - mis + 1, mb_row - 1, mb_col + 1, 12) + get_nzc_4x4_y(cm, cur - mis + 1, mb_row - 1, mb_col + 1, 13) + get_nzc_4x4_y(cm, cur - mis + 1, mb_row - 1, mb_col + 1, 14) + get_nzc_4x4_y(cm, cur - mis + 1, mb_row - 1, mb_col + 1, 15) + get_nzc_4x4_y(cm, cur - 1, mb_row, mb_col - 1, 3) + get_nzc_4x4_y(cm, cur - 1, mb_row, mb_col - 1, 7) + get_nzc_4x4_y(cm, cur - 1, mb_row, mb_col - 1, 11) + get_nzc_4x4_y(cm, cur - 1, mb_row, mb_col - 1, 15) + get_nzc_4x4_y(cm, cur - 1 + mis, mb_row + 1, mb_col - 1, 3) + get_nzc_4x4_y(cm, cur - 1 + mis, mb_row + 1, mb_col - 1, 7) + get_nzc_4x4_y(cm, cur - 1 + mis, mb_row + 1, mb_col - 1, 11) + get_nzc_4x4_y(cm, cur - 1 + mis, mb_row + 1, mb_col - 1, 15)) << 2; // Note nzc_exp is 64 times the average value expected at 32x32 scale return choose_nzc_context(nzc_exp, NZC_T2_32X32, NZC_T1_32X32); break; case TX_16X16: assert((block & 15) == 0); if (block < 32) { int o = (block >> 4) & 1; nzc_exp = get_nzc_4x4_y(cm, cur - mis + o, mb_row - 1, mb_col + o, 12) + get_nzc_4x4_y(cm, cur - mis + o, mb_row - 1, mb_col + o, 13) + get_nzc_4x4_y(cm, cur - mis + o, mb_row - 1, mb_col + o, 14) + get_nzc_4x4_y(cm, cur - mis + o, mb_row - 1, mb_col + o, 15); } else { nzc_exp = cur->mbmi.nzcs[block - 32] << 4; } if ((block & 31) == 0) { int o = block >> 5; nzc_exp += get_nzc_4x4_y(cm, cur - 1 + o * mis, mb_row + o, mb_col - 1, 3) + get_nzc_4x4_y(cm, cur - 1 + o * mis, mb_row + o, mb_col - 1, 7) + get_nzc_4x4_y(cm, cur - 1 + o * mis, mb_row + o, mb_col - 1, 11) + get_nzc_4x4_y(cm, cur - 1 + o * mis, mb_row + o, mb_col - 1, 15); } else { nzc_exp += cur->mbmi.nzcs[block - 16] << 4; } nzc_exp <<= 1; // Note nzc_exp is 64 times the average value expected at 16x16 scale return choose_nzc_context(nzc_exp, NZC_T2_16X16, NZC_T1_16X16); break; case TX_8X8: assert((block & 3) == 0); if (block < 16) { int o = block >> 3; int p = ((block >> 2) & 1) ? 14 : 12; nzc_exp = get_nzc_4x4_y(cm, cur - mis + o, mb_row - 1, mb_col + o, p) + get_nzc_4x4_y(cm, cur - mis + o, mb_row - 1, mb_col + o, p + 1); } else { nzc_exp = cur->mbmi.nzcs[block - 16] << 5; } if ((block & 15) == 0) { int o = block >> 5; int p = ((block >> 4) & 1) ? 11 : 3; nzc_exp += get_nzc_4x4_y(cm, cur - 1 + o * mis, mb_row + o, mb_col - 1, p) + get_nzc_4x4_y(cm, cur - 1 + o * mis, mb_row + o, mb_col - 1, p + 4); } else { nzc_exp += cur->mbmi.nzcs[block - 4] << 5; } // Note nzc_exp is 64 times the average value expected at 8x8 scale return choose_nzc_context(nzc_exp, NZC_T2_8X8, NZC_T1_8X8); break; case TX_4X4: if (block < 8) { int o = block >> 2; int p = block & 3; nzc_exp = get_nzc_4x4_y(cm, cur - mis + o, mb_row - 1, mb_col + o, 12 + p); } else { nzc_exp = (cur->mbmi.nzcs[block - 8] << 6); } if ((block & 7) == 0) { int o = block >> 5; int p = (block >> 3) & 3; nzc_exp += get_nzc_4x4_y(cm, cur - 1 + o * mis, mb_row + o, mb_col - 1, 3 + 4 * p); } else { nzc_exp += (cur->mbmi.nzcs[block - 1] << 6); } nzc_exp >>= 1; // Note nzc_exp is 64 times the average value expected at 4x4 scale return choose_nzc_context(nzc_exp, NZC_T2_4X4, NZC_T1_4X4); break; default: return 0; break; } } int vp9_get_nzc_context_y_mb16(VP9_COMMON *cm, MODE_INFO *cur, int mb_row, int mb_col, int block) { // returns an index in [0, MAX_NZC_CONTEXTS - 1] to reflect how busy // neighboring blocks are int mis = cm->mode_info_stride; int nzc_exp = 0; TX_SIZE txfm_size = cur->mbmi.txfm_size; assert(block < 16); switch (txfm_size) { case TX_16X16: assert(block == 0); nzc_exp = get_nzc_4x4_y(cm, cur - mis, mb_row - 1, mb_col, 12) + get_nzc_4x4_y(cm, cur - mis, mb_row - 1, mb_col, 13) + get_nzc_4x4_y(cm, cur - mis, mb_row - 1, mb_col, 14) + get_nzc_4x4_y(cm, cur - mis, mb_row - 1, mb_col, 15) + get_nzc_4x4_y(cm, cur - 1, mb_row, mb_col - 1, 3) + get_nzc_4x4_y(cm, cur - 1, mb_row, mb_col - 1, 7) + get_nzc_4x4_y(cm, cur - 1, mb_row, mb_col - 1, 11) + get_nzc_4x4_y(cm, cur - 1, mb_row, mb_col - 1, 15); nzc_exp <<= 1; // Note nzc_exp is 64 times the average value expected at 16x16 scale return choose_nzc_context(nzc_exp, NZC_T2_16X16, NZC_T1_16X16); case TX_8X8: assert((block & 3) == 0); if (block < 8) { int p = ((block >> 2) & 1) ? 14 : 12; nzc_exp = get_nzc_4x4_y(cm, cur - mis, mb_row - 1, mb_col, p) + get_nzc_4x4_y(cm, cur - mis, mb_row - 1, mb_col, p + 1); } else { nzc_exp = cur->mbmi.nzcs[block - 8] << 5; } if ((block & 7) == 0) { int p = ((block >> 3) & 1) ? 11 : 3; nzc_exp += get_nzc_4x4_y(cm, cur - 1, mb_row, mb_col - 1, p) + get_nzc_4x4_y(cm, cur - 1, mb_row, mb_col - 1, p + 4); } else { nzc_exp += cur->mbmi.nzcs[block - 4] << 5; } // Note nzc_exp is 64 times the average value expected at 8x8 scale return choose_nzc_context(nzc_exp, NZC_T2_8X8, NZC_T1_8X8); case TX_4X4: if (block < 4) { int p = block & 3; nzc_exp = get_nzc_4x4_y(cm, cur - mis, mb_row - 1, mb_col, 12 + p); } else { nzc_exp = (cur->mbmi.nzcs[block - 4] << 6); } if ((block & 3) == 0) { int p = (block >> 2) & 3; nzc_exp += get_nzc_4x4_y(cm, cur - 1, mb_row, mb_col - 1, 3 + 4 * p); } else { nzc_exp += (cur->mbmi.nzcs[block - 1] << 6); } nzc_exp >>= 1; // Note nzc_exp is 64 times the average value expected at 4x4 scale return choose_nzc_context(nzc_exp, NZC_T2_4X4, NZC_T1_4X4); default: return 0; break; } } int vp9_get_nzc_context_uv_sb64(VP9_COMMON *cm, MODE_INFO *cur, int mb_row, int mb_col, int block) { // returns an index in [0, MAX_NZC_CONTEXTS - 1] to reflect how busy // neighboring blocks are int mis = cm->mode_info_stride; int nzc_exp = 0; const int base = block - (block & 63); const int boff = (block & 63); const int base_mb16 = base >> 4; TX_SIZE txfm_size = cur->mbmi.txfm_size; TX_SIZE txfm_size_uv; assert(block >= 256 && block < 384); txfm_size_uv = txfm_size; switch (txfm_size_uv) { case TX_32X32: assert(block == 256 || block == 320); nzc_exp = get_nzc_4x4_uv(cm, cur - mis, mb_row - 1, mb_col, base_mb16 + 2) + get_nzc_4x4_uv(cm, cur - mis, mb_row - 1, mb_col, base_mb16 + 3) + get_nzc_4x4_uv(cm, cur - mis + 1, mb_row - 1, mb_col + 1, base_mb16 + 2) + get_nzc_4x4_uv(cm, cur - mis + 1, mb_row - 1, mb_col + 1, base_mb16 + 3) + get_nzc_4x4_uv(cm, cur - mis + 2, mb_row - 1, mb_col + 2, base_mb16 + 2) + get_nzc_4x4_uv(cm, cur - mis + 2, mb_row - 1, mb_col + 2, base_mb16 + 3) + get_nzc_4x4_uv(cm, cur - mis + 3, mb_row - 1, mb_col + 3, base_mb16 + 2) + get_nzc_4x4_uv(cm, cur - mis + 3, mb_row - 1, mb_col + 3, base_mb16 + 3) + get_nzc_4x4_uv(cm, cur - 1, mb_row, mb_col - 1, base_mb16 + 1) + get_nzc_4x4_uv(cm, cur - 1, mb_row, mb_col - 1, base_mb16 + 3) + get_nzc_4x4_uv(cm, cur - 1 + mis, mb_row + 1, mb_col - 1, base_mb16 + 1) + get_nzc_4x4_uv(cm, cur - 1 + mis, mb_row + 1, mb_col - 1, base_mb16 + 3) + get_nzc_4x4_uv(cm, cur - 1 + 2 * mis, mb_row + 2, mb_col - 1, base_mb16 + 1) + get_nzc_4x4_uv(cm, cur - 1 + 2 * mis, mb_row + 2, mb_col - 1, base_mb16 + 3) + get_nzc_4x4_uv(cm, cur - 1 + 3 * mis, mb_row + 3, mb_col - 1, base_mb16 + 1) + get_nzc_4x4_uv(cm, cur - 1 + 3 * mis, mb_row + 3, mb_col - 1, base_mb16 + 3); nzc_exp <<= 2; // Note nzc_exp is 64 times the average value expected at 32x32 scale return choose_nzc_context(nzc_exp, NZC_T2_32X32, NZC_T1_32X32); case TX_16X16: // uv txfm_size 16x16 assert((block & 15) == 0); if (boff < 32) { int o = (boff >> 4) & 1; nzc_exp = get_nzc_4x4_uv(cm, cur - mis + o, mb_row - 1, mb_col + o, base_mb16 + 2) + get_nzc_4x4_uv(cm, cur - mis + o, mb_row - 1, mb_col + o, base_mb16 + 3) + get_nzc_4x4_uv(cm, cur - mis + o + 1, mb_row - 1, mb_col + o + 1, base_mb16 + 2) + get_nzc_4x4_uv(cm, cur - mis + o + 1, mb_row - 1, mb_col + o + 1, base_mb16 + 3); } else { nzc_exp = cur->mbmi.nzcs[block - 32] << 4; } if ((boff & 31) == 0) { int o = boff >> 5; nzc_exp += get_nzc_4x4_uv(cm, cur - 1 + o * mis, mb_row + o, mb_col - 1, base_mb16 + 1) + get_nzc_4x4_uv(cm, cur - 1 + o * mis, mb_row + o, mb_col - 1, base_mb16 + 3) + get_nzc_4x4_uv(cm, cur - 1 + o * mis + mis, mb_row + o + 1, mb_col - 1, base_mb16 + 1) + get_nzc_4x4_uv(cm, cur - 1 + o * mis + mis, mb_row + o + 1, mb_col - 1, base_mb16 + 3); } else { nzc_exp += cur->mbmi.nzcs[block - 16] << 4; } nzc_exp <<= 1; // Note nzc_exp is 64 times the average value expected at 16x16 scale return choose_nzc_context(nzc_exp, NZC_T2_16X16, NZC_T1_16X16); case TX_8X8: assert((block & 3) == 0); if (boff < 16) { int o = boff >> 2; nzc_exp = get_nzc_4x4_uv(cm, cur - mis + o, mb_row - 1, mb_col + o, base_mb16 + 2) + get_nzc_4x4_uv(cm, cur - mis + o, mb_row - 1, mb_col + o, base_mb16 + 3); } else { nzc_exp = cur->mbmi.nzcs[block - 16] << 5; } if ((boff & 15) == 0) { int o = boff >> 4; nzc_exp += get_nzc_4x4_uv(cm, cur - 1 + o * mis, mb_row + o, mb_col - 1, base_mb16 + 1) + get_nzc_4x4_uv(cm, cur - 1 + o * mis, mb_row + o, mb_col - 1, base_mb16 + 3); } else { nzc_exp += cur->mbmi.nzcs[block - 4] << 5; } // Note nzc_exp is 64 times the average value expected at 8x8 scale return choose_nzc_context(nzc_exp, NZC_T2_8X8, NZC_T1_8X8); case TX_4X4: if (boff < 8) { int o = boff >> 1; int p = boff & 1; nzc_exp = get_nzc_4x4_uv(cm, cur - mis + o, mb_row - 1, mb_col + o, base_mb16 + 2 + p); } else { nzc_exp = (cur->mbmi.nzcs[block - 8] << 6); } if ((boff & 7) == 0) { int o = boff >> 4; int p = (boff >> 3) & 1; nzc_exp += get_nzc_4x4_uv(cm, cur - 1 + o * mis, mb_row + o, mb_col - 1, base_mb16 + 1 + 2 * p); } else { nzc_exp += (cur->mbmi.nzcs[block - 1] << 6); } nzc_exp >>= 1; // Note nzc_exp is 64 times the average value expected at 4x4 scale return choose_nzc_context(nzc_exp, NZC_T2_4X4, NZC_T1_4X4); default: return 0; } } int vp9_get_nzc_context_uv_sb32(VP9_COMMON *cm, MODE_INFO *cur, int mb_row, int mb_col, int block) { // returns an index in [0, MAX_NZC_CONTEXTS - 1] to reflect how busy // neighboring blocks are int mis = cm->mode_info_stride; int nzc_exp = 0; const int base = block - (block & 15); const int boff = (block & 15); const int base_mb16 = base >> 2; TX_SIZE txfm_size = cur->mbmi.txfm_size; TX_SIZE txfm_size_uv; assert(block >= 64 && block < 96); if (txfm_size == TX_32X32) txfm_size_uv = TX_16X16; else txfm_size_uv = txfm_size; switch (txfm_size_uv) { case TX_16X16: // uv txfm_size 16x16 assert(block == 64 || block == 80); nzc_exp = get_nzc_4x4_uv(cm, cur - mis, mb_row - 1, mb_col, base_mb16 + 2) + get_nzc_4x4_uv(cm, cur - mis, mb_row - 1, mb_col, base_mb16 + 3) + get_nzc_4x4_uv(cm, cur - mis + 1, mb_row - 1, mb_col + 1, base_mb16 + 2) + get_nzc_4x4_uv(cm, cur - mis + 1, mb_row - 1, mb_col + 1, base_mb16 + 3) + get_nzc_4x4_uv(cm, cur - 1 + mis, mb_row, mb_col - 1, base_mb16 + 1) + get_nzc_4x4_uv(cm, cur - 1 + mis, mb_row, mb_col - 1, base_mb16 + 3) + get_nzc_4x4_uv(cm, cur - 1 + mis, mb_row + 1, mb_col - 1, base_mb16 + 1) + get_nzc_4x4_uv(cm, cur - 1 + mis, mb_row + 1, mb_col - 1, base_mb16 + 3); nzc_exp <<= 1; // Note nzc_exp is 64 times the average value expected at 16x16 scale return choose_nzc_context(nzc_exp, NZC_T2_16X16, NZC_T1_16X16); break; case TX_8X8: assert((block & 3) == 0); if (boff < 8) { int o = boff >> 2; nzc_exp = get_nzc_4x4_uv(cm, cur - mis + o, mb_row - 1, mb_col + o, base_mb16 + 2) + get_nzc_4x4_uv(cm, cur - mis + o, mb_row - 1, mb_col + o, base_mb16 + 3); } else { nzc_exp = cur->mbmi.nzcs[block - 8] << 5; } if ((boff & 7) == 0) { int o = boff >> 3; nzc_exp += get_nzc_4x4_uv(cm, cur - 1 + o * mis, mb_row + o, mb_col - 1, base_mb16 + 1) + get_nzc_4x4_uv(cm, cur - 1 + o * mis, mb_row + o, mb_col - 1, base_mb16 + 3); } else { nzc_exp += cur->mbmi.nzcs[block - 4] << 5; } // Note nzc_exp is 64 times the average value expected at 8x8 scale return choose_nzc_context(nzc_exp, NZC_T2_8X8, NZC_T1_8X8); case TX_4X4: if (boff < 4) { int o = boff >> 1; int p = boff & 1; nzc_exp = get_nzc_4x4_uv(cm, cur - mis + o, mb_row - 1, mb_col + o, base_mb16 + 2 + p); } else { nzc_exp = (cur->mbmi.nzcs[block - 4] << 6); } if ((boff & 3) == 0) { int o = boff >> 3; int p = (boff >> 2) & 1; nzc_exp += get_nzc_4x4_uv(cm, cur - 1 + o * mis, mb_row + o, mb_col - 1, base_mb16 + 1 + 2 * p); } else { nzc_exp += (cur->mbmi.nzcs[block - 1] << 6); } nzc_exp >>= 1; // Note nzc_exp is 64 times the average value expected at 4x4 scale return choose_nzc_context(nzc_exp, NZC_T2_4X4, NZC_T1_4X4); default: return 0; } } int vp9_get_nzc_context_uv_mb16(VP9_COMMON *cm, MODE_INFO *cur, int mb_row, int mb_col, int block) { // returns an index in [0, MAX_NZC_CONTEXTS - 1] to reflect how busy // neighboring blocks are int mis = cm->mode_info_stride; int nzc_exp = 0; const int base = block - (block & 3); const int boff = (block & 3); const int base_mb16 = base; TX_SIZE txfm_size = cur->mbmi.txfm_size; TX_SIZE txfm_size_uv; assert(block >= 16 && block < 24); if (txfm_size == TX_16X16) txfm_size_uv = TX_8X8; else if (txfm_size == TX_8X8 && (cur->mbmi.mode == I8X8_PRED || cur->mbmi.mode == SPLITMV)) txfm_size_uv = TX_4X4; else txfm_size_uv = txfm_size; switch (txfm_size_uv) { case TX_8X8: assert((block & 3) == 0); nzc_exp = get_nzc_4x4_uv(cm, cur - mis, mb_row - 1, mb_col, base_mb16 + 2) + get_nzc_4x4_uv(cm, cur - mis, mb_row - 1, mb_col, base_mb16 + 3) + get_nzc_4x4_uv(cm, cur - 1, mb_row, mb_col - 1, base_mb16 + 1) + get_nzc_4x4_uv(cm, cur - 1, mb_row, mb_col - 1, base_mb16 + 3); // Note nzc_exp is 64 times the average value expected at 8x8 scale return choose_nzc_context(nzc_exp, NZC_T2_8X8, NZC_T1_8X8); case TX_4X4: if (boff < 2) { int p = boff & 1; nzc_exp = get_nzc_4x4_uv(cm, cur - mis, mb_row - 1, mb_col, base_mb16 + 2 + p); } else { nzc_exp = (cur->mbmi.nzcs[block - 2] << 6); } if ((boff & 1) == 0) { int p = (boff >> 1) & 1; nzc_exp += get_nzc_4x4_uv(cm, cur - 1, mb_row, mb_col - 1, base_mb16 + 1 + 2 * p); } else { nzc_exp += (cur->mbmi.nzcs[block - 1] << 6); } nzc_exp >>= 1; // Note nzc_exp is 64 times the average value expected at 4x4 scale return choose_nzc_context(nzc_exp, NZC_T2_4X4, NZC_T1_4X4); default: return 0; } } int vp9_get_nzc_context(VP9_COMMON *cm, MACROBLOCKD *xd, int block) { if (xd->mode_info_context->mbmi.sb_type == BLOCK_SIZE_SB64X64) { assert(block < 384); if (block < 256) return vp9_get_nzc_context_y_sb64(cm, xd->mode_info_context, get_mb_row(xd), get_mb_col(xd), block); else return vp9_get_nzc_context_uv_sb64(cm, xd->mode_info_context, get_mb_row(xd), get_mb_col(xd), block); } else if (xd->mode_info_context->mbmi.sb_type == BLOCK_SIZE_SB32X32) { assert(block < 96); if (block < 64) return vp9_get_nzc_context_y_sb32(cm, xd->mode_info_context, get_mb_row(xd), get_mb_col(xd), block); else return vp9_get_nzc_context_uv_sb32(cm, xd->mode_info_context, get_mb_row(xd), get_mb_col(xd), block); } else { assert(block < 64); if (block < 16) return vp9_get_nzc_context_y_mb16(cm, xd->mode_info_context, get_mb_row(xd), get_mb_col(xd), block); else return vp9_get_nzc_context_uv_mb16(cm, xd->mode_info_context, get_mb_row(xd), get_mb_col(xd), block); } } static void update_nzc(VP9_COMMON *cm, uint16_t nzc, int nzc_context, TX_SIZE tx_size, int ref, int type) { int e, c; c = codenzc(nzc); if (tx_size == TX_32X32) cm->fc.nzc_counts_32x32[nzc_context][ref][type][c]++; else if (tx_size == TX_16X16) cm->fc.nzc_counts_16x16[nzc_context][ref][type][c]++; else if (tx_size == TX_8X8) cm->fc.nzc_counts_8x8[nzc_context][ref][type][c]++; else if (tx_size == TX_4X4) cm->fc.nzc_counts_4x4[nzc_context][ref][type][c]++; else assert(0); if ((e = vp9_extranzcbits[c])) { int x = nzc - vp9_basenzcvalue[c]; while (e--) { int b = (x >> e) & 1; cm->fc.nzc_pcat_counts[nzc_context][c - NZC_TOKENS_NOEXTRA][e][b]++; } } } static void update_nzcs_sb64(VP9_COMMON *cm, MACROBLOCKD *xd, int mb_row, int mb_col) { MODE_INFO *m = xd->mode_info_context; MB_MODE_INFO *const mi = &m->mbmi; int j, nzc_context; const int ref = m->mbmi.ref_frame != INTRA_FRAME; assert(mb_col == get_mb_col(xd)); assert(mb_row == get_mb_row(xd)); if (mi->mb_skip_coeff) return; switch (mi->txfm_size) { case TX_32X32: for (j = 0; j < 256; j += 64) { nzc_context = vp9_get_nzc_context_y_sb64(cm, m, mb_row, mb_col, j); update_nzc(cm, m->mbmi.nzcs[j], nzc_context, TX_32X32, ref, 0); } for (j = 256; j < 384; j += 64) { nzc_context = vp9_get_nzc_context_uv_sb64(cm, m, mb_row, mb_col, j); update_nzc(cm, m->mbmi.nzcs[j], nzc_context, TX_32X32, ref, 1); } break; case TX_16X16: for (j = 0; j < 256; j += 16) { nzc_context = vp9_get_nzc_context_y_sb64(cm, m, mb_row, mb_col, j); update_nzc(cm, m->mbmi.nzcs[j], nzc_context, TX_16X16, ref, 0); } for (j = 256; j < 384; j += 16) { nzc_context = vp9_get_nzc_context_uv_sb64(cm, m, mb_row, mb_col, j); update_nzc(cm, m->mbmi.nzcs[j], nzc_context, TX_16X16, ref, 1); } break; case TX_8X8: for (j = 0; j < 256; j += 4) { nzc_context = vp9_get_nzc_context_y_sb64(cm, m, mb_row, mb_col, j); update_nzc(cm, m->mbmi.nzcs[j], nzc_context, TX_8X8, ref, 0); } for (j = 256; j < 384; j += 4) { nzc_context = vp9_get_nzc_context_uv_sb64(cm, m, mb_row, mb_col, j); update_nzc(cm, m->mbmi.nzcs[j], nzc_context, TX_8X8, ref, 1); } break; case TX_4X4: for (j = 0; j < 256; ++j) { nzc_context = vp9_get_nzc_context_y_sb64(cm, m, mb_row, mb_col, j); update_nzc(cm, m->mbmi.nzcs[j], nzc_context, TX_4X4, ref, 0); } for (j = 256; j < 384; ++j) { nzc_context = vp9_get_nzc_context_uv_sb64(cm, m, mb_row, mb_col, j); update_nzc(cm, m->mbmi.nzcs[j], nzc_context, TX_4X4, ref, 1); } break; default: break; } } static void update_nzcs_sb32(VP9_COMMON *cm, MACROBLOCKD *xd, int mb_row, int mb_col) { MODE_INFO *m = xd->mode_info_context; MB_MODE_INFO *const mi = &m->mbmi; int j, nzc_context; const int ref = m->mbmi.ref_frame != INTRA_FRAME; assert(mb_col == get_mb_col(xd)); assert(mb_row == get_mb_row(xd)); if (mi->mb_skip_coeff) return; switch (mi->txfm_size) { case TX_32X32: for (j = 0; j < 64; j += 64) { nzc_context = vp9_get_nzc_context_y_sb32(cm, m, mb_row, mb_col, j); update_nzc(cm, m->mbmi.nzcs[j], nzc_context, TX_32X32, ref, 0); } for (j = 64; j < 96; j += 16) { nzc_context = vp9_get_nzc_context_uv_sb32(cm, m, mb_row, mb_col, j); update_nzc(cm, m->mbmi.nzcs[j], nzc_context, TX_16X16, ref, 1); } break; case TX_16X16: for (j = 0; j < 64; j += 16) { nzc_context = vp9_get_nzc_context_y_sb32(cm, m, mb_row, mb_col, j); update_nzc(cm, m->mbmi.nzcs[j], nzc_context, TX_16X16, ref, 0); } for (j = 64; j < 96; j += 16) { nzc_context = vp9_get_nzc_context_uv_sb32(cm, m, mb_row, mb_col, j); update_nzc(cm, m->mbmi.nzcs[j], nzc_context, TX_16X16, ref, 1); } break; case TX_8X8: for (j = 0; j < 64; j += 4) { nzc_context = vp9_get_nzc_context_y_sb32(cm, m, mb_row, mb_col, j); update_nzc(cm, m->mbmi.nzcs[j], nzc_context, TX_8X8, ref, 0); } for (j = 64; j < 96; j += 4) { nzc_context = vp9_get_nzc_context_uv_sb32(cm, m, mb_row, mb_col, j); update_nzc(cm, m->mbmi.nzcs[j], nzc_context, TX_8X8, ref, 1); } break; case TX_4X4: for (j = 0; j < 64; ++j) { nzc_context = vp9_get_nzc_context_y_sb32(cm, m, mb_row, mb_col, j); update_nzc(cm, m->mbmi.nzcs[j], nzc_context, TX_4X4, ref, 0); } for (j = 64; j < 96; ++j) { nzc_context = vp9_get_nzc_context_uv_sb32(cm, m, mb_row, mb_col, j); update_nzc(cm, m->mbmi.nzcs[j], nzc_context, TX_4X4, ref, 1); } break; default: break; } } static void update_nzcs_mb16(VP9_COMMON *cm, MACROBLOCKD *xd, int mb_row, int mb_col) { MODE_INFO *m = xd->mode_info_context; MB_MODE_INFO *const mi = &m->mbmi; int j, nzc_context; const int ref = m->mbmi.ref_frame != INTRA_FRAME; assert(mb_col == get_mb_col(xd)); assert(mb_row == get_mb_row(xd)); if (mi->mb_skip_coeff) return; switch (mi->txfm_size) { case TX_16X16: for (j = 0; j < 16; j += 16) { nzc_context = vp9_get_nzc_context_y_mb16(cm, m, mb_row, mb_col, j); update_nzc(cm, m->mbmi.nzcs[j], nzc_context, TX_16X16, ref, 0); } for (j = 16; j < 24; j += 4) { nzc_context = vp9_get_nzc_context_uv_mb16(cm, m, mb_row, mb_col, j); update_nzc(cm, m->mbmi.nzcs[j], nzc_context, TX_8X8, ref, 1); } break; case TX_8X8: for (j = 0; j < 16; j += 4) { nzc_context = vp9_get_nzc_context_y_mb16(cm, m, mb_row, mb_col, j); update_nzc(cm, m->mbmi.nzcs[j], nzc_context, TX_8X8, ref, 0); } if (mi->mode == I8X8_PRED || mi->mode == SPLITMV) { for (j = 16; j < 24; ++j) { nzc_context = vp9_get_nzc_context_uv_mb16(cm, m, mb_row, mb_col, j); update_nzc(cm, m->mbmi.nzcs[j], nzc_context, TX_4X4, ref, 1); } } else { for (j = 16; j < 24; j += 4) { nzc_context = vp9_get_nzc_context_uv_mb16(cm, m, mb_row, mb_col, j); update_nzc(cm, m->mbmi.nzcs[j], nzc_context, TX_8X8, ref, 1); } } break; case TX_4X4: for (j = 0; j < 16; ++j) { nzc_context = vp9_get_nzc_context_y_mb16(cm, m, mb_row, mb_col, j); update_nzc(cm, m->mbmi.nzcs[j], nzc_context, TX_4X4, ref, 0); } for (j = 16; j < 24; ++j) { nzc_context = vp9_get_nzc_context_uv_mb16(cm, m, mb_row, mb_col, j); update_nzc(cm, m->mbmi.nzcs[j], nzc_context, TX_4X4, ref, 1); } break; default: break; } } void vp9_update_nzc_counts(VP9_COMMON *cm, MACROBLOCKD *xd, int mb_row, int mb_col) { if (xd->mode_info_context->mbmi.sb_type == BLOCK_SIZE_SB64X64) update_nzcs_sb64(cm, xd, mb_row, mb_col); else if (xd->mode_info_context->mbmi.sb_type == BLOCK_SIZE_SB32X32) update_nzcs_sb32(cm, xd, mb_row, mb_col); else update_nzcs_mb16(cm, xd, mb_row, mb_col); } #endif // CONFIG_CODE_NONZEROCOUNT // #define COEF_COUNT_TESTING #define COEF_COUNT_SAT 24 #define COEF_MAX_UPDATE_FACTOR 112 #define COEF_COUNT_SAT_KEY 24 #define COEF_MAX_UPDATE_FACTOR_KEY 112 #define COEF_COUNT_SAT_AFTER_KEY 24 #define COEF_MAX_UPDATE_FACTOR_AFTER_KEY 128 static void adapt_coef_probs(vp9_coeff_probs *dst_coef_probs, vp9_coeff_probs *pre_coef_probs, int block_types, vp9_coeff_count *coef_counts, int count_sat, int update_factor) { int t, i, j, k, l, count; unsigned int branch_ct[ENTROPY_NODES][2]; vp9_prob coef_probs[ENTROPY_NODES]; int factor; for (i = 0; i < block_types; ++i) for (j = 0; j < REF_TYPES; ++j) for (k = 0; k < COEF_BANDS; ++k) for (l = 0; l < PREV_COEF_CONTEXTS; ++l) { if (l >= 3 && k == 0) continue; vp9_tree_probs_from_distribution(vp9_coef_tree, coef_probs, branch_ct, coef_counts[i][j][k][l], 0); for (t = 0; t < ENTROPY_NODES; ++t) { count = branch_ct[t][0] + branch_ct[t][1]; count = count > count_sat ? count_sat : count; factor = (update_factor * count / count_sat); dst_coef_probs[i][j][k][l][t] = weighted_prob(pre_coef_probs[i][j][k][l][t], coef_probs[t], factor); } } } void vp9_adapt_coef_probs(VP9_COMMON *cm) { int count_sat; int update_factor; /* denominator 256 */ // printf("Frame type: %d\n", cm->frame_type); if (cm->frame_type == KEY_FRAME) { update_factor = COEF_MAX_UPDATE_FACTOR_KEY; count_sat = COEF_COUNT_SAT_KEY; } else if (cm->last_frame_type == KEY_FRAME) { update_factor = COEF_MAX_UPDATE_FACTOR_AFTER_KEY; /* adapt quickly */ count_sat = COEF_COUNT_SAT_AFTER_KEY; } else { update_factor = COEF_MAX_UPDATE_FACTOR; count_sat = COEF_COUNT_SAT; } adapt_coef_probs(cm->fc.coef_probs_4x4, cm->fc.pre_coef_probs_4x4, BLOCK_TYPES, cm->fc.coef_counts_4x4, count_sat, update_factor); adapt_coef_probs(cm->fc.coef_probs_8x8, cm->fc.pre_coef_probs_8x8, BLOCK_TYPES, cm->fc.coef_counts_8x8, count_sat, update_factor); adapt_coef_probs(cm->fc.coef_probs_16x16, cm->fc.pre_coef_probs_16x16, BLOCK_TYPES, cm->fc.coef_counts_16x16, count_sat, update_factor); adapt_coef_probs(cm->fc.coef_probs_32x32, cm->fc.pre_coef_probs_32x32, BLOCK_TYPES, cm->fc.coef_counts_32x32, count_sat, update_factor); } #if CONFIG_CODE_NONZEROCOUNT static void adapt_nzc_probs(VP9_COMMON *cm, int block_size, int count_sat, int update_factor) { int c, r, b, n; int count, factor; unsigned int nzc_branch_ct[NZC32X32_NODES][2]; vp9_prob nzc_probs[NZC32X32_NODES]; int tokens, nodes; const vp9_tree_index *nzc_tree; vp9_prob *dst_nzc_probs; vp9_prob *pre_nzc_probs; unsigned int *nzc_counts; if (block_size == 32) { tokens = NZC32X32_TOKENS; nzc_tree = vp9_nzc32x32_tree; dst_nzc_probs = cm->fc.nzc_probs_32x32[0][0][0]; pre_nzc_probs = cm->fc.pre_nzc_probs_32x32[0][0][0]; nzc_counts = cm->fc.nzc_counts_32x32[0][0][0]; } else if (block_size == 16) { tokens = NZC16X16_TOKENS; nzc_tree = vp9_nzc16x16_tree; dst_nzc_probs = cm->fc.nzc_probs_16x16[0][0][0]; pre_nzc_probs = cm->fc.pre_nzc_probs_16x16[0][0][0]; nzc_counts = cm->fc.nzc_counts_16x16[0][0][0]; } else if (block_size == 8) { tokens = NZC8X8_TOKENS; nzc_tree = vp9_nzc8x8_tree; dst_nzc_probs = cm->fc.nzc_probs_8x8[0][0][0]; pre_nzc_probs = cm->fc.pre_nzc_probs_8x8[0][0][0]; nzc_counts = cm->fc.nzc_counts_8x8[0][0][0]; } else { nzc_tree = vp9_nzc4x4_tree; tokens = NZC4X4_TOKENS; dst_nzc_probs = cm->fc.nzc_probs_4x4[0][0][0]; pre_nzc_probs = cm->fc.pre_nzc_probs_4x4[0][0][0]; nzc_counts = cm->fc.nzc_counts_4x4[0][0][0]; } nodes = tokens - 1; for (c = 0; c < MAX_NZC_CONTEXTS; ++c) for (r = 0; r < REF_TYPES; ++r) for (b = 0; b < BLOCK_TYPES; ++b) { int offset = c * REF_TYPES * BLOCK_TYPES + r * BLOCK_TYPES + b; int offset_nodes = offset * nodes; int offset_tokens = offset * tokens; vp9_tree_probs_from_distribution(nzc_tree, nzc_probs, nzc_branch_ct, nzc_counts + offset_tokens, 0); for (n = 0; n < nodes; ++n) { count = nzc_branch_ct[n][0] + nzc_branch_ct[n][1]; count = count > count_sat ? count_sat : count; factor = (update_factor * count / count_sat); dst_nzc_probs[offset_nodes + n] = weighted_prob(pre_nzc_probs[offset_nodes + n], nzc_probs[n], factor); } } } static void adapt_nzc_pcat(VP9_COMMON *cm, int count_sat, int update_factor) { int c, t; int count, factor; for (c = 0; c < MAX_NZC_CONTEXTS; ++c) { for (t = 0; t < NZC_TOKENS_EXTRA; ++t) { int bits = vp9_extranzcbits[t + NZC_TOKENS_NOEXTRA]; int b; for (b = 0; b < bits; ++b) { vp9_prob prob = get_binary_prob(cm->fc.nzc_pcat_counts[c][t][b][0], cm->fc.nzc_pcat_counts[c][t][b][1]); count = cm->fc.nzc_pcat_counts[c][t][b][0] + cm->fc.nzc_pcat_counts[c][t][b][1]; count = count > count_sat ? count_sat : count; factor = (update_factor * count / count_sat); cm->fc.nzc_pcat_probs[c][t][b] = weighted_prob( cm->fc.pre_nzc_pcat_probs[c][t][b], prob, factor); } } } } // #define NZC_COUNT_TESTING void vp9_adapt_nzc_probs(VP9_COMMON *cm) { int count_sat; int update_factor; /* denominator 256 */ #ifdef NZC_COUNT_TESTING int c, r, b, t; printf("\n"); for (c = 0; c < MAX_NZC_CONTEXTS; ++c) for (r = 0; r < REF_TYPES; ++r) { for (b = 0; b < BLOCK_TYPES; ++b) { printf(" {"); for (t = 0; t < NZC4X4_TOKENS; ++t) { printf(" %d,", cm->fc.nzc_counts_4x4[c][r][b][t]); } printf("}\n"); } printf("\n"); } #endif if (cm->frame_type == KEY_FRAME) { update_factor = COEF_MAX_UPDATE_FACTOR_KEY; count_sat = COEF_COUNT_SAT_KEY; } else if (cm->last_frame_type == KEY_FRAME) { update_factor = COEF_MAX_UPDATE_FACTOR_AFTER_KEY; /* adapt quickly */ count_sat = COEF_COUNT_SAT_AFTER_KEY; } else { update_factor = COEF_MAX_UPDATE_FACTOR; count_sat = COEF_COUNT_SAT; } adapt_nzc_probs(cm, 4, count_sat, update_factor); adapt_nzc_probs(cm, 8, count_sat, update_factor); adapt_nzc_probs(cm, 16, count_sat, update_factor); adapt_nzc_probs(cm, 32, count_sat, update_factor); adapt_nzc_pcat(cm, count_sat, update_factor); } #endif // CONFIG_CODE_NONZEROCOUNT