ref: 42ce25821d628560dd7b5c7af875697ab5111cf0
dir: /vpx_dsp/ppc/sad_vsx.c/
/* * Copyright (c) 2017 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 <stdlib.h> #include "vpx_dsp/ppc/types_vsx.h" #include "vpx/vpx_integer.h" #define PROCESS16(offset) \ v_a = vec_vsx_ld(offset, a); \ v_b = vec_vsx_ld(offset, b); \ v_ah = unpack_to_s16_h(v_a); \ v_al = unpack_to_s16_l(v_a); \ v_bh = unpack_to_s16_h(v_b); \ v_bl = unpack_to_s16_l(v_b); \ v_subh = vec_sub(v_ah, v_bh); \ v_subl = vec_sub(v_al, v_bl); \ v_absh = vec_abs(v_subh); \ v_absl = vec_abs(v_subl); \ v_sad = vec_sum4s(v_absh, v_sad); \ v_sad = vec_sum4s(v_absl, v_sad); #define SAD16(height) \ unsigned int vpx_sad16x##height##_vsx(const uint8_t *a, int a_stride, \ const uint8_t *b, int b_stride) { \ int y; \ unsigned int sad[4]; \ uint8x16_t v_a, v_b; \ int16x8_t v_ah, v_al, v_bh, v_bl, v_absh, v_absl, v_subh, v_subl; \ int32x4_t v_sad = vec_splat_s32(0); \ \ for (y = 0; y < height; y++) { \ PROCESS16(0); \ \ a += a_stride; \ b += b_stride; \ } \ vec_vsx_st((uint32x4_t)v_sad, 0, sad); \ \ return sad[3] + sad[2] + sad[1] + sad[0]; \ } #define SAD32(height) \ unsigned int vpx_sad32x##height##_vsx(const uint8_t *a, int a_stride, \ const uint8_t *b, int b_stride) { \ int y; \ unsigned int sad[4]; \ uint8x16_t v_a, v_b; \ int16x8_t v_ah, v_al, v_bh, v_bl, v_absh, v_absl, v_subh, v_subl; \ int32x4_t v_sad = vec_splat_s32(0); \ \ for (y = 0; y < height; y++) { \ PROCESS16(0); \ PROCESS16(16); \ \ a += a_stride; \ b += b_stride; \ } \ vec_vsx_st((uint32x4_t)v_sad, 0, sad); \ \ return sad[3] + sad[2] + sad[1] + sad[0]; \ } #define SAD64(height) \ unsigned int vpx_sad64x##height##_vsx(const uint8_t *a, int a_stride, \ const uint8_t *b, int b_stride) { \ int y; \ unsigned int sad[4]; \ uint8x16_t v_a, v_b; \ int16x8_t v_ah, v_al, v_bh, v_bl, v_absh, v_absl, v_subh, v_subl; \ int32x4_t v_sad = vec_splat_s32(0); \ \ for (y = 0; y < height; y++) { \ PROCESS16(0); \ PROCESS16(16); \ PROCESS16(32); \ PROCESS16(48); \ \ a += a_stride; \ b += b_stride; \ } \ vec_vsx_st((uint32x4_t)v_sad, 0, sad); \ \ return sad[3] + sad[2] + sad[1] + sad[0]; \ } SAD16(8); SAD16(16); SAD16(32); SAD32(16); SAD32(32); SAD32(64); SAD64(32); SAD64(64);