1 /* 2 * Copyright (c) 2010 The WebM project authors. All Rights Reserved. 3 * 4 * Use of this source code is governed by a BSD-style license 5 * that can be found in the LICENSE file in the root of the source 6 * tree. An additional intellectual property rights grant can be found 7 * in the file PATENTS. All contributing project authors may 8 * be found in the AUTHORS file in the root of the source tree. 9 */ 10 11 12 #include <math.h> 13 #include "vpx_mem/vpx_mem.h" 14 15 #include "onyx_int.h" 16 #include "quantize.h" 17 #include "vp8/common/quant_common.h" 18 19 #define EXACT_QUANT 20 21 #ifdef EXACT_FASTQUANT 22 void vp8_fast_quantize_b_c(BLOCK *b, BLOCKD *d) 23 { 24 int i, rc, eob; 25 int zbin; 26 int x, y, z, sz; 27 short *coeff_ptr = b->coeff; 28 short *zbin_ptr = b->zbin; 29 short *round_ptr = b->round; 30 short *quant_ptr = b->quant_fast; 31 unsigned char *quant_shift_ptr = b->quant_shift; 32 short *qcoeff_ptr = d->qcoeff; 33 short *dqcoeff_ptr = d->dqcoeff; 34 short *dequant_ptr = d->dequant; 35 36 vpx_memset(qcoeff_ptr, 0, 32); 37 vpx_memset(dqcoeff_ptr, 0, 32); 38 39 eob = -1; 40 41 for (i = 0; i < 16; i++) 42 { 43 rc = vp8_default_zig_zag1d[i]; 44 z = coeff_ptr[rc]; 45 zbin = zbin_ptr[rc] ; 46 47 sz = (z >> 31); /* sign of z */ 48 x = (z ^ sz) - sz; /* x = abs(z) */ 49 50 if (x >= zbin) 51 { 52 x += round_ptr[rc]; 53 y = ((((x * quant_ptr[rc]) >> 16) + x) 54 * quant_shift_ptr[rc]) >> 16; /* quantize (x) */ 55 x = (y ^ sz) - sz; /* get the sign back */ 56 qcoeff_ptr[rc] = x; /* write to destination */ 57 dqcoeff_ptr[rc] = x * dequant_ptr[rc]; /* dequantized value */ 58 59 if (y) 60 { 61 eob = i; /* last nonzero coeffs */ 62 } 63 } 64 } 65 *d->eob = (char)(eob + 1); 66 } 67 68 #else 69 70 void vp8_fast_quantize_b_c(BLOCK *b, BLOCKD *d) 71 { 72 int i, rc, eob; 73 int x, y, z, sz; 74 short *coeff_ptr = b->coeff; 75 short *round_ptr = b->round; 76 short *quant_ptr = b->quant_fast; 77 short *qcoeff_ptr = d->qcoeff; 78 short *dqcoeff_ptr = d->dqcoeff; 79 short *dequant_ptr = d->dequant; 80 81 eob = -1; 82 for (i = 0; i < 16; i++) 83 { 84 rc = vp8_default_zig_zag1d[i]; 85 z = coeff_ptr[rc]; 86 87 sz = (z >> 31); /* sign of z */ 88 x = (z ^ sz) - sz; /* x = abs(z) */ 89 90 y = ((x + round_ptr[rc]) * quant_ptr[rc]) >> 16; /* quantize (x) */ 91 x = (y ^ sz) - sz; /* get the sign back */ 92 qcoeff_ptr[rc] = x; /* write to destination */ 93 dqcoeff_ptr[rc] = x * dequant_ptr[rc]; /* dequantized value */ 94 95 if (y) 96 { 97 eob = i; /* last nonzero coeffs */ 98 } 99 } 100 *d->eob = (char)(eob + 1); 101 } 102 103 #endif 104 105 #ifdef EXACT_QUANT 106 void vp8_regular_quantize_b_c(BLOCK *b, BLOCKD *d) 107 { 108 int i, rc, eob; 109 int zbin; 110 int x, y, z, sz; 111 short *zbin_boost_ptr = b->zrun_zbin_boost; 112 short *coeff_ptr = b->coeff; 113 short *zbin_ptr = b->zbin; 114 short *round_ptr = b->round; 115 short *quant_ptr = b->quant; 116 short *quant_shift_ptr = b->quant_shift; 117 short *qcoeff_ptr = d->qcoeff; 118 short *dqcoeff_ptr = d->dqcoeff; 119 short *dequant_ptr = d->dequant; 120 short zbin_oq_value = b->zbin_extra; 121 122 vpx_memset(qcoeff_ptr, 0, 32); 123 vpx_memset(dqcoeff_ptr, 0, 32); 124 125 eob = -1; 126 127 for (i = 0; i < 16; i++) 128 { 129 rc = vp8_default_zig_zag1d[i]; 130 z = coeff_ptr[rc]; 131 132 zbin = zbin_ptr[rc] + *zbin_boost_ptr + zbin_oq_value; 133 134 zbin_boost_ptr ++; 135 sz = (z >> 31); /* sign of z */ 136 x = (z ^ sz) - sz; /* x = abs(z) */ 137 138 if (x >= zbin) 139 { 140 x += round_ptr[rc]; 141 y = ((((x * quant_ptr[rc]) >> 16) + x) 142 * quant_shift_ptr[rc]) >> 16; /* quantize (x) */ 143 x = (y ^ sz) - sz; /* get the sign back */ 144 qcoeff_ptr[rc] = x; /* write to destination */ 145 dqcoeff_ptr[rc] = x * dequant_ptr[rc]; /* dequantized value */ 146 147 if (y) 148 { 149 eob = i; /* last nonzero coeffs */ 150 zbin_boost_ptr = b->zrun_zbin_boost; /* reset zero runlength */ 151 } 152 } 153 } 154 155 *d->eob = (char)(eob + 1); 156 } 157 158 /* Perform regular quantization, with unbiased rounding and no zero bin. */ 159 void vp8_strict_quantize_b_c(BLOCK *b, BLOCKD *d) 160 { 161 int i; 162 int rc; 163 int eob; 164 int x; 165 int y; 166 int z; 167 int sz; 168 short *coeff_ptr; 169 short *quant_ptr; 170 short *quant_shift_ptr; 171 short *qcoeff_ptr; 172 short *dqcoeff_ptr; 173 short *dequant_ptr; 174 175 coeff_ptr = b->coeff; 176 quant_ptr = b->quant; 177 quant_shift_ptr = b->quant_shift; 178 qcoeff_ptr = d->qcoeff; 179 dqcoeff_ptr = d->dqcoeff; 180 dequant_ptr = d->dequant; 181 eob = - 1; 182 vpx_memset(qcoeff_ptr, 0, 32); 183 vpx_memset(dqcoeff_ptr, 0, 32); 184 for (i = 0; i < 16; i++) 185 { 186 int dq; 187 int rounding; 188 189 /*TODO: These arrays should be stored in zig-zag order.*/ 190 rc = vp8_default_zig_zag1d[i]; 191 z = coeff_ptr[rc]; 192 dq = dequant_ptr[rc]; 193 rounding = dq >> 1; 194 /* Sign of z. */ 195 sz = -(z < 0); 196 x = (z + sz) ^ sz; 197 x += rounding; 198 if (x >= dq) 199 { 200 /* Quantize x. */ 201 y = ((((x * quant_ptr[rc]) >> 16) + x) * quant_shift_ptr[rc]) >> 16; 202 /* Put the sign back. */ 203 x = (y + sz) ^ sz; 204 /* Save the coefficient and its dequantized value. */ 205 qcoeff_ptr[rc] = x; 206 dqcoeff_ptr[rc] = x * dq; 207 /* Remember the last non-zero coefficient. */ 208 if (y) 209 eob = i; 210 } 211 } 212 213 *d->eob = (char)(eob + 1); 214 } 215 216 #else 217 218 void vp8_regular_quantize_b_c(BLOCK *b, BLOCKD *d) 219 { 220 int i, rc, eob; 221 int zbin; 222 int x, y, z, sz; 223 short *zbin_boost_ptr = b->zrun_zbin_boost; 224 short *coeff_ptr = b->coeff; 225 short *zbin_ptr = b->zbin; 226 short *round_ptr = b->round; 227 short *quant_ptr = b->quant; 228 short *qcoeff_ptr = d->qcoeff; 229 short *dqcoeff_ptr = d->dqcoeff; 230 short *dequant_ptr = d->dequant; 231 short zbin_oq_value = b->zbin_extra; 232 233 vpx_memset(qcoeff_ptr, 0, 32); 234 vpx_memset(dqcoeff_ptr, 0, 32); 235 236 eob = -1; 237 238 for (i = 0; i < 16; i++) 239 { 240 rc = vp8_default_zig_zag1d[i]; 241 z = coeff_ptr[rc]; 242 243 zbin = zbin_ptr[rc] + *zbin_boost_ptr + zbin_oq_value; 244 245 zbin_boost_ptr ++; 246 sz = (z >> 31); /* sign of z */ 247 x = (z ^ sz) - sz; /* x = abs(z) */ 248 249 if (x >= zbin) 250 { 251 y = ((x + round_ptr[rc]) * quant_ptr[rc]) >> 16; /* quantize (x) */ 252 x = (y ^ sz) - sz; /* get the sign back */ 253 qcoeff_ptr[rc] = x; /* write to destination */ 254 dqcoeff_ptr[rc] = x * dequant_ptr[rc]; /* dequantized value */ 255 256 if (y) 257 { 258 eob = i; /* last nonzero coeffs */ 259 zbin_boost_ptr = &b->zrun_zbin_boost[0]; /* reset zrl */ 260 } 261 } 262 } 263 264 *d->eob = (char)(eob + 1); 265 } 266 267 #endif 268 269 void vp8_quantize_mby_c(MACROBLOCK *x) 270 { 271 int i; 272 int has_2nd_order = (x->e_mbd.mode_info_context->mbmi.mode != B_PRED 273 && x->e_mbd.mode_info_context->mbmi.mode != SPLITMV); 274 275 for (i = 0; i < 16; i++) 276 x->quantize_b(&x->block[i], &x->e_mbd.block[i]); 277 278 if(has_2nd_order) 279 x->quantize_b(&x->block[24], &x->e_mbd.block[24]); 280 } 281 282 void vp8_quantize_mb_c(MACROBLOCK *x) 283 { 284 int i; 285 int has_2nd_order=(x->e_mbd.mode_info_context->mbmi.mode != B_PRED 286 && x->e_mbd.mode_info_context->mbmi.mode != SPLITMV); 287 288 for (i = 0; i < 24+has_2nd_order; i++) 289 x->quantize_b(&x->block[i], &x->e_mbd.block[i]); 290 } 291 292 293 void vp8_quantize_mbuv_c(MACROBLOCK *x) 294 { 295 int i; 296 297 for (i = 16; i < 24; i++) 298 x->quantize_b(&x->block[i], &x->e_mbd.block[i]); 299 } 300 301 /* quantize_b_pair function pointer in MACROBLOCK structure is set to one of 302 * these two C functions if corresponding optimized routine is not available. 303 * NEON optimized version implements currently the fast quantization for pair 304 * of blocks. */ 305 void vp8_regular_quantize_b_pair(BLOCK *b1, BLOCK *b2, BLOCKD *d1, BLOCKD *d2) 306 { 307 vp8_regular_quantize_b(b1, d1); 308 vp8_regular_quantize_b(b2, d2); 309 } 310 311 void vp8_fast_quantize_b_pair_c(BLOCK *b1, BLOCK *b2, BLOCKD *d1, BLOCKD *d2) 312 { 313 vp8_fast_quantize_b_c(b1, d1); 314 vp8_fast_quantize_b_c(b2, d2); 315 } 316 317 318 static const int qrounding_factors[129] = 319 { 320 48, 48, 48, 48, 48, 48, 48, 48, 321 48, 48, 48, 48, 48, 48, 48, 48, 322 48, 48, 48, 48, 48, 48, 48, 48, 323 48, 48, 48, 48, 48, 48, 48, 48, 324 48, 48, 48, 48, 48, 48, 48, 48, 325 48, 48, 48, 48, 48, 48, 48, 48, 326 48, 48, 48, 48, 48, 48, 48, 48, 327 48, 48, 48, 48, 48, 48, 48, 48, 328 48, 48, 48, 48, 48, 48, 48, 48, 329 48, 48, 48, 48, 48, 48, 48, 48, 330 48, 48, 48, 48, 48, 48, 48, 48, 331 48, 48, 48, 48, 48, 48, 48, 48, 332 48, 48, 48, 48, 48, 48, 48, 48, 333 48, 48, 48, 48, 48, 48, 48, 48, 334 48, 48, 48, 48, 48, 48, 48, 48, 335 48, 48, 48, 48, 48, 48, 48, 48, 336 48 337 }; 338 339 340 static const int qzbin_factors[129] = 341 { 342 84, 84, 84, 84, 84, 84, 84, 84, 343 84, 84, 84, 84, 84, 84, 84, 84, 344 84, 84, 84, 84, 84, 84, 84, 84, 345 84, 84, 84, 84, 84, 84, 84, 84, 346 84, 84, 84, 84, 84, 84, 84, 84, 347 84, 84, 84, 84, 84, 84, 84, 84, 348 80, 80, 80, 80, 80, 80, 80, 80, 349 80, 80, 80, 80, 80, 80, 80, 80, 350 80, 80, 80, 80, 80, 80, 80, 80, 351 80, 80, 80, 80, 80, 80, 80, 80, 352 80, 80, 80, 80, 80, 80, 80, 80, 353 80, 80, 80, 80, 80, 80, 80, 80, 354 80, 80, 80, 80, 80, 80, 80, 80, 355 80, 80, 80, 80, 80, 80, 80, 80, 356 80, 80, 80, 80, 80, 80, 80, 80, 357 80, 80, 80, 80, 80, 80, 80, 80, 358 80 359 }; 360 361 362 static const int qrounding_factors_y2[129] = 363 { 364 48, 48, 48, 48, 48, 48, 48, 48, 365 48, 48, 48, 48, 48, 48, 48, 48, 366 48, 48, 48, 48, 48, 48, 48, 48, 367 48, 48, 48, 48, 48, 48, 48, 48, 368 48, 48, 48, 48, 48, 48, 48, 48, 369 48, 48, 48, 48, 48, 48, 48, 48, 370 48, 48, 48, 48, 48, 48, 48, 48, 371 48, 48, 48, 48, 48, 48, 48, 48, 372 48, 48, 48, 48, 48, 48, 48, 48, 373 48, 48, 48, 48, 48, 48, 48, 48, 374 48, 48, 48, 48, 48, 48, 48, 48, 375 48, 48, 48, 48, 48, 48, 48, 48, 376 48, 48, 48, 48, 48, 48, 48, 48, 377 48, 48, 48, 48, 48, 48, 48, 48, 378 48, 48, 48, 48, 48, 48, 48, 48, 379 48, 48, 48, 48, 48, 48, 48, 48, 380 48 381 }; 382 383 384 static const int qzbin_factors_y2[129] = 385 { 386 84, 84, 84, 84, 84, 84, 84, 84, 387 84, 84, 84, 84, 84, 84, 84, 84, 388 84, 84, 84, 84, 84, 84, 84, 84, 389 84, 84, 84, 84, 84, 84, 84, 84, 390 84, 84, 84, 84, 84, 84, 84, 84, 391 84, 84, 84, 84, 84, 84, 84, 84, 392 80, 80, 80, 80, 80, 80, 80, 80, 393 80, 80, 80, 80, 80, 80, 80, 80, 394 80, 80, 80, 80, 80, 80, 80, 80, 395 80, 80, 80, 80, 80, 80, 80, 80, 396 80, 80, 80, 80, 80, 80, 80, 80, 397 80, 80, 80, 80, 80, 80, 80, 80, 398 80, 80, 80, 80, 80, 80, 80, 80, 399 80, 80, 80, 80, 80, 80, 80, 80, 400 80, 80, 80, 80, 80, 80, 80, 80, 401 80, 80, 80, 80, 80, 80, 80, 80, 402 80 403 }; 404 405 406 #define EXACT_QUANT 407 #ifdef EXACT_QUANT 408 static void invert_quant(int improved_quant, short *quant, 409 short *shift, short d) 410 { 411 if(improved_quant) 412 { 413 unsigned t; 414 int l; 415 t = d; 416 for(l = 0; t > 1; l++) 417 t>>=1; 418 t = 1 + (1<<(16+l))/d; 419 *quant = (short)(t - (1<<16)); 420 *shift = l; 421 /* use multiplication and constant shift by 16 */ 422 *shift = 1 << (16 - *shift); 423 } 424 else 425 { 426 *quant = (1 << 16) / d; 427 *shift = 0; 428 /* use multiplication and constant shift by 16 */ 429 *shift = 1 << (16 - *shift); 430 } 431 } 432 433 434 void vp8cx_init_quantizer(VP8_COMP *cpi) 435 { 436 int i; 437 int quant_val; 438 int Q; 439 440 int zbin_boost[16] = {0, 0, 8, 10, 12, 14, 16, 20, 24, 28, 32, 36, 40, 44, 441 44, 44}; 442 443 for (Q = 0; Q < QINDEX_RANGE; Q++) 444 { 445 /* dc values */ 446 quant_val = vp8_dc_quant(Q, cpi->common.y1dc_delta_q); 447 cpi->Y1quant_fast[Q][0] = (1 << 16) / quant_val; 448 invert_quant(cpi->sf.improved_quant, cpi->Y1quant[Q] + 0, 449 cpi->Y1quant_shift[Q] + 0, quant_val); 450 cpi->Y1zbin[Q][0] = ((qzbin_factors[Q] * quant_val) + 64) >> 7; 451 cpi->Y1round[Q][0] = (qrounding_factors[Q] * quant_val) >> 7; 452 cpi->common.Y1dequant[Q][0] = quant_val; 453 cpi->zrun_zbin_boost_y1[Q][0] = (quant_val * zbin_boost[0]) >> 7; 454 455 quant_val = vp8_dc2quant(Q, cpi->common.y2dc_delta_q); 456 cpi->Y2quant_fast[Q][0] = (1 << 16) / quant_val; 457 invert_quant(cpi->sf.improved_quant, cpi->Y2quant[Q] + 0, 458 cpi->Y2quant_shift[Q] + 0, quant_val); 459 cpi->Y2zbin[Q][0] = ((qzbin_factors_y2[Q] * quant_val) + 64) >> 7; 460 cpi->Y2round[Q][0] = (qrounding_factors_y2[Q] * quant_val) >> 7; 461 cpi->common.Y2dequant[Q][0] = quant_val; 462 cpi->zrun_zbin_boost_y2[Q][0] = (quant_val * zbin_boost[0]) >> 7; 463 464 quant_val = vp8_dc_uv_quant(Q, cpi->common.uvdc_delta_q); 465 cpi->UVquant_fast[Q][0] = (1 << 16) / quant_val; 466 invert_quant(cpi->sf.improved_quant, cpi->UVquant[Q] + 0, 467 cpi->UVquant_shift[Q] + 0, quant_val); 468 cpi->UVzbin[Q][0] = ((qzbin_factors[Q] * quant_val) + 64) >> 7;; 469 cpi->UVround[Q][0] = (qrounding_factors[Q] * quant_val) >> 7; 470 cpi->common.UVdequant[Q][0] = quant_val; 471 cpi->zrun_zbin_boost_uv[Q][0] = (quant_val * zbin_boost[0]) >> 7; 472 473 /* all the ac values = ; */ 474 quant_val = vp8_ac_yquant(Q); 475 cpi->Y1quant_fast[Q][1] = (1 << 16) / quant_val; 476 invert_quant(cpi->sf.improved_quant, cpi->Y1quant[Q] + 1, 477 cpi->Y1quant_shift[Q] + 1, quant_val); 478 cpi->Y1zbin[Q][1] = ((qzbin_factors[Q] * quant_val) + 64) >> 7; 479 cpi->Y1round[Q][1] = (qrounding_factors[Q] * quant_val) >> 7; 480 cpi->common.Y1dequant[Q][1] = quant_val; 481 cpi->zrun_zbin_boost_y1[Q][1] = (quant_val * zbin_boost[1]) >> 7; 482 483 quant_val = vp8_ac2quant(Q, cpi->common.y2ac_delta_q); 484 cpi->Y2quant_fast[Q][1] = (1 << 16) / quant_val; 485 invert_quant(cpi->sf.improved_quant, cpi->Y2quant[Q] + 1, 486 cpi->Y2quant_shift[Q] + 1, quant_val); 487 cpi->Y2zbin[Q][1] = ((qzbin_factors_y2[Q] * quant_val) + 64) >> 7; 488 cpi->Y2round[Q][1] = (qrounding_factors_y2[Q] * quant_val) >> 7; 489 cpi->common.Y2dequant[Q][1] = quant_val; 490 cpi->zrun_zbin_boost_y2[Q][1] = (quant_val * zbin_boost[1]) >> 7; 491 492 quant_val = vp8_ac_uv_quant(Q, cpi->common.uvac_delta_q); 493 cpi->UVquant_fast[Q][1] = (1 << 16) / quant_val; 494 invert_quant(cpi->sf.improved_quant, cpi->UVquant[Q] + 1, 495 cpi->UVquant_shift[Q] + 1, quant_val); 496 cpi->UVzbin[Q][1] = ((qzbin_factors[Q] * quant_val) + 64) >> 7; 497 cpi->UVround[Q][1] = (qrounding_factors[Q] * quant_val) >> 7; 498 cpi->common.UVdequant[Q][1] = quant_val; 499 cpi->zrun_zbin_boost_uv[Q][1] = (quant_val * zbin_boost[1]) >> 7; 500 501 for (i = 2; i < 16; i++) 502 { 503 cpi->Y1quant_fast[Q][i] = cpi->Y1quant_fast[Q][1]; 504 cpi->Y1quant[Q][i] = cpi->Y1quant[Q][1]; 505 cpi->Y1quant_shift[Q][i] = cpi->Y1quant_shift[Q][1]; 506 cpi->Y1zbin[Q][i] = cpi->Y1zbin[Q][1]; 507 cpi->Y1round[Q][i] = cpi->Y1round[Q][1]; 508 cpi->zrun_zbin_boost_y1[Q][i] = (cpi->common.Y1dequant[Q][1] * 509 zbin_boost[i]) >> 7; 510 511 cpi->Y2quant_fast[Q][i] = cpi->Y2quant_fast[Q][1]; 512 cpi->Y2quant[Q][i] = cpi->Y2quant[Q][1]; 513 cpi->Y2quant_shift[Q][i] = cpi->Y2quant_shift[Q][1]; 514 cpi->Y2zbin[Q][i] = cpi->Y2zbin[Q][1]; 515 cpi->Y2round[Q][i] = cpi->Y2round[Q][1]; 516 cpi->zrun_zbin_boost_y2[Q][i] = (cpi->common.Y2dequant[Q][1] * 517 zbin_boost[i]) >> 7; 518 519 cpi->UVquant_fast[Q][i] = cpi->UVquant_fast[Q][1]; 520 cpi->UVquant[Q][i] = cpi->UVquant[Q][1]; 521 cpi->UVquant_shift[Q][i] = cpi->UVquant_shift[Q][1]; 522 cpi->UVzbin[Q][i] = cpi->UVzbin[Q][1]; 523 cpi->UVround[Q][i] = cpi->UVround[Q][1]; 524 cpi->zrun_zbin_boost_uv[Q][i] = (cpi->common.UVdequant[Q][1] * 525 zbin_boost[i]) >> 7; 526 } 527 } 528 } 529 #else 530 void vp8cx_init_quantizer(VP8_COMP *cpi) 531 { 532 int i; 533 int quant_val; 534 int Q; 535 536 int zbin_boost[16] = {0, 0, 8, 10, 12, 14, 16, 20, 24, 28, 32, 36, 40, 44, 44, 44}; 537 538 for (Q = 0; Q < QINDEX_RANGE; Q++) 539 { 540 /* dc values */ 541 quant_val = vp8_dc_quant(Q, cpi->common.y1dc_delta_q); 542 cpi->Y1quant[Q][0] = (1 << 16) / quant_val; 543 cpi->Y1zbin[Q][0] = ((qzbin_factors[Q] * quant_val) + 64) >> 7; 544 cpi->Y1round[Q][0] = (qrounding_factors[Q] * quant_val) >> 7; 545 cpi->common.Y1dequant[Q][0] = quant_val; 546 cpi->zrun_zbin_boost_y1[Q][0] = (quant_val * zbin_boost[0]) >> 7; 547 548 quant_val = vp8_dc2quant(Q, cpi->common.y2dc_delta_q); 549 cpi->Y2quant[Q][0] = (1 << 16) / quant_val; 550 cpi->Y2zbin[Q][0] = ((qzbin_factors_y2[Q] * quant_val) + 64) >> 7; 551 cpi->Y2round[Q][0] = (qrounding_factors_y2[Q] * quant_val) >> 7; 552 cpi->common.Y2dequant[Q][0] = quant_val; 553 cpi->zrun_zbin_boost_y2[Q][0] = (quant_val * zbin_boost[0]) >> 7; 554 555 quant_val = vp8_dc_uv_quant(Q, cpi->common.uvdc_delta_q); 556 cpi->UVquant[Q][0] = (1 << 16) / quant_val; 557 cpi->UVzbin[Q][0] = ((qzbin_factors[Q] * quant_val) + 64) >> 7;; 558 cpi->UVround[Q][0] = (qrounding_factors[Q] * quant_val) >> 7; 559 cpi->common.UVdequant[Q][0] = quant_val; 560 cpi->zrun_zbin_boost_uv[Q][0] = (quant_val * zbin_boost[0]) >> 7; 561 562 /* all the ac values = ; */ 563 for (i = 1; i < 16; i++) 564 { 565 int rc = vp8_default_zig_zag1d[i]; 566 567 quant_val = vp8_ac_yquant(Q); 568 cpi->Y1quant[Q][rc] = (1 << 16) / quant_val; 569 cpi->Y1zbin[Q][rc] = ((qzbin_factors[Q] * quant_val) + 64) >> 7; 570 cpi->Y1round[Q][rc] = (qrounding_factors[Q] * quant_val) >> 7; 571 cpi->common.Y1dequant[Q][rc] = quant_val; 572 cpi->zrun_zbin_boost_y1[Q][i] = (quant_val * zbin_boost[i]) >> 7; 573 574 quant_val = vp8_ac2quant(Q, cpi->common.y2ac_delta_q); 575 cpi->Y2quant[Q][rc] = (1 << 16) / quant_val; 576 cpi->Y2zbin[Q][rc] = ((qzbin_factors_y2[Q] * quant_val) + 64) >> 7; 577 cpi->Y2round[Q][rc] = (qrounding_factors_y2[Q] * quant_val) >> 7; 578 cpi->common.Y2dequant[Q][rc] = quant_val; 579 cpi->zrun_zbin_boost_y2[Q][i] = (quant_val * zbin_boost[i]) >> 7; 580 581 quant_val = vp8_ac_uv_quant(Q, cpi->common.uvac_delta_q); 582 cpi->UVquant[Q][rc] = (1 << 16) / quant_val; 583 cpi->UVzbin[Q][rc] = ((qzbin_factors[Q] * quant_val) + 64) >> 7; 584 cpi->UVround[Q][rc] = (qrounding_factors[Q] * quant_val) >> 7; 585 cpi->common.UVdequant[Q][rc] = quant_val; 586 cpi->zrun_zbin_boost_uv[Q][i] = (quant_val * zbin_boost[i]) >> 7; 587 } 588 } 589 } 590 #endif 591 592 #define ZBIN_EXTRA_Y \ 593 (( cpi->common.Y1dequant[QIndex][1] * \ 594 ( x->zbin_over_quant + \ 595 x->zbin_mode_boost + \ 596 x->act_zbin_adj ) ) >> 7) 597 598 #define ZBIN_EXTRA_UV \ 599 (( cpi->common.UVdequant[QIndex][1] * \ 600 ( x->zbin_over_quant + \ 601 x->zbin_mode_boost + \ 602 x->act_zbin_adj ) ) >> 7) 603 604 #define ZBIN_EXTRA_Y2 \ 605 (( cpi->common.Y2dequant[QIndex][1] * \ 606 ( (x->zbin_over_quant / 2) + \ 607 x->zbin_mode_boost + \ 608 x->act_zbin_adj ) ) >> 7) 609 610 void vp8cx_mb_init_quantizer(VP8_COMP *cpi, MACROBLOCK *x, int ok_to_skip) 611 { 612 int i; 613 int QIndex; 614 MACROBLOCKD *xd = &x->e_mbd; 615 int zbin_extra; 616 617 /* Select the baseline MB Q index. */ 618 if (xd->segmentation_enabled) 619 { 620 /* Abs Value */ 621 if (xd->mb_segement_abs_delta == SEGMENT_ABSDATA) 622 QIndex = xd->segment_feature_data[MB_LVL_ALT_Q][xd->mode_info_context->mbmi.segment_id]; 623 /* Delta Value */ 624 else 625 { 626 QIndex = cpi->common.base_qindex + xd->segment_feature_data[MB_LVL_ALT_Q][xd->mode_info_context->mbmi.segment_id]; 627 /* Clamp to valid range */ 628 QIndex = (QIndex >= 0) ? ((QIndex <= MAXQ) ? QIndex : MAXQ) : 0; 629 } 630 } 631 else 632 QIndex = cpi->common.base_qindex; 633 634 /* This initialization should be called at least once. Use ok_to_skip to 635 * decide if it is ok to skip. 636 * Before encoding a frame, this function is always called with ok_to_skip 637 * =0, which means no skiping of calculations. The "last" values are 638 * initialized at that time. 639 */ 640 if (!ok_to_skip || QIndex != x->q_index) 641 { 642 643 xd->dequant_y1_dc[0] = 1; 644 xd->dequant_y1[0] = cpi->common.Y1dequant[QIndex][0]; 645 xd->dequant_y2[0] = cpi->common.Y2dequant[QIndex][0]; 646 xd->dequant_uv[0] = cpi->common.UVdequant[QIndex][0]; 647 648 for (i = 1; i < 16; i++) 649 { 650 xd->dequant_y1_dc[i] = 651 xd->dequant_y1[i] = cpi->common.Y1dequant[QIndex][1]; 652 xd->dequant_y2[i] = cpi->common.Y2dequant[QIndex][1]; 653 xd->dequant_uv[i] = cpi->common.UVdequant[QIndex][1]; 654 } 655 #if 1 656 /*TODO: Remove dequant from BLOCKD. This is a temporary solution until 657 * the quantizer code uses a passed in pointer to the dequant constants. 658 * This will also require modifications to the x86 and neon assembly. 659 * */ 660 for (i = 0; i < 16; i++) 661 x->e_mbd.block[i].dequant = xd->dequant_y1; 662 for (i = 16; i < 24; i++) 663 x->e_mbd.block[i].dequant = xd->dequant_uv; 664 x->e_mbd.block[24].dequant = xd->dequant_y2; 665 #endif 666 667 /* Y */ 668 zbin_extra = ZBIN_EXTRA_Y; 669 670 for (i = 0; i < 16; i++) 671 { 672 x->block[i].quant = cpi->Y1quant[QIndex]; 673 x->block[i].quant_fast = cpi->Y1quant_fast[QIndex]; 674 x->block[i].quant_shift = cpi->Y1quant_shift[QIndex]; 675 x->block[i].zbin = cpi->Y1zbin[QIndex]; 676 x->block[i].round = cpi->Y1round[QIndex]; 677 x->block[i].zrun_zbin_boost = cpi->zrun_zbin_boost_y1[QIndex]; 678 x->block[i].zbin_extra = (short)zbin_extra; 679 } 680 681 /* UV */ 682 zbin_extra = ZBIN_EXTRA_UV; 683 684 for (i = 16; i < 24; i++) 685 { 686 x->block[i].quant = cpi->UVquant[QIndex]; 687 x->block[i].quant_fast = cpi->UVquant_fast[QIndex]; 688 x->block[i].quant_shift = cpi->UVquant_shift[QIndex]; 689 x->block[i].zbin = cpi->UVzbin[QIndex]; 690 x->block[i].round = cpi->UVround[QIndex]; 691 x->block[i].zrun_zbin_boost = cpi->zrun_zbin_boost_uv[QIndex]; 692 x->block[i].zbin_extra = (short)zbin_extra; 693 } 694 695 /* Y2 */ 696 zbin_extra = ZBIN_EXTRA_Y2; 697 698 x->block[24].quant_fast = cpi->Y2quant_fast[QIndex]; 699 x->block[24].quant = cpi->Y2quant[QIndex]; 700 x->block[24].quant_shift = cpi->Y2quant_shift[QIndex]; 701 x->block[24].zbin = cpi->Y2zbin[QIndex]; 702 x->block[24].round = cpi->Y2round[QIndex]; 703 x->block[24].zrun_zbin_boost = cpi->zrun_zbin_boost_y2[QIndex]; 704 x->block[24].zbin_extra = (short)zbin_extra; 705 706 /* save this macroblock QIndex for vp8_update_zbin_extra() */ 707 x->q_index = QIndex; 708 709 x->last_zbin_over_quant = x->zbin_over_quant; 710 x->last_zbin_mode_boost = x->zbin_mode_boost; 711 x->last_act_zbin_adj = x->act_zbin_adj; 712 713 714 715 } 716 else if(x->last_zbin_over_quant != x->zbin_over_quant 717 || x->last_zbin_mode_boost != x->zbin_mode_boost 718 || x->last_act_zbin_adj != x->act_zbin_adj) 719 { 720 /* Y */ 721 zbin_extra = ZBIN_EXTRA_Y; 722 723 for (i = 0; i < 16; i++) 724 x->block[i].zbin_extra = (short)zbin_extra; 725 726 /* UV */ 727 zbin_extra = ZBIN_EXTRA_UV; 728 729 for (i = 16; i < 24; i++) 730 x->block[i].zbin_extra = (short)zbin_extra; 731 732 /* Y2 */ 733 zbin_extra = ZBIN_EXTRA_Y2; 734 x->block[24].zbin_extra = (short)zbin_extra; 735 736 x->last_zbin_over_quant = x->zbin_over_quant; 737 x->last_zbin_mode_boost = x->zbin_mode_boost; 738 x->last_act_zbin_adj = x->act_zbin_adj; 739 } 740 } 741 742 void vp8_update_zbin_extra(VP8_COMP *cpi, MACROBLOCK *x) 743 { 744 int i; 745 int QIndex = x->q_index; 746 int zbin_extra; 747 748 /* Y */ 749 zbin_extra = ZBIN_EXTRA_Y; 750 751 for (i = 0; i < 16; i++) 752 x->block[i].zbin_extra = (short)zbin_extra; 753 754 /* UV */ 755 zbin_extra = ZBIN_EXTRA_UV; 756 757 for (i = 16; i < 24; i++) 758 x->block[i].zbin_extra = (short)zbin_extra; 759 760 /* Y2 */ 761 zbin_extra = ZBIN_EXTRA_Y2; 762 x->block[24].zbin_extra = (short)zbin_extra; 763 } 764 #undef ZBIN_EXTRA_Y 765 #undef ZBIN_EXTRA_UV 766 #undef ZBIN_EXTRA_Y2 767 768 void vp8cx_frame_init_quantizer(VP8_COMP *cpi) 769 { 770 /* Clear Zbin mode boost for default case */ 771 cpi->mb.zbin_mode_boost = 0; 772 773 /* MB level quantizer setup */ 774 vp8cx_mb_init_quantizer(cpi, &cpi->mb, 0); 775 } 776 777 778 void vp8_set_quantizer(struct VP8_COMP *cpi, int Q) 779 { 780 VP8_COMMON *cm = &cpi->common; 781 MACROBLOCKD *mbd = &cpi->mb.e_mbd; 782 int update = 0; 783 int new_delta_q; 784 cm->base_qindex = Q; 785 786 /* if any of the delta_q values are changing update flag has to be set */ 787 /* currently only y2dc_delta_q may change */ 788 789 cm->y1dc_delta_q = 0; 790 cm->y2ac_delta_q = 0; 791 cm->uvdc_delta_q = 0; 792 cm->uvac_delta_q = 0; 793 794 if (Q < 4) 795 { 796 new_delta_q = 4-Q; 797 } 798 else 799 new_delta_q = 0; 800 801 update |= cm->y2dc_delta_q != new_delta_q; 802 cm->y2dc_delta_q = new_delta_q; 803 804 805 /* Set Segment specific quatizers */ 806 mbd->segment_feature_data[MB_LVL_ALT_Q][0] = cpi->segment_feature_data[MB_LVL_ALT_Q][0]; 807 mbd->segment_feature_data[MB_LVL_ALT_Q][1] = cpi->segment_feature_data[MB_LVL_ALT_Q][1]; 808 mbd->segment_feature_data[MB_LVL_ALT_Q][2] = cpi->segment_feature_data[MB_LVL_ALT_Q][2]; 809 mbd->segment_feature_data[MB_LVL_ALT_Q][3] = cpi->segment_feature_data[MB_LVL_ALT_Q][3]; 810 811 /* quantizer has to be reinitialized for any delta_q changes */ 812 if(update) 813 vp8cx_init_quantizer(cpi); 814 815 } 816
