1 // Auto-generated file. Do not edit!
2 // Template: src/s16-window/neon.c.in
3 // Generator: tools/xngen
4 //
5 // Copyright 2022 Google LLC
6 //
7 // This source code is licensed under the BSD-style license found in the
8 // LICENSE file in the root directory of this source tree.
9
10 #include <assert.h>
11 #include <stddef.h>
12 #include <stdint.h>
13
14 #include <arm_neon.h>
15
16 #include <xnnpack/math.h>
17 #include <xnnpack/window.h>
18
19
xnn_s16_window_ukernel__neon_x8(size_t rows,size_t batch_size,const int16_t * input,const int16_t * weights,int16_t * output,uint32_t shift)20 void xnn_s16_window_ukernel__neon_x8(
21 size_t rows,
22 size_t batch_size,
23 const int16_t* input,
24 const int16_t* weights,
25 int16_t* output,
26 uint32_t shift)
27 {
28 assert(rows != 0);
29 assert(batch_size != 0);
30 assert(input != NULL);
31 assert(weights != NULL);
32 assert(output != NULL);
33 assert(shift < 32);
34
35 const int32x4_t vshift = vdupq_n_s32(-(int32_t)shift); // negative to shift right.
36
37 do {
38 const int16_t* w = weights;
39 size_t n = batch_size * sizeof(int16_t);
40
41 // Remainder of full vectors
42 for (; n >= 8 * sizeof(int16_t); n -= 8 * sizeof(int16_t)) {
43 const int16x8_t vi = vld1q_s16(input); input += 8;
44 const int16x8_t vw = vld1q_s16(w); w += 8;
45 int32x4_t vacc_lo = vmull_s16(vget_low_s16(vi), vget_low_s16(vw));
46 int32x4_t vacc_hi = vmull_s16(vget_high_s16(vi), vget_high_s16(vw));
47 vacc_lo = vshlq_s32(vacc_lo, vshift);
48 vacc_hi = vshlq_s32(vacc_hi, vshift);
49 const int16x8_t vout = vcombine_s16(vqmovn_s32(vacc_lo), vqmovn_s32(vacc_hi));
50 vst1q_s16(output, vout); output += 8;
51 }
52
53 assert(n % 2 == 0);
54 // Remainder of 1 to 7 batch_size
55 if XNN_UNLIKELY(n != 0) {
56 const int16x8_t vi = vld1q_s16(input); input = (const int16_t*) ((uintptr_t) input + n);
57 const int16x8_t vw = vld1q_s16(w);
58 int32x4_t vacc = vmull_s16(vget_low_s16(vi), vget_low_s16(vw));
59 vacc = vshlq_s32(vacc, vshift);
60 int16x4_t vout = vqmovn_s32(vacc);
61 if (n & (4 * sizeof(int16_t))) {
62 vst1_s16(output, vout); output += 4;
63 vacc = vmull_s16(vget_high_s16(vi), vget_high_s16(vw));
64 vacc = vshlq_s32(vacc, vshift);
65 vout = vqmovn_s32(vacc);
66 }
67 if (n & (2 * sizeof(int16_t))) {
68 vst1_lane_u32((void*) output, vreinterpret_u32_s16(vout), 0); output += 2;
69 vout = vext_s16(vout, vout, 2);
70 }
71 if (n & (1 * sizeof(int16_t))) {
72 vst1_lane_s16(output, vout, 0); output += 1;
73 }
74 }
75
76 } while (--rows != 0);
77 }
78