/* * Copyright (c) 2014 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 "./vpx_dsp_rtcd.h" #include "./vpx_config.h" #include "vpx/vpx_integer.h" #include "vpx_dsp/variance.h" #include "vpx_dsp/arm/mem_neon.h" // Process a block exactly 4 wide and a multiple of 2 high. static void var_filter_block2d_bil_w4(const uint8_t *src_ptr, uint8_t *dst_ptr, int src_stride, int pixel_step, int dst_height, int filter_offset) { const uint8x8_t f0 = vdup_n_u8(8 - filter_offset); const uint8x8_t f1 = vdup_n_u8(filter_offset); int i = dst_height; do { uint8x8_t s0 = load_unaligned_u8(src_ptr, src_stride); uint8x8_t s1 = load_unaligned_u8(src_ptr + pixel_step, src_stride); uint16x8_t blend = vmlal_u8(vmull_u8(s0, f0), s1, f1); uint8x8_t blend_u8 = vrshrn_n_u16(blend, 3); vst1_u8(dst_ptr, blend_u8); src_ptr += 2 * src_stride; dst_ptr += 2 * 4; i -= 2; } while (i != 0); } // Process a block exactly 8 wide and any height. static void var_filter_block2d_bil_w8(const uint8_t *src_ptr, uint8_t *dst_ptr, int src_stride, int pixel_step, int dst_height, int filter_offset) { const uint8x8_t f0 = vdup_n_u8(8 - filter_offset); const uint8x8_t f1 = vdup_n_u8(filter_offset); int i = dst_height; do { uint8x8_t s0 = vld1_u8(src_ptr); uint8x8_t s1 = vld1_u8(src_ptr + pixel_step); uint16x8_t blend = vmlal_u8(vmull_u8(s0, f0), s1, f1); uint8x8_t blend_u8 = vrshrn_n_u16(blend, 3); vst1_u8(dst_ptr, blend_u8); src_ptr += src_stride; dst_ptr += 8; } while (--i != 0); } // Process a block which is a mutiple of 16 wide and any height. static void var_filter_block2d_bil_large(const uint8_t *src_ptr, uint8_t *dst_ptr, int src_stride, int pixel_step, int dst_width, int dst_height, int filter_offset) { const uint8x8_t f0 = vdup_n_u8(8 - filter_offset); const uint8x8_t f1 = vdup_n_u8(filter_offset); int i = dst_height; do { int j = 0; do { uint8x16_t s0 = vld1q_u8(src_ptr + j); uint8x16_t s1 = vld1q_u8(src_ptr + j + pixel_step); uint16x8_t blend_l = vmlal_u8(vmull_u8(vget_low_u8(s0), f0), vget_low_u8(s1), f1); uint16x8_t blend_h = vmlal_u8(vmull_u8(vget_high_u8(s0), f0), vget_high_u8(s1), f1); uint8x8_t out_lo = vrshrn_n_u16(blend_l, 3); uint8x8_t out_hi = vrshrn_n_u16(blend_h, 3); vst1q_u8(dst_ptr + j, vcombine_u8(out_lo, out_hi)); j += 16; } while (j < dst_width); src_ptr += src_stride; dst_ptr += dst_width; } while (--i != 0); } static void var_filter_block2d_bil_w16(const uint8_t *src_ptr, uint8_t *dst_ptr, int src_stride, int pixel_step, int dst_height, int filter_offset) { var_filter_block2d_bil_large(src_ptr, dst_ptr, src_stride, pixel_step, 16, dst_height, filter_offset); } static void var_filter_block2d_bil_w32(const uint8_t *src_ptr, uint8_t *dst_ptr, int src_stride, int pixel_step, int dst_height, int filter_offset) { var_filter_block2d_bil_large(src_ptr, dst_ptr, src_stride, pixel_step, 32, dst_height, filter_offset); } static void var_filter_block2d_bil_w64(const uint8_t *src_ptr, uint8_t *dst_ptr, int src_stride, int pixel_step, int dst_height, int filter_offset) { var_filter_block2d_bil_large(src_ptr, dst_ptr, src_stride, pixel_step, 64, dst_height, filter_offset); } static void var_filter_block2d_avg(const uint8_t *src_ptr, uint8_t *dst_ptr, int src_stride, int pixel_step, int dst_width, int dst_height) { int i = dst_height; // We only specialize on the filter values for large block sizes (>= 16x16.) assert(dst_width >= 16 && dst_width % 16 == 0); do { int j = 0; do { uint8x16_t s0 = vld1q_u8(src_ptr + j); uint8x16_t s1 = vld1q_u8(src_ptr + j + pixel_step); uint8x16_t avg = vrhaddq_u8(s0, s1); vst1q_u8(dst_ptr + j, avg); j += 16; } while (j < dst_width); src_ptr += src_stride; dst_ptr += dst_width; } while (--i != 0); } #define SUBPEL_VARIANCE_WXH_NEON(w, h, padding) \ unsigned int vpx_sub_pixel_variance##w##x##h##_neon( \ const uint8_t *src, int src_stride, int xoffset, int yoffset, \ const uint8_t *ref, int ref_stride, uint32_t *sse) { \ uint8_t tmp0[w * (h + padding)]; \ uint8_t tmp1[w * h]; \ var_filter_block2d_bil_w##w(src, tmp0, src_stride, 1, (h + padding), \ xoffset); \ var_filter_block2d_bil_w##w(tmp0, tmp1, w, w, h, yoffset); \ return vpx_variance##w##x##h(tmp1, w, ref, ref_stride, sse); \ } #define SPECIALIZED_SUBPEL_VARIANCE_WXH_NEON(w, h, padding) \ unsigned int vpx_sub_pixel_variance##w##x##h##_neon( \ const uint8_t *src, int src_stride, int xoffset, int yoffset, \ const uint8_t *ref, int ref_stride, unsigned int *sse) { \ if (xoffset == 0) { \ if (yoffset == 0) { \ return vpx_variance##w##x##h(src, src_stride, ref, ref_stride, sse); \ } else if (yoffset == 4) { \ uint8_t tmp[w * h]; \ var_filter_block2d_avg(src, tmp, src_stride, src_stride, w, h); \ return vpx_variance##w##x##h(tmp, w, ref, ref_stride, sse); \ } else { \ uint8_t tmp[w * h]; \ var_filter_block2d_bil_w##w(src, tmp, src_stride, src_stride, h, \ yoffset); \ return vpx_variance##w##x##h(tmp, w, ref, ref_stride, sse); \ } \ } else if (xoffset == 4) { \ uint8_t tmp0[w * (h + padding)]; \ if (yoffset == 0) { \ var_filter_block2d_avg(src, tmp0, src_stride, 1, w, h); \ return vpx_variance##w##x##h(tmp0, w, ref, ref_stride, sse); \ } else if (yoffset == 4) { \ uint8_t tmp1[w * (h + padding)]; \ var_filter_block2d_avg(src, tmp0, src_stride, 1, w, (h + padding)); \ var_filter_block2d_avg(tmp0, tmp1, w, w, w, h); \ return vpx_variance##w##x##h(tmp1, w, ref, ref_stride, sse); \ } else { \ uint8_t tmp1[w * (h + padding)]; \ var_filter_block2d_avg(src, tmp0, src_stride, 1, w, (h + padding)); \ var_filter_block2d_bil_w##w(tmp0, tmp1, w, w, h, yoffset); \ return vpx_variance##w##x##h(tmp1, w, ref, ref_stride, sse); \ } \ } else { \ uint8_t tmp0[w * (h + padding)]; \ if (yoffset == 0) { \ var_filter_block2d_bil_w##w(src, tmp0, src_stride, 1, h, xoffset); \ return vpx_variance##w##x##h(tmp0, w, ref, ref_stride, sse); \ } else if (yoffset == 4) { \ uint8_t tmp1[w * h]; \ var_filter_block2d_bil_w##w(src, tmp0, src_stride, 1, (h + padding), \ xoffset); \ var_filter_block2d_avg(tmp0, tmp1, w, w, w, h); \ return vpx_variance##w##x##h(tmp1, w, ref, ref_stride, sse); \ } else { \ uint8_t tmp1[w * h]; \ var_filter_block2d_bil_w##w(src, tmp0, src_stride, 1, (h + padding), \ xoffset); \ var_filter_block2d_bil_w##w(tmp0, tmp1, w, w, h, yoffset); \ return vpx_variance##w##x##h(tmp1, w, ref, ref_stride, sse); \ } \ } \ } // 4x blocks are processed two rows at a time, so require an extra row of // padding. SUBPEL_VARIANCE_WXH_NEON(4, 4, 2) SUBPEL_VARIANCE_WXH_NEON(4, 8, 2) SUBPEL_VARIANCE_WXH_NEON(8, 4, 1) SUBPEL_VARIANCE_WXH_NEON(8, 8, 1) SUBPEL_VARIANCE_WXH_NEON(8, 16, 1) SPECIALIZED_SUBPEL_VARIANCE_WXH_NEON(16, 8, 1) SPECIALIZED_SUBPEL_VARIANCE_WXH_NEON(16, 16, 1) SPECIALIZED_SUBPEL_VARIANCE_WXH_NEON(16, 32, 1) SPECIALIZED_SUBPEL_VARIANCE_WXH_NEON(32, 16, 1) SPECIALIZED_SUBPEL_VARIANCE_WXH_NEON(32, 32, 1) SPECIALIZED_SUBPEL_VARIANCE_WXH_NEON(32, 64, 1) SPECIALIZED_SUBPEL_VARIANCE_WXH_NEON(64, 32, 1) SPECIALIZED_SUBPEL_VARIANCE_WXH_NEON(64, 64, 1) // Combine bilinear filter with vpx_comp_avg_pred for blocks having width 4. static void avg_pred_var_filter_block2d_bil_w4(const uint8_t *src_ptr, uint8_t *dst_ptr, int src_stride, int pixel_step, int dst_height, int filter_offset, const uint8_t *second_pred) { const uint8x8_t f0 = vdup_n_u8(8 - filter_offset); const uint8x8_t f1 = vdup_n_u8(filter_offset); int i = dst_height; do { uint8x8_t s0 = load_unaligned_u8(src_ptr, src_stride); uint8x8_t s1 = load_unaligned_u8(src_ptr + pixel_step, src_stride); uint16x8_t blend = vmlal_u8(vmull_u8(s0, f0), s1, f1); uint8x8_t blend_u8 = vrshrn_n_u16(blend, 3); uint8x8_t p = vld1_u8(second_pred); uint8x8_t avg = vrhadd_u8(blend_u8, p); vst1_u8(dst_ptr, avg); src_ptr += 2 * src_stride; dst_ptr += 2 * 4; second_pred += 2 * 4; i -= 2; } while (i != 0); } // Combine bilinear filter with vpx_comp_avg_pred for blocks having width 8. static void avg_pred_var_filter_block2d_bil_w8(const uint8_t *src_ptr, uint8_t *dst_ptr, int src_stride, int pixel_step, int dst_height, int filter_offset, const uint8_t *second_pred) { const uint8x8_t f0 = vdup_n_u8(8 - filter_offset); const uint8x8_t f1 = vdup_n_u8(filter_offset); int i = dst_height; do { uint8x8_t s0 = vld1_u8(src_ptr); uint8x8_t s1 = vld1_u8(src_ptr + pixel_step); uint16x8_t blend = vmlal_u8(vmull_u8(s0, f0), s1, f1); uint8x8_t blend_u8 = vrshrn_n_u16(blend, 3); uint8x8_t p = vld1_u8(second_pred); uint8x8_t avg = vrhadd_u8(blend_u8, p); vst1_u8(dst_ptr, avg); src_ptr += src_stride; dst_ptr += 8; second_pred += 8; } while (--i > 0); } // Combine bilinear filter with vpx_comp_avg_pred for large blocks. static void avg_pred_var_filter_block2d_bil_large( const uint8_t *src_ptr, uint8_t *dst_ptr, int src_stride, int pixel_step, int dst_width, int dst_height, int filter_offset, const uint8_t *second_pred) { const uint8x8_t f0 = vdup_n_u8(8 - filter_offset); const uint8x8_t f1 = vdup_n_u8(filter_offset); int i = dst_height; do { int j = 0; do { uint8x16_t s0 = vld1q_u8(src_ptr + j); uint8x16_t s1 = vld1q_u8(src_ptr + j + pixel_step); uint16x8_t blend_l = vmlal_u8(vmull_u8(vget_low_u8(s0), f0), vget_low_u8(s1), f1); uint16x8_t blend_h = vmlal_u8(vmull_u8(vget_high_u8(s0), f0), vget_high_u8(s1), f1); uint8x16_t blend_u8 = vcombine_u8(vrshrn_n_u16(blend_l, 3), vrshrn_n_u16(blend_h, 3)); uint8x16_t p = vld1q_u8(second_pred); uint8x16_t avg = vrhaddq_u8(blend_u8, p); vst1q_u8(dst_ptr + j, avg); j += 16; second_pred += 16; } while (j < dst_width); src_ptr += src_stride; dst_ptr += dst_width; } while (--i != 0); } // Combine bilinear filter with vpx_comp_avg_pred for blocks having width 16. static void avg_pred_var_filter_block2d_bil_w16( const uint8_t *src_ptr, uint8_t *dst_ptr, int src_stride, int pixel_step, int dst_height, int filter_offset, const uint8_t *second_pred) { avg_pred_var_filter_block2d_bil_large(src_ptr, dst_ptr, src_stride, pixel_step, 16, dst_height, filter_offset, second_pred); } // Combine bilinear filter with vpx_comp_avg_pred for blocks having width 32. static void avg_pred_var_filter_block2d_bil_w32( const uint8_t *src_ptr, uint8_t *dst_ptr, int src_stride, int pixel_step, int dst_height, int filter_offset, const uint8_t *second_pred) { avg_pred_var_filter_block2d_bil_large(src_ptr, dst_ptr, src_stride, pixel_step, 32, dst_height, filter_offset, second_pred); } // Combine bilinear filter with vpx_comp_avg_pred for blocks having width 64. static void avg_pred_var_filter_block2d_bil_w64( const uint8_t *src_ptr, uint8_t *dst_ptr, int src_stride, int pixel_step, int dst_height, int filter_offset, const uint8_t *second_pred) { avg_pred_var_filter_block2d_bil_large(src_ptr, dst_ptr, src_stride, pixel_step, 64, dst_height, filter_offset, second_pred); } // Combine averaging subpel filter with vpx_comp_avg_pred. static void avg_pred_var_filter_block2d_avg(const uint8_t *src_ptr, uint8_t *dst_ptr, int src_stride, int pixel_step, int dst_width, int dst_height, const uint8_t *second_pred) { int i = dst_height; // We only specialize on the filter values for large block sizes (>= 16x16.) assert(dst_width >= 16 && dst_width % 16 == 0); do { int j = 0; do { uint8x16_t s0 = vld1q_u8(src_ptr + j); uint8x16_t s1 = vld1q_u8(src_ptr + j + pixel_step); uint8x16_t avg = vrhaddq_u8(s0, s1); uint8x16_t p = vld1q_u8(second_pred); avg = vrhaddq_u8(avg, p); vst1q_u8(dst_ptr + j, avg); j += 16; second_pred += 16; } while (j < dst_width); src_ptr += src_stride; dst_ptr += dst_width; } while (--i != 0); } // Implementation of vpx_comp_avg_pred for blocks having width >= 16. static void avg_pred(const uint8_t *src_ptr, uint8_t *dst_ptr, int src_stride, int dst_width, int dst_height, const uint8_t *second_pred) { int i = dst_height; // We only specialize on the filter values for large block sizes (>= 16x16.) assert(dst_width >= 16 && dst_width % 16 == 0); do { int j = 0; do { uint8x16_t s = vld1q_u8(src_ptr + j); uint8x16_t p = vld1q_u8(second_pred); uint8x16_t avg = vrhaddq_u8(s, p); vst1q_u8(dst_ptr + j, avg); j += 16; second_pred += 16; } while (j < dst_width); src_ptr += src_stride; dst_ptr += dst_width; } while (--i != 0); } #define SUBPEL_AVG_VARIANCE_WXH_NEON(w, h, padding) \ unsigned int vpx_sub_pixel_avg_variance##w##x##h##_neon( \ const uint8_t *src, int source_stride, int xoffset, int yoffset, \ const uint8_t *ref, int ref_stride, uint32_t *sse, \ const uint8_t *second_pred) { \ uint8_t tmp0[w * (h + padding)]; \ uint8_t tmp1[w * h]; \ var_filter_block2d_bil_w##w(src, tmp0, source_stride, 1, (h + padding), \ xoffset); \ avg_pred_var_filter_block2d_bil_w##w(tmp0, tmp1, w, w, h, yoffset, \ second_pred); \ return vpx_variance##w##x##h(tmp1, w, ref, ref_stride, sse); \ } #define SPECIALIZED_SUBPEL_AVG_VARIANCE_WXH_NEON(w, h, padding) \ unsigned int vpx_sub_pixel_avg_variance##w##x##h##_neon( \ const uint8_t *src, int source_stride, int xoffset, int yoffset, \ const uint8_t *ref, int ref_stride, unsigned int *sse, \ const uint8_t *second_pred) { \ if (xoffset == 0) { \ uint8_t tmp[w * h]; \ if (yoffset == 0) { \ avg_pred(src, tmp, source_stride, w, h, second_pred); \ return vpx_variance##w##x##h(tmp, w, ref, ref_stride, sse); \ } else if (yoffset == 4) { \ avg_pred_var_filter_block2d_avg(src, tmp, source_stride, \ source_stride, w, h, second_pred); \ return vpx_variance##w##x##h(tmp, w, ref, ref_stride, sse); \ } else { \ avg_pred_var_filter_block2d_bil_w##w( \ src, tmp, source_stride, source_stride, h, yoffset, second_pred); \ return vpx_variance##w##x##h(tmp, w, ref, ref_stride, sse); \ } \ } else if (xoffset == 4) { \ uint8_t tmp0[w * (h + padding)]; \ if (yoffset == 0) { \ avg_pred_var_filter_block2d_avg(src, tmp0, source_stride, 1, w, h, \ second_pred); \ return vpx_variance##w##x##h(tmp0, w, ref, ref_stride, sse); \ } else if (yoffset == 4) { \ uint8_t tmp1[w * (h + padding)]; \ var_filter_block2d_avg(src, tmp0, source_stride, 1, w, (h + padding)); \ avg_pred_var_filter_block2d_avg(tmp0, tmp1, w, w, w, h, second_pred); \ return vpx_variance##w##x##h(tmp1, w, ref, ref_stride, sse); \ } else { \ uint8_t tmp1[w * (h + padding)]; \ var_filter_block2d_avg(src, tmp0, source_stride, 1, w, (h + padding)); \ avg_pred_var_filter_block2d_bil_w##w(tmp0, tmp1, w, w, h, yoffset, \ second_pred); \ return vpx_variance##w##x##h(tmp1, w, ref, ref_stride, sse); \ } \ } else { \ uint8_t tmp0[w * (h + padding)]; \ if (yoffset == 0) { \ avg_pred_var_filter_block2d_bil_w##w(src, tmp0, source_stride, 1, h, \ xoffset, second_pred); \ return vpx_variance##w##x##h(tmp0, w, ref, ref_stride, sse); \ } else if (yoffset == 4) { \ uint8_t tmp1[w * h]; \ var_filter_block2d_bil_w##w(src, tmp0, source_stride, 1, \ (h + padding), xoffset); \ avg_pred_var_filter_block2d_avg(tmp0, tmp1, w, w, w, h, second_pred); \ return vpx_variance##w##x##h(tmp1, w, ref, ref_stride, sse); \ } else { \ uint8_t tmp1[w * h]; \ var_filter_block2d_bil_w##w(src, tmp0, source_stride, 1, \ (h + padding), xoffset); \ avg_pred_var_filter_block2d_bil_w##w(tmp0, tmp1, w, w, h, yoffset, \ second_pred); \ return vpx_variance##w##x##h(tmp1, w, ref, ref_stride, sse); \ } \ } \ } // 4x blocks are processed two rows at a time, so require an extra row of // padding. SUBPEL_AVG_VARIANCE_WXH_NEON(4, 4, 2) SUBPEL_AVG_VARIANCE_WXH_NEON(4, 8, 2) SUBPEL_AVG_VARIANCE_WXH_NEON(8, 4, 1) SUBPEL_AVG_VARIANCE_WXH_NEON(8, 8, 1) SUBPEL_AVG_VARIANCE_WXH_NEON(8, 16, 1) SPECIALIZED_SUBPEL_AVG_VARIANCE_WXH_NEON(16, 8, 1) SPECIALIZED_SUBPEL_AVG_VARIANCE_WXH_NEON(16, 16, 1) SPECIALIZED_SUBPEL_AVG_VARIANCE_WXH_NEON(16, 32, 1) SPECIALIZED_SUBPEL_AVG_VARIANCE_WXH_NEON(32, 16, 1) SPECIALIZED_SUBPEL_AVG_VARIANCE_WXH_NEON(32, 32, 1) SPECIALIZED_SUBPEL_AVG_VARIANCE_WXH_NEON(32, 64, 1) SPECIALIZED_SUBPEL_AVG_VARIANCE_WXH_NEON(64, 32, 1) SPECIALIZED_SUBPEL_AVG_VARIANCE_WXH_NEON(64, 64, 1)