blob: f505a669bd73259bc69f115c7c79f882fd5285bc [file] [log] [blame]
Marat Dukhan98042f22021-06-15 00:43:13 -07001// Auto-generated file. Do not edit!
2// Template: src/qs8-dwconv/unipass-sse-mul32.c.in
3// Generator: tools/xngen
4//
5// Copyright 2020 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
12#include <immintrin.h>
13
14#include <xnnpack/dwconv.h>
15#include <xnnpack/intrinsics-polyfill.h>
16
17
18void xnn_qc8_dwconv_minmax_fp32_ukernel_up24x25__sse41_mul32(
19 size_t channels,
20 size_t output_width,
21 const int8_t** input,
22 const void* weights,
23 int8_t* output,
24 size_t input_stride,
25 size_t output_increment,
26 size_t input_offset,
27 const int8_t* zero,
28 const union xnn_qs8_minmax_params params[restrict XNN_MIN_ELEMENTS(1)]) XNN_DISABLE_TSAN
29{
30 assert(channels != 0);
31 assert(output_width != 0);
32
33 do {
34 const int8_t* i0 = input[0];
35 assert(i0 != NULL);
36 if XNN_UNPREDICTABLE(i0 != zero) {
37 i0 = (const int8_t*) ((uintptr_t) i0 + input_offset);
38 }
39 const int8_t* i1 = input[1];
40 assert(i1 != NULL);
41 if XNN_UNPREDICTABLE(i1 != zero) {
42 i1 = (const int8_t*) ((uintptr_t) i1 + input_offset);
43 }
44 const int8_t* i2 = input[2];
45 assert(i2 != NULL);
46 if XNN_UNPREDICTABLE(i2 != zero) {
47 i2 = (const int8_t*) ((uintptr_t) i2 + input_offset);
48 }
49 const int8_t* i3 = input[3];
50 assert(i3 != NULL);
51 if XNN_UNPREDICTABLE(i3 != zero) {
52 i3 = (const int8_t*) ((uintptr_t) i3 + input_offset);
53 }
54 const int8_t* i4 = input[4];
55 assert(i4 != NULL);
56 if XNN_UNPREDICTABLE(i4 != zero) {
57 i4 = (const int8_t*) ((uintptr_t) i4 + input_offset);
58 }
59 const int8_t* i5 = input[5];
60 assert(i5 != NULL);
61 if XNN_UNPREDICTABLE(i5 != zero) {
62 i5 = (const int8_t*) ((uintptr_t) i5 + input_offset);
63 }
64 const int8_t* i6 = input[6];
65 assert(i6 != NULL);
66 if XNN_UNPREDICTABLE(i6 != zero) {
67 i6 = (const int8_t*) ((uintptr_t) i6 + input_offset);
68 }
69 const int8_t* i7 = input[7];
70 assert(i7 != NULL);
71 if XNN_UNPREDICTABLE(i7 != zero) {
72 i7 = (const int8_t*) ((uintptr_t) i7 + input_offset);
73 }
74 const int8_t* i8 = input[8];
75 assert(i8 != NULL);
76 if XNN_UNPREDICTABLE(i8 != zero) {
77 i8 = (const int8_t*) ((uintptr_t) i8 + input_offset);
78 }
79 const int8_t* i9 = input[9];
80 assert(i9 != NULL);
81 if XNN_UNPREDICTABLE(i9 != zero) {
82 i9 = (const int8_t*) ((uintptr_t) i9 + input_offset);
83 }
84 const int8_t* i10 = input[10];
85 assert(i10 != NULL);
86 if XNN_UNPREDICTABLE(i10 != zero) {
87 i10 = (const int8_t*) ((uintptr_t) i10 + input_offset);
88 }
89 const int8_t* i11 = input[11];
90 assert(i11 != NULL);
91 if XNN_UNPREDICTABLE(i11 != zero) {
92 i11 = (const int8_t*) ((uintptr_t) i11 + input_offset);
93 }
94 const int8_t* i12 = input[12];
95 assert(i12 != NULL);
96 if XNN_UNPREDICTABLE(i12 != zero) {
97 i12 = (const int8_t*) ((uintptr_t) i12 + input_offset);
98 }
99 const int8_t* i13 = input[13];
100 assert(i13 != NULL);
101 if XNN_UNPREDICTABLE(i13 != zero) {
102 i13 = (const int8_t*) ((uintptr_t) i13 + input_offset);
103 }
104 const int8_t* i14 = input[14];
105 assert(i14 != NULL);
106 if XNN_UNPREDICTABLE(i14 != zero) {
107 i14 = (const int8_t*) ((uintptr_t) i14 + input_offset);
108 }
109 const int8_t* i15 = input[15];
110 assert(i15 != NULL);
111 if XNN_UNPREDICTABLE(i15 != zero) {
112 i15 = (const int8_t*) ((uintptr_t) i15 + input_offset);
113 }
114 const int8_t* i16 = input[16];
115 assert(i16 != NULL);
116 if XNN_UNPREDICTABLE(i16 != zero) {
117 i16 = (const int8_t*) ((uintptr_t) i16 + input_offset);
118 }
119 const int8_t* i17 = input[17];
120 assert(i17 != NULL);
121 if XNN_UNPREDICTABLE(i17 != zero) {
122 i17 = (const int8_t*) ((uintptr_t) i17 + input_offset);
123 }
124 const int8_t* i18 = input[18];
125 assert(i18 != NULL);
126 if XNN_UNPREDICTABLE(i18 != zero) {
127 i18 = (const int8_t*) ((uintptr_t) i18 + input_offset);
128 }
129 const int8_t* i19 = input[19];
130 assert(i19 != NULL);
131 if XNN_UNPREDICTABLE(i19 != zero) {
132 i19 = (const int8_t*) ((uintptr_t) i19 + input_offset);
133 }
134 const int8_t* i20 = input[20];
135 assert(i20 != NULL);
136 if XNN_UNPREDICTABLE(i20 != zero) {
137 i20 = (const int8_t*) ((uintptr_t) i20 + input_offset);
138 }
139 const int8_t* i21 = input[21];
140 assert(i21 != NULL);
141 if XNN_UNPREDICTABLE(i21 != zero) {
142 i21 = (const int8_t*) ((uintptr_t) i21 + input_offset);
143 }
144 const int8_t* i22 = input[22];
145 assert(i22 != NULL);
146 if XNN_UNPREDICTABLE(i22 != zero) {
147 i22 = (const int8_t*) ((uintptr_t) i22 + input_offset);
148 }
149 const int8_t* i23 = input[23];
150 assert(i23 != NULL);
151 if XNN_UNPREDICTABLE(i23 != zero) {
152 i23 = (const int8_t*) ((uintptr_t) i23 + input_offset);
153 }
154 const int8_t* i24 = input[24];
155 assert(i24 != NULL);
156 if XNN_UNPREDICTABLE(i24 != zero) {
157 i24 = (const int8_t*) ((uintptr_t) i24 + input_offset);
158 }
159 input = (const int8_t**) ((uintptr_t) input + input_stride);
160
161 size_t c = channels;
162 const void* w = weights;
163 for (; c >= 24; c -= 24) {
164 __m128i vacc0123 = _mm_loadu_si128((const __m128i*) w);
Marat Dukhane5eee462021-07-01 19:34:39 -0700165 __m128i vacc4567 = _mm_loadu_si128((const __m128i*) ((const int32_t*) w + 4));
166 __m128i vacc89AB = _mm_loadu_si128((const __m128i*) ((const int32_t*) w + 8));
167 __m128i vaccCDEF = _mm_loadu_si128((const __m128i*) ((const int32_t*) w + 12));
168 __m128i vaccGHIJ = _mm_loadu_si128((const __m128i*) ((const int32_t*) w + 16));
169 __m128i vaccKLMN = _mm_loadu_si128((const __m128i*) ((const int32_t*) w + 20));
Marat Dukhan98042f22021-06-15 00:43:13 -0700170
171
172 const __m128i vi0x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32(i0));
173 const __m128i vk0x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 0 * sizeof(int8_t))));
174 const __m128i vi0x4567 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i0 + 4)));
175 const __m128i vk0x4567 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 4 * sizeof(int8_t))));
176 const __m128i vi0x89AB = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i0 + 8)));
177 const __m128i vk0x89AB = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 8 * sizeof(int8_t))));
178 const __m128i vi0xCDEF = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i0 + 12)));
179 const __m128i vk0xCDEF = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 12 * sizeof(int8_t))));
180 const __m128i vi0xGHIJ = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i0 + 16)));
181 const __m128i vk0xGHIJ = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 16 * sizeof(int8_t))));
182 const __m128i vi0xKLMN = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i0 + 20)));
183 const __m128i vk0xKLMN = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 20 * sizeof(int8_t))));
184 i0 += 24;
185
186 vacc0123 = _mm_add_epi32(vacc0123, _mm_mullo_epi32(vi0x0123, vk0x0123));
187 vacc4567 = _mm_add_epi32(vacc4567, _mm_mullo_epi32(vi0x4567, vk0x4567));
188 vacc89AB = _mm_add_epi32(vacc89AB, _mm_mullo_epi32(vi0x89AB, vk0x89AB));
189 vaccCDEF = _mm_add_epi32(vaccCDEF, _mm_mullo_epi32(vi0xCDEF, vk0xCDEF));
190 vaccGHIJ = _mm_add_epi32(vaccGHIJ, _mm_mullo_epi32(vi0xGHIJ, vk0xGHIJ));
191 vaccKLMN = _mm_add_epi32(vaccKLMN, _mm_mullo_epi32(vi0xKLMN, vk0xKLMN));
192
193 const __m128i vi1x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32(i1));
194 const __m128i vk1x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 24 * sizeof(int8_t))));
195 const __m128i vi1x4567 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i1 + 4)));
196 const __m128i vk1x4567 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 28 * sizeof(int8_t))));
197 const __m128i vi1x89AB = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i1 + 8)));
198 const __m128i vk1x89AB = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 32 * sizeof(int8_t))));
199 const __m128i vi1xCDEF = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i1 + 12)));
200 const __m128i vk1xCDEF = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 36 * sizeof(int8_t))));
201 const __m128i vi1xGHIJ = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i1 + 16)));
202 const __m128i vk1xGHIJ = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 40 * sizeof(int8_t))));
203 const __m128i vi1xKLMN = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i1 + 20)));
204 const __m128i vk1xKLMN = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 44 * sizeof(int8_t))));
205 i1 += 24;
206
207 vacc0123 = _mm_add_epi32(vacc0123, _mm_mullo_epi32(vi1x0123, vk1x0123));
208 vacc4567 = _mm_add_epi32(vacc4567, _mm_mullo_epi32(vi1x4567, vk1x4567));
209 vacc89AB = _mm_add_epi32(vacc89AB, _mm_mullo_epi32(vi1x89AB, vk1x89AB));
210 vaccCDEF = _mm_add_epi32(vaccCDEF, _mm_mullo_epi32(vi1xCDEF, vk1xCDEF));
211 vaccGHIJ = _mm_add_epi32(vaccGHIJ, _mm_mullo_epi32(vi1xGHIJ, vk1xGHIJ));
212 vaccKLMN = _mm_add_epi32(vaccKLMN, _mm_mullo_epi32(vi1xKLMN, vk1xKLMN));
213
214 const __m128i vi2x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32(i2));
215 const __m128i vk2x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 48 * sizeof(int8_t))));
216 const __m128i vi2x4567 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i2 + 4)));
217 const __m128i vk2x4567 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 52 * sizeof(int8_t))));
218 const __m128i vi2x89AB = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i2 + 8)));
219 const __m128i vk2x89AB = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 56 * sizeof(int8_t))));
220 const __m128i vi2xCDEF = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i2 + 12)));
221 const __m128i vk2xCDEF = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 60 * sizeof(int8_t))));
222 const __m128i vi2xGHIJ = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i2 + 16)));
223 const __m128i vk2xGHIJ = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 64 * sizeof(int8_t))));
224 const __m128i vi2xKLMN = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i2 + 20)));
225 const __m128i vk2xKLMN = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 68 * sizeof(int8_t))));
226 i2 += 24;
227
228 vacc0123 = _mm_add_epi32(vacc0123, _mm_mullo_epi32(vi2x0123, vk2x0123));
229 vacc4567 = _mm_add_epi32(vacc4567, _mm_mullo_epi32(vi2x4567, vk2x4567));
230 vacc89AB = _mm_add_epi32(vacc89AB, _mm_mullo_epi32(vi2x89AB, vk2x89AB));
231 vaccCDEF = _mm_add_epi32(vaccCDEF, _mm_mullo_epi32(vi2xCDEF, vk2xCDEF));
232 vaccGHIJ = _mm_add_epi32(vaccGHIJ, _mm_mullo_epi32(vi2xGHIJ, vk2xGHIJ));
233 vaccKLMN = _mm_add_epi32(vaccKLMN, _mm_mullo_epi32(vi2xKLMN, vk2xKLMN));
234
235 const __m128i vi3x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32(i3));
236 const __m128i vk3x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 72 * sizeof(int8_t))));
237 const __m128i vi3x4567 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i3 + 4)));
238 const __m128i vk3x4567 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 76 * sizeof(int8_t))));
239 const __m128i vi3x89AB = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i3 + 8)));
240 const __m128i vk3x89AB = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 80 * sizeof(int8_t))));
241 const __m128i vi3xCDEF = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i3 + 12)));
242 const __m128i vk3xCDEF = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 84 * sizeof(int8_t))));
243 const __m128i vi3xGHIJ = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i3 + 16)));
244 const __m128i vk3xGHIJ = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 88 * sizeof(int8_t))));
245 const __m128i vi3xKLMN = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i3 + 20)));
246 const __m128i vk3xKLMN = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 92 * sizeof(int8_t))));
247 i3 += 24;
248
249 vacc0123 = _mm_add_epi32(vacc0123, _mm_mullo_epi32(vi3x0123, vk3x0123));
250 vacc4567 = _mm_add_epi32(vacc4567, _mm_mullo_epi32(vi3x4567, vk3x4567));
251 vacc89AB = _mm_add_epi32(vacc89AB, _mm_mullo_epi32(vi3x89AB, vk3x89AB));
252 vaccCDEF = _mm_add_epi32(vaccCDEF, _mm_mullo_epi32(vi3xCDEF, vk3xCDEF));
253 vaccGHIJ = _mm_add_epi32(vaccGHIJ, _mm_mullo_epi32(vi3xGHIJ, vk3xGHIJ));
254 vaccKLMN = _mm_add_epi32(vaccKLMN, _mm_mullo_epi32(vi3xKLMN, vk3xKLMN));
255
256 const __m128i vi4x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32(i4));
257 const __m128i vk4x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 96 * sizeof(int8_t))));
258 const __m128i vi4x4567 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i4 + 4)));
259 const __m128i vk4x4567 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 100 * sizeof(int8_t))));
260 const __m128i vi4x89AB = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i4 + 8)));
261 const __m128i vk4x89AB = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 104 * sizeof(int8_t))));
262 const __m128i vi4xCDEF = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i4 + 12)));
263 const __m128i vk4xCDEF = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 108 * sizeof(int8_t))));
264 const __m128i vi4xGHIJ = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i4 + 16)));
265 const __m128i vk4xGHIJ = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 112 * sizeof(int8_t))));
266 const __m128i vi4xKLMN = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i4 + 20)));
267 const __m128i vk4xKLMN = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 116 * sizeof(int8_t))));
268 i4 += 24;
269
270 vacc0123 = _mm_add_epi32(vacc0123, _mm_mullo_epi32(vi4x0123, vk4x0123));
271 vacc4567 = _mm_add_epi32(vacc4567, _mm_mullo_epi32(vi4x4567, vk4x4567));
272 vacc89AB = _mm_add_epi32(vacc89AB, _mm_mullo_epi32(vi4x89AB, vk4x89AB));
273 vaccCDEF = _mm_add_epi32(vaccCDEF, _mm_mullo_epi32(vi4xCDEF, vk4xCDEF));
274 vaccGHIJ = _mm_add_epi32(vaccGHIJ, _mm_mullo_epi32(vi4xGHIJ, vk4xGHIJ));
275 vaccKLMN = _mm_add_epi32(vaccKLMN, _mm_mullo_epi32(vi4xKLMN, vk4xKLMN));
276
277 const __m128i vi5x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32(i5));
278 const __m128i vk5x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 120 * sizeof(int8_t))));
279 const __m128i vi5x4567 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i5 + 4)));
280 const __m128i vk5x4567 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 124 * sizeof(int8_t))));
281 const __m128i vi5x89AB = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i5 + 8)));
282 const __m128i vk5x89AB = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 128 * sizeof(int8_t))));
283 const __m128i vi5xCDEF = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i5 + 12)));
284 const __m128i vk5xCDEF = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 132 * sizeof(int8_t))));
285 const __m128i vi5xGHIJ = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i5 + 16)));
286 const __m128i vk5xGHIJ = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 136 * sizeof(int8_t))));
287 const __m128i vi5xKLMN = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i5 + 20)));
288 const __m128i vk5xKLMN = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 140 * sizeof(int8_t))));
289 i5 += 24;
290
291 vacc0123 = _mm_add_epi32(vacc0123, _mm_mullo_epi32(vi5x0123, vk5x0123));
292 vacc4567 = _mm_add_epi32(vacc4567, _mm_mullo_epi32(vi5x4567, vk5x4567));
293 vacc89AB = _mm_add_epi32(vacc89AB, _mm_mullo_epi32(vi5x89AB, vk5x89AB));
294 vaccCDEF = _mm_add_epi32(vaccCDEF, _mm_mullo_epi32(vi5xCDEF, vk5xCDEF));
295 vaccGHIJ = _mm_add_epi32(vaccGHIJ, _mm_mullo_epi32(vi5xGHIJ, vk5xGHIJ));
296 vaccKLMN = _mm_add_epi32(vaccKLMN, _mm_mullo_epi32(vi5xKLMN, vk5xKLMN));
297
298 const __m128i vi6x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32(i6));
299 const __m128i vk6x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 144 * sizeof(int8_t))));
300 const __m128i vi6x4567 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i6 + 4)));
301 const __m128i vk6x4567 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 148 * sizeof(int8_t))));
302 const __m128i vi6x89AB = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i6 + 8)));
303 const __m128i vk6x89AB = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 152 * sizeof(int8_t))));
304 const __m128i vi6xCDEF = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i6 + 12)));
305 const __m128i vk6xCDEF = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 156 * sizeof(int8_t))));
306 const __m128i vi6xGHIJ = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i6 + 16)));
307 const __m128i vk6xGHIJ = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 160 * sizeof(int8_t))));
308 const __m128i vi6xKLMN = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i6 + 20)));
309 const __m128i vk6xKLMN = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 164 * sizeof(int8_t))));
310 i6 += 24;
311
312 vacc0123 = _mm_add_epi32(vacc0123, _mm_mullo_epi32(vi6x0123, vk6x0123));
313 vacc4567 = _mm_add_epi32(vacc4567, _mm_mullo_epi32(vi6x4567, vk6x4567));
314 vacc89AB = _mm_add_epi32(vacc89AB, _mm_mullo_epi32(vi6x89AB, vk6x89AB));
315 vaccCDEF = _mm_add_epi32(vaccCDEF, _mm_mullo_epi32(vi6xCDEF, vk6xCDEF));
316 vaccGHIJ = _mm_add_epi32(vaccGHIJ, _mm_mullo_epi32(vi6xGHIJ, vk6xGHIJ));
317 vaccKLMN = _mm_add_epi32(vaccKLMN, _mm_mullo_epi32(vi6xKLMN, vk6xKLMN));
318
319 const __m128i vi7x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32(i7));
320 const __m128i vk7x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 168 * sizeof(int8_t))));
321 const __m128i vi7x4567 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i7 + 4)));
322 const __m128i vk7x4567 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 172 * sizeof(int8_t))));
323 const __m128i vi7x89AB = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i7 + 8)));
324 const __m128i vk7x89AB = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 176 * sizeof(int8_t))));
325 const __m128i vi7xCDEF = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i7 + 12)));
326 const __m128i vk7xCDEF = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 180 * sizeof(int8_t))));
327 const __m128i vi7xGHIJ = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i7 + 16)));
328 const __m128i vk7xGHIJ = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 184 * sizeof(int8_t))));
329 const __m128i vi7xKLMN = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i7 + 20)));
330 const __m128i vk7xKLMN = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 188 * sizeof(int8_t))));
331 i7 += 24;
332
333 vacc0123 = _mm_add_epi32(vacc0123, _mm_mullo_epi32(vi7x0123, vk7x0123));
334 vacc4567 = _mm_add_epi32(vacc4567, _mm_mullo_epi32(vi7x4567, vk7x4567));
335 vacc89AB = _mm_add_epi32(vacc89AB, _mm_mullo_epi32(vi7x89AB, vk7x89AB));
336 vaccCDEF = _mm_add_epi32(vaccCDEF, _mm_mullo_epi32(vi7xCDEF, vk7xCDEF));
337 vaccGHIJ = _mm_add_epi32(vaccGHIJ, _mm_mullo_epi32(vi7xGHIJ, vk7xGHIJ));
338 vaccKLMN = _mm_add_epi32(vaccKLMN, _mm_mullo_epi32(vi7xKLMN, vk7xKLMN));
339
340 const __m128i vi8x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32(i8));
341 const __m128i vk8x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 192 * sizeof(int8_t))));
342 const __m128i vi8x4567 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i8 + 4)));
343 const __m128i vk8x4567 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 196 * sizeof(int8_t))));
344 const __m128i vi8x89AB = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i8 + 8)));
345 const __m128i vk8x89AB = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 200 * sizeof(int8_t))));
346 const __m128i vi8xCDEF = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i8 + 12)));
347 const __m128i vk8xCDEF = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 204 * sizeof(int8_t))));
348 const __m128i vi8xGHIJ = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i8 + 16)));
349 const __m128i vk8xGHIJ = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 208 * sizeof(int8_t))));
350 const __m128i vi8xKLMN = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i8 + 20)));
351 const __m128i vk8xKLMN = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 212 * sizeof(int8_t))));
352 i8 += 24;
353
354 vacc0123 = _mm_add_epi32(vacc0123, _mm_mullo_epi32(vi8x0123, vk8x0123));
355 vacc4567 = _mm_add_epi32(vacc4567, _mm_mullo_epi32(vi8x4567, vk8x4567));
356 vacc89AB = _mm_add_epi32(vacc89AB, _mm_mullo_epi32(vi8x89AB, vk8x89AB));
357 vaccCDEF = _mm_add_epi32(vaccCDEF, _mm_mullo_epi32(vi8xCDEF, vk8xCDEF));
358 vaccGHIJ = _mm_add_epi32(vaccGHIJ, _mm_mullo_epi32(vi8xGHIJ, vk8xGHIJ));
359 vaccKLMN = _mm_add_epi32(vaccKLMN, _mm_mullo_epi32(vi8xKLMN, vk8xKLMN));
360
361 const __m128i vi9x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32(i9));
362 const __m128i vk9x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 216 * sizeof(int8_t))));
363 const __m128i vi9x4567 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i9 + 4)));
364 const __m128i vk9x4567 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 220 * sizeof(int8_t))));
365 const __m128i vi9x89AB = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i9 + 8)));
366 const __m128i vk9x89AB = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 224 * sizeof(int8_t))));
367 const __m128i vi9xCDEF = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i9 + 12)));
368 const __m128i vk9xCDEF = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 228 * sizeof(int8_t))));
369 const __m128i vi9xGHIJ = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i9 + 16)));
370 const __m128i vk9xGHIJ = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 232 * sizeof(int8_t))));
371 const __m128i vi9xKLMN = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i9 + 20)));
372 const __m128i vk9xKLMN = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 236 * sizeof(int8_t))));
373 i9 += 24;
374
375 vacc0123 = _mm_add_epi32(vacc0123, _mm_mullo_epi32(vi9x0123, vk9x0123));
376 vacc4567 = _mm_add_epi32(vacc4567, _mm_mullo_epi32(vi9x4567, vk9x4567));
377 vacc89AB = _mm_add_epi32(vacc89AB, _mm_mullo_epi32(vi9x89AB, vk9x89AB));
378 vaccCDEF = _mm_add_epi32(vaccCDEF, _mm_mullo_epi32(vi9xCDEF, vk9xCDEF));
379 vaccGHIJ = _mm_add_epi32(vaccGHIJ, _mm_mullo_epi32(vi9xGHIJ, vk9xGHIJ));
380 vaccKLMN = _mm_add_epi32(vaccKLMN, _mm_mullo_epi32(vi9xKLMN, vk9xKLMN));
381
382 const __m128i vi10x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32(i10));
383 const __m128i vk10x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 240 * sizeof(int8_t))));
384 const __m128i vi10x4567 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i10 + 4)));
385 const __m128i vk10x4567 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 244 * sizeof(int8_t))));
386 const __m128i vi10x89AB = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i10 + 8)));
387 const __m128i vk10x89AB = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 248 * sizeof(int8_t))));
388 const __m128i vi10xCDEF = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i10 + 12)));
389 const __m128i vk10xCDEF = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 252 * sizeof(int8_t))));
390 const __m128i vi10xGHIJ = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i10 + 16)));
391 const __m128i vk10xGHIJ = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 256 * sizeof(int8_t))));
392 const __m128i vi10xKLMN = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i10 + 20)));
393 const __m128i vk10xKLMN = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 260 * sizeof(int8_t))));
394 i10 += 24;
395
396 vacc0123 = _mm_add_epi32(vacc0123, _mm_mullo_epi32(vi10x0123, vk10x0123));
397 vacc4567 = _mm_add_epi32(vacc4567, _mm_mullo_epi32(vi10x4567, vk10x4567));
398 vacc89AB = _mm_add_epi32(vacc89AB, _mm_mullo_epi32(vi10x89AB, vk10x89AB));
399 vaccCDEF = _mm_add_epi32(vaccCDEF, _mm_mullo_epi32(vi10xCDEF, vk10xCDEF));
400 vaccGHIJ = _mm_add_epi32(vaccGHIJ, _mm_mullo_epi32(vi10xGHIJ, vk10xGHIJ));
401 vaccKLMN = _mm_add_epi32(vaccKLMN, _mm_mullo_epi32(vi10xKLMN, vk10xKLMN));
402
403 const __m128i vi11x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32(i11));
404 const __m128i vk11x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 264 * sizeof(int8_t))));
405 const __m128i vi11x4567 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i11 + 4)));
406 const __m128i vk11x4567 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 268 * sizeof(int8_t))));
407 const __m128i vi11x89AB = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i11 + 8)));
408 const __m128i vk11x89AB = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 272 * sizeof(int8_t))));
409 const __m128i vi11xCDEF = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i11 + 12)));
410 const __m128i vk11xCDEF = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 276 * sizeof(int8_t))));
411 const __m128i vi11xGHIJ = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i11 + 16)));
412 const __m128i vk11xGHIJ = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 280 * sizeof(int8_t))));
413 const __m128i vi11xKLMN = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i11 + 20)));
414 const __m128i vk11xKLMN = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 284 * sizeof(int8_t))));
415 i11 += 24;
416
417 vacc0123 = _mm_add_epi32(vacc0123, _mm_mullo_epi32(vi11x0123, vk11x0123));
418 vacc4567 = _mm_add_epi32(vacc4567, _mm_mullo_epi32(vi11x4567, vk11x4567));
419 vacc89AB = _mm_add_epi32(vacc89AB, _mm_mullo_epi32(vi11x89AB, vk11x89AB));
420 vaccCDEF = _mm_add_epi32(vaccCDEF, _mm_mullo_epi32(vi11xCDEF, vk11xCDEF));
421 vaccGHIJ = _mm_add_epi32(vaccGHIJ, _mm_mullo_epi32(vi11xGHIJ, vk11xGHIJ));
422 vaccKLMN = _mm_add_epi32(vaccKLMN, _mm_mullo_epi32(vi11xKLMN, vk11xKLMN));
423
424 const __m128i vi12x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32(i12));
425 const __m128i vk12x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 288 * sizeof(int8_t))));
426 const __m128i vi12x4567 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i12 + 4)));
427 const __m128i vk12x4567 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 292 * sizeof(int8_t))));
428 const __m128i vi12x89AB = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i12 + 8)));
429 const __m128i vk12x89AB = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 296 * sizeof(int8_t))));
430 const __m128i vi12xCDEF = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i12 + 12)));
431 const __m128i vk12xCDEF = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 300 * sizeof(int8_t))));
432 const __m128i vi12xGHIJ = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i12 + 16)));
433 const __m128i vk12xGHIJ = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 304 * sizeof(int8_t))));
434 const __m128i vi12xKLMN = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i12 + 20)));
435 const __m128i vk12xKLMN = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 308 * sizeof(int8_t))));
436 i12 += 24;
437
438 vacc0123 = _mm_add_epi32(vacc0123, _mm_mullo_epi32(vi12x0123, vk12x0123));
439 vacc4567 = _mm_add_epi32(vacc4567, _mm_mullo_epi32(vi12x4567, vk12x4567));
440 vacc89AB = _mm_add_epi32(vacc89AB, _mm_mullo_epi32(vi12x89AB, vk12x89AB));
441 vaccCDEF = _mm_add_epi32(vaccCDEF, _mm_mullo_epi32(vi12xCDEF, vk12xCDEF));
442 vaccGHIJ = _mm_add_epi32(vaccGHIJ, _mm_mullo_epi32(vi12xGHIJ, vk12xGHIJ));
443 vaccKLMN = _mm_add_epi32(vaccKLMN, _mm_mullo_epi32(vi12xKLMN, vk12xKLMN));
444
445 const __m128i vi13x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32(i13));
446 const __m128i vk13x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 312 * sizeof(int8_t))));
447 const __m128i vi13x4567 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i13 + 4)));
448 const __m128i vk13x4567 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 316 * sizeof(int8_t))));
449 const __m128i vi13x89AB = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i13 + 8)));
450 const __m128i vk13x89AB = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 320 * sizeof(int8_t))));
451 const __m128i vi13xCDEF = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i13 + 12)));
452 const __m128i vk13xCDEF = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 324 * sizeof(int8_t))));
453 const __m128i vi13xGHIJ = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i13 + 16)));
454 const __m128i vk13xGHIJ = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 328 * sizeof(int8_t))));
455 const __m128i vi13xKLMN = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i13 + 20)));
456 const __m128i vk13xKLMN = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 332 * sizeof(int8_t))));
457 i13 += 24;
458
459 vacc0123 = _mm_add_epi32(vacc0123, _mm_mullo_epi32(vi13x0123, vk13x0123));
460 vacc4567 = _mm_add_epi32(vacc4567, _mm_mullo_epi32(vi13x4567, vk13x4567));
461 vacc89AB = _mm_add_epi32(vacc89AB, _mm_mullo_epi32(vi13x89AB, vk13x89AB));
462 vaccCDEF = _mm_add_epi32(vaccCDEF, _mm_mullo_epi32(vi13xCDEF, vk13xCDEF));
463 vaccGHIJ = _mm_add_epi32(vaccGHIJ, _mm_mullo_epi32(vi13xGHIJ, vk13xGHIJ));
464 vaccKLMN = _mm_add_epi32(vaccKLMN, _mm_mullo_epi32(vi13xKLMN, vk13xKLMN));
465
466 const __m128i vi14x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32(i14));
467 const __m128i vk14x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 336 * sizeof(int8_t))));
468 const __m128i vi14x4567 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i14 + 4)));
469 const __m128i vk14x4567 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 340 * sizeof(int8_t))));
470 const __m128i vi14x89AB = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i14 + 8)));
471 const __m128i vk14x89AB = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 344 * sizeof(int8_t))));
472 const __m128i vi14xCDEF = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i14 + 12)));
473 const __m128i vk14xCDEF = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 348 * sizeof(int8_t))));
474 const __m128i vi14xGHIJ = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i14 + 16)));
475 const __m128i vk14xGHIJ = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 352 * sizeof(int8_t))));
476 const __m128i vi14xKLMN = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i14 + 20)));
477 const __m128i vk14xKLMN = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 356 * sizeof(int8_t))));
478 i14 += 24;
479
480 vacc0123 = _mm_add_epi32(vacc0123, _mm_mullo_epi32(vi14x0123, vk14x0123));
481 vacc4567 = _mm_add_epi32(vacc4567, _mm_mullo_epi32(vi14x4567, vk14x4567));
482 vacc89AB = _mm_add_epi32(vacc89AB, _mm_mullo_epi32(vi14x89AB, vk14x89AB));
483 vaccCDEF = _mm_add_epi32(vaccCDEF, _mm_mullo_epi32(vi14xCDEF, vk14xCDEF));
484 vaccGHIJ = _mm_add_epi32(vaccGHIJ, _mm_mullo_epi32(vi14xGHIJ, vk14xGHIJ));
485 vaccKLMN = _mm_add_epi32(vaccKLMN, _mm_mullo_epi32(vi14xKLMN, vk14xKLMN));
486
487 const __m128i vi15x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32(i15));
488 const __m128i vk15x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 360 * sizeof(int8_t))));
489 const __m128i vi15x4567 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i15 + 4)));
490 const __m128i vk15x4567 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 364 * sizeof(int8_t))));
491 const __m128i vi15x89AB = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i15 + 8)));
492 const __m128i vk15x89AB = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 368 * sizeof(int8_t))));
493 const __m128i vi15xCDEF = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i15 + 12)));
494 const __m128i vk15xCDEF = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 372 * sizeof(int8_t))));
495 const __m128i vi15xGHIJ = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i15 + 16)));
496 const __m128i vk15xGHIJ = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 376 * sizeof(int8_t))));
497 const __m128i vi15xKLMN = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i15 + 20)));
498 const __m128i vk15xKLMN = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 380 * sizeof(int8_t))));
499 i15 += 24;
500
501 vacc0123 = _mm_add_epi32(vacc0123, _mm_mullo_epi32(vi15x0123, vk15x0123));
502 vacc4567 = _mm_add_epi32(vacc4567, _mm_mullo_epi32(vi15x4567, vk15x4567));
503 vacc89AB = _mm_add_epi32(vacc89AB, _mm_mullo_epi32(vi15x89AB, vk15x89AB));
504 vaccCDEF = _mm_add_epi32(vaccCDEF, _mm_mullo_epi32(vi15xCDEF, vk15xCDEF));
505 vaccGHIJ = _mm_add_epi32(vaccGHIJ, _mm_mullo_epi32(vi15xGHIJ, vk15xGHIJ));
506 vaccKLMN = _mm_add_epi32(vaccKLMN, _mm_mullo_epi32(vi15xKLMN, vk15xKLMN));
507
508 const __m128i vi16x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32(i16));
509 const __m128i vk16x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 384 * sizeof(int8_t))));
510 const __m128i vi16x4567 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i16 + 4)));
511 const __m128i vk16x4567 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 388 * sizeof(int8_t))));
512 const __m128i vi16x89AB = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i16 + 8)));
513 const __m128i vk16x89AB = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 392 * sizeof(int8_t))));
514 const __m128i vi16xCDEF = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i16 + 12)));
515 const __m128i vk16xCDEF = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 396 * sizeof(int8_t))));
516 const __m128i vi16xGHIJ = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i16 + 16)));
517 const __m128i vk16xGHIJ = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 400 * sizeof(int8_t))));
518 const __m128i vi16xKLMN = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i16 + 20)));
519 const __m128i vk16xKLMN = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 404 * sizeof(int8_t))));
520 i16 += 24;
521
522 vacc0123 = _mm_add_epi32(vacc0123, _mm_mullo_epi32(vi16x0123, vk16x0123));
523 vacc4567 = _mm_add_epi32(vacc4567, _mm_mullo_epi32(vi16x4567, vk16x4567));
524 vacc89AB = _mm_add_epi32(vacc89AB, _mm_mullo_epi32(vi16x89AB, vk16x89AB));
525 vaccCDEF = _mm_add_epi32(vaccCDEF, _mm_mullo_epi32(vi16xCDEF, vk16xCDEF));
526 vaccGHIJ = _mm_add_epi32(vaccGHIJ, _mm_mullo_epi32(vi16xGHIJ, vk16xGHIJ));
527 vaccKLMN = _mm_add_epi32(vaccKLMN, _mm_mullo_epi32(vi16xKLMN, vk16xKLMN));
528
529 const __m128i vi17x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32(i17));
530 const __m128i vk17x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 408 * sizeof(int8_t))));
531 const __m128i vi17x4567 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i17 + 4)));
532 const __m128i vk17x4567 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 412 * sizeof(int8_t))));
533 const __m128i vi17x89AB = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i17 + 8)));
534 const __m128i vk17x89AB = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 416 * sizeof(int8_t))));
535 const __m128i vi17xCDEF = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i17 + 12)));
536 const __m128i vk17xCDEF = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 420 * sizeof(int8_t))));
537 const __m128i vi17xGHIJ = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i17 + 16)));
538 const __m128i vk17xGHIJ = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 424 * sizeof(int8_t))));
539 const __m128i vi17xKLMN = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i17 + 20)));
540 const __m128i vk17xKLMN = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 428 * sizeof(int8_t))));
541 i17 += 24;
542
543 vacc0123 = _mm_add_epi32(vacc0123, _mm_mullo_epi32(vi17x0123, vk17x0123));
544 vacc4567 = _mm_add_epi32(vacc4567, _mm_mullo_epi32(vi17x4567, vk17x4567));
545 vacc89AB = _mm_add_epi32(vacc89AB, _mm_mullo_epi32(vi17x89AB, vk17x89AB));
546 vaccCDEF = _mm_add_epi32(vaccCDEF, _mm_mullo_epi32(vi17xCDEF, vk17xCDEF));
547 vaccGHIJ = _mm_add_epi32(vaccGHIJ, _mm_mullo_epi32(vi17xGHIJ, vk17xGHIJ));
548 vaccKLMN = _mm_add_epi32(vaccKLMN, _mm_mullo_epi32(vi17xKLMN, vk17xKLMN));
549
550 const __m128i vi18x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32(i18));
551 const __m128i vk18x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 432 * sizeof(int8_t))));
552 const __m128i vi18x4567 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i18 + 4)));
553 const __m128i vk18x4567 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 436 * sizeof(int8_t))));
554 const __m128i vi18x89AB = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i18 + 8)));
555 const __m128i vk18x89AB = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 440 * sizeof(int8_t))));
556 const __m128i vi18xCDEF = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i18 + 12)));
557 const __m128i vk18xCDEF = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 444 * sizeof(int8_t))));
558 const __m128i vi18xGHIJ = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i18 + 16)));
559 const __m128i vk18xGHIJ = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 448 * sizeof(int8_t))));
560 const __m128i vi18xKLMN = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i18 + 20)));
561 const __m128i vk18xKLMN = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 452 * sizeof(int8_t))));
562 i18 += 24;
563
564 vacc0123 = _mm_add_epi32(vacc0123, _mm_mullo_epi32(vi18x0123, vk18x0123));
565 vacc4567 = _mm_add_epi32(vacc4567, _mm_mullo_epi32(vi18x4567, vk18x4567));
566 vacc89AB = _mm_add_epi32(vacc89AB, _mm_mullo_epi32(vi18x89AB, vk18x89AB));
567 vaccCDEF = _mm_add_epi32(vaccCDEF, _mm_mullo_epi32(vi18xCDEF, vk18xCDEF));
568 vaccGHIJ = _mm_add_epi32(vaccGHIJ, _mm_mullo_epi32(vi18xGHIJ, vk18xGHIJ));
569 vaccKLMN = _mm_add_epi32(vaccKLMN, _mm_mullo_epi32(vi18xKLMN, vk18xKLMN));
570
571 const __m128i vi19x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32(i19));
572 const __m128i vk19x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 456 * sizeof(int8_t))));
573 const __m128i vi19x4567 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i19 + 4)));
574 const __m128i vk19x4567 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 460 * sizeof(int8_t))));
575 const __m128i vi19x89AB = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i19 + 8)));
576 const __m128i vk19x89AB = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 464 * sizeof(int8_t))));
577 const __m128i vi19xCDEF = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i19 + 12)));
578 const __m128i vk19xCDEF = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 468 * sizeof(int8_t))));
579 const __m128i vi19xGHIJ = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i19 + 16)));
580 const __m128i vk19xGHIJ = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 472 * sizeof(int8_t))));
581 const __m128i vi19xKLMN = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i19 + 20)));
582 const __m128i vk19xKLMN = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 476 * sizeof(int8_t))));
583 i19 += 24;
584
585 vacc0123 = _mm_add_epi32(vacc0123, _mm_mullo_epi32(vi19x0123, vk19x0123));
586 vacc4567 = _mm_add_epi32(vacc4567, _mm_mullo_epi32(vi19x4567, vk19x4567));
587 vacc89AB = _mm_add_epi32(vacc89AB, _mm_mullo_epi32(vi19x89AB, vk19x89AB));
588 vaccCDEF = _mm_add_epi32(vaccCDEF, _mm_mullo_epi32(vi19xCDEF, vk19xCDEF));
589 vaccGHIJ = _mm_add_epi32(vaccGHIJ, _mm_mullo_epi32(vi19xGHIJ, vk19xGHIJ));
590 vaccKLMN = _mm_add_epi32(vaccKLMN, _mm_mullo_epi32(vi19xKLMN, vk19xKLMN));
591
592 const __m128i vi20x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32(i20));
593 const __m128i vk20x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 480 * sizeof(int8_t))));
594 const __m128i vi20x4567 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i20 + 4)));
595 const __m128i vk20x4567 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 484 * sizeof(int8_t))));
596 const __m128i vi20x89AB = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i20 + 8)));
597 const __m128i vk20x89AB = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 488 * sizeof(int8_t))));
598 const __m128i vi20xCDEF = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i20 + 12)));
599 const __m128i vk20xCDEF = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 492 * sizeof(int8_t))));
600 const __m128i vi20xGHIJ = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i20 + 16)));
601 const __m128i vk20xGHIJ = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 496 * sizeof(int8_t))));
602 const __m128i vi20xKLMN = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i20 + 20)));
603 const __m128i vk20xKLMN = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 500 * sizeof(int8_t))));
604 i20 += 24;
605
606 vacc0123 = _mm_add_epi32(vacc0123, _mm_mullo_epi32(vi20x0123, vk20x0123));
607 vacc4567 = _mm_add_epi32(vacc4567, _mm_mullo_epi32(vi20x4567, vk20x4567));
608 vacc89AB = _mm_add_epi32(vacc89AB, _mm_mullo_epi32(vi20x89AB, vk20x89AB));
609 vaccCDEF = _mm_add_epi32(vaccCDEF, _mm_mullo_epi32(vi20xCDEF, vk20xCDEF));
610 vaccGHIJ = _mm_add_epi32(vaccGHIJ, _mm_mullo_epi32(vi20xGHIJ, vk20xGHIJ));
611 vaccKLMN = _mm_add_epi32(vaccKLMN, _mm_mullo_epi32(vi20xKLMN, vk20xKLMN));
612
613 const __m128i vi21x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32(i21));
614 const __m128i vk21x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 504 * sizeof(int8_t))));
615 const __m128i vi21x4567 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i21 + 4)));
616 const __m128i vk21x4567 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 508 * sizeof(int8_t))));
617 const __m128i vi21x89AB = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i21 + 8)));
618 const __m128i vk21x89AB = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 512 * sizeof(int8_t))));
619 const __m128i vi21xCDEF = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i21 + 12)));
620 const __m128i vk21xCDEF = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 516 * sizeof(int8_t))));
621 const __m128i vi21xGHIJ = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i21 + 16)));
622 const __m128i vk21xGHIJ = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 520 * sizeof(int8_t))));
623 const __m128i vi21xKLMN = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i21 + 20)));
624 const __m128i vk21xKLMN = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 524 * sizeof(int8_t))));
625 i21 += 24;
626
627 vacc0123 = _mm_add_epi32(vacc0123, _mm_mullo_epi32(vi21x0123, vk21x0123));
628 vacc4567 = _mm_add_epi32(vacc4567, _mm_mullo_epi32(vi21x4567, vk21x4567));
629 vacc89AB = _mm_add_epi32(vacc89AB, _mm_mullo_epi32(vi21x89AB, vk21x89AB));
630 vaccCDEF = _mm_add_epi32(vaccCDEF, _mm_mullo_epi32(vi21xCDEF, vk21xCDEF));
631 vaccGHIJ = _mm_add_epi32(vaccGHIJ, _mm_mullo_epi32(vi21xGHIJ, vk21xGHIJ));
632 vaccKLMN = _mm_add_epi32(vaccKLMN, _mm_mullo_epi32(vi21xKLMN, vk21xKLMN));
633
634 const __m128i vi22x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32(i22));
635 const __m128i vk22x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 528 * sizeof(int8_t))));
636 const __m128i vi22x4567 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i22 + 4)));
637 const __m128i vk22x4567 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 532 * sizeof(int8_t))));
638 const __m128i vi22x89AB = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i22 + 8)));
639 const __m128i vk22x89AB = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 536 * sizeof(int8_t))));
640 const __m128i vi22xCDEF = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i22 + 12)));
641 const __m128i vk22xCDEF = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 540 * sizeof(int8_t))));
642 const __m128i vi22xGHIJ = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i22 + 16)));
643 const __m128i vk22xGHIJ = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 544 * sizeof(int8_t))));
644 const __m128i vi22xKLMN = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i22 + 20)));
645 const __m128i vk22xKLMN = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 548 * sizeof(int8_t))));
646 i22 += 24;
647
648 vacc0123 = _mm_add_epi32(vacc0123, _mm_mullo_epi32(vi22x0123, vk22x0123));
649 vacc4567 = _mm_add_epi32(vacc4567, _mm_mullo_epi32(vi22x4567, vk22x4567));
650 vacc89AB = _mm_add_epi32(vacc89AB, _mm_mullo_epi32(vi22x89AB, vk22x89AB));
651 vaccCDEF = _mm_add_epi32(vaccCDEF, _mm_mullo_epi32(vi22xCDEF, vk22xCDEF));
652 vaccGHIJ = _mm_add_epi32(vaccGHIJ, _mm_mullo_epi32(vi22xGHIJ, vk22xGHIJ));
653 vaccKLMN = _mm_add_epi32(vaccKLMN, _mm_mullo_epi32(vi22xKLMN, vk22xKLMN));
654
655 const __m128i vi23x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32(i23));
656 const __m128i vk23x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 552 * sizeof(int8_t))));
657 const __m128i vi23x4567 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i23 + 4)));
658 const __m128i vk23x4567 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 556 * sizeof(int8_t))));
659 const __m128i vi23x89AB = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i23 + 8)));
660 const __m128i vk23x89AB = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 560 * sizeof(int8_t))));
661 const __m128i vi23xCDEF = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i23 + 12)));
662 const __m128i vk23xCDEF = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 564 * sizeof(int8_t))));
663 const __m128i vi23xGHIJ = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i23 + 16)));
664 const __m128i vk23xGHIJ = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 568 * sizeof(int8_t))));
665 const __m128i vi23xKLMN = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i23 + 20)));
666 const __m128i vk23xKLMN = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 572 * sizeof(int8_t))));
667 i23 += 24;
668
669 vacc0123 = _mm_add_epi32(vacc0123, _mm_mullo_epi32(vi23x0123, vk23x0123));
670 vacc4567 = _mm_add_epi32(vacc4567, _mm_mullo_epi32(vi23x4567, vk23x4567));
671 vacc89AB = _mm_add_epi32(vacc89AB, _mm_mullo_epi32(vi23x89AB, vk23x89AB));
672 vaccCDEF = _mm_add_epi32(vaccCDEF, _mm_mullo_epi32(vi23xCDEF, vk23xCDEF));
673 vaccGHIJ = _mm_add_epi32(vaccGHIJ, _mm_mullo_epi32(vi23xGHIJ, vk23xGHIJ));
674 vaccKLMN = _mm_add_epi32(vaccKLMN, _mm_mullo_epi32(vi23xKLMN, vk23xKLMN));
675
676 const __m128i vi24x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32(i24));
677 const __m128i vk24x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 576 * sizeof(int8_t))));
678 const __m128i vi24x4567 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i24 + 4)));
679 const __m128i vk24x4567 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 580 * sizeof(int8_t))));
680 const __m128i vi24x89AB = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i24 + 8)));
681 const __m128i vk24x89AB = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 584 * sizeof(int8_t))));
682 const __m128i vi24xCDEF = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i24 + 12)));
683 const __m128i vk24xCDEF = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 588 * sizeof(int8_t))));
684 const __m128i vi24xGHIJ = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i24 + 16)));
685 const __m128i vk24xGHIJ = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 592 * sizeof(int8_t))));
686 const __m128i vi24xKLMN = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (i24 + 20)));
687 const __m128i vk24xKLMN = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 596 * sizeof(int8_t))));
688 i24 += 24;
689
690 vacc0123 = _mm_add_epi32(vacc0123, _mm_mullo_epi32(vi24x0123, vk24x0123));
691 vacc4567 = _mm_add_epi32(vacc4567, _mm_mullo_epi32(vi24x4567, vk24x4567));
692 vacc89AB = _mm_add_epi32(vacc89AB, _mm_mullo_epi32(vi24x89AB, vk24x89AB));
693 vaccCDEF = _mm_add_epi32(vaccCDEF, _mm_mullo_epi32(vi24xCDEF, vk24xCDEF));
694 vaccGHIJ = _mm_add_epi32(vaccGHIJ, _mm_mullo_epi32(vi24xGHIJ, vk24xGHIJ));
695 vaccKLMN = _mm_add_epi32(vaccKLMN, _mm_mullo_epi32(vi24xKLMN, vk24xKLMN));
696
697 w = (const void*) ((uintptr_t) w + 24 * sizeof(int32_t) + 600 * sizeof(int8_t));
698
699 __m128 vscaled0123 = _mm_cvtepi32_ps(vacc0123);
700 __m128 vscaled4567 = _mm_cvtepi32_ps(vacc4567);
701 __m128 vscaled89AB = _mm_cvtepi32_ps(vacc89AB);
702 __m128 vscaledCDEF = _mm_cvtepi32_ps(vaccCDEF);
703 __m128 vscaledGHIJ = _mm_cvtepi32_ps(vaccGHIJ);
704 __m128 vscaledKLMN = _mm_cvtepi32_ps(vaccKLMN);
705
706 const __m128 vscale0123 = _mm_loadu_ps((const float*) w);
Marat Dukhane5eee462021-07-01 19:34:39 -0700707 const __m128 vscale4567 = _mm_loadu_ps((const float*) w + 4);
708 const __m128 vscale89AB = _mm_loadu_ps((const float*) w + 8);
709 const __m128 vscaleCDEF = _mm_loadu_ps((const float*) w + 12);
710 const __m128 vscaleGHIJ = _mm_loadu_ps((const float*) w + 16);
711 const __m128 vscaleKLMN = _mm_loadu_ps((const float*) w + 20);
712 w = (const void*) ((const float*) w + 24);
Marat Dukhan98042f22021-06-15 00:43:13 -0700713 vscaled0123 = _mm_mul_ps(vscaled0123, vscale0123);
714 vscaled4567 = _mm_mul_ps(vscaled4567, vscale4567);
715 vscaled89AB = _mm_mul_ps(vscaled89AB, vscale89AB);
716 vscaledCDEF = _mm_mul_ps(vscaledCDEF, vscaleCDEF);
717 vscaledGHIJ = _mm_mul_ps(vscaledGHIJ, vscaleGHIJ);
718 vscaledKLMN = _mm_mul_ps(vscaledKLMN, vscaleKLMN);
719
720 vacc0123 = _mm_cvtps_epi32(vscaled0123);
721 vacc4567 = _mm_cvtps_epi32(vscaled4567);
722 vacc89AB = _mm_cvtps_epi32(vscaled89AB);
723 vaccCDEF = _mm_cvtps_epi32(vscaledCDEF);
724 vaccGHIJ = _mm_cvtps_epi32(vscaledGHIJ);
725 vaccKLMN = _mm_cvtps_epi32(vscaledKLMN);
726
727 const __m128i voutput_zero_point = _mm_load_si128((const __m128i*) params->sse4.output_zero_point);
728 __m128i vout01234567 = _mm_adds_epi16(_mm_packs_epi32(vacc0123, vacc4567), voutput_zero_point);
729 __m128i vout89ABCDEF = _mm_adds_epi16(_mm_packs_epi32(vacc89AB, vaccCDEF), voutput_zero_point);
730 __m128i voutGHIJKLMN = _mm_adds_epi16(_mm_packs_epi32(vaccGHIJ, vaccKLMN), voutput_zero_point);
731
732 const __m128i voutput_min = _mm_load_si128((const __m128i*) params->sse4.output_min);
733 const __m128i voutput_max = _mm_load_si128((const __m128i*) params->sse4.output_max);
734 __m128i vout0123456789ABCDEF = _mm_packs_epi16(vout01234567, vout89ABCDEF);
735 vout0123456789ABCDEF = _mm_max_epi8(vout0123456789ABCDEF, voutput_min);
736 vout0123456789ABCDEF = _mm_min_epi8(vout0123456789ABCDEF, voutput_max);
737 __m128i voutGHIJKLMNGHIJKLMN = _mm_packs_epi16(voutGHIJKLMN, voutGHIJKLMN);
738 voutGHIJKLMNGHIJKLMN = _mm_max_epi8(voutGHIJKLMNGHIJKLMN, voutput_min);
739 voutGHIJKLMNGHIJKLMN = _mm_min_epi8(voutGHIJKLMNGHIJKLMN, voutput_max);
740
741 _mm_storeu_si128((__m128i*) output, vout0123456789ABCDEF);
742 _mm_storel_epi64((__m128i*) (output + 16), voutGHIJKLMNGHIJKLMN);
743 output += 24;
744 }
745 if XNN_UNLIKELY(c != 0) {
Marat Dukhane5eee462021-07-01 19:34:39 -0700746 const int8_t* k = (const int8_t*) ((const int32_t*) w + 24);
Marat Dukhan98042f22021-06-15 00:43:13 -0700747 do {
748 __m128i vacc0123 = _mm_loadu_si128((const __m128i*) w);
749
750
751 const __m128i vi0x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32(i0));
752 const __m128i vk0x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32(k));
753 i0 += 4;
754
755 vacc0123 = _mm_add_epi32(vacc0123, _mm_mullo_epi32(vi0x0123, vk0x0123));
756
757 const __m128i vi1x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32(i1));
758 const __m128i vk1x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (k + 24)));
759 i1 += 4;
760
761 vacc0123 = _mm_add_epi32(vacc0123, _mm_mullo_epi32(vi1x0123, vk1x0123));
762
763 const __m128i vi2x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32(i2));
764 const __m128i vk2x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (k + 48)));
765 i2 += 4;
766
767 vacc0123 = _mm_add_epi32(vacc0123, _mm_mullo_epi32(vi2x0123, vk2x0123));
768
769 const __m128i vi3x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32(i3));
770 const __m128i vk3x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (k + 72)));
771 i3 += 4;
772
773 vacc0123 = _mm_add_epi32(vacc0123, _mm_mullo_epi32(vi3x0123, vk3x0123));
774
775 const __m128i vi4x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32(i4));
776 const __m128i vk4x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (k + 96)));
777 i4 += 4;
778
779 vacc0123 = _mm_add_epi32(vacc0123, _mm_mullo_epi32(vi4x0123, vk4x0123));
780
781 const __m128i vi5x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32(i5));
782 const __m128i vk5x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (k + 120)));
783 i5 += 4;
784
785 vacc0123 = _mm_add_epi32(vacc0123, _mm_mullo_epi32(vi5x0123, vk5x0123));
786
787 const __m128i vi6x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32(i6));
788 const __m128i vk6x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (k + 144)));
789 i6 += 4;
790
791 vacc0123 = _mm_add_epi32(vacc0123, _mm_mullo_epi32(vi6x0123, vk6x0123));
792
793 const __m128i vi7x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32(i7));
794 const __m128i vk7x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (k + 168)));
795 i7 += 4;
796
797 vacc0123 = _mm_add_epi32(vacc0123, _mm_mullo_epi32(vi7x0123, vk7x0123));
798
799 const __m128i vi8x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32(i8));
800 const __m128i vk8x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (k + 192)));
801 i8 += 4;
802
803 vacc0123 = _mm_add_epi32(vacc0123, _mm_mullo_epi32(vi8x0123, vk8x0123));
804
805 const __m128i vi9x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32(i9));
806 const __m128i vk9x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (k + 216)));
807 i9 += 4;
808
809 vacc0123 = _mm_add_epi32(vacc0123, _mm_mullo_epi32(vi9x0123, vk9x0123));
810
811 const __m128i vi10x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32(i10));
812 const __m128i vk10x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (k + 240)));
813 i10 += 4;
814
815 vacc0123 = _mm_add_epi32(vacc0123, _mm_mullo_epi32(vi10x0123, vk10x0123));
816
817 const __m128i vi11x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32(i11));
818 const __m128i vk11x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (k + 264)));
819 i11 += 4;
820
821 vacc0123 = _mm_add_epi32(vacc0123, _mm_mullo_epi32(vi11x0123, vk11x0123));
822
823 const __m128i vi12x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32(i12));
824 const __m128i vk12x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (k + 288)));
825 i12 += 4;
826
827 vacc0123 = _mm_add_epi32(vacc0123, _mm_mullo_epi32(vi12x0123, vk12x0123));
828
829 const __m128i vi13x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32(i13));
830 const __m128i vk13x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (k + 312)));
831 i13 += 4;
832
833 vacc0123 = _mm_add_epi32(vacc0123, _mm_mullo_epi32(vi13x0123, vk13x0123));
834
835 const __m128i vi14x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32(i14));
836 const __m128i vk14x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (k + 336)));
837 i14 += 4;
838
839 vacc0123 = _mm_add_epi32(vacc0123, _mm_mullo_epi32(vi14x0123, vk14x0123));
840
841 const __m128i vi15x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32(i15));
842 const __m128i vk15x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (k + 360)));
843 i15 += 4;
844
845 vacc0123 = _mm_add_epi32(vacc0123, _mm_mullo_epi32(vi15x0123, vk15x0123));
846
847 const __m128i vi16x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32(i16));
848 const __m128i vk16x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (k + 384)));
849 i16 += 4;
850
851 vacc0123 = _mm_add_epi32(vacc0123, _mm_mullo_epi32(vi16x0123, vk16x0123));
852
853 const __m128i vi17x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32(i17));
854 const __m128i vk17x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (k + 408)));
855 i17 += 4;
856
857 vacc0123 = _mm_add_epi32(vacc0123, _mm_mullo_epi32(vi17x0123, vk17x0123));
858
859 const __m128i vi18x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32(i18));
860 const __m128i vk18x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (k + 432)));
861 i18 += 4;
862
863 vacc0123 = _mm_add_epi32(vacc0123, _mm_mullo_epi32(vi18x0123, vk18x0123));
864
865 const __m128i vi19x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32(i19));
866 const __m128i vk19x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (k + 456)));
867 i19 += 4;
868
869 vacc0123 = _mm_add_epi32(vacc0123, _mm_mullo_epi32(vi19x0123, vk19x0123));
870
871 const __m128i vi20x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32(i20));
872 const __m128i vk20x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (k + 480)));
873 i20 += 4;
874
875 vacc0123 = _mm_add_epi32(vacc0123, _mm_mullo_epi32(vi20x0123, vk20x0123));
876
877 const __m128i vi21x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32(i21));
878 const __m128i vk21x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (k + 504)));
879 i21 += 4;
880
881 vacc0123 = _mm_add_epi32(vacc0123, _mm_mullo_epi32(vi21x0123, vk21x0123));
882
883 const __m128i vi22x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32(i22));
884 const __m128i vk22x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (k + 528)));
885 i22 += 4;
886
887 vacc0123 = _mm_add_epi32(vacc0123, _mm_mullo_epi32(vi22x0123, vk22x0123));
888
889 const __m128i vi23x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32(i23));
890 const __m128i vk23x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (k + 552)));
891 i23 += 4;
892
893 vacc0123 = _mm_add_epi32(vacc0123, _mm_mullo_epi32(vi23x0123, vk23x0123));
894
895 const __m128i vi24x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32(i24));
896 const __m128i vk24x0123 = _mm_cvtepi8_epi32(_mm_loadu_si32((const void*) (k + 576)));
897 i24 += 4;
898
899 vacc0123 = _mm_add_epi32(vacc0123, _mm_mullo_epi32(vi24x0123, vk24x0123));
900
901 k += 4;
902
903 __m128 vscaled0123 = _mm_cvtepi32_ps(vacc0123);
904 const __m128 vscale0123 = _mm_loadu_ps((const float*) ((uintptr_t) w + 24 * sizeof(int32_t) + 600 * sizeof(int8_t)));
905 vscaled0123 = _mm_mul_ps(vscaled0123, vscale0123);
906 vacc0123 = _mm_cvtps_epi32(vscaled0123);
907
Marat Dukhane5eee462021-07-01 19:34:39 -0700908 w = (const void*) ((const int32_t*) w + 4);
Marat Dukhan98042f22021-06-15 00:43:13 -0700909
910 const __m128i voutput_zero_point = _mm_load_si128((const __m128i*) params->sse4.output_zero_point);
911 __m128i vout0123 = _mm_adds_epi16(_mm_packs_epi32(vacc0123, vacc0123), voutput_zero_point);
912
913 vout0123 = _mm_packs_epi16(vout0123, vout0123);
914 vout0123 = _mm_max_epi8(vout0123, _mm_load_si128((const __m128i*) params->sse4.output_min));
915 vout0123 = _mm_min_epi8(vout0123, _mm_load_si128((const __m128i*) params->sse4.output_max));
916
917 if XNN_LIKELY(c >= 4) {
918 _mm_storeu_si32(output, vout0123);
919 output += 4;
920 c -= 4;
921 } else {
922 if (c & 2) {
923 *((uint16_t*) output) = (uint16_t) _mm_extract_epi16(vout0123, 0);
924 vout0123 = _mm_srli_epi32(vout0123, 16);
925 output += 2;
926 }
927 if (c & 1) {
928 *output = (int8_t) _mm_extract_epi8(vout0123, 0);
929 output += 1;
930 }
931 c = 0;
932 }
933 } while (c != 0);
934 }
935
936 output = (int8_t*) ((uintptr_t) output + output_increment);
937 } while (--output_width != 0);
938}