736
|
1 ///
|
|
2 /// Replacement of synth_1to1() with AMD's 3DNowEx(DSP)! SIMD operations support
|
|
3 ///
|
|
4 /// This code based 'decode_3dnow.s' by Syuuhei Kashiyama
|
|
5 /// <squash@mb.kcom.ne.jp>,only some types of changes have been made:
|
|
6 ///
|
|
7 /// - change function name for support 3DNowEx! automatic detect
|
|
8 ///
|
|
9 /// note: because K7 processors are an aggresive out-of-order three-way
|
|
10 /// superscalar ones instruction order is not significand for them.
|
|
11 ///
|
|
12 /// Modified by Nick Kurshev <nickols_k@mail.ru>
|
|
13 ///
|
|
14 / synth_1to1_3dnow works the same way as the c version of
|
|
15 / synth_1to1. this assembler code based 'decode-i586.s'
|
|
16 / (by Stefan Bieschewski <stb@acm.org>), two types of changes
|
|
17 / have been made:
|
|
18 / - use {MMX,3DNow!} instruction for reduce cpu
|
|
19 / - remove unused(?) local symbols
|
|
20 /
|
|
21 / useful sources of information on optimizing 3DNow! code include:
|
|
22 / AMD 3DNow! Technology Manual (Publication #21928)
|
|
23 / English: http://www.amd.com/K6/k6docs/pdf/21928d.pdf
|
|
24 / (Japanese: http://www.amd.com/japan/K6/k6docs/j21928c.pdf)
|
|
25 / AMD-K6-2 Processor Code Optimization Application Note (Publication #21924)
|
|
26 / English: http://www.amd.com/K6/k6docs/pdf/21924b.pdf
|
|
27 /
|
|
28 / This code was tested only AMD-K6-2 processor Linux systems,
|
|
29 / please tell me:
|
|
30 / - whether this code works on other 3DNow! capable processors
|
|
31 / (ex.IDT-C6-2) or not
|
|
32 / - whether this code works on other OSes or not
|
|
33 /
|
|
34 / by KIMURA Takuhiro <kim@hannah.ipc.miyakyo-u.ac.jp> - until 31.Mar.1998
|
|
35 / <kim@comtec.co.jp> - after 1.Apr.1998
|
|
36
|
|
37 / Enhancments for q-word operation by Michael Hipp
|
|
38
|
|
39 .bss
|
|
40 .comm buffs,4352,4
|
|
41 .data
|
|
42 .align 4
|
|
43 bo:
|
|
44 .long 1
|
|
45 .text
|
|
46 .globl synth_1to1_3dnowex
|
|
47 synth_1to1_3dnowex:
|
|
48 subl $12,%esp
|
|
49 pushl %ebp
|
|
50 pushl %edi
|
|
51 pushl %esi
|
|
52 pushl %ebx
|
|
53
|
|
54 movl 32(%esp),%eax
|
|
55 movl 40(%esp),%esi
|
|
56 movl $0,%edi
|
|
57 movl bo,%ebp
|
|
58 cmpl %edi,36(%esp)
|
|
59 jne .L48
|
|
60 decl %ebp
|
|
61 andl $15,%ebp
|
|
62 movl %ebp,bo
|
|
63 movl $buffs,%ecx
|
|
64 jmp .L49
|
|
65 .L48:
|
|
66 addl $2,%esi
|
|
67 movl $buffs+2176,%ecx
|
|
68 .L49:
|
|
69 testl $1,%ebp
|
|
70 je .L50
|
|
71 movl %ecx,%ebx
|
|
72 movl %ebp,16(%esp)
|
|
73 pushl %eax
|
|
74 movl 20(%esp),%edx
|
|
75 leal (%ebx,%edx,4),%eax
|
|
76 pushl %eax
|
|
77 movl 24(%esp),%eax
|
|
78 incl %eax
|
|
79 andl $15,%eax
|
|
80 leal 1088(,%eax,4),%eax
|
|
81 addl %ebx,%eax
|
|
82 jmp .L74
|
|
83 .L50:
|
|
84 leal 1088(%ecx),%ebx
|
|
85 leal 1(%ebp),%edx
|
|
86 movl %edx,16(%esp)
|
|
87 pushl %eax
|
|
88 leal 1092(%ecx,%ebp,4),%eax
|
|
89 pushl %eax
|
|
90 leal (%ecx,%ebp,4),%eax
|
|
91 .L74:
|
|
92 pushl %eax
|
|
93 call dct64_3dnowex
|
|
94 addl $12,%esp
|
|
95 movl 16(%esp),%edx
|
|
96 leal 0(,%edx,4),%edx
|
|
97 movl $decwin+64,%eax
|
|
98 movl %eax,%ecx
|
|
99 subl %edx,%ecx
|
|
100 movl $16,%ebp
|
|
101
|
|
102 .L55:
|
|
103 movq (%ecx),%mm4
|
|
104 movq (%ebx),%mm3
|
|
105 movq 8(%ecx),%mm0
|
|
106 movq 8(%ebx),%mm1
|
|
107 pfmul %mm3,%mm4
|
|
108
|
|
109 movq 16(%ecx),%mm2
|
|
110 pfmul %mm1,%mm0
|
|
111 movq 16(%ebx),%mm3
|
|
112 pfadd %mm0,%mm4
|
|
113
|
|
114 movq 24(%ecx),%mm0
|
|
115 pfmul %mm2,%mm3
|
|
116 movq 24(%ebx),%mm1
|
|
117 pfadd %mm3,%mm4
|
|
118
|
|
119 movq 32(%ecx),%mm2
|
|
120 pfmul %mm1,%mm0
|
|
121 movq 32(%ebx),%mm3
|
|
122 pfadd %mm0,%mm4
|
|
123
|
|
124 movq 40(%ecx),%mm0
|
|
125 pfmul %mm2,%mm3
|
|
126 movq 40(%ebx),%mm1
|
|
127 pfadd %mm3,%mm4
|
|
128
|
|
129 movq 48(%ecx),%mm2
|
|
130 pfmul %mm1,%mm0
|
|
131 movq 48(%ebx),%mm3
|
|
132 pfadd %mm0,%mm4
|
|
133
|
|
134 movq 56(%ecx),%mm0
|
|
135 pfmul %mm2,%mm3
|
|
136 movq 56(%ebx),%mm1
|
|
137 pfadd %mm3,%mm4
|
|
138
|
|
139 pfmul %mm1,%mm0
|
|
140 pfadd %mm0,%mm4
|
|
141
|
|
142 movq %mm4,%mm0
|
|
143 psrlq $32,%mm0
|
|
144 pfsub %mm0,%mm4
|
|
145
|
|
146 pf2id %mm4,%mm4
|
|
147 movd %mm4,%eax
|
|
148
|
|
149 sar $16,%eax
|
|
150 movw %ax,(%esi)
|
|
151
|
|
152 addl $64,%ebx
|
|
153 subl $-128,%ecx
|
|
154 addl $4,%esi
|
|
155 decl %ebp
|
|
156 jnz .L55
|
|
157
|
|
158 / --- end of loop 1 ---
|
|
159
|
|
160 movd (%ecx),%mm2
|
|
161 movd (%ebx),%mm1
|
|
162 pfmul %mm1,%mm2
|
|
163
|
|
164 movd 8(%ecx),%mm0
|
|
165 movd 8(%ebx),%mm1
|
|
166 pfmul %mm0,%mm1
|
|
167 pfadd %mm1,%mm2
|
|
168
|
|
169 movd 16(%ecx),%mm0
|
|
170 movd 16(%ebx),%mm1
|
|
171 pfmul %mm0,%mm1
|
|
172 pfadd %mm1,%mm2
|
|
173
|
|
174 movd 24(%ecx),%mm0
|
|
175 movd 24(%ebx),%mm1
|
|
176 pfmul %mm0,%mm1
|
|
177 pfadd %mm1,%mm2
|
|
178
|
|
179 movd 32(%ecx),%mm0
|
|
180 movd 32(%ebx),%mm1
|
|
181 pfmul %mm0,%mm1
|
|
182 pfadd %mm1,%mm2
|
|
183
|
|
184 movd 40(%ecx),%mm0
|
|
185 movd 40(%ebx),%mm1
|
|
186 pfmul %mm0,%mm1
|
|
187 pfadd %mm1,%mm2
|
|
188
|
|
189 movd 48(%ecx),%mm0
|
|
190 movd 48(%ebx),%mm1
|
|
191 pfmul %mm0,%mm1
|
|
192 pfadd %mm1,%mm2
|
|
193
|
|
194 movd 56(%ecx),%mm0
|
|
195 movd 56(%ebx),%mm1
|
|
196 pfmul %mm0,%mm1
|
|
197 pfadd %mm1,%mm2
|
|
198
|
|
199 pf2id %mm2,%mm2
|
|
200 movd %mm2,%eax
|
|
201
|
|
202 sar $16,%eax
|
|
203
|
|
204 movw %ax,(%esi)
|
|
205
|
|
206 addl $-64,%ebx
|
|
207 addl $4,%esi
|
|
208 addl $256,%ecx
|
|
209 movl $15,%ebp
|
|
210
|
|
211 .L68:
|
|
212 psubd %mm0,%mm0
|
|
213
|
|
214 movq (%ebx),%mm1
|
|
215 movq (%ecx),%mm2
|
|
216 pfmul %mm1,%mm2
|
|
217 pfsub %mm2,%mm0
|
|
218
|
|
219 movq 8(%ebx),%mm3
|
|
220 movq 8(%ecx),%mm4
|
|
221 pfmul %mm3,%mm4
|
|
222 pfsub %mm4,%mm0
|
|
223
|
|
224 movq 16(%ebx),%mm1
|
|
225 movq 16(%ecx),%mm2
|
|
226 pfmul %mm1,%mm2
|
|
227 pfsub %mm2,%mm0
|
|
228
|
|
229 movq 24(%ebx),%mm3
|
|
230 movq 24(%ecx),%mm4
|
|
231 pfmul %mm3,%mm4
|
|
232 pfsub %mm4,%mm0
|
|
233
|
|
234 movq 32(%ebx),%mm1
|
|
235 movq 32(%ecx),%mm2
|
|
236 pfmul %mm1,%mm2
|
|
237 pfsub %mm2,%mm0
|
|
238
|
|
239 movq 40(%ebx),%mm3
|
|
240 movq 40(%ecx),%mm4
|
|
241 pfmul %mm3,%mm4
|
|
242 pfsub %mm4,%mm0
|
|
243
|
|
244 movq 48(%ebx),%mm1
|
|
245 movq 48(%ecx),%mm2
|
|
246 pfmul %mm1,%mm2
|
|
247 pfsub %mm2,%mm0
|
|
248
|
|
249 movq 56(%ebx),%mm3
|
|
250 movq 56(%ecx),%mm4
|
|
251 pfmul %mm3,%mm4
|
|
252 pfsub %mm4,%mm0
|
|
253
|
|
254 pfacc %mm0,%mm0
|
|
255
|
|
256 pf2id %mm0,%mm0
|
|
257 movd %mm0,%eax
|
|
258
|
|
259 sar $16,%eax
|
|
260
|
|
261 movw %ax,(%esi)
|
|
262
|
|
263 addl $-64,%ebx
|
|
264 subl $-128,%ecx
|
|
265 addl $4,%esi
|
|
266 decl %ebp
|
|
267 jnz .L68
|
|
268
|
|
269 / --- end of loop 2
|
|
270
|
|
271 femms
|
|
272
|
|
273 movl %edi,%eax
|
|
274 popl %ebx
|
|
275 popl %esi
|
|
276 popl %edi
|
|
277 popl %ebp
|
|
278 addl $12,%esp
|
|
279 ret
|