/
synth_sse_s32.asm
299 lines (266 loc) · 5.16 KB
/
synth_sse_s32.asm
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
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
; 1 "synth_sse_s32.S"
; 1 "<built-in>"
; 1 "<command line>"
; 1 "synth_sse_s32.S"
; 9 "synth_sse_s32.S"
; 1 "mangle.h" 1
; 13 "mangle.h"
; 1 "config.h" 1
; 14 "mangle.h" 2
; 1 "intsym.h" 1
; 15 "mangle.h" 2
; 10 "synth_sse_s32.S" 2
; 28 "synth_sse_s32.S"
%include "asm_nasm.inc"
_sym_mangle synth_1to1_s32_sse_asm
SECTION .data
align 32
scale_s32:
dd 1199570944
dd 1199570944
dd 1199570944
dd 1199570944
align 16
maxmin_s32:
dd 1191182335
dd 1191182335
dd 1191182335
dd 1191182335
dd -956301312
dd -956301312
dd -956301312
dd -956301312
SECTION .text
align 16
GLOBAL synth_1to1_s32_sse_asm
synth_1to1_s32_sse_asm:
push dword ebp
mov dword ebp, esp
push dword ebx
push dword esi
pxor mm7, mm7
mov dword ebx, [ebp+8]
mov dword edx, [ebp+12]
mov dword esi, [ebp+16]
mov dword eax, [ebp+20]
shl dword eax, 2
lea dword ebx, [ebx+64]
sub dword ebx, eax
mov dword ecx, 4
align 16
Loop_start_1:
movups xmm0, [ebx]
movups xmm1, [ebx+16]
movups xmm2, [ebx+32]
movups xmm3, [ebx+48]
movups xmm4, [ebx+128]
movups xmm5, [ebx+144]
movups xmm6, [ebx+160]
movups xmm7, [ebx+176]
mulps xmm0, [edx+0]
mulps xmm1, [edx+16]
mulps xmm2, [edx+32]
mulps xmm3, [edx+48]
mulps xmm4, [edx+64]
mulps xmm5, [edx+80]
mulps xmm6, [edx+96]
mulps xmm7, [edx+112]
addps xmm0, xmm1
addps xmm2, xmm3
addps xmm4, xmm5
addps xmm6, xmm7
addps xmm0, xmm2
addps xmm4, xmm6
movaps xmm5, xmm4
movaps xmm4, xmm0
lea dword ebx, [ebx+256]
lea dword edx, [edx+128]
movups xmm0, [ebx]
movups xmm1, [ebx+16]
movups xmm2, [ebx+32]
movups xmm3, [ebx+48]
movups xmm6, [ebx+128]
movups xmm7, [ebx+144]
mulps xmm0, [edx]
mulps xmm1, [edx+16]
mulps xmm2, [edx+32]
mulps xmm3, [edx+48]
mulps xmm6, [edx+64]
mulps xmm7, [edx+80]
addps xmm0, xmm1
addps xmm2, xmm3
addps xmm6, xmm7
movups xmm1, [ebx+160]
movups xmm3, [ebx+176]
mulps xmm1, [edx+96]
mulps xmm3, [edx+112]
addps xmm0, xmm2
addps xmm1, xmm3
addps xmm6, xmm1
movaps xmm7, xmm6
movaps xmm6, xmm0
lea dword ebx, [ebx+256]
lea dword edx, [edx+128]
movaps xmm0, xmm4
movaps xmm1, xmm6
unpcklps xmm4, xmm5
unpcklps xmm6, xmm7
unpckhps xmm0, xmm5
unpckhps xmm1, xmm7
movaps xmm2, xmm4
movaps xmm3, xmm0
movlhps xmm4, xmm6
movhlps xmm6, xmm2
movlhps xmm0, xmm1
movhlps xmm1, xmm3
subps xmm4, xmm6
subps xmm0, xmm1
addps xmm0, xmm4
movaps xmm1, xmm0
movaps xmm2, xmm0
mulps xmm0, [scale_s32]
cmpnleps xmm1, [maxmin_s32]
cmpltps xmm2, [maxmin_s32+16]
cvtps2pi mm0, xmm0
movhlps xmm0, xmm0
cvtps2pi mm1, xmm0
cvtps2pi mm2, xmm1
movhlps xmm1, xmm1
cvtps2pi mm3, xmm1
psrad mm2, 31
psrad mm3, 31
pxor mm0, mm2
pxor mm1, mm3
movd [esi], mm0
psrlq mm0, 32
movd [esi+8], mm0
movd [esi+16], mm1
psrlq mm1, 32
movd [esi+24], mm1
cvtps2pi mm0, xmm2
movhlps xmm2, xmm2
cvtps2pi mm1, xmm2
packssdw mm2, mm3
packssdw mm0, mm1
psrlw mm2, 15
psrlw mm0, 15
paddw mm0, mm2
paddw mm7, mm0
lea dword esi, [esi+32]
dec dword ecx
jnz Loop_start_1
mov dword ecx, 4
align 16
Loop_start_2:
movups xmm0, [ebx]
movups xmm1, [ebx+16]
movups xmm2, [ebx+32]
movups xmm3, [ebx+48]
movups xmm4, [ebx+128]
movups xmm5, [ebx+144]
movups xmm6, [ebx+160]
movups xmm7, [ebx+176]
mulps xmm0, [edx+0]
mulps xmm1, [edx+16]
mulps xmm2, [edx+32]
mulps xmm3, [edx+48]
mulps xmm4, [edx+-64]
mulps xmm5, [edx+-48]
mulps xmm6, [edx+-32]
mulps xmm7, [edx+-16]
addps xmm0, xmm1
addps xmm2, xmm3
addps xmm4, xmm5
addps xmm6, xmm7
addps xmm0, xmm2
addps xmm4, xmm6
movaps xmm5, xmm4
movaps xmm4, xmm0
lea dword ebx, [ebx+256]
lea dword edx, [edx+-128]
movups xmm0, [ebx]
movups xmm1, [ebx+16]
movups xmm2, [ebx+32]
movups xmm3, [ebx+48]
movups xmm6, [ebx+128]
movups xmm7, [ebx+144]
mulps xmm0, [edx]
mulps xmm1, [edx+16]
mulps xmm2, [edx+32]
mulps xmm3, [edx+48]
mulps xmm6, [edx+-64]
mulps xmm7, [edx+-48]
addps xmm0, xmm1
addps xmm2, xmm3
addps xmm6, xmm7
movups xmm1, [ebx+160]
movups xmm3, [ebx+176]
mulps xmm1, [edx+-32]
mulps xmm3, [edx+-16]
addps xmm0, xmm2
addps xmm1, xmm3
addps xmm6, xmm1
movaps xmm7, xmm6
movaps xmm6, xmm0
lea dword ebx, [ebx+256]
lea dword edx, [edx+-128]
movaps xmm0, xmm4
movaps xmm1, xmm6
unpcklps xmm4, xmm5
unpcklps xmm6, xmm7
unpckhps xmm0, xmm5
unpckhps xmm1, xmm7
movaps xmm2, xmm4
movaps xmm3, xmm0
movlhps xmm4, xmm6
movhlps xmm6, xmm2
movlhps xmm0, xmm1
movhlps xmm1, xmm3
addps xmm4, xmm6
addps xmm0, xmm1
addps xmm0, xmm4
movaps xmm1, xmm0
movaps xmm2, xmm0
mulps xmm0, [scale_s32]
cmpnleps xmm1, [maxmin_s32]
cmpltps xmm2, [maxmin_s32+16]
cvtps2pi mm0, xmm0
movhlps xmm0, xmm0
cvtps2pi mm1, xmm0
cvtps2pi mm2, xmm1
movhlps xmm1, xmm1
cvtps2pi mm3, xmm1
psrad mm2, 31
psrad mm3, 31
pxor mm0, mm2
pxor mm1, mm3
movd [esi], mm0
psrlq mm0, 32
movd [esi+8], mm0
movd [esi+16], mm1
psrlq mm1, 32
movd [esi+24], mm1
cvtps2pi mm0, xmm2
movhlps xmm2, xmm2
cvtps2pi mm1, xmm2
packssdw mm2, mm3
packssdw mm0, mm1
psrlw mm2, 15
psrlw mm0, 15
paddw mm0, mm2
paddw mm7, mm0
lea dword esi, [esi+32]
dec dword ecx
jnz Loop_start_2
pshufw mm0, mm7, 0xee
paddw mm0, mm7
pshufw mm1, mm0, 0x55
paddw mm0, mm1
movd eax, mm0
and dword eax, 0xffff
pop dword esi
pop dword ebx
mov dword esp, ebp
pop dword ebp
emms
ret