/* * Copyright (c) 2023, Alliance for Open Media. All rights reserved. * * This source code is subject to the terms of the BSD 2 Clause License and * the Alliance for Open Media Patent License 1.0. If the BSD 2 Clause License * was not distributed with this source code in the LICENSE file, you can * obtain it at www.aomedia.org/license/software. If the Alliance for Open * Media Patent License 1.0 was not distributed with this source code in the * PATENTS file, you can obtain it at www.aomedia.org/license/patent. */ #include #include #include "config/aom_config.h" #include "config/aom_dsp_rtcd.h" #include "aom_dsp/arm/blend_neon.h" #include "aom_dsp/arm/mem_neon.h" #include "aom_dsp/blend.h" #define HBD_BLEND_A64_D16_MASK(bd, round0_bits) \ static inline uint16x8_t alpha_##bd##_blend_a64_d16_u16x8( \ uint16x8_t m, uint16x8_t a, uint16x8_t b, int32x4_t round_offset) { \ const uint16x8_t m_inv = \ vsubq_u16(vdupq_n_u16(AOM_BLEND_A64_MAX_ALPHA), m); \ \ uint32x4_t blend_u32_lo = vmlal_u16(vreinterpretq_u32_s32(round_offset), \ vget_low_u16(m), vget_low_u16(a)); \ uint32x4_t blend_u32_hi = vmlal_u16(vreinterpretq_u32_s32(round_offset), \ vget_high_u16(m), vget_high_u16(a)); \ \ blend_u32_lo = \ vmlal_u16(blend_u32_lo, vget_low_u16(m_inv), vget_low_u16(b)); \ blend_u32_hi = \ vmlal_u16(blend_u32_hi, vget_high_u16(m_inv), vget_high_u16(b)); \ \ uint16x4_t blend_u16_lo = \ vqrshrun_n_s32(vreinterpretq_s32_u32(blend_u32_lo), \ AOM_BLEND_A64_ROUND_BITS + 2 * FILTER_BITS - \ round0_bits - COMPOUND_ROUND1_BITS); \ uint16x4_t blend_u16_hi = \ vqrshrun_n_s32(vreinterpretq_s32_u32(blend_u32_hi), \ AOM_BLEND_A64_ROUND_BITS + 2 * FILTER_BITS - \ round0_bits - COMPOUND_ROUND1_BITS); \ \ uint16x8_t blend_u16 = vcombine_u16(blend_u16_lo, blend_u16_hi); \ blend_u16 = vminq_u16(blend_u16, vdupq_n_u16((1 << bd) - 1)); \ \ return blend_u16; \ } \ \ static inline void highbd_##bd##_blend_a64_d16_mask_neon( \ uint16_t *dst, uint32_t dst_stride, const CONV_BUF_TYPE *src0, \ uint32_t src0_stride, const CONV_BUF_TYPE *src1, uint32_t src1_stride, \ const uint8_t *mask, uint32_t mask_stride, int w, int h, int subw, \ int subh) { \ const int offset_bits = bd + 2 * FILTER_BITS - round0_bits; \ int32_t round_offset = (1 << (offset_bits - COMPOUND_ROUND1_BITS)) + \ (1 << (offset_bits - COMPOUND_ROUND1_BITS - 1)); \ int32x4_t offset = \ vdupq_n_s32(-(round_offset << AOM_BLEND_A64_ROUND_BITS)); \ \ if ((subw | subh) == 0) { \ if (w >= 8) { \ do { \ int i = 0; \ do { \ uint16x8_t m0 = vmovl_u8(vld1_u8(mask + i)); \ uint16x8_t s0 = vld1q_u16(src0 + i); \ uint16x8_t s1 = vld1q_u16(src1 + i); \ \ uint16x8_t blend = \ alpha_##bd##_blend_a64_d16_u16x8(m0, s0, s1, offset); \ \ vst1q_u16(dst + i, blend); \ i += 8; \ } while (i < w); \ \ mask += mask_stride; \ src0 += src0_stride; \ src1 += src1_stride; \ dst += dst_stride; \ } while (--h != 0); \ } else { \ do { \ uint16x8_t m0 = vmovl_u8(load_unaligned_u8_4x2(mask, mask_stride)); \ uint16x8_t s0 = load_unaligned_u16_4x2(src0, src0_stride); \ uint16x8_t s1 = load_unaligned_u16_4x2(src1, src1_stride); \ \ uint16x8_t blend = \ alpha_##bd##_blend_a64_d16_u16x8(m0, s0, s1, offset); \ \ store_u16x4_strided_x2(dst, dst_stride, blend); \ \ mask += 2 * mask_stride; \ src0 += 2 * src0_stride; \ src1 += 2 * src1_stride; \ dst += 2 * dst_stride; \ h -= 2; \ } while (h != 0); \ } \ } else if ((subw & subh) == 1) { \ if (w >= 8) { \ do { \ int i = 0; \ do { \ uint8x16_t m0 = vld1q_u8(mask + 0 * mask_stride + 2 * i); \ uint8x16_t m1 = vld1q_u8(mask + 1 * mask_stride + 2 * i); \ uint16x8_t s0 = vld1q_u16(src0 + i); \ uint16x8_t s1 = vld1q_u16(src1 + i); \ \ uint16x8_t m_avg = vmovl_u8(avg_blend_pairwise_u8x8_4( \ vget_low_u8(m0), vget_low_u8(m1), vget_high_u8(m0), \ vget_high_u8(m1))); \ uint16x8_t blend = \ alpha_##bd##_blend_a64_d16_u16x8(m_avg, s0, s1, offset); \ \ vst1q_u16(dst + i, blend); \ i += 8; \ } while (i < w); \ \ mask += 2 * mask_stride; \ src0 += src0_stride; \ src1 += src1_stride; \ dst += dst_stride; \ } while (--h != 0); \ } else { \ do { \ uint8x8_t m0 = vld1_u8(mask + 0 * mask_stride); \ uint8x8_t m1 = vld1_u8(mask + 1 * mask_stride); \ uint8x8_t m2 = vld1_u8(mask + 2 * mask_stride); \ uint8x8_t m3 = vld1_u8(mask + 3 * mask_stride); \ uint16x8_t s0 = load_unaligned_u16_4x2(src0, src0_stride); \ uint16x8_t s1 = load_unaligned_u16_4x2(src1, src1_stride); \ \ uint16x8_t m_avg = \ vmovl_u8(avg_blend_pairwise_u8x8_4(m0, m1, m2, m3)); \ uint16x8_t blend = \ alpha_##bd##_blend_a64_d16_u16x8(m_avg, s0, s1, offset); \ \ store_u16x4_strided_x2(dst, dst_stride, blend); \ \ mask += 4 * mask_stride; \ src0 += 2 * src0_stride; \ src1 += 2 * src1_stride; \ dst += 2 * dst_stride; \ h -= 2; \ } while (h != 0); \ } \ } else if (subw == 1 && subh == 0) { \ if (w >= 8) { \ do { \ int i = 0; \ do { \ uint8x8_t m0 = vld1_u8(mask + 2 * i); \ uint8x8_t m1 = vld1_u8(mask + 2 * i + 8); \ uint16x8_t s0 = vld1q_u16(src0 + i); \ uint16x8_t s1 = vld1q_u16(src1 + i); \ \ uint16x8_t m_avg = vmovl_u8(avg_blend_pairwise_u8x8(m0, m1)); \ uint16x8_t blend = \ alpha_##bd##_blend_a64_d16_u16x8(m_avg, s0, s1, offset); \ \ vst1q_u16(dst + i, blend); \ i += 8; \ } while (i < w); \ \ mask += mask_stride; \ src0 += src0_stride; \ src1 += src1_stride; \ dst += dst_stride; \ } while (--h != 0); \ } else { \ do { \ uint8x8_t m0 = vld1_u8(mask + 0 * mask_stride); \ uint8x8_t m1 = vld1_u8(mask + 1 * mask_stride); \ uint16x8_t s0 = load_unaligned_u16_4x2(src0, src0_stride); \ uint16x8_t s1 = load_unaligned_u16_4x2(src1, src1_stride); \ \ uint16x8_t m_avg = vmovl_u8(avg_blend_pairwise_u8x8(m0, m1)); \ uint16x8_t blend = \ alpha_##bd##_blend_a64_d16_u16x8(m_avg, s0, s1, offset); \ \ store_u16x4_strided_x2(dst, dst_stride, blend); \ \ mask += 2 * mask_stride; \ src0 += 2 * src0_stride; \ src1 += 2 * src1_stride; \ dst += 2 * dst_stride; \ h -= 2; \ } while (h != 0); \ } \ } else { \ if (w >= 8) { \ do { \ int i = 0; \ do { \ uint8x8_t m0 = vld1_u8(mask + 0 * mask_stride + i); \ uint8x8_t m1 = vld1_u8(mask + 1 * mask_stride + i); \ uint16x8_t s0 = vld1q_u16(src0 + i); \ uint16x8_t s1 = vld1q_u16(src1 + i); \ \ uint16x8_t m_avg = vmovl_u8(avg_blend_u8x8(m0, m1)); \ uint16x8_t blend = \ alpha_##bd##_blend_a64_d16_u16x8(m_avg, s0, s1, offset); \ \ vst1q_u16(dst + i, blend); \ i += 8; \ } while (i < w); \ \ mask += 2 * mask_stride; \ src0 += src0_stride; \ src1 += src1_stride; \ dst += dst_stride; \ } while (--h != 0); \ } else { \ do { \ uint8x8_t m0_2 = \ load_unaligned_u8_4x2(mask + 0 * mask_stride, 2 * mask_stride); \ uint8x8_t m1_3 = \ load_unaligned_u8_4x2(mask + 1 * mask_stride, 2 * mask_stride); \ uint16x8_t s0 = load_unaligned_u16_4x2(src0, src0_stride); \ uint16x8_t s1 = load_unaligned_u16_4x2(src1, src1_stride); \ \ uint16x8_t m_avg = vmovl_u8(avg_blend_u8x8(m0_2, m1_3)); \ uint16x8_t blend = \ alpha_##bd##_blend_a64_d16_u16x8(m_avg, s0, s1, offset); \ \ store_u16x4_strided_x2(dst, dst_stride, blend); \ \ mask += 4 * mask_stride; \ src0 += 2 * src0_stride; \ src1 += 2 * src1_stride; \ dst += 2 * dst_stride; \ h -= 2; \ } while (h != 0); \ } \ } \ } // 12 bitdepth HBD_BLEND_A64_D16_MASK(12, (ROUND0_BITS + 2)) // 10 bitdepth HBD_BLEND_A64_D16_MASK(10, ROUND0_BITS) // 8 bitdepth HBD_BLEND_A64_D16_MASK(8, ROUND0_BITS) void aom_highbd_blend_a64_d16_mask_neon( uint8_t *dst_8, uint32_t dst_stride, const CONV_BUF_TYPE *src0, uint32_t src0_stride, const CONV_BUF_TYPE *src1, uint32_t src1_stride, const uint8_t *mask, uint32_t mask_stride, int w, int h, int subw, int subh, ConvolveParams *conv_params, const int bd) { (void)conv_params; assert(h >= 1); assert(w >= 1); assert(IS_POWER_OF_TWO(h)); assert(IS_POWER_OF_TWO(w)); uint16_t *dst = CONVERT_TO_SHORTPTR(dst_8); assert(IMPLIES(src0 == dst, src0_stride == dst_stride)); assert(IMPLIES(src1 == dst, src1_stride == dst_stride)); if (bd == 12) { highbd_12_blend_a64_d16_mask_neon(dst, dst_stride, src0, src0_stride, src1, src1_stride, mask, mask_stride, w, h, subw, subh); } else if (bd == 10) { highbd_10_blend_a64_d16_mask_neon(dst, dst_stride, src0, src0_stride, src1, src1_stride, mask, mask_stride, w, h, subw, subh); } else { highbd_8_blend_a64_d16_mask_neon(dst, dst_stride, src0, src0_stride, src1, src1_stride, mask, mask_stride, w, h, subw, subh); } } void aom_highbd_blend_a64_mask_neon(uint8_t *dst_8, uint32_t dst_stride, const uint8_t *src0_8, uint32_t src0_stride, const uint8_t *src1_8, uint32_t src1_stride, const uint8_t *mask, uint32_t mask_stride, int w, int h, int subw, int subh, int bd) { (void)bd; const uint16_t *src0 = CONVERT_TO_SHORTPTR(src0_8); const uint16_t *src1 = CONVERT_TO_SHORTPTR(src1_8); uint16_t *dst = CONVERT_TO_SHORTPTR(dst_8); assert(IMPLIES(src0 == dst, src0_stride == dst_stride)); assert(IMPLIES(src1 == dst, src1_stride == dst_stride)); assert(h >= 1); assert(w >= 1); assert(IS_POWER_OF_TWO(h)); assert(IS_POWER_OF_TWO(w)); assert(bd == 8 || bd == 10 || bd == 12); if ((subw | subh) == 0) { if (w >= 8) { do { int i = 0; do { uint16x8_t m0 = vmovl_u8(vld1_u8(mask + i)); uint16x8_t s0 = vld1q_u16(src0 + i); uint16x8_t s1 = vld1q_u16(src1 + i); uint16x8_t blend = alpha_blend_a64_u16x8(m0, s0, s1); vst1q_u16(dst + i, blend); i += 8; } while (i < w); mask += mask_stride; src0 += src0_stride; src1 += src1_stride; dst += dst_stride; } while (--h != 0); } else { do { uint16x8_t m0 = vmovl_u8(load_unaligned_u8_4x2(mask, mask_stride)); uint16x8_t s0 = load_unaligned_u16_4x2(src0, src0_stride); uint16x8_t s1 = load_unaligned_u16_4x2(src1, src1_stride); uint16x8_t blend = alpha_blend_a64_u16x8(m0, s0, s1); store_u16x4_strided_x2(dst, dst_stride, blend); mask += 2 * mask_stride; src0 += 2 * src0_stride; src1 += 2 * src1_stride; dst += 2 * dst_stride; h -= 2; } while (h != 0); } } else if ((subw & subh) == 1) { if (w >= 8) { do { int i = 0; do { uint8x8_t m0 = vld1_u8(mask + 0 * mask_stride + 2 * i); uint8x8_t m1 = vld1_u8(mask + 1 * mask_stride + 2 * i); uint8x8_t m2 = vld1_u8(mask + 0 * mask_stride + 2 * i + 8); uint8x8_t m3 = vld1_u8(mask + 1 * mask_stride + 2 * i + 8); uint16x8_t s0 = vld1q_u16(src0 + i); uint16x8_t s1 = vld1q_u16(src1 + i); uint16x8_t m_avg = vmovl_u8(avg_blend_pairwise_u8x8_4(m0, m1, m2, m3)); uint16x8_t blend = alpha_blend_a64_u16x8(m_avg, s0, s1); vst1q_u16(dst + i, blend); i += 8; } while (i < w); mask += 2 * mask_stride; src0 += src0_stride; src1 += src1_stride; dst += dst_stride; } while (--h != 0); } else { do { uint8x8_t m0 = vld1_u8(mask + 0 * mask_stride); uint8x8_t m1 = vld1_u8(mask + 1 * mask_stride); uint8x8_t m2 = vld1_u8(mask + 2 * mask_stride); uint8x8_t m3 = vld1_u8(mask + 3 * mask_stride); uint16x8_t s0 = load_unaligned_u16_4x2(src0, src0_stride); uint16x8_t s1 = load_unaligned_u16_4x2(src1, src1_stride); uint16x8_t m_avg = vmovl_u8(avg_blend_pairwise_u8x8_4(m0, m1, m2, m3)); uint16x8_t blend = alpha_blend_a64_u16x8(m_avg, s0, s1); store_u16x4_strided_x2(dst, dst_stride, blend); mask += 4 * mask_stride; src0 += 2 * src0_stride; src1 += 2 * src1_stride; dst += 2 * dst_stride; h -= 2; } while (h != 0); } } else if (subw == 1 && subh == 0) { if (w >= 8) { do { int i = 0; do { uint8x8_t m0 = vld1_u8(mask + 2 * i); uint8x8_t m1 = vld1_u8(mask + 2 * i + 8); uint16x8_t s0 = vld1q_u16(src0 + i); uint16x8_t s1 = vld1q_u16(src1 + i); uint16x8_t m_avg = vmovl_u8(avg_blend_pairwise_u8x8(m0, m1)); uint16x8_t blend = alpha_blend_a64_u16x8(m_avg, s0, s1); vst1q_u16(dst + i, blend); i += 8; } while (i < w); mask += mask_stride; src0 += src0_stride; src1 += src1_stride; dst += dst_stride; } while (--h != 0); } else { do { uint8x8_t m0 = vld1_u8(mask + 0 * mask_stride); uint8x8_t m1 = vld1_u8(mask + 1 * mask_stride); uint16x8_t s0 = load_unaligned_u16_4x2(src0, src0_stride); uint16x8_t s1 = load_unaligned_u16_4x2(src1, src1_stride); uint16x8_t m_avg = vmovl_u8(avg_blend_pairwise_u8x8(m0, m1)); uint16x8_t blend = alpha_blend_a64_u16x8(m_avg, s0, s1); store_u16x4_strided_x2(dst, dst_stride, blend); mask += 2 * mask_stride; src0 += 2 * src0_stride; src1 += 2 * src1_stride; dst += 2 * dst_stride; h -= 2; } while (h != 0); } } else { if (w >= 8) { do { int i = 0; do { uint8x8_t m0 = vld1_u8(mask + 0 * mask_stride + i); uint8x8_t m1 = vld1_u8(mask + 1 * mask_stride + i); uint16x8_t s0 = vld1q_u16(src0 + i); uint16x8_t s1 = vld1q_u16(src1 + i); uint16x8_t m_avg = vmovl_u8(avg_blend_u8x8(m0, m1)); uint16x8_t blend = alpha_blend_a64_u16x8(m_avg, s0, s1); vst1q_u16(dst + i, blend); i += 8; } while (i < w); mask += 2 * mask_stride; src0 += src0_stride; src1 += src1_stride; dst += dst_stride; } while (--h != 0); } else { do { uint8x8_t m0_2 = load_unaligned_u8_4x2(mask + 0 * mask_stride, 2 * mask_stride); uint8x8_t m1_3 = load_unaligned_u8_4x2(mask + 1 * mask_stride, 2 * mask_stride); uint16x8_t s0 = load_unaligned_u16_4x2(src0, src0_stride); uint16x8_t s1 = load_unaligned_u16_4x2(src1, src1_stride); uint16x8_t m_avg = vmovl_u8(avg_blend_u8x8(m0_2, m1_3)); uint16x8_t blend = alpha_blend_a64_u16x8(m_avg, s0, s1); store_u16x4_strided_x2(dst, dst_stride, blend); mask += 4 * mask_stride; src0 += 2 * src0_stride; src1 += 2 * src1_stride; dst += 2 * dst_stride; h -= 2; } while (h != 0); } } }