-
Notifications
You must be signed in to change notification settings - Fork 326
/
avx.c.in
143 lines (129 loc) · 4.58 KB
/
avx.c.in
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
// Copyright 2024 Google LLC
//
// This source code is licensed under the BSD-style license found in the
// LICENSE file in the root directory of this source tree.
$ABC = "0123456789ABCDEFGHIJKLMNOPQRSTUVWXYZ"
#include <assert.h>
#include <immintrin.h>
#include <xnnpack/common.h>
#include <xnnpack/reduce.h>
#include <xnnpack/math.h>
$UNROLL = CHANNELS >> 3
void xnn_f32_rdsum_ukernel_${ACCUMULATORS}p${ACCUMULATORS}x__avx_c${CHANNELS}(
size_t rows,
size_t channels,
const float* input,
size_t input_stride,
const float* zero,
float* output,
const union xnn_f32_scale_params params[restrict XNN_MIN_ELEMENTS(1)])
{
assert(rows != 0);
assert(channels != 0);
assert(input != NULL);
assert(output != NULL);
const __m256 vscale = _mm256_set1_ps(params->avx.scale);
size_t input_increment = ${ACCUMULATORS} * input_stride;
for (; channels >= ${CHANNELS}; channels -= ${CHANNELS}) {
const float* i0 = input;
$for i in range(1, ACCUMULATORS):
const float* i${i} = (const float*) ((uintptr_t) input + ${i} * input_stride);
$for i in range(UNROLL):
__m256 vacc${i} = _mm256_setzero_ps();
for (int r = rows; r > 0; r -= ${ACCUMULATORS}) {
$for N in range(1, ACCUMULATORS, 2):
if XNN_UNPREDICTABLE(r < ${N+1}) {
i${N} = zero;
}
if XNN_UNPREDICTABLE(r <= ${N+1}) {
i${N+1} = zero;
}
$for c in range(UNROLL):
__m256 vin${c};
$for j in range(ACCUMULATORS):
$for c in range(UNROLL):
vin${c} = _mm256_loadu_ps(&i${j}[${c*8}]);
$for c in range(UNROLL):
vacc${c} = _mm256_add_ps(vin${c}, vacc${c});
$for N in range(0, ACCUMULATORS):
i${N} = (const float*) ((uintptr_t) i${N} + input_increment);
}
$for i in range(UNROLL):
vacc${i} = _mm256_mul_ps(vacc${i}, vscale);
const float* o = output;
$for i in range(0, UNROLL):
__m256 vo${i} = _mm256_loadu_ps(o); o += 8;
$for i in range(0, UNROLL):
vacc${i} = _mm256_add_ps(vo${i}, vacc${i});
$for i in range(0, UNROLL):
_mm256_storeu_ps(output, vacc${i}); output += 8;
input = (const float*) ((uintptr_t) input + ${CHANNELS} * sizeof(float));
}
__m256i vmask;
if (channels != 0) {
input_increment = ${ACCUMULATORS} * input_stride;
const float* i0 = input;
$for i in range(1, ACCUMULATORS):
const float* i${i} = (const float*) ((uintptr_t) input + ${i} * input_stride);
__m256 vacc[${UNROLL}];
$for i in range(UNROLL):
vacc[${i}] = _mm256_setzero_ps();
const size_t num_full_chunks = channels >> 3;
const size_t num_chunks = round_up_po2(channels, 8) >> 3;
const size_t remainder = channels & 0x7;
for (int r = rows; r > 0; r -= ${ACCUMULATORS}) {
$for N in range(1, ACCUMULATORS, 2):
if XNN_UNPREDICTABLE(r < ${N+1}) {
i${N} = zero;
}
if XNN_UNPREDICTABLE(r <= ${N+1}) {
i${N+1} = zero;
}
for (int i = 0; i < num_full_chunks; ++i) {
$for c in range(ACCUMULATORS):
vacc[i] = _mm256_add_ps(_mm256_loadu_ps(&i${c}[i*8]), vacc[i]);
}
if (remainder) {
vmask = _mm256_loadu_si256((const __m256i*) ((uintptr_t) ¶ms->avx.mask_table[7] - (channels & 0x7) * sizeof(float)));
$for c in range(ACCUMULATORS):
vacc[num_full_chunks] = _mm256_add_ps(_mm256_maskload_ps(&i${c}[num_full_chunks*8], vmask), vacc[num_full_chunks]);
}
$for N in range(ACCUMULATORS):
i${N} = (const float*) ((uintptr_t) i${N} + input_increment);
}
for (size_t i = 0; i < num_chunks; ++i) {
vacc[i] = _mm256_mul_ps(vacc[i], vscale);
}
__m256 vo[${UNROLL}];
const float* o = output;
for (int i = 0; i < channels >> 3; ++i) {
vo[i] = _mm256_loadu_ps(o); o += 8;
}
for (int i = 0; i < channels >> 3; ++i) {
vacc[i] = _mm256_add_ps(vo[i], vacc[i]);
}
for (int i = 0; i < channels >> 3; ++i) {
_mm256_storeu_ps(output, vacc[i]); output += 8;
}
if (remainder) {
const size_t pos = num_full_chunks;
__m256 vout = vacc[pos];
const __m256 vdata = _mm256_maskload_ps(output, vmask);
vout = _mm256_add_ps(vout, vdata);
__m128 vout_lo = _mm256_castps256_ps128(vout);
if (channels & 4) {
_mm_storeu_ps(output, vout_lo);
vout_lo = _mm256_extractf128_ps(vout, 1);
output += 4;
}
if (channels & 2) {
_mm_storel_pi((__m64*) output, vout_lo);
vout_lo = _mm_movehl_ps(vout_lo, vout_lo);
output += 2;
}
if (channels & 1) {
_mm_store_ss(output, vout_lo);
}
}
}
}