test_vreinterpret_f64_s8 ()
{
int8x8_t a;
float64x1_t b;
int8_t c[8] = { 0x18, 0x2D, 0x44, 0x54, 0xFB, 0x21, 0x09, 0x40 };
float64_t d[1] = { PI_F64 };
float64_t e[1];
int i;
a = vld1_s8 (c);
b = wrap_vreinterpret_f64_s8 (a);
vst1_f64 (e, b);
if (!DOUBLE_EQUALS (d[0], e[0], __DBL_EPSILON__))
return 1;
return 0;
};
test_vdupb_lane_s8 ()
{
int8x8_t a;
int8_t b;
int8_t c[8] = { 0, 1, 2, 3, 4, 5, 6, 7 };
a = vld1_s8 (c);
b = wrap_vdupb_lane_s8_0 (a, a);
if (c[0] != b)
return 1;
b = wrap_vdupb_lane_s8_1 (a);
if (c[1] != b)
return 1;
return 0;
}
test_vld1_vst1 ()
{
int8x8_t a;
int8x8_t b;
int i = 0;
int8_t c[8] = { 0, 1, 2, 3, 4, 5, 6, 7 };
int8_t d[8];
a = vld1_s8 (c);
asm volatile ("":::"memory");
vst1_s8 (d, a);
asm volatile ("":::"memory");
for (; i < 8; i++)
if (c[i] != d[i])
return 1;
return 0;
}
test_vdup_lane_s8 ()
{
int8x8_t a;
int8x8_t b;
int i;
/* Only two first cases are interesting. */
int8_t c[8] = { 0, 1, 2, 3, 4, 5, 6, 7 };
int8_t d[8];
a = vld1_s8 (c);
b = wrap_vdup_lane_s8_0 (a);
vst1_s8 (d, b);
for (i = 0; i < 8; i++)
if (c[0] != d[i])
return 1;
b = wrap_vdup_lane_s8_1 (a);
vst1_s8 (d, b);
for (i = 0; i < 8; i++)
if (c[1] != d[i])
return 1;
return 0;
}
inline int8x8_t vld1(const s8 * ptr) { return vld1_s8(ptr); }
void test_vld1s8 (void)
{
int8x8_t out_int8x8_t;
out_int8x8_t = vld1_s8 (0);
}
void vp8_sixtap_predict16x16_neon(
unsigned char *src_ptr,
int src_pixels_per_line,
int xoffset,
int yoffset,
unsigned char *dst_ptr,
int dst_pitch) {
unsigned char *src, *src_tmp, *dst, *tmpp;
unsigned char tmp[336];
int i, j;
uint8x8_t d0u8, d1u8, d2u8, d3u8, d4u8, d5u8, d6u8, d7u8, d8u8, d9u8;
uint8x8_t d10u8, d11u8, d12u8, d13u8, d14u8, d15u8, d18u8, d19u8;
uint8x8_t d20u8, d21u8, d22u8, d23u8, d24u8, d25u8, d26u8, d27u8;
uint8x8_t d28u8, d29u8, d30u8, d31u8;
int8x8_t dtmps8, d0s8, d1s8, d2s8, d3s8, d4s8, d5s8;
uint8x16_t q3u8, q4u8;
uint16x8_t q3u16, q4u16, q5u16, q6u16, q7u16, q8u16, q9u16, q10u16;
uint16x8_t q11u16, q12u16, q13u16, q15u16;
int16x8_t q3s16, q4s16, q5s16, q6s16, q7s16, q8s16, q9s16, q10s16;
int16x8_t q11s16, q12s16, q13s16, q15s16;
if (xoffset == 0) { // secondpass_filter8x8_only
// load second_pass filter
dtmps8 = vld1_s8(vp8_sub_pel_filters[yoffset]);
d0s8 = vdup_lane_s8(dtmps8, 0);
d1s8 = vdup_lane_s8(dtmps8, 1);
d2s8 = vdup_lane_s8(dtmps8, 2);
d3s8 = vdup_lane_s8(dtmps8, 3);
d4s8 = vdup_lane_s8(dtmps8, 4);
d5s8 = vdup_lane_s8(dtmps8, 5);
d0u8 = vreinterpret_u8_s8(vabs_s8(d0s8));
d1u8 = vreinterpret_u8_s8(vabs_s8(d1s8));
d2u8 = vreinterpret_u8_s8(vabs_s8(d2s8));
d3u8 = vreinterpret_u8_s8(vabs_s8(d3s8));
d4u8 = vreinterpret_u8_s8(vabs_s8(d4s8));
d5u8 = vreinterpret_u8_s8(vabs_s8(d5s8));
// load src data
src_tmp = src_ptr - src_pixels_per_line * 2;
for (i = 0; i < 2; i++) {
src = src_tmp + i * 8;
dst = dst_ptr + i * 8;
d18u8 = vld1_u8(src);
src += src_pixels_per_line;
d19u8 = vld1_u8(src);
src += src_pixels_per_line;
d20u8 = vld1_u8(src);
src += src_pixels_per_line;
d21u8 = vld1_u8(src);
src += src_pixels_per_line;
d22u8 = vld1_u8(src);
src += src_pixels_per_line;
for (j = 0; j < 4; j++) {
d23u8 = vld1_u8(src);
src += src_pixels_per_line;
d24u8 = vld1_u8(src);
src += src_pixels_per_line;
d25u8 = vld1_u8(src);
src += src_pixels_per_line;
d26u8 = vld1_u8(src);
src += src_pixels_per_line;
q3u16 = vmull_u8(d18u8, d0u8);
q4u16 = vmull_u8(d19u8, d0u8);
q5u16 = vmull_u8(d20u8, d0u8);
q6u16 = vmull_u8(d21u8, d0u8);
q3u16 = vmlsl_u8(q3u16, d19u8, d1u8);
q4u16 = vmlsl_u8(q4u16, d20u8, d1u8);
q5u16 = vmlsl_u8(q5u16, d21u8, d1u8);
q6u16 = vmlsl_u8(q6u16, d22u8, d1u8);
q3u16 = vmlsl_u8(q3u16, d22u8, d4u8);
q4u16 = vmlsl_u8(q4u16, d23u8, d4u8);
q5u16 = vmlsl_u8(q5u16, d24u8, d4u8);
q6u16 = vmlsl_u8(q6u16, d25u8, d4u8);
q3u16 = vmlal_u8(q3u16, d20u8, d2u8);
q4u16 = vmlal_u8(q4u16, d21u8, d2u8);
q5u16 = vmlal_u8(q5u16, d22u8, d2u8);
q6u16 = vmlal_u8(q6u16, d23u8, d2u8);
q3u16 = vmlal_u8(q3u16, d23u8, d5u8);
q4u16 = vmlal_u8(q4u16, d24u8, d5u8);
q5u16 = vmlal_u8(q5u16, d25u8, d5u8);
q6u16 = vmlal_u8(q6u16, d26u8, d5u8);
q7u16 = vmull_u8(d21u8, d3u8);
q8u16 = vmull_u8(d22u8, d3u8);
q9u16 = vmull_u8(d23u8, d3u8);
q10u16 = vmull_u8(d24u8, d3u8);
q3s16 = vreinterpretq_s16_u16(q3u16);
q4s16 = vreinterpretq_s16_u16(q4u16);
q5s16 = vreinterpretq_s16_u16(q5u16);
q6s16 = vreinterpretq_s16_u16(q6u16);
q7s16 = vreinterpretq_s16_u16(q7u16);
q8s16 = vreinterpretq_s16_u16(q8u16);
q9s16 = vreinterpretq_s16_u16(q9u16);
q10s16 = vreinterpretq_s16_u16(q10u16);
q7s16 = vqaddq_s16(q7s16, q3s16);
q8s16 = vqaddq_s16(q8s16, q4s16);
q9s16 = vqaddq_s16(q9s16, q5s16);
q10s16 = vqaddq_s16(q10s16, q6s16);
d6u8 = vqrshrun_n_s16(q7s16, 7);
d7u8 = vqrshrun_n_s16(q8s16, 7);
d8u8 = vqrshrun_n_s16(q9s16, 7);
d9u8 = vqrshrun_n_s16(q10s16, 7);
d18u8 = d22u8;
d19u8 = d23u8;
d20u8 = d24u8;
d21u8 = d25u8;
d22u8 = d26u8;
vst1_u8(dst, d6u8);
dst += dst_pitch;
vst1_u8(dst, d7u8);
dst += dst_pitch;
vst1_u8(dst, d8u8);
dst += dst_pitch;
vst1_u8(dst, d9u8);
dst += dst_pitch;
}
}
return;
}
// load first_pass filter
dtmps8 = vld1_s8(vp8_sub_pel_filters[xoffset]);
d0s8 = vdup_lane_s8(dtmps8, 0);
d1s8 = vdup_lane_s8(dtmps8, 1);
d2s8 = vdup_lane_s8(dtmps8, 2);
d3s8 = vdup_lane_s8(dtmps8, 3);
d4s8 = vdup_lane_s8(dtmps8, 4);
d5s8 = vdup_lane_s8(dtmps8, 5);
d0u8 = vreinterpret_u8_s8(vabs_s8(d0s8));
d1u8 = vreinterpret_u8_s8(vabs_s8(d1s8));
d2u8 = vreinterpret_u8_s8(vabs_s8(d2s8));
d3u8 = vreinterpret_u8_s8(vabs_s8(d3s8));
d4u8 = vreinterpret_u8_s8(vabs_s8(d4s8));
d5u8 = vreinterpret_u8_s8(vabs_s8(d5s8));
// First pass: output_height lines x output_width columns (9x4)
if (yoffset == 0) { // firstpass_filter4x4_only
src = src_ptr - 2;
dst = dst_ptr;
for (i = 0; i < 8; i++) {
d6u8 = vld1_u8(src);
d7u8 = vld1_u8(src + 8);
d8u8 = vld1_u8(src + 16);
src += src_pixels_per_line;
d9u8 = vld1_u8(src);
d10u8 = vld1_u8(src + 8);
d11u8 = vld1_u8(src + 16);
src += src_pixels_per_line;
__builtin_prefetch(src);
__builtin_prefetch(src + src_pixels_per_line);
q6u16 = vmull_u8(d6u8, d0u8);
q7u16 = vmull_u8(d7u8, d0u8);
q8u16 = vmull_u8(d9u8, d0u8);
q9u16 = vmull_u8(d10u8, d0u8);
d20u8 = vext_u8(d6u8, d7u8, 1);
d21u8 = vext_u8(d9u8, d10u8, 1);
d22u8 = vext_u8(d7u8, d8u8, 1);
d23u8 = vext_u8(d10u8, d11u8, 1);
d24u8 = vext_u8(d6u8, d7u8, 4);
d25u8 = vext_u8(d9u8, d10u8, 4);
d26u8 = vext_u8(d7u8, d8u8, 4);
d27u8 = vext_u8(d10u8, d11u8, 4);
d28u8 = vext_u8(d6u8, d7u8, 5);
d29u8 = vext_u8(d9u8, d10u8, 5);
q6u16 = vmlsl_u8(q6u16, d20u8, d1u8);
q8u16 = vmlsl_u8(q8u16, d21u8, d1u8);
q7u16 = vmlsl_u8(q7u16, d22u8, d1u8);
q9u16 = vmlsl_u8(q9u16, d23u8, d1u8);
q6u16 = vmlsl_u8(q6u16, d24u8, d4u8);
q8u16 = vmlsl_u8(q8u16, d25u8, d4u8);
q7u16 = vmlsl_u8(q7u16, d26u8, d4u8);
q9u16 = vmlsl_u8(q9u16, d27u8, d4u8);
q6u16 = vmlal_u8(q6u16, d28u8, d5u8);
q8u16 = vmlal_u8(q8u16, d29u8, d5u8);
d20u8 = vext_u8(d7u8, d8u8, 5);
d21u8 = vext_u8(d10u8, d11u8, 5);
d22u8 = vext_u8(d6u8, d7u8, 2);
d23u8 = vext_u8(d9u8, d10u8, 2);
d24u8 = vext_u8(d7u8, d8u8, 2);
d25u8 = vext_u8(d10u8, d11u8, 2);
d26u8 = vext_u8(d6u8, d7u8, 3);
d27u8 = vext_u8(d9u8, d10u8, 3);
d28u8 = vext_u8(d7u8, d8u8, 3);
d29u8 = vext_u8(d10u8, d11u8, 3);
q7u16 = vmlal_u8(q7u16, d20u8, d5u8);
q9u16 = vmlal_u8(q9u16, d21u8, d5u8);
q6u16 = vmlal_u8(q6u16, d22u8, d2u8);
q8u16 = vmlal_u8(q8u16, d23u8, d2u8);
q7u16 = vmlal_u8(q7u16, d24u8, d2u8);
q9u16 = vmlal_u8(q9u16, d25u8, d2u8);
q10u16 = vmull_u8(d26u8, d3u8);
q11u16 = vmull_u8(d27u8, d3u8);
q12u16 = vmull_u8(d28u8, d3u8);
q15u16 = vmull_u8(d29u8, d3u8);
q6s16 = vreinterpretq_s16_u16(q6u16);
q7s16 = vreinterpretq_s16_u16(q7u16);
q8s16 = vreinterpretq_s16_u16(q8u16);
q9s16 = vreinterpretq_s16_u16(q9u16);
q10s16 = vreinterpretq_s16_u16(q10u16);
q11s16 = vreinterpretq_s16_u16(q11u16);
q12s16 = vreinterpretq_s16_u16(q12u16);
q15s16 = vreinterpretq_s16_u16(q15u16);
q6s16 = vqaddq_s16(q6s16, q10s16);
q8s16 = vqaddq_s16(q8s16, q11s16);
q7s16 = vqaddq_s16(q7s16, q12s16);
q9s16 = vqaddq_s16(q9s16, q15s16);
d6u8 = vqrshrun_n_s16(q6s16, 7);
d7u8 = vqrshrun_n_s16(q7s16, 7);
d8u8 = vqrshrun_n_s16(q8s16, 7);
d9u8 = vqrshrun_n_s16(q9s16, 7);
q3u8 = vcombine_u8(d6u8, d7u8);
q4u8 = vcombine_u8(d8u8, d9u8);
vst1q_u8(dst, q3u8);
dst += dst_pitch;
vst1q_u8(dst, q4u8);
dst += dst_pitch;
}
return;
}
src = src_ptr - 2 - src_pixels_per_line * 2;
tmpp = tmp;
for (i = 0; i < 7; i++) {
d6u8 = vld1_u8(src);
d7u8 = vld1_u8(src + 8);
d8u8 = vld1_u8(src + 16);
src += src_pixels_per_line;
d9u8 = vld1_u8(src);
d10u8 = vld1_u8(src + 8);
d11u8 = vld1_u8(src + 16);
src += src_pixels_per_line;
d12u8 = vld1_u8(src);
d13u8 = vld1_u8(src + 8);
d14u8 = vld1_u8(src + 16);
src += src_pixels_per_line;
__builtin_prefetch(src);
__builtin_prefetch(src + src_pixels_per_line);
__builtin_prefetch(src + src_pixels_per_line * 2);
q8u16 = vmull_u8(d6u8, d0u8);
q9u16 = vmull_u8(d7u8, d0u8);
q10u16 = vmull_u8(d9u8, d0u8);
q11u16 = vmull_u8(d10u8, d0u8);
q12u16 = vmull_u8(d12u8, d0u8);
q13u16 = vmull_u8(d13u8, d0u8);
d28u8 = vext_u8(d6u8, d7u8, 1);
d29u8 = vext_u8(d9u8, d10u8, 1);
d30u8 = vext_u8(d12u8, d13u8, 1);
q8u16 = vmlsl_u8(q8u16, d28u8, d1u8);
q10u16 = vmlsl_u8(q10u16, d29u8, d1u8);
q12u16 = vmlsl_u8(q12u16, d30u8, d1u8);
d28u8 = vext_u8(d7u8, d8u8, 1);
d29u8 = vext_u8(d10u8, d11u8, 1);
d30u8 = vext_u8(d13u8, d14u8, 1);
q9u16 = vmlsl_u8(q9u16, d28u8, d1u8);
q11u16 = vmlsl_u8(q11u16, d29u8, d1u8);
q13u16 = vmlsl_u8(q13u16, d30u8, d1u8);
d28u8 = vext_u8(d6u8, d7u8, 4);
d29u8 = vext_u8(d9u8, d10u8, 4);
d30u8 = vext_u8(d12u8, d13u8, 4);
q8u16 = vmlsl_u8(q8u16, d28u8, d4u8);
q10u16 = vmlsl_u8(q10u16, d29u8, d4u8);
q12u16 = vmlsl_u8(q12u16, d30u8, d4u8);
d28u8 = vext_u8(d7u8, d8u8, 4);
d29u8 = vext_u8(d10u8, d11u8, 4);
d30u8 = vext_u8(d13u8, d14u8, 4);
q9u16 = vmlsl_u8(q9u16, d28u8, d4u8);
q11u16 = vmlsl_u8(q11u16, d29u8, d4u8);
q13u16 = vmlsl_u8(q13u16, d30u8, d4u8);
d28u8 = vext_u8(d6u8, d7u8, 5);
d29u8 = vext_u8(d9u8, d10u8, 5);
d30u8 = vext_u8(d12u8, d13u8, 5);
q8u16 = vmlal_u8(q8u16, d28u8, d5u8);
q10u16 = vmlal_u8(q10u16, d29u8, d5u8);
q12u16 = vmlal_u8(q12u16, d30u8, d5u8);
d28u8 = vext_u8(d7u8, d8u8, 5);
d29u8 = vext_u8(d10u8, d11u8, 5);
d30u8 = vext_u8(d13u8, d14u8, 5);
q9u16 = vmlal_u8(q9u16, d28u8, d5u8);
q11u16 = vmlal_u8(q11u16, d29u8, d5u8);
q13u16 = vmlal_u8(q13u16, d30u8, d5u8);
d28u8 = vext_u8(d6u8, d7u8, 2);
d29u8 = vext_u8(d9u8, d10u8, 2);
d30u8 = vext_u8(d12u8, d13u8, 2);
q8u16 = vmlal_u8(q8u16, d28u8, d2u8);
q10u16 = vmlal_u8(q10u16, d29u8, d2u8);
q12u16 = vmlal_u8(q12u16, d30u8, d2u8);
d28u8 = vext_u8(d7u8, d8u8, 2);
d29u8 = vext_u8(d10u8, d11u8, 2);
d30u8 = vext_u8(d13u8, d14u8, 2);
q9u16 = vmlal_u8(q9u16, d28u8, d2u8);
q11u16 = vmlal_u8(q11u16, d29u8, d2u8);
q13u16 = vmlal_u8(q13u16, d30u8, d2u8);
d28u8 = vext_u8(d6u8, d7u8, 3);
d29u8 = vext_u8(d9u8, d10u8, 3);
d30u8 = vext_u8(d12u8, d13u8, 3);
d15u8 = vext_u8(d7u8, d8u8, 3);
d31u8 = vext_u8(d10u8, d11u8, 3);
d6u8 = vext_u8(d13u8, d14u8, 3);
q4u16 = vmull_u8(d28u8, d3u8);
q5u16 = vmull_u8(d29u8, d3u8);
q6u16 = vmull_u8(d30u8, d3u8);
q4s16 = vreinterpretq_s16_u16(q4u16);
q5s16 = vreinterpretq_s16_u16(q5u16);
q6s16 = vreinterpretq_s16_u16(q6u16);
q8s16 = vreinterpretq_s16_u16(q8u16);
q10s16 = vreinterpretq_s16_u16(q10u16);
q12s16 = vreinterpretq_s16_u16(q12u16);
q8s16 = vqaddq_s16(q8s16, q4s16);
q10s16 = vqaddq_s16(q10s16, q5s16);
q12s16 = vqaddq_s16(q12s16, q6s16);
q6u16 = vmull_u8(d15u8, d3u8);
q7u16 = vmull_u8(d31u8, d3u8);
q3u16 = vmull_u8(d6u8, d3u8);
q3s16 = vreinterpretq_s16_u16(q3u16);
q6s16 = vreinterpretq_s16_u16(q6u16);
q7s16 = vreinterpretq_s16_u16(q7u16);
q9s16 = vreinterpretq_s16_u16(q9u16);
q11s16 = vreinterpretq_s16_u16(q11u16);
q13s16 = vreinterpretq_s16_u16(q13u16);
q9s16 = vqaddq_s16(q9s16, q6s16);
q11s16 = vqaddq_s16(q11s16, q7s16);
q13s16 = vqaddq_s16(q13s16, q3s16);
d6u8 = vqrshrun_n_s16(q8s16, 7);
d7u8 = vqrshrun_n_s16(q9s16, 7);
d8u8 = vqrshrun_n_s16(q10s16, 7);
d9u8 = vqrshrun_n_s16(q11s16, 7);
d10u8 = vqrshrun_n_s16(q12s16, 7);
d11u8 = vqrshrun_n_s16(q13s16, 7);
vst1_u8(tmpp, d6u8);
tmpp += 8;
vst1_u8(tmpp, d7u8);
tmpp += 8;
vst1_u8(tmpp, d8u8);
tmpp += 8;
vst1_u8(tmpp, d9u8);
tmpp += 8;
vst1_u8(tmpp, d10u8);
tmpp += 8;
vst1_u8(tmpp, d11u8);
tmpp += 8;
}
// Second pass: 16x16
dtmps8 = vld1_s8(vp8_sub_pel_filters[yoffset]);
d0s8 = vdup_lane_s8(dtmps8, 0);
d1s8 = vdup_lane_s8(dtmps8, 1);
d2s8 = vdup_lane_s8(dtmps8, 2);
d3s8 = vdup_lane_s8(dtmps8, 3);
d4s8 = vdup_lane_s8(dtmps8, 4);
d5s8 = vdup_lane_s8(dtmps8, 5);
d0u8 = vreinterpret_u8_s8(vabs_s8(d0s8));
d1u8 = vreinterpret_u8_s8(vabs_s8(d1s8));
d2u8 = vreinterpret_u8_s8(vabs_s8(d2s8));
d3u8 = vreinterpret_u8_s8(vabs_s8(d3s8));
d4u8 = vreinterpret_u8_s8(vabs_s8(d4s8));
d5u8 = vreinterpret_u8_s8(vabs_s8(d5s8));
for (i = 0; i < 2; i++) {
dst = dst_ptr + 8 * i;
tmpp = tmp + 8 * i;
d18u8 = vld1_u8(tmpp);
tmpp += 16;
d19u8 = vld1_u8(tmpp);
tmpp += 16;
d20u8 = vld1_u8(tmpp);
tmpp += 16;
d21u8 = vld1_u8(tmpp);
tmpp += 16;
d22u8 = vld1_u8(tmpp);
tmpp += 16;
for (j = 0; j < 4; j++) {
d23u8 = vld1_u8(tmpp);
tmpp += 16;
d24u8 = vld1_u8(tmpp);
tmpp += 16;
d25u8 = vld1_u8(tmpp);
tmpp += 16;
d26u8 = vld1_u8(tmpp);
tmpp += 16;
q3u16 = vmull_u8(d18u8, d0u8);
q4u16 = vmull_u8(d19u8, d0u8);
q5u16 = vmull_u8(d20u8, d0u8);
q6u16 = vmull_u8(d21u8, d0u8);
q3u16 = vmlsl_u8(q3u16, d19u8, d1u8);
q4u16 = vmlsl_u8(q4u16, d20u8, d1u8);
q5u16 = vmlsl_u8(q5u16, d21u8, d1u8);
q6u16 = vmlsl_u8(q6u16, d22u8, d1u8);
q3u16 = vmlsl_u8(q3u16, d22u8, d4u8);
q4u16 = vmlsl_u8(q4u16, d23u8, d4u8);
q5u16 = vmlsl_u8(q5u16, d24u8, d4u8);
q6u16 = vmlsl_u8(q6u16, d25u8, d4u8);
q3u16 = vmlal_u8(q3u16, d20u8, d2u8);
q4u16 = vmlal_u8(q4u16, d21u8, d2u8);
q5u16 = vmlal_u8(q5u16, d22u8, d2u8);
q6u16 = vmlal_u8(q6u16, d23u8, d2u8);
q3u16 = vmlal_u8(q3u16, d23u8, d5u8);
q4u16 = vmlal_u8(q4u16, d24u8, d5u8);
q5u16 = vmlal_u8(q5u16, d25u8, d5u8);
q6u16 = vmlal_u8(q6u16, d26u8, d5u8);
q7u16 = vmull_u8(d21u8, d3u8);
q8u16 = vmull_u8(d22u8, d3u8);
q9u16 = vmull_u8(d23u8, d3u8);
q10u16 = vmull_u8(d24u8, d3u8);
q3s16 = vreinterpretq_s16_u16(q3u16);
q4s16 = vreinterpretq_s16_u16(q4u16);
q5s16 = vreinterpretq_s16_u16(q5u16);
q6s16 = vreinterpretq_s16_u16(q6u16);
q7s16 = vreinterpretq_s16_u16(q7u16);
q8s16 = vreinterpretq_s16_u16(q8u16);
q9s16 = vreinterpretq_s16_u16(q9u16);
q10s16 = vreinterpretq_s16_u16(q10u16);
q7s16 = vqaddq_s16(q7s16, q3s16);
q8s16 = vqaddq_s16(q8s16, q4s16);
q9s16 = vqaddq_s16(q9s16, q5s16);
q10s16 = vqaddq_s16(q10s16, q6s16);
d6u8 = vqrshrun_n_s16(q7s16, 7);
d7u8 = vqrshrun_n_s16(q8s16, 7);
d8u8 = vqrshrun_n_s16(q9s16, 7);
d9u8 = vqrshrun_n_s16(q10s16, 7);
d18u8 = d22u8;
d19u8 = d23u8;
d20u8 = d24u8;
d21u8 = d25u8;
d22u8 = d26u8;
vst1_u8(dst, d6u8);
dst += dst_pitch;
vst1_u8(dst, d7u8);
dst += dst_pitch;
vst1_u8(dst, d8u8);
dst += dst_pitch;
vst1_u8(dst, d9u8);
dst += dst_pitch;
}
}
return;
}
void vp8_sixtap_predict8x8_neon(
unsigned char *src_ptr,
int src_pixels_per_line,
int xoffset,
int yoffset,
unsigned char *dst_ptr,
int dst_pitch) {
unsigned char *src, *tmpp;
unsigned char tmp[64];
int i;
uint8x8_t d0u8, d1u8, d2u8, d3u8, d4u8, d5u8, d6u8, d7u8, d8u8, d9u8;
uint8x8_t d18u8, d19u8, d20u8, d21u8, d22u8, d23u8, d24u8, d25u8;
uint8x8_t d26u8, d27u8, d28u8, d29u8, d30u8, d31u8;
int8x8_t dtmps8, d0s8, d1s8, d2s8, d3s8, d4s8, d5s8;
uint16x8_t q3u16, q4u16, q5u16, q6u16, q7u16;
uint16x8_t q8u16, q9u16, q10u16, q11u16, q12u16;
int16x8_t q3s16, q4s16, q5s16, q6s16, q7s16;
int16x8_t q8s16, q9s16, q10s16, q11s16, q12s16;
uint8x16_t q3u8, q4u8, q5u8, q6u8, q7u8, q9u8, q10u8, q11u8, q12u8;
if (xoffset == 0) { // secondpass_filter8x8_only
// load second_pass filter
dtmps8 = vld1_s8(vp8_sub_pel_filters[yoffset]);
d0s8 = vdup_lane_s8(dtmps8, 0);
d1s8 = vdup_lane_s8(dtmps8, 1);
d2s8 = vdup_lane_s8(dtmps8, 2);
d3s8 = vdup_lane_s8(dtmps8, 3);
d4s8 = vdup_lane_s8(dtmps8, 4);
d5s8 = vdup_lane_s8(dtmps8, 5);
d0u8 = vreinterpret_u8_s8(vabs_s8(d0s8));
d1u8 = vreinterpret_u8_s8(vabs_s8(d1s8));
d2u8 = vreinterpret_u8_s8(vabs_s8(d2s8));
d3u8 = vreinterpret_u8_s8(vabs_s8(d3s8));
d4u8 = vreinterpret_u8_s8(vabs_s8(d4s8));
d5u8 = vreinterpret_u8_s8(vabs_s8(d5s8));
// load src data
src = src_ptr - src_pixels_per_line * 2;
d18u8 = vld1_u8(src);
src += src_pixels_per_line;
d19u8 = vld1_u8(src);
src += src_pixels_per_line;
d20u8 = vld1_u8(src);
src += src_pixels_per_line;
d21u8 = vld1_u8(src);
src += src_pixels_per_line;
d22u8 = vld1_u8(src);
src += src_pixels_per_line;
d23u8 = vld1_u8(src);
src += src_pixels_per_line;
d24u8 = vld1_u8(src);
src += src_pixels_per_line;
d25u8 = vld1_u8(src);
src += src_pixels_per_line;
d26u8 = vld1_u8(src);
src += src_pixels_per_line;
d27u8 = vld1_u8(src);
src += src_pixels_per_line;
d28u8 = vld1_u8(src);
src += src_pixels_per_line;
d29u8 = vld1_u8(src);
src += src_pixels_per_line;
d30u8 = vld1_u8(src);
for (i = 2; i > 0; i--) {
q3u16 = vmull_u8(d18u8, d0u8);
q4u16 = vmull_u8(d19u8, d0u8);
q5u16 = vmull_u8(d20u8, d0u8);
q6u16 = vmull_u8(d21u8, d0u8);
q3u16 = vmlsl_u8(q3u16, d19u8, d1u8);
q4u16 = vmlsl_u8(q4u16, d20u8, d1u8);
q5u16 = vmlsl_u8(q5u16, d21u8, d1u8);
q6u16 = vmlsl_u8(q6u16, d22u8, d1u8);
q3u16 = vmlsl_u8(q3u16, d22u8, d4u8);
q4u16 = vmlsl_u8(q4u16, d23u8, d4u8);
q5u16 = vmlsl_u8(q5u16, d24u8, d4u8);
q6u16 = vmlsl_u8(q6u16, d25u8, d4u8);
q3u16 = vmlal_u8(q3u16, d20u8, d2u8);
q4u16 = vmlal_u8(q4u16, d21u8, d2u8);
q5u16 = vmlal_u8(q5u16, d22u8, d2u8);
q6u16 = vmlal_u8(q6u16, d23u8, d2u8);
q3u16 = vmlal_u8(q3u16, d23u8, d5u8);
q4u16 = vmlal_u8(q4u16, d24u8, d5u8);
q5u16 = vmlal_u8(q5u16, d25u8, d5u8);
q6u16 = vmlal_u8(q6u16, d26u8, d5u8);
q7u16 = vmull_u8(d21u8, d3u8);
q8u16 = vmull_u8(d22u8, d3u8);
q9u16 = vmull_u8(d23u8, d3u8);
q10u16 = vmull_u8(d24u8, d3u8);
q3s16 = vreinterpretq_s16_u16(q3u16);
q4s16 = vreinterpretq_s16_u16(q4u16);
q5s16 = vreinterpretq_s16_u16(q5u16);
q6s16 = vreinterpretq_s16_u16(q6u16);
q7s16 = vreinterpretq_s16_u16(q7u16);
q8s16 = vreinterpretq_s16_u16(q8u16);
q9s16 = vreinterpretq_s16_u16(q9u16);
q10s16 = vreinterpretq_s16_u16(q10u16);
q7s16 = vqaddq_s16(q7s16, q3s16);
q8s16 = vqaddq_s16(q8s16, q4s16);
q9s16 = vqaddq_s16(q9s16, q5s16);
q10s16 = vqaddq_s16(q10s16, q6s16);
d6u8 = vqrshrun_n_s16(q7s16, 7);
d7u8 = vqrshrun_n_s16(q8s16, 7);
d8u8 = vqrshrun_n_s16(q9s16, 7);
d9u8 = vqrshrun_n_s16(q10s16, 7);
d18u8 = d22u8;
d19u8 = d23u8;
d20u8 = d24u8;
d21u8 = d25u8;
d22u8 = d26u8;
d23u8 = d27u8;
d24u8 = d28u8;
d25u8 = d29u8;
d26u8 = d30u8;
vst1_u8(dst_ptr, d6u8);
dst_ptr += dst_pitch;
vst1_u8(dst_ptr, d7u8);
dst_ptr += dst_pitch;
vst1_u8(dst_ptr, d8u8);
dst_ptr += dst_pitch;
vst1_u8(dst_ptr, d9u8);
dst_ptr += dst_pitch;
}
return;
}
// load first_pass filter
dtmps8 = vld1_s8(vp8_sub_pel_filters[xoffset]);
d0s8 = vdup_lane_s8(dtmps8, 0);
d1s8 = vdup_lane_s8(dtmps8, 1);
d2s8 = vdup_lane_s8(dtmps8, 2);
d3s8 = vdup_lane_s8(dtmps8, 3);
d4s8 = vdup_lane_s8(dtmps8, 4);
d5s8 = vdup_lane_s8(dtmps8, 5);
d0u8 = vreinterpret_u8_s8(vabs_s8(d0s8));
d1u8 = vreinterpret_u8_s8(vabs_s8(d1s8));
d2u8 = vreinterpret_u8_s8(vabs_s8(d2s8));
d3u8 = vreinterpret_u8_s8(vabs_s8(d3s8));
d4u8 = vreinterpret_u8_s8(vabs_s8(d4s8));
d5u8 = vreinterpret_u8_s8(vabs_s8(d5s8));
// First pass: output_height lines x output_width columns (9x4)
if (yoffset == 0) // firstpass_filter4x4_only
src = src_ptr - 2;
else
src = src_ptr - 2 - (src_pixels_per_line * 2);
tmpp = tmp;
for (i = 2; i > 0; i--) {
q3u8 = vld1q_u8(src);
src += src_pixels_per_line;
q4u8 = vld1q_u8(src);
src += src_pixels_per_line;
q5u8 = vld1q_u8(src);
src += src_pixels_per_line;
q6u8 = vld1q_u8(src);
src += src_pixels_per_line;
__builtin_prefetch(src);
__builtin_prefetch(src + src_pixels_per_line);
__builtin_prefetch(src + src_pixels_per_line * 2);
q7u16 = vmull_u8(vget_low_u8(q3u8), d0u8);
q8u16 = vmull_u8(vget_low_u8(q4u8), d0u8);
q9u16 = vmull_u8(vget_low_u8(q5u8), d0u8);
q10u16 = vmull_u8(vget_low_u8(q6u8), d0u8);
d28u8 = vext_u8(vget_low_u8(q3u8), vget_high_u8(q3u8), 1);
d29u8 = vext_u8(vget_low_u8(q4u8), vget_high_u8(q4u8), 1);
d30u8 = vext_u8(vget_low_u8(q5u8), vget_high_u8(q5u8), 1);
d31u8 = vext_u8(vget_low_u8(q6u8), vget_high_u8(q6u8), 1);
q7u16 = vmlsl_u8(q7u16, d28u8, d1u8);
q8u16 = vmlsl_u8(q8u16, d29u8, d1u8);
q9u16 = vmlsl_u8(q9u16, d30u8, d1u8);
q10u16 = vmlsl_u8(q10u16, d31u8, d1u8);
d28u8 = vext_u8(vget_low_u8(q3u8), vget_high_u8(q3u8), 4);
d29u8 = vext_u8(vget_low_u8(q4u8), vget_high_u8(q4u8), 4);
d30u8 = vext_u8(vget_low_u8(q5u8), vget_high_u8(q5u8), 4);
d31u8 = vext_u8(vget_low_u8(q6u8), vget_high_u8(q6u8), 4);
q7u16 = vmlsl_u8(q7u16, d28u8, d4u8);
q8u16 = vmlsl_u8(q8u16, d29u8, d4u8);
q9u16 = vmlsl_u8(q9u16, d30u8, d4u8);
q10u16 = vmlsl_u8(q10u16, d31u8, d4u8);
d28u8 = vext_u8(vget_low_u8(q3u8), vget_high_u8(q3u8), 2);
d29u8 = vext_u8(vget_low_u8(q4u8), vget_high_u8(q4u8), 2);
d30u8 = vext_u8(vget_low_u8(q5u8), vget_high_u8(q5u8), 2);
d31u8 = vext_u8(vget_low_u8(q6u8), vget_high_u8(q6u8), 2);
q7u16 = vmlal_u8(q7u16, d28u8, d2u8);
q8u16 = vmlal_u8(q8u16, d29u8, d2u8);
q9u16 = vmlal_u8(q9u16, d30u8, d2u8);
q10u16 = vmlal_u8(q10u16, d31u8, d2u8);
d28u8 = vext_u8(vget_low_u8(q3u8), vget_high_u8(q3u8), 5);
d29u8 = vext_u8(vget_low_u8(q4u8), vget_high_u8(q4u8), 5);
d30u8 = vext_u8(vget_low_u8(q5u8), vget_high_u8(q5u8), 5);
d31u8 = vext_u8(vget_low_u8(q6u8), vget_high_u8(q6u8), 5);
q7u16 = vmlal_u8(q7u16, d28u8, d5u8);
q8u16 = vmlal_u8(q8u16, d29u8, d5u8);
q9u16 = vmlal_u8(q9u16, d30u8, d5u8);
q10u16 = vmlal_u8(q10u16, d31u8, d5u8);
d28u8 = vext_u8(vget_low_u8(q3u8), vget_high_u8(q3u8), 3);
d29u8 = vext_u8(vget_low_u8(q4u8), vget_high_u8(q4u8), 3);
d30u8 = vext_u8(vget_low_u8(q5u8), vget_high_u8(q5u8), 3);
d31u8 = vext_u8(vget_low_u8(q6u8), vget_high_u8(q6u8), 3);
q3u16 = vmull_u8(d28u8, d3u8);
q4u16 = vmull_u8(d29u8, d3u8);
q5u16 = vmull_u8(d30u8, d3u8);
q6u16 = vmull_u8(d31u8, d3u8);
q3s16 = vreinterpretq_s16_u16(q3u16);
q4s16 = vreinterpretq_s16_u16(q4u16);
q5s16 = vreinterpretq_s16_u16(q5u16);
q6s16 = vreinterpretq_s16_u16(q6u16);
q7s16 = vreinterpretq_s16_u16(q7u16);
q8s16 = vreinterpretq_s16_u16(q8u16);
q9s16 = vreinterpretq_s16_u16(q9u16);
q10s16 = vreinterpretq_s16_u16(q10u16);
q7s16 = vqaddq_s16(q7s16, q3s16);
q8s16 = vqaddq_s16(q8s16, q4s16);
q9s16 = vqaddq_s16(q9s16, q5s16);
q10s16 = vqaddq_s16(q10s16, q6s16);
d22u8 = vqrshrun_n_s16(q7s16, 7);
d23u8 = vqrshrun_n_s16(q8s16, 7);
d24u8 = vqrshrun_n_s16(q9s16, 7);
d25u8 = vqrshrun_n_s16(q10s16, 7);
if (yoffset == 0) { // firstpass_filter8x4_only
vst1_u8(dst_ptr, d22u8);
dst_ptr += dst_pitch;
vst1_u8(dst_ptr, d23u8);
dst_ptr += dst_pitch;
vst1_u8(dst_ptr, d24u8);
dst_ptr += dst_pitch;
vst1_u8(dst_ptr, d25u8);
dst_ptr += dst_pitch;
} else {
vst1_u8(tmpp, d22u8);
tmpp += 8;
vst1_u8(tmpp, d23u8);
tmpp += 8;
vst1_u8(tmpp, d24u8);
tmpp += 8;
vst1_u8(tmpp, d25u8);
tmpp += 8;
}
}
if (yoffset == 0)
return;
// First Pass on rest 5-line data
q3u8 = vld1q_u8(src);
src += src_pixels_per_line;
q4u8 = vld1q_u8(src);
src += src_pixels_per_line;
q5u8 = vld1q_u8(src);
src += src_pixels_per_line;
q6u8 = vld1q_u8(src);
src += src_pixels_per_line;
q7u8 = vld1q_u8(src);
q8u16 = vmull_u8(vget_low_u8(q3u8), d0u8);
q9u16 = vmull_u8(vget_low_u8(q4u8), d0u8);
q10u16 = vmull_u8(vget_low_u8(q5u8), d0u8);
q11u16 = vmull_u8(vget_low_u8(q6u8), d0u8);
q12u16 = vmull_u8(vget_low_u8(q7u8), d0u8);
d27u8 = vext_u8(vget_low_u8(q3u8), vget_high_u8(q3u8), 1);
d28u8 = vext_u8(vget_low_u8(q4u8), vget_high_u8(q4u8), 1);
d29u8 = vext_u8(vget_low_u8(q5u8), vget_high_u8(q5u8), 1);
d30u8 = vext_u8(vget_low_u8(q6u8), vget_high_u8(q6u8), 1);
d31u8 = vext_u8(vget_low_u8(q7u8), vget_high_u8(q7u8), 1);
q8u16 = vmlsl_u8(q8u16, d27u8, d1u8);
q9u16 = vmlsl_u8(q9u16, d28u8, d1u8);
q10u16 = vmlsl_u8(q10u16, d29u8, d1u8);
q11u16 = vmlsl_u8(q11u16, d30u8, d1u8);
q12u16 = vmlsl_u8(q12u16, d31u8, d1u8);
d27u8 = vext_u8(vget_low_u8(q3u8), vget_high_u8(q3u8), 4);
d28u8 = vext_u8(vget_low_u8(q4u8), vget_high_u8(q4u8), 4);
d29u8 = vext_u8(vget_low_u8(q5u8), vget_high_u8(q5u8), 4);
d30u8 = vext_u8(vget_low_u8(q6u8), vget_high_u8(q6u8), 4);
d31u8 = vext_u8(vget_low_u8(q7u8), vget_high_u8(q7u8), 4);
q8u16 = vmlsl_u8(q8u16, d27u8, d4u8);
q9u16 = vmlsl_u8(q9u16, d28u8, d4u8);
q10u16 = vmlsl_u8(q10u16, d29u8, d4u8);
q11u16 = vmlsl_u8(q11u16, d30u8, d4u8);
q12u16 = vmlsl_u8(q12u16, d31u8, d4u8);
d27u8 = vext_u8(vget_low_u8(q3u8), vget_high_u8(q3u8), 2);
d28u8 = vext_u8(vget_low_u8(q4u8), vget_high_u8(q4u8), 2);
d29u8 = vext_u8(vget_low_u8(q5u8), vget_high_u8(q5u8), 2);
d30u8 = vext_u8(vget_low_u8(q6u8), vget_high_u8(q6u8), 2);
d31u8 = vext_u8(vget_low_u8(q7u8), vget_high_u8(q7u8), 2);
q8u16 = vmlal_u8(q8u16, d27u8, d2u8);
q9u16 = vmlal_u8(q9u16, d28u8, d2u8);
q10u16 = vmlal_u8(q10u16, d29u8, d2u8);
q11u16 = vmlal_u8(q11u16, d30u8, d2u8);
q12u16 = vmlal_u8(q12u16, d31u8, d2u8);
d27u8 = vext_u8(vget_low_u8(q3u8), vget_high_u8(q3u8), 5);
d28u8 = vext_u8(vget_low_u8(q4u8), vget_high_u8(q4u8), 5);
d29u8 = vext_u8(vget_low_u8(q5u8), vget_high_u8(q5u8), 5);
d30u8 = vext_u8(vget_low_u8(q6u8), vget_high_u8(q6u8), 5);
d31u8 = vext_u8(vget_low_u8(q7u8), vget_high_u8(q7u8), 5);
q8u16 = vmlal_u8(q8u16, d27u8, d5u8);
q9u16 = vmlal_u8(q9u16, d28u8, d5u8);
q10u16 = vmlal_u8(q10u16, d29u8, d5u8);
q11u16 = vmlal_u8(q11u16, d30u8, d5u8);
q12u16 = vmlal_u8(q12u16, d31u8, d5u8);
d27u8 = vext_u8(vget_low_u8(q3u8), vget_high_u8(q3u8), 3);
d28u8 = vext_u8(vget_low_u8(q4u8), vget_high_u8(q4u8), 3);
d29u8 = vext_u8(vget_low_u8(q5u8), vget_high_u8(q5u8), 3);
d30u8 = vext_u8(vget_low_u8(q6u8), vget_high_u8(q6u8), 3);
d31u8 = vext_u8(vget_low_u8(q7u8), vget_high_u8(q7u8), 3);
q3u16 = vmull_u8(d27u8, d3u8);
q4u16 = vmull_u8(d28u8, d3u8);
q5u16 = vmull_u8(d29u8, d3u8);
q6u16 = vmull_u8(d30u8, d3u8);
q7u16 = vmull_u8(d31u8, d3u8);
q3s16 = vreinterpretq_s16_u16(q3u16);
q4s16 = vreinterpretq_s16_u16(q4u16);
q5s16 = vreinterpretq_s16_u16(q5u16);
q6s16 = vreinterpretq_s16_u16(q6u16);
q7s16 = vreinterpretq_s16_u16(q7u16);
q8s16 = vreinterpretq_s16_u16(q8u16);
q9s16 = vreinterpretq_s16_u16(q9u16);
q10s16 = vreinterpretq_s16_u16(q10u16);
q11s16 = vreinterpretq_s16_u16(q11u16);
q12s16 = vreinterpretq_s16_u16(q12u16);
q8s16 = vqaddq_s16(q8s16, q3s16);
q9s16 = vqaddq_s16(q9s16, q4s16);
q10s16 = vqaddq_s16(q10s16, q5s16);
q11s16 = vqaddq_s16(q11s16, q6s16);
q12s16 = vqaddq_s16(q12s16, q7s16);
d26u8 = vqrshrun_n_s16(q8s16, 7);
d27u8 = vqrshrun_n_s16(q9s16, 7);
d28u8 = vqrshrun_n_s16(q10s16, 7);
d29u8 = vqrshrun_n_s16(q11s16, 7);
d30u8 = vqrshrun_n_s16(q12s16, 7);
// Second pass: 8x8
dtmps8 = vld1_s8(vp8_sub_pel_filters[yoffset]);
d0s8 = vdup_lane_s8(dtmps8, 0);
d1s8 = vdup_lane_s8(dtmps8, 1);
d2s8 = vdup_lane_s8(dtmps8, 2);
d3s8 = vdup_lane_s8(dtmps8, 3);
d4s8 = vdup_lane_s8(dtmps8, 4);
d5s8 = vdup_lane_s8(dtmps8, 5);
d0u8 = vreinterpret_u8_s8(vabs_s8(d0s8));
d1u8 = vreinterpret_u8_s8(vabs_s8(d1s8));
d2u8 = vreinterpret_u8_s8(vabs_s8(d2s8));
d3u8 = vreinterpret_u8_s8(vabs_s8(d3s8));
d4u8 = vreinterpret_u8_s8(vabs_s8(d4s8));
d5u8 = vreinterpret_u8_s8(vabs_s8(d5s8));
tmpp = tmp;
q9u8 = vld1q_u8(tmpp);
tmpp += 16;
q10u8 = vld1q_u8(tmpp);
tmpp += 16;
q11u8 = vld1q_u8(tmpp);
tmpp += 16;
q12u8 = vld1q_u8(tmpp);
d18u8 = vget_low_u8(q9u8);
d19u8 = vget_high_u8(q9u8);
d20u8 = vget_low_u8(q10u8);
d21u8 = vget_high_u8(q10u8);
d22u8 = vget_low_u8(q11u8);
d23u8 = vget_high_u8(q11u8);
d24u8 = vget_low_u8(q12u8);
d25u8 = vget_high_u8(q12u8);
for (i = 2; i > 0; i--) {
q3u16 = vmull_u8(d18u8, d0u8);
q4u16 = vmull_u8(d19u8, d0u8);
q5u16 = vmull_u8(d20u8, d0u8);
q6u16 = vmull_u8(d21u8, d0u8);
q3u16 = vmlsl_u8(q3u16, d19u8, d1u8);
q4u16 = vmlsl_u8(q4u16, d20u8, d1u8);
q5u16 = vmlsl_u8(q5u16, d21u8, d1u8);
q6u16 = vmlsl_u8(q6u16, d22u8, d1u8);
q3u16 = vmlsl_u8(q3u16, d22u8, d4u8);
q4u16 = vmlsl_u8(q4u16, d23u8, d4u8);
q5u16 = vmlsl_u8(q5u16, d24u8, d4u8);
q6u16 = vmlsl_u8(q6u16, d25u8, d4u8);
q3u16 = vmlal_u8(q3u16, d20u8, d2u8);
q4u16 = vmlal_u8(q4u16, d21u8, d2u8);
q5u16 = vmlal_u8(q5u16, d22u8, d2u8);
q6u16 = vmlal_u8(q6u16, d23u8, d2u8);
q3u16 = vmlal_u8(q3u16, d23u8, d5u8);
q4u16 = vmlal_u8(q4u16, d24u8, d5u8);
q5u16 = vmlal_u8(q5u16, d25u8, d5u8);
q6u16 = vmlal_u8(q6u16, d26u8, d5u8);
q7u16 = vmull_u8(d21u8, d3u8);
q8u16 = vmull_u8(d22u8, d3u8);
q9u16 = vmull_u8(d23u8, d3u8);
q10u16 = vmull_u8(d24u8, d3u8);
q3s16 = vreinterpretq_s16_u16(q3u16);
q4s16 = vreinterpretq_s16_u16(q4u16);
q5s16 = vreinterpretq_s16_u16(q5u16);
q6s16 = vreinterpretq_s16_u16(q6u16);
q7s16 = vreinterpretq_s16_u16(q7u16);
q8s16 = vreinterpretq_s16_u16(q8u16);
q9s16 = vreinterpretq_s16_u16(q9u16);
q10s16 = vreinterpretq_s16_u16(q10u16);
q7s16 = vqaddq_s16(q7s16, q3s16);
q8s16 = vqaddq_s16(q8s16, q4s16);
q9s16 = vqaddq_s16(q9s16, q5s16);
q10s16 = vqaddq_s16(q10s16, q6s16);
d6u8 = vqrshrun_n_s16(q7s16, 7);
d7u8 = vqrshrun_n_s16(q8s16, 7);
d8u8 = vqrshrun_n_s16(q9s16, 7);
d9u8 = vqrshrun_n_s16(q10s16, 7);
d18u8 = d22u8;
d19u8 = d23u8;
d20u8 = d24u8;
d21u8 = d25u8;
d22u8 = d26u8;
d23u8 = d27u8;
d24u8 = d28u8;
d25u8 = d29u8;
d26u8 = d30u8;
vst1_u8(dst_ptr, d6u8);
dst_ptr += dst_pitch;
vst1_u8(dst_ptr, d7u8);
dst_ptr += dst_pitch;
vst1_u8(dst_ptr, d8u8);
dst_ptr += dst_pitch;
vst1_u8(dst_ptr, d9u8);
dst_ptr += dst_pitch;
}
return;
}
f64 dotProduct(const Size2D &_size,
const s8 * src0Base, ptrdiff_t src0Stride,
const s8 * src1Base, ptrdiff_t src1Stride)
{
internal::assertSupportedConfiguration();
#ifdef CAROTENE_NEON
Size2D size(_size);
if (src0Stride == src1Stride &&
src0Stride == (ptrdiff_t)(size.width))
{
size.width *= size.height;
size.height = 1;
}
// It is possible to accumulate up to 131071 schar multiplication results in sint32 without overflow
// We process 16 elements and accumulate two new elements per step. So we could handle 131071/2*16 elements
#define DOT_INT_BLOCKSIZE 131070*8
f64 result = 0.0;
for (size_t row = 0; row < size.height; ++row)
{
const s8 * src0 = internal::getRowPtr(src0Base, src0Stride, row);
const s8 * src1 = internal::getRowPtr(src1Base, src1Stride, row);
size_t i = 0;
int64x2_t ws = vmovq_n_s64(0);
while(i + 16 <= size.width)
{
size_t lim = std::min(i + DOT_UINT_BLOCKSIZE, size.width) - 16;
int32x4_t s1 = vmovq_n_s32(0);
int32x4_t s2 = vmovq_n_s32(0);
for (; i <= lim; i += 16)
{
internal::prefetch(src0 + i);
internal::prefetch(src1 + i);
int8x16_t vs1 = vld1q_s8(src0 + i);
int8x16_t vs2 = vld1q_s8(src1 + i);
int16x8_t vdot1 = vmull_s8(vget_low_s8(vs1), vget_low_s8(vs2));
int16x8_t vdot2 = vmull_s8(vget_high_s8(vs1), vget_high_s8(vs2));
s1 = vpadalq_s16(s1, vdot1);
s2 = vpadalq_s16(s2, vdot2);
}
ws = vpadalq_s32(ws, s1);
ws = vpadalq_s32(ws, s2);
}
if(i + 8 <= size.width)
{
int8x8_t vs1 = vld1_s8(src0 + i);
int8x8_t vs2 = vld1_s8(src1 + i);
ws = vpadalq_s32(ws, vpaddlq_s16(vmull_s8(vs1, vs2)));
i += 8;
}
result += (double)vget_lane_s64(vadd_s64(vget_low_s64(ws), vget_high_s64(ws)), 0);
for (; i < size.width; ++i)
result += s32(src0[i]) * s32(src1[i]);
}
return result;
#else
(void)_size;
(void)src0Base;
(void)src0Stride;
(void)src1Base;
(void)src1Stride;
return 0;
#endif
}
