| /* |
| * Copyright (c) 2015 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. |
| */ |
| #ifndef VPX_VPX_DSP_X86_CONVOLVE_H_ |
| #define VPX_VPX_DSP_X86_CONVOLVE_H_ |
| |
| #include <assert.h> |
| |
| #include "./vpx_config.h" |
| #include "vpx/vpx_integer.h" |
| #include "vpx_ports/compiler_attributes.h" |
| |
| // TODO(chiyotsai@google.com): Refactor the code here. Currently this is pretty |
| // hacky and awful to read. Note that there is a filter_x[3] == 128 check in |
| // HIGHBD_FUN_CONV_2D to avoid seg fault due to the fact that the c function |
| // assumes the filter is always 8 tap. |
| typedef void filter8_1dfunction(const uint8_t *src_ptr, ptrdiff_t src_pitch, |
| uint8_t *output_ptr, ptrdiff_t out_pitch, |
| uint32_t output_height, const int16_t *filter); |
| |
| // TODO(chiyotsai@google.com): Remove the is_avg argument to the MACROS once we |
| // have 4-tap vert avg filter. |
| #define FUN_CONV_1D(name, offset, step_q4, dir, src_start, avg, opt, is_avg) \ |
| void vpx_convolve8_##name##_##opt( \ |
| const uint8_t *src, ptrdiff_t src_stride, uint8_t *dst, \ |
| ptrdiff_t dst_stride, const InterpKernel *filter, int x0_q4, \ |
| int x_step_q4, int y0_q4, int y_step_q4, int w, int h) { \ |
| const int16_t *filter_row = filter[offset]; \ |
| (void)x0_q4; \ |
| (void)x_step_q4; \ |
| (void)y0_q4; \ |
| (void)y_step_q4; \ |
| assert(filter_row[3] != 128); \ |
| assert(step_q4 == 16); \ |
| if (filter_row[0] | filter_row[1] | filter_row[6] | filter_row[7]) { \ |
| const int num_taps = 8; \ |
| while (w >= 16) { \ |
| vpx_filter_block1d16_##dir##8_##avg##opt(src_start, src_stride, dst, \ |
| dst_stride, h, filter_row); \ |
| src += 16; \ |
| dst += 16; \ |
| w -= 16; \ |
| } \ |
| if (w == 8) { \ |
| vpx_filter_block1d8_##dir##8_##avg##opt(src_start, src_stride, dst, \ |
| dst_stride, h, filter_row); \ |
| } else if (w == 4) { \ |
| vpx_filter_block1d4_##dir##8_##avg##opt(src_start, src_stride, dst, \ |
| dst_stride, h, filter_row); \ |
| } \ |
| (void)num_taps; \ |
| } else if (filter_row[2] | filter_row[5]) { \ |
| const int num_taps = is_avg ? 8 : 4; \ |
| while (w >= 16) { \ |
| vpx_filter_block1d16_##dir##4_##avg##opt(src_start, src_stride, dst, \ |
| dst_stride, h, filter_row); \ |
| src += 16; \ |
| dst += 16; \ |
| w -= 16; \ |
| } \ |
| if (w == 8) { \ |
| vpx_filter_block1d8_##dir##4_##avg##opt(src_start, src_stride, dst, \ |
| dst_stride, h, filter_row); \ |
| } else if (w == 4) { \ |
| vpx_filter_block1d4_##dir##4_##avg##opt(src_start, src_stride, dst, \ |
| dst_stride, h, filter_row); \ |
| } \ |
| (void)num_taps; \ |
| } else { \ |
| const int num_taps = 2; \ |
| while (w >= 16) { \ |
| vpx_filter_block1d16_##dir##2_##avg##opt(src_start, src_stride, dst, \ |
| dst_stride, h, filter_row); \ |
| src += 16; \ |
| dst += 16; \ |
| w -= 16; \ |
| } \ |
| if (w == 8) { \ |
| vpx_filter_block1d8_##dir##2_##avg##opt(src_start, src_stride, dst, \ |
| dst_stride, h, filter_row); \ |
| } else if (w == 4) { \ |
| vpx_filter_block1d4_##dir##2_##avg##opt(src_start, src_stride, dst, \ |
| dst_stride, h, filter_row); \ |
| } \ |
| (void)num_taps; \ |
| } \ |
| } |
| |
| #define FUN_CONV_2D(avg, opt, is_avg) \ |
| void vpx_convolve8_##avg##opt( \ |
| const uint8_t *src, ptrdiff_t src_stride, uint8_t *dst, \ |
| ptrdiff_t dst_stride, const InterpKernel *filter, int x0_q4, \ |
| int x_step_q4, int y0_q4, int y_step_q4, int w, int h) { \ |
| const int16_t *filter_x = filter[x0_q4]; \ |
| const int16_t *filter_y = filter[y0_q4]; \ |
| (void)filter_y; \ |
| assert(filter_x[3] != 128); \ |
| assert(filter_y[3] != 128); \ |
| assert(w <= 64); \ |
| assert(h <= 64); \ |
| assert(x_step_q4 == 16); \ |
| assert(y_step_q4 == 16); \ |
| if (filter_x[0] | filter_x[1] | filter_x[6] | filter_x[7]) { \ |
| DECLARE_ALIGNED(16, uint8_t, fdata2[64 * 71] VPX_UNINITIALIZED); \ |
| vpx_convolve8_horiz_##opt(src - 3 * src_stride, src_stride, fdata2, 64, \ |
| filter, x0_q4, x_step_q4, y0_q4, y_step_q4, w, \ |
| h + 7); \ |
| vpx_convolve8_##avg##vert_##opt(fdata2 + 3 * 64, 64, dst, dst_stride, \ |
| filter, x0_q4, x_step_q4, y0_q4, \ |
| y_step_q4, w, h); \ |
| } else if (filter_x[2] | filter_x[5]) { \ |
| const int num_taps = is_avg ? 8 : 4; \ |
| DECLARE_ALIGNED(16, uint8_t, fdata2[64 * 71] VPX_UNINITIALIZED); \ |
| vpx_convolve8_horiz_##opt( \ |
| src - (num_taps / 2 - 1) * src_stride, src_stride, fdata2, 64, \ |
| filter, x0_q4, x_step_q4, y0_q4, y_step_q4, w, h + num_taps - 1); \ |
| vpx_convolve8_##avg##vert_##opt(fdata2 + 64 * (num_taps / 2 - 1), 64, \ |
| dst, dst_stride, filter, x0_q4, \ |
| x_step_q4, y0_q4, y_step_q4, w, h); \ |
| } else { \ |
| DECLARE_ALIGNED(16, uint8_t, fdata2[64 * 65] VPX_UNINITIALIZED); \ |
| vpx_convolve8_horiz_##opt(src, src_stride, fdata2, 64, filter, x0_q4, \ |
| x_step_q4, y0_q4, y_step_q4, w, h + 1); \ |
| vpx_convolve8_##avg##vert_##opt(fdata2, 64, dst, dst_stride, filter, \ |
| x0_q4, x_step_q4, y0_q4, y_step_q4, w, \ |
| h); \ |
| } \ |
| } |
| |
| #if CONFIG_VP9_HIGHBITDEPTH |
| |
| typedef void highbd_filter8_1dfunction(const uint16_t *src_ptr, |
| const ptrdiff_t src_pitch, |
| uint16_t *output_ptr, |
| ptrdiff_t out_pitch, |
| unsigned int output_height, |
| const int16_t *filter, int bd); |
| |
| #define HIGH_FUN_CONV_1D(name, offset, step_q4, dir, src_start, avg, opt, \ |
| is_avg) \ |
| void vpx_highbd_convolve8_##name##_##opt( \ |
| const uint16_t *src, ptrdiff_t src_stride, uint16_t *dst, \ |
| ptrdiff_t dst_stride, const InterpKernel *filter_kernel, int x0_q4, \ |
| int x_step_q4, int y0_q4, int y_step_q4, int w, int h, int bd) { \ |
| const int16_t *filter_row = filter_kernel[offset]; \ |
| if (step_q4 == 16 && filter_row[3] != 128) { \ |
| if (filter_row[0] | filter_row[1] | filter_row[6] | filter_row[7]) { \ |
| const int num_taps = 8; \ |
| while (w >= 16) { \ |
| vpx_highbd_filter_block1d16_##dir##8_##avg##opt( \ |
| src_start, src_stride, dst, dst_stride, h, filter_row, bd); \ |
| src += 16; \ |
| dst += 16; \ |
| w -= 16; \ |
| } \ |
| while (w >= 8) { \ |
| vpx_highbd_filter_block1d8_##dir##8_##avg##opt( \ |
| src_start, src_stride, dst, dst_stride, h, filter_row, bd); \ |
| src += 8; \ |
| dst += 8; \ |
| w -= 8; \ |
| } \ |
| while (w >= 4) { \ |
| vpx_highbd_filter_block1d4_##dir##8_##avg##opt( \ |
| src_start, src_stride, dst, dst_stride, h, filter_row, bd); \ |
| src += 4; \ |
| dst += 4; \ |
| w -= 4; \ |
| } \ |
| (void)num_taps; \ |
| } else if (filter_row[2] | filter_row[5]) { \ |
| const int num_taps = is_avg ? 8 : 4; \ |
| while (w >= 16) { \ |
| vpx_highbd_filter_block1d16_##dir##4_##avg##opt( \ |
| src_start, src_stride, dst, dst_stride, h, filter_row, bd); \ |
| src += 16; \ |
| dst += 16; \ |
| w -= 16; \ |
| } \ |
| while (w >= 8) { \ |
| vpx_highbd_filter_block1d8_##dir##4_##avg##opt( \ |
| src_start, src_stride, dst, dst_stride, h, filter_row, bd); \ |
| src += 8; \ |
| dst += 8; \ |
| w -= 8; \ |
| } \ |
| while (w >= 4) { \ |
| vpx_highbd_filter_block1d4_##dir##4_##avg##opt( \ |
| src_start, src_stride, dst, dst_stride, h, filter_row, bd); \ |
| src += 4; \ |
| dst += 4; \ |
| w -= 4; \ |
| } \ |
| (void)num_taps; \ |
| } else { \ |
| const int num_taps = 2; \ |
| while (w >= 16) { \ |
| vpx_highbd_filter_block1d16_##dir##2_##avg##opt( \ |
| src_start, src_stride, dst, dst_stride, h, filter_row, bd); \ |
| src += 16; \ |
| dst += 16; \ |
| w -= 16; \ |
| } \ |
| while (w >= 8) { \ |
| vpx_highbd_filter_block1d8_##dir##2_##avg##opt( \ |
| src_start, src_stride, dst, dst_stride, h, filter_row, bd); \ |
| src += 8; \ |
| dst += 8; \ |
| w -= 8; \ |
| } \ |
| while (w >= 4) { \ |
| vpx_highbd_filter_block1d4_##dir##2_##avg##opt( \ |
| src_start, src_stride, dst, dst_stride, h, filter_row, bd); \ |
| src += 4; \ |
| dst += 4; \ |
| w -= 4; \ |
| } \ |
| (void)num_taps; \ |
| } \ |
| } \ |
| if (w) { \ |
| vpx_highbd_convolve8_##name##_c(src, src_stride, dst, dst_stride, \ |
| filter_kernel, x0_q4, x_step_q4, y0_q4, \ |
| y_step_q4, w, h, bd); \ |
| } \ |
| } |
| |
| #define HIGH_FUN_CONV_2D(avg, opt, is_avg) \ |
| void vpx_highbd_convolve8_##avg##opt( \ |
| const uint16_t *src, ptrdiff_t src_stride, uint16_t *dst, \ |
| ptrdiff_t dst_stride, const InterpKernel *filter, int x0_q4, \ |
| int x_step_q4, int y0_q4, int y_step_q4, int w, int h, int bd) { \ |
| const int16_t *filter_x = filter[x0_q4]; \ |
| assert(w <= 64); \ |
| assert(h <= 64); \ |
| if (x_step_q4 == 16 && y_step_q4 == 16) { \ |
| if ((filter_x[0] | filter_x[1] | filter_x[6] | filter_x[7]) || \ |
| filter_x[3] == 128) { \ |
| DECLARE_ALIGNED(16, uint16_t, fdata2[64 * 71] VPX_UNINITIALIZED); \ |
| vpx_highbd_convolve8_horiz_##opt(src - 3 * src_stride, src_stride, \ |
| fdata2, 64, filter, x0_q4, x_step_q4, \ |
| y0_q4, y_step_q4, w, h + 7, bd); \ |
| vpx_highbd_convolve8_##avg##vert_##opt( \ |
| fdata2 + 192, 64, dst, dst_stride, filter, x0_q4, x_step_q4, \ |
| y0_q4, y_step_q4, w, h, bd); \ |
| } else if (filter_x[2] | filter_x[5]) { \ |
| const int num_taps = is_avg ? 8 : 4; \ |
| DECLARE_ALIGNED(16, uint16_t, fdata2[64 * 71] VPX_UNINITIALIZED); \ |
| vpx_highbd_convolve8_horiz_##opt( \ |
| src - (num_taps / 2 - 1) * src_stride, src_stride, fdata2, 64, \ |
| filter, x0_q4, x_step_q4, y0_q4, y_step_q4, w, h + num_taps - 1, \ |
| bd); \ |
| vpx_highbd_convolve8_##avg##vert_##opt( \ |
| fdata2 + 64 * (num_taps / 2 - 1), 64, dst, dst_stride, filter, \ |
| x0_q4, x_step_q4, y0_q4, y_step_q4, w, h, bd); \ |
| } else { \ |
| DECLARE_ALIGNED(16, uint16_t, fdata2[64 * 65] VPX_UNINITIALIZED); \ |
| vpx_highbd_convolve8_horiz_##opt(src, src_stride, fdata2, 64, filter, \ |
| x0_q4, x_step_q4, y0_q4, y_step_q4, \ |
| w, h + 1, bd); \ |
| vpx_highbd_convolve8_##avg##vert_##opt(fdata2, 64, dst, dst_stride, \ |
| filter, x0_q4, x_step_q4, \ |
| y0_q4, y_step_q4, w, h, bd); \ |
| } \ |
| } else { \ |
| vpx_highbd_convolve8_##avg##c(src, src_stride, dst, dst_stride, filter, \ |
| x0_q4, x_step_q4, y0_q4, y_step_q4, w, h, \ |
| bd); \ |
| } \ |
| } |
| |
| #endif // CONFIG_VP9_HIGHBITDEPTH |
| #endif // VPX_VPX_DSP_X86_CONVOLVE_H_ |