Index: source/libvpx/vp9/encoder/vp9_encodemb.c |
=================================================================== |
--- source/libvpx/vp9/encoder/vp9_encodemb.c (revision 232232) |
+++ source/libvpx/vp9/encoder/vp9_encodemb.c (working copy) |
@@ -8,19 +8,22 @@ |
* be found in the AUTHORS file in the root of the source tree. |
*/ |
+ |
+#include "./vp9_rtcd.h" |
#include "./vpx_config.h" |
-#include "vp9/encoder/vp9_encodemb.h" |
+ |
+#include "vpx_mem/vpx_mem.h" |
+ |
+#include "vp9/common/vp9_idct.h" |
#include "vp9/common/vp9_reconinter.h" |
-#include "vp9/encoder/vp9_quantize.h" |
-#include "vp9/encoder/vp9_tokenize.h" |
#include "vp9/common/vp9_reconintra.h" |
-#include "vpx_mem/vpx_mem.h" |
-#include "vp9/encoder/vp9_rdopt.h" |
#include "vp9/common/vp9_systemdependent.h" |
-#include "vp9_rtcd.h" |
-DECLARE_ALIGNED(16, extern const uint8_t, |
- vp9_pt_energy_class[MAX_ENTROPY_TOKENS]); |
+#include "vp9/encoder/vp9_dct.h" |
+#include "vp9/encoder/vp9_encodemb.h" |
+#include "vp9/encoder/vp9_quantize.h" |
+#include "vp9/encoder/vp9_rdopt.h" |
+#include "vp9/encoder/vp9_tokenize.h" |
void vp9_subtract_block_c(int rows, int cols, |
int16_t *diff_ptr, ptrdiff_t diff_stride, |
@@ -38,37 +41,6 @@ |
} |
} |
-static void inverse_transform_b_4x4_add(MACROBLOCKD *xd, int eob, |
- int16_t *dqcoeff, uint8_t *dest, |
- int stride) { |
- if (eob <= 1) |
- xd->inv_txm4x4_1_add(dqcoeff, dest, stride); |
- else |
- xd->inv_txm4x4_add(dqcoeff, dest, stride); |
-} |
- |
-static void inverse_transform_b_8x8_add(int eob, |
- int16_t *dqcoeff, uint8_t *dest, |
- int stride) { |
- if (eob <= 1) |
- vp9_short_idct8x8_1_add(dqcoeff, dest, stride); |
- else if (eob <= 10) |
- vp9_short_idct10_8x8_add(dqcoeff, dest, stride); |
- else |
- vp9_short_idct8x8_add(dqcoeff, dest, stride); |
-} |
- |
-static void inverse_transform_b_16x16_add(int eob, |
- int16_t *dqcoeff, uint8_t *dest, |
- int stride) { |
- if (eob <= 1) |
- vp9_short_idct16x16_1_add(dqcoeff, dest, stride); |
- else if (eob <= 10) |
- vp9_short_idct10_16x16_add(dqcoeff, dest, stride); |
- else |
- vp9_short_idct16x16_add(dqcoeff, dest, stride); |
-} |
- |
static void subtract_plane(MACROBLOCK *x, BLOCK_SIZE bsize, int plane) { |
struct macroblock_plane *const p = &x->plane[plane]; |
const MACROBLOCKD *const xd = &x->e_mbd; |
@@ -97,8 +69,7 @@ |
vp9_subtract_sbuv(x, bsize); |
} |
- |
-#define RDTRUNC(RM,DM,R,D) ( (128+(R)*(RM)) & 0xFF ) |
+#define RDTRUNC(RM, DM, R, D) ((128 + (R) * (RM)) & 0xFF) |
typedef struct vp9_token_state vp9_token_state; |
struct vp9_token_state { |
@@ -109,7 +80,7 @@ |
short qc; |
}; |
-// TODO: experiments to find optimal multiple numbers |
+// TODO(jimbankoski): experiment to find optimal RD numbers. |
#define Y1_RD_MULT 4 |
#define UV_RD_MULT 2 |
@@ -147,7 +118,7 @@ |
TX_SIZE tx_size) { |
MACROBLOCKD *const xd = &mb->e_mbd; |
struct macroblockd_plane *pd = &xd->plane[plane]; |
- const int ref = is_inter_block(&xd->this_mi->mbmi); |
+ const int ref = is_inter_block(&xd->mi_8x8[0]->mbmi); |
vp9_token_state tokens[1025][2]; |
unsigned best_index[1025][2]; |
const int16_t *coeff_ptr = BLOCK_OFFSET(mb->plane[plane].coeff, block); |
@@ -161,40 +132,18 @@ |
int best, band, pt; |
PLANE_TYPE type = pd->plane_type; |
int err_mult = plane_rd_mult[type]; |
- int default_eob; |
+ const int default_eob = 16 << (tx_size << 1); |
const int16_t *scan, *nb; |
const int mul = 1 + (tx_size == TX_32X32); |
uint8_t token_cache[1024]; |
const int ib = txfrm_block_to_raster_block(plane_bsize, tx_size, block); |
const int16_t *dequant_ptr = pd->dequant; |
- const uint8_t * band_translate; |
+ const uint8_t *const band_translate = get_band_translate(tx_size); |
assert((!type && !plane) || (type && plane)); |
dqcoeff_ptr = BLOCK_OFFSET(pd->dqcoeff, block); |
qcoeff_ptr = BLOCK_OFFSET(pd->qcoeff, block); |
- switch (tx_size) { |
- default: |
- case TX_4X4: |
- default_eob = 16; |
- scan = get_scan_4x4(get_tx_type_4x4(type, xd, ib)); |
- band_translate = vp9_coefband_trans_4x4; |
- break; |
- case TX_8X8: |
- scan = get_scan_8x8(get_tx_type_8x8(type, xd)); |
- default_eob = 64; |
- band_translate = vp9_coefband_trans_8x8plus; |
- break; |
- case TX_16X16: |
- scan = get_scan_16x16(get_tx_type_16x16(type, xd)); |
- default_eob = 256; |
- band_translate = vp9_coefband_trans_8x8plus; |
- break; |
- case TX_32X32: |
- scan = vp9_default_scan_32x32; |
- default_eob = 1024; |
- band_translate = vp9_coefband_trans_8x8plus; |
- break; |
- } |
+ get_scan(xd, tx_size, type, ib, &scan, &nb); |
assert(eob <= default_eob); |
/* Now set up a Viterbi trellis to evaluate alternative roundings. */ |
@@ -213,7 +162,6 @@ |
for (i = 0; i < eob; i++) |
token_cache[scan[i]] = vp9_pt_energy_class[vp9_dct_value_tokens_ptr[ |
qcoeff_ptr[scan[i]]].token]; |
- nb = vp9_get_coef_neighbors_handle(scan); |
for (i = eob; i-- > i0;) { |
int base_bits, d2, dx; |
@@ -312,11 +260,10 @@ |
best_index[i][1] = best; |
/* Finally, make this the new head of the trellis. */ |
next = i; |
- } |
- /* There's no choice to make for a zero coefficient, so we don't |
- * add a new trellis node, but we do need to update the costs. |
- */ |
- else { |
+ } else { |
+ /* There's no choice to make for a zero coefficient, so we don't |
+ * add a new trellis node, but we do need to update the costs. |
+ */ |
band = get_coef_band(band_translate, i + 1); |
t0 = tokens[next][0].token; |
t1 = tokens[next][1].token; |
@@ -385,38 +332,12 @@ |
const BLOCK_SIZE plane_bsize = get_plane_block_size(bsize, pd); |
const int num_4x4_w = num_4x4_blocks_wide_lookup[plane_bsize]; |
const int num_4x4_h = num_4x4_blocks_high_lookup[plane_bsize]; |
- const MB_MODE_INFO *mbmi = &xd->this_mi->mbmi; |
+ const MB_MODE_INFO *mbmi = &xd->mi_8x8[0]->mbmi; |
const TX_SIZE tx_size = plane ? get_uv_tx_size(mbmi) : mbmi->tx_size; |
- int i; |
- switch (tx_size) { |
- case TX_4X4: |
- vpx_memcpy(args->ctx->ta[plane], pd->above_context, |
- sizeof(ENTROPY_CONTEXT) * num_4x4_w); |
- vpx_memcpy(args->ctx->tl[plane], pd->left_context, |
- sizeof(ENTROPY_CONTEXT) * num_4x4_h); |
- break; |
- case TX_8X8: |
- for (i = 0; i < num_4x4_w; i += 2) |
- args->ctx->ta[plane][i] = !!*(uint16_t *)&pd->above_context[i]; |
- for (i = 0; i < num_4x4_h; i += 2) |
- args->ctx->tl[plane][i] = !!*(uint16_t *)&pd->left_context[i]; |
- break; |
- case TX_16X16: |
- for (i = 0; i < num_4x4_w; i += 4) |
- args->ctx->ta[plane][i] = !!*(uint32_t *)&pd->above_context[i]; |
- for (i = 0; i < num_4x4_h; i += 4) |
- args->ctx->tl[plane][i] = !!*(uint32_t *)&pd->left_context[i]; |
- break; |
- case TX_32X32: |
- for (i = 0; i < num_4x4_w; i += 8) |
- args->ctx->ta[plane][i] = !!*(uint64_t *)&pd->above_context[i]; |
- for (i = 0; i < num_4x4_h; i += 8) |
- args->ctx->tl[plane][i] = !!*(uint64_t *)&pd->left_context[i]; |
- break; |
- default: |
- assert(0); |
- } |
+ vp9_get_entropy_contexts(tx_size, args->ctx->ta[plane], args->ctx->tl[plane], |
+ pd->above_context, pd->left_context, |
+ num_4x4_w, num_4x4_h); |
} |
void vp9_xform_quant(int plane, int block, BLOCK_SIZE plane_bsize, |
@@ -445,9 +366,9 @@ |
yoff = 32 * (block >> twl); |
src_diff = p->src_diff + 4 * bw * yoff + xoff; |
if (x->use_lp32x32fdct) |
- vp9_short_fdct32x32_rd(src_diff, coeff, bw * 8); |
+ vp9_fdct32x32_rd(src_diff, coeff, bw * 4); |
else |
- vp9_short_fdct32x32(src_diff, coeff, bw * 8); |
+ vp9_fdct32x32(src_diff, coeff, bw * 4); |
vp9_quantize_b_32x32(coeff, 1024, x->skip_block, p->zbin, p->round, |
p->quant, p->quant_shift, qcoeff, dqcoeff, |
pd->dequant, p->zbin_extra, eob, scan, iscan); |
@@ -459,7 +380,7 @@ |
xoff = 16 * (block & twmask); |
yoff = 16 * (block >> twl); |
src_diff = p->src_diff + 4 * bw * yoff + xoff; |
- x->fwd_txm16x16(src_diff, coeff, bw * 8); |
+ vp9_fdct16x16(src_diff, coeff, bw * 4); |
vp9_quantize_b(coeff, 256, x->skip_block, p->zbin, p->round, |
p->quant, p->quant_shift, qcoeff, dqcoeff, |
pd->dequant, p->zbin_extra, eob, scan, iscan); |
@@ -471,7 +392,7 @@ |
xoff = 8 * (block & twmask); |
yoff = 8 * (block >> twl); |
src_diff = p->src_diff + 4 * bw * yoff + xoff; |
- x->fwd_txm8x8(src_diff, coeff, bw * 8); |
+ vp9_fdct8x8(src_diff, coeff, bw * 4); |
vp9_quantize_b(coeff, 64, x->skip_block, p->zbin, p->round, |
p->quant, p->quant_shift, qcoeff, dqcoeff, |
pd->dequant, p->zbin_extra, eob, scan, iscan); |
@@ -482,7 +403,7 @@ |
xoff = 4 * (block & twmask); |
yoff = 4 * (block >> twl); |
src_diff = p->src_diff + 4 * bw * yoff + xoff; |
- x->fwd_txm4x4(src_diff, coeff, bw * 8); |
+ x->fwd_txm4x4(src_diff, coeff, bw * 4); |
vp9_quantize_b(coeff, 16, x->skip_block, p->zbin, p->round, |
p->quant, p->quant_shift, qcoeff, dqcoeff, |
pd->dequant, p->zbin_extra, eob, scan, iscan); |
@@ -497,6 +418,7 @@ |
struct encode_b_args *const args = arg; |
MACROBLOCK *const x = args->x; |
MACROBLOCKD *const xd = &x->e_mbd; |
+ struct optimize_ctx *const ctx = args->ctx; |
struct macroblockd_plane *const pd = &xd->plane[plane]; |
const int raster_block = txfrm_block_to_raster_block(plane_bsize, tx_size, |
block); |
@@ -504,38 +426,68 @@ |
int16_t *const dqcoeff = BLOCK_OFFSET(pd->dqcoeff, block); |
uint8_t *const dst = raster_block_offset_uint8(plane_bsize, raster_block, |
pd->dst.buf, pd->dst.stride); |
+ |
+ // TODO(jingning): per transformed block zero forcing only enabled for |
+ // luma component. will integrate chroma components as well. |
+ if (x->zcoeff_blk[tx_size][block] && plane == 0) { |
+ int x, y; |
+ pd->eobs[block] = 0; |
+ txfrm_block_to_raster_xy(plane_bsize, tx_size, block, &x, &y); |
+ ctx->ta[plane][x] = 0; |
+ ctx->tl[plane][y] = 0; |
+ return; |
+ } |
+ |
vp9_xform_quant(plane, block, plane_bsize, tx_size, arg); |
if (x->optimize) |
- vp9_optimize_b(plane, block, plane_bsize, tx_size, x, args->ctx); |
+ vp9_optimize_b(plane, block, plane_bsize, tx_size, x, ctx); |
if (x->skip_encode || pd->eobs[block] == 0) |
return; |
switch (tx_size) { |
case TX_32X32: |
- vp9_short_idct32x32_add(dqcoeff, dst, pd->dst.stride); |
+ vp9_idct32x32_add(dqcoeff, dst, pd->dst.stride, pd->eobs[block]); |
break; |
case TX_16X16: |
- inverse_transform_b_16x16_add(pd->eobs[block], dqcoeff, dst, |
- pd->dst.stride); |
+ vp9_idct16x16_add(dqcoeff, dst, pd->dst.stride, pd->eobs[block]); |
break; |
case TX_8X8: |
- inverse_transform_b_8x8_add(pd->eobs[block], dqcoeff, dst, |
- pd->dst.stride); |
+ vp9_idct8x8_add(dqcoeff, dst, pd->dst.stride, pd->eobs[block]); |
break; |
case TX_4X4: |
// this is like vp9_short_idct4x4 but has a special case around eob<=1 |
// which is significant (not just an optimization) for the lossless |
// case. |
- inverse_transform_b_4x4_add(xd, pd->eobs[block], dqcoeff, |
- dst, pd->dst.stride); |
+ xd->itxm_add(dqcoeff, dst, pd->dst.stride, pd->eobs[block]); |
break; |
default: |
assert(!"Invalid transform size"); |
} |
} |
+static void encode_block_pass1(int plane, int block, BLOCK_SIZE plane_bsize, |
+ TX_SIZE tx_size, void *arg) { |
+ struct encode_b_args *const args = arg; |
+ MACROBLOCK *const x = args->x; |
+ MACROBLOCKD *const xd = &x->e_mbd; |
+ struct macroblockd_plane *const pd = &xd->plane[plane]; |
+ const int raster_block = txfrm_block_to_raster_block(plane_bsize, tx_size, |
+ block); |
+ |
+ int16_t *const dqcoeff = BLOCK_OFFSET(pd->dqcoeff, block); |
+ uint8_t *const dst = raster_block_offset_uint8(plane_bsize, raster_block, |
+ pd->dst.buf, pd->dst.stride); |
+ |
+ vp9_xform_quant(plane, block, plane_bsize, tx_size, arg); |
+ |
+ if (pd->eobs[block] == 0) |
+ return; |
+ |
+ xd->itxm_add(dqcoeff, dst, pd->dst.stride, pd->eobs[block]); |
+} |
+ |
void vp9_encode_sby(MACROBLOCK *x, BLOCK_SIZE bsize) { |
MACROBLOCKD *const xd = &x->e_mbd; |
struct optimize_ctx ctx; |
@@ -545,7 +497,7 @@ |
if (x->optimize) |
optimize_init_b(0, bsize, &arg); |
- foreach_transformed_block_in_plane(xd, bsize, 0, encode_block, &arg); |
+ foreach_transformed_block_in_plane(xd, bsize, 0, encode_block_pass1, &arg); |
} |
void vp9_encode_sb(MACROBLOCK *x, BLOCK_SIZE bsize) { |
@@ -569,7 +521,7 @@ |
struct encode_b_args* const args = arg; |
MACROBLOCK *const x = args->x; |
MACROBLOCKD *const xd = &x->e_mbd; |
- MB_MODE_INFO *mbmi = &xd->this_mi->mbmi; |
+ MB_MODE_INFO *mbmi = &xd->mi_8x8[0]->mbmi; |
struct macroblock_plane *const p = &x->plane[plane]; |
struct macroblockd_plane *const pd = &xd->plane[plane]; |
int16_t *coeff = BLOCK_OFFSET(p->coeff, block); |
@@ -607,14 +559,14 @@ |
vp9_subtract_block(32, 32, src_diff, bw * 4, |
src, p->src.stride, dst, pd->dst.stride); |
if (x->use_lp32x32fdct) |
- vp9_short_fdct32x32_rd(src_diff, coeff, bw * 8); |
+ vp9_fdct32x32_rd(src_diff, coeff, bw * 4); |
else |
- vp9_short_fdct32x32(src_diff, coeff, bw * 8); |
+ vp9_fdct32x32(src_diff, coeff, bw * 4); |
vp9_quantize_b_32x32(coeff, 1024, x->skip_block, p->zbin, p->round, |
p->quant, p->quant_shift, qcoeff, dqcoeff, |
pd->dequant, p->zbin_extra, eob, scan, iscan); |
if (!x->skip_encode && *eob) |
- vp9_short_idct32x32_add(dqcoeff, dst, pd->dst.stride); |
+ vp9_idct32x32_add(dqcoeff, dst, pd->dst.stride, *eob); |
break; |
case TX_16X16: |
tx_type = get_tx_type_16x16(pd->plane_type, xd); |
@@ -631,19 +583,12 @@ |
dst, pd->dst.stride, dst, pd->dst.stride); |
vp9_subtract_block(16, 16, src_diff, bw * 4, |
src, p->src.stride, dst, pd->dst.stride); |
- if (tx_type != DCT_DCT) |
- vp9_short_fht16x16(src_diff, coeff, bw * 4, tx_type); |
- else |
- x->fwd_txm16x16(src_diff, coeff, bw * 8); |
+ vp9_fht16x16(tx_type, src_diff, coeff, bw * 4); |
vp9_quantize_b(coeff, 256, x->skip_block, p->zbin, p->round, |
p->quant, p->quant_shift, qcoeff, dqcoeff, |
pd->dequant, p->zbin_extra, eob, scan, iscan); |
- if (!x->skip_encode && *eob) { |
- if (tx_type == DCT_DCT) |
- inverse_transform_b_16x16_add(*eob, dqcoeff, dst, pd->dst.stride); |
- else |
- vp9_short_iht16x16_add(dqcoeff, dst, pd->dst.stride, tx_type); |
- } |
+ if (!x->skip_encode && *eob) |
+ vp9_iht16x16_add(tx_type, dqcoeff, dst, pd->dst.stride, *eob); |
break; |
case TX_8X8: |
tx_type = get_tx_type_8x8(pd->plane_type, xd); |
@@ -660,26 +605,19 @@ |
dst, pd->dst.stride, dst, pd->dst.stride); |
vp9_subtract_block(8, 8, src_diff, bw * 4, |
src, p->src.stride, dst, pd->dst.stride); |
- if (tx_type != DCT_DCT) |
- vp9_short_fht8x8(src_diff, coeff, bw * 4, tx_type); |
- else |
- x->fwd_txm8x8(src_diff, coeff, bw * 8); |
+ vp9_fht8x8(tx_type, src_diff, coeff, bw * 4); |
vp9_quantize_b(coeff, 64, x->skip_block, p->zbin, p->round, p->quant, |
p->quant_shift, qcoeff, dqcoeff, |
pd->dequant, p->zbin_extra, eob, scan, iscan); |
- if (!x->skip_encode && *eob) { |
- if (tx_type == DCT_DCT) |
- inverse_transform_b_8x8_add(*eob, dqcoeff, dst, pd->dst.stride); |
- else |
- vp9_short_iht8x8_add(dqcoeff, dst, pd->dst.stride, tx_type); |
- } |
+ if (!x->skip_encode && *eob) |
+ vp9_iht8x8_add(tx_type, dqcoeff, dst, pd->dst.stride, *eob); |
break; |
case TX_4X4: |
tx_type = get_tx_type_4x4(pd->plane_type, xd, block); |
scan = get_scan_4x4(tx_type); |
iscan = get_iscan_4x4(tx_type); |
if (mbmi->sb_type < BLOCK_8X8 && plane == 0) |
- mode = xd->this_mi->bmi[block].as_mode; |
+ mode = xd->mi_8x8[0]->bmi[block].as_mode; |
else |
mode = plane == 0 ? mbmi->mode : mbmi->uv_mode; |
@@ -695,7 +633,7 @@ |
if (tx_type != DCT_DCT) |
vp9_short_fht4x4(src_diff, coeff, bw * 4, tx_type); |
else |
- x->fwd_txm4x4(src_diff, coeff, bw * 8); |
+ x->fwd_txm4x4(src_diff, coeff, bw * 4); |
vp9_quantize_b(coeff, 16, x->skip_block, p->zbin, p->round, p->quant, |
p->quant_shift, qcoeff, dqcoeff, |
pd->dequant, p->zbin_extra, eob, scan, iscan); |
@@ -704,9 +642,9 @@ |
// this is like vp9_short_idct4x4 but has a special case around eob<=1 |
// which is significant (not just an optimization) for the lossless |
// case. |
- inverse_transform_b_4x4_add(xd, *eob, dqcoeff, dst, pd->dst.stride); |
+ xd->itxm_add(dqcoeff, dst, pd->dst.stride, *eob); |
else |
- vp9_short_iht4x4_add(dqcoeff, dst, pd->dst.stride, tx_type); |
+ vp9_iht4x4_16_add(dqcoeff, dst, pd->dst.stride, tx_type); |
} |
break; |
default: |