blob: 4639d4c41690eec22297481258b5a294b62d367f [file] [log] [blame]
/*
*
* Copyright (c) 2018, 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 <arm_neon.h>
#include <assert.h>
#include "aom/aom_integer.h"
#include "aom_dsp/aom_dsp_common.h"
#include "aom_dsp/blend.h"
#include "aom_dsp/arm/mem_neon.h"
#include "aom_ports/mem.h"
#include "config/aom_dsp_rtcd.h"
void aom_blend_a64_hmask_neon(uint8_t *dst, uint32_t dst_stride,
const uint8_t *src0, uint32_t src0_stride,
const uint8_t *src1, uint32_t src1_stride,
const uint8_t *mask, int w, int h) {
assert(IMPLIES(src0 == dst, src0_stride == dst_stride));
assert(IMPLIES(src1 == dst, src1_stride == dst_stride));
assert(h >= 2);
assert(w >= 2);
assert(IS_POWER_OF_TWO(h));
assert(IS_POWER_OF_TWO(w));
uint8x8_t tmp0, tmp1;
uint8x16_t res_q;
uint16x8_t res, res_low, res_high;
uint32x2_t tmp0_32 = vdup_n_u32(0), tmp1_32 = vdup_n_u32(0);
uint16x4_t tmp0_16 = vdup_n_u16(0), tmp1_16 = vdup_n_u16(0);
const uint8x8_t vdup_64 = vdup_n_u8((uint8_t)64);
if (w >= 16) {
const uint8x16_t vdup_64_q = vdupq_n_u8((uint8_t)64);
for (int i = 0; i < h; ++i) {
for (int j = 0; j < w; j += 16) {
__builtin_prefetch(src0);
__builtin_prefetch(src1);
const uint8x16_t tmp0_q = vld1q_u8(src0);
const uint8x16_t tmp1_q = vld1q_u8(src1);
const uint8x16_t m_q = vld1q_u8(mask);
const uint8x16_t max_minus_m_q = vsubq_u8(vdup_64_q, m_q);
res_low = vmull_u8(vget_low_u8(m_q), vget_low_u8(tmp0_q));
res_low =
vmlal_u8(res_low, vget_low_u8(max_minus_m_q), vget_low_u8(tmp1_q));
res_high = vmull_u8(vget_high_u8(m_q), vget_high_u8(tmp0_q));
res_high = vmlal_u8(res_high, vget_high_u8(max_minus_m_q),
vget_high_u8(tmp1_q));
res_q = vcombine_u8(vrshrn_n_u16(res_low, AOM_BLEND_A64_ROUND_BITS),
vrshrn_n_u16(res_high, AOM_BLEND_A64_ROUND_BITS));
vst1q_u8(dst, res_q);
src0 += 16;
src1 += 16;
dst += 16;
mask += 16;
}
src0 += src0_stride - w;
src1 += src1_stride - w;
dst += dst_stride - w;
mask -= w;
}
} else if (w == 8) {
const uint8x8_t m = vld1_u8(mask);
const uint8x8_t max_minus_m = vsub_u8(vdup_64, m);
for (int i = 0; i < h; ++i) {
__builtin_prefetch(src0);
__builtin_prefetch(src1);
tmp0 = vld1_u8(src0);
tmp1 = vld1_u8(src1);
res = vmull_u8(m, tmp0);
res = vmlal_u8(res, max_minus_m, tmp1);
vst1_u8(dst, vrshrn_n_u16(res, AOM_BLEND_A64_ROUND_BITS));
src0 += src0_stride;
src1 += src1_stride;
dst += dst_stride;
}
} else if (w == 4) {
const uint8x8_t m = vreinterpret_u8_u32(vld1_dup_u32((uint32_t *)mask));
const uint8x8_t max_minus_m = vsub_u8(vdup_64, m);
for (int i = 0; i < h; i += 2) {
__builtin_prefetch(src0 + 0 * src0_stride);
__builtin_prefetch(src0 + 1 * src0_stride);
__builtin_prefetch(src1 + 0 * src1_stride);
__builtin_prefetch(src1 + 1 * src1_stride);
load_unaligned_u8_4x2(src0, src0_stride, &tmp0_32);
tmp0 = vreinterpret_u8_u32(tmp0_32);
load_unaligned_u8_4x2(src1, src1_stride, &tmp1_32);
tmp1 = vreinterpret_u8_u32(tmp1_32);
res = vmull_u8(m, tmp0);
res = vmlal_u8(res, max_minus_m, tmp1);
vst1_lane_u32(
(uint32_t *)(dst + (0 * dst_stride)),
vreinterpret_u32_u8(vrshrn_n_u16(res, AOM_BLEND_A64_ROUND_BITS)), 0);
vst1_lane_u32(
(uint32_t *)(dst + (1 * dst_stride)),
vreinterpret_u32_u8(vrshrn_n_u16(res, AOM_BLEND_A64_ROUND_BITS)), 1);
src0 += (2 * src0_stride);
src1 += (2 * src1_stride);
dst += (2 * dst_stride);
}
} else if (w == 2) {
const uint8x8_t m = vreinterpret_u8_u16(vld1_dup_u16((uint16_t *)mask));
const uint8x8_t max_minus_m = vsub_u8(vdup_64, m);
for (int i = 0; i < h; i += 2) {
__builtin_prefetch(src0 + 0 * src0_stride);
__builtin_prefetch(src0 + 1 * src0_stride);
__builtin_prefetch(src1 + 0 * src1_stride);
__builtin_prefetch(src1 + 1 * src1_stride);
load_unaligned_u8_2x2(src0, src0_stride, &tmp0_16);
tmp0 = vreinterpret_u8_u16(tmp0_16);
load_unaligned_u8_2x2(src1, src1_stride, &tmp1_16);
tmp1 = vreinterpret_u8_u16(tmp1_16);
res = vmull_u8(m, tmp0);
res = vmlal_u8(res, max_minus_m, tmp1);
vst1_lane_u16(
(uint16_t *)(dst + (0 * dst_stride)),
vreinterpret_u16_u8(vrshrn_n_u16(res, AOM_BLEND_A64_ROUND_BITS)), 0);
vst1_lane_u16(
(uint16_t *)(dst + (1 * dst_stride)),
vreinterpret_u16_u8(vrshrn_n_u16(res, AOM_BLEND_A64_ROUND_BITS)), 1);
src0 += (2 * src0_stride);
src1 += (2 * src1_stride);
dst += (2 * dst_stride);
}
}
}