[svn] / branches / dev-api-4 / xvidcore / src / dct / x86_asm / idct_3dne.asm Repository:
ViewVC logotype

Annotation of /branches/dev-api-4/xvidcore/src/dct/x86_asm/idct_3dne.asm

Parent Directory Parent Directory | Revision Log Revision Log

Revision 1192 - (view) (download)

1 : edgomez 1190 ;/****************************************************************************
2 :     ; *
3 :     ; * XVID MPEG-4 VIDEO CODEC
4 :     ; * - MMX and XMM forward discrete cosine transform -
5 : edgomez 1192 ; *
6 : edgomez 1190 ; * Copyright(C) 2001 Peter Ross <pross@xvid.org>
7 :     ; * 2002 Jaan Kalda
8 :     ; *
9 :     ; * This program is free software; you can redistribute it and/or modify it
10 :     ; * under the terms of the GNU General Public License as published by
11 :     ; * the Free Software Foundation; either version 2 of the License, or
12 :     ; * (at your option) any later version.
13 :     ; *
14 :     ; * This program is distributed in the hope that it will be useful,
15 :     ; * but WITHOUT ANY WARRANTY; without even the implied warranty of
16 :     ; * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
17 :     ; * GNU General Public License for more details.
18 :     ; *
19 :     ; * You should have received a copy of the GNU General Public License
20 :     ; * along with this program; if not, write to the Free Software
21 :     ; * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
22 :     ; *
23 : edgomez 1192 ; * $Id: idct_3dne.asm,v 1.2.2.2 2003-10-28 22:23:03 edgomez Exp $
24 : edgomez 1190 ; *
25 :     ; ***************************************************************************/
26 :    
27 :     ; ****************************************************************************
28 :     ;
29 : edgomez 851 ; Originally provided by Intel at AP-922
30 :     ; http://developer.intel.com/vtune/cbts/strmsimd/922down.htm
31 :     ; (See more app notes at http://developer.intel.com/vtune/cbts/strmsimd/appnotes.htm)
32 :     ; but in a limited edition.
33 :     ; New macro implements a column part for precise iDCT
34 : edgomez 1192 ; The routine precision now satisfies IEEE standard 1180-1990.
35 : edgomez 851 ;
36 : edgomez 1190 ; Copyright(C) 2000-2001 Peter Gubanov <peter@elecard.net.ru>
37 :     ; Rounding trick Copyright(C) 2000 Michel Lespinasse <walken@zoy.org>
38 : edgomez 851 ;
39 :     ; http://www.elecard.com/peter/idct.html
40 :     ; http://www.linuxvideo.org/mpeg2dec/
41 :     ;
42 : edgomez 1190 ; ***************************************************************************/
43 : edgomez 851 ;
44 :     ; These examples contain code fragments for first stage iDCT 8x8
45 :     ; (for rows) and first stage DCT 8x8 (for columns)
46 :     ;
47 :    
48 : edgomez 1190 ; ***************************************************************************/
49 : edgomez 1192 ; this 3dne function is compatible with iSSE, but is optimized specifically for
50 : edgomez 1190 ; K7 pipelines (ca 5% gain), for implementation details see the idct_mmx.asm
51 :     ; file
52 : edgomez 851 ;
53 : edgomez 1190 ; ----------------------------------------------------------------------------
54 : edgomez 1192 ; Athlon optimizations contributed by Jaan Kalda
55 : edgomez 1190 ;-----------------------------------------------------------------------------
56 : edgomez 851
57 : edgomez 1190 BITS 32
58 : edgomez 851
59 : edgomez 1190 ;=============================================================================
60 :     ; Macros and other preprocessor constants
61 :     ;=============================================================================
62 :    
63 : edgomez 1192 %macro cglobal 1
64 : edgomez 851 %ifdef PREFIX
65 : edgomez 1192 global _%1
66 : edgomez 851 %define %1 _%1
67 :     %else
68 :     global %1
69 :     %endif
70 :     %endmacro
71 :    
72 : edgomez 1190 %define BITS_INV_ACC 5 ; 4 or 5 for IEEE
73 : edgomez 851 %define SHIFT_INV_ROW 16 - BITS_INV_ACC
74 :     %define SHIFT_INV_COL 1 + BITS_INV_ACC
75 :     %define RND_INV_ROW 1024 * (6 - BITS_INV_ACC) ; 1 << (SHIFT_INV_ROW-1)
76 :     %define RND_INV_COL 16 * (BITS_INV_ACC - 3) ; 1 << (SHIFT_INV_COL-1)
77 : edgomez 1190 %define RND_INV_CORR RND_INV_COL - 1 ; correction -1.0 and round
78 : edgomez 851
79 : edgomez 1190 %define BITS_FRW_ACC 3 ; 2 or 3 for accuracy
80 : edgomez 851 %define SHIFT_FRW_COL BITS_FRW_ACC
81 :     %define SHIFT_FRW_ROW BITS_FRW_ACC + 17
82 : edgomez 1190 %define RND_FRW_ROW 262144*(BITS_FRW_ACC - 1) ; 1 << (SHIFT_FRW_ROW-1)
83 : edgomez 851
84 : edgomez 1190 ;=============================================================================
85 :     ; Local Data (Read Only)
86 :     ;=============================================================================
87 : edgomez 851
88 : edgomez 1190 SECTION .rodata
89 : edgomez 851
90 : edgomez 1190 ;-----------------------------------------------------------------------------
91 :     ; Various memory constants (trigonometric values or rounding values)
92 :     ;-----------------------------------------------------------------------------
93 : edgomez 851
94 : edgomez 1190 ALIGN 16
95 :     one_corr:
96 : edgomez 1192 dw 1, 1, 1, 1
97 : edgomez 1190 round_inv_row:
98 : edgomez 1192 dd RND_INV_ROW, RND_INV_ROW
99 : edgomez 1190 round_inv_col:
100 : edgomez 1192 dw RND_INV_COL, RND_INV_COL, RND_INV_COL, RND_INV_COL
101 : edgomez 1190 round_inv_corr:
102 : edgomez 1192 dw RND_INV_CORR, RND_INV_CORR, RND_INV_CORR, RND_INV_CORR
103 : edgomez 1190 round_frw_row:
104 : edgomez 1192 dd RND_FRW_ROW, RND_FRW_ROW
105 : edgomez 1190 tg_1_16:
106 : edgomez 1192 dw 13036, 13036, 13036, 13036 ; tg * (2<<16) + 0.5
107 : edgomez 1190 tg_2_16:
108 : edgomez 1192 dw 27146, 27146, 27146, 27146 ; tg * (2<<16) + 0.5
109 : edgomez 1190 tg_3_16:
110 : edgomez 1192 dw -21746, -21746, -21746, -21746 ; tg * (2<<16) + 0.5
111 : edgomez 1190 cos_4_16:
112 : edgomez 1192 dw -19195, -19195, -19195, -19195 ; cos * (2<<16) + 0.5
113 : edgomez 1190 ocos_4_16:
114 : edgomez 1192 dw 23170, 23170, 23170, 23170 ; cos * (2<<15) + 0.5
115 : edgomez 1190 otg_3_16:
116 : edgomez 1192 dw 21895, 21895, 21895, 21895 ; tg * (2<<16) + 0.5
117 : edgomez 851
118 :     %if SHIFT_INV_ROW == 12 ; assume SHIFT_INV_ROW == 12
119 : edgomez 1190 rounder_0:
120 : edgomez 1192 dd 65536, 65536
121 : edgomez 1190 rounder_4:
122 : edgomez 1192 dd 0, 0
123 : edgomez 1190 rounder_1:
124 : edgomez 1192 dd 7195, 7195
125 : edgomez 1190 rounder_7
126 : edgomez 1192 dd 1024, 1024
127 : edgomez 1190 rounder_2:
128 : edgomez 1192 dd 4520, 4520
129 : edgomez 1190 rounder_6:
130 : edgomez 1192 dd 1024, 1024
131 : edgomez 1190 rounder_3:
132 : edgomez 1192 dd 2407, 2407
133 : edgomez 1190 rounder_5:
134 : edgomez 1192 dd 240, 240
135 : edgomez 851
136 :     %elif SHIFT_INV_ROW == 11 ; assume SHIFT_INV_ROW == 11
137 : edgomez 1190 rounder_0:
138 : edgomez 1192 dd 65536, 65536
139 : edgomez 1190 rounder_4:
140 : edgomez 1192 dd 0, 0
141 : edgomez 1190 rounder_1:
142 : edgomez 1192 dd 3597, 3597
143 : edgomez 1190 rounder_7:
144 : edgomez 1192 dd 512, 512
145 : edgomez 1190 rounder_2:
146 : edgomez 1192 dd 2260, 2260
147 : edgomez 1190 rounder_6:
148 : edgomez 1192 dd 512, 512
149 : edgomez 1190 rounder_3:
150 : edgomez 1192 dd 1203, 1203
151 : edgomez 1190 rounder_5:
152 : edgomez 1192 dd 120, 120
153 : edgomez 851 %else
154 :    
155 : edgomez 1190 %error invalid SHIFT_INV_ROW
156 : edgomez 851
157 :     %endif
158 :    
159 :     ;-----------------------------------------------------------------------------
160 : edgomez 1190 ; Tables for xmm processors
161 : edgomez 851 ;-----------------------------------------------------------------------------
162 :    
163 :     ; %3 for rows 0,4 - constants are multiplied by cos_4_16
164 : edgomez 1190 tab_i_04_xmm:
165 : edgomez 1192 dw 16384, 21407, 16384, 8867 ; movq-> w05 w04 w01 w00
166 :     dw 16384, 8867, -16384, -21407 ; w07 w06 w03 w02
167 :     dw 16384, -8867, 16384, -21407 ; w13 w12 w09 w08
168 :     dw -16384, 21407, 16384, -8867 ; w15 w14 w11 w10
169 :     dw 22725, 19266, 19266, -4520 ; w21 w20 w17 w16
170 :     dw 12873, 4520, -22725, -12873 ; w23 w22 w19 w18
171 :     dw 12873, -22725, 4520, -12873 ; w29 w28 w25 w24
172 :     dw 4520, 19266, 19266, -22725 ; w31 w30 w27 w26
173 : edgomez 851
174 :     ; %3 for rows 1,7 - constants are multiplied by cos_1_16
175 : edgomez 1190 tab_i_17_xmm:
176 : edgomez 1192 dw 22725, 29692, 22725, 12299 ; movq-> w05 w04 w01 w00
177 :     dw 22725, 12299, -22725, -29692 ; w07 w06 w03 w02
178 :     dw 22725, -12299, 22725, -29692 ; w13 w12 w09 w08
179 :     dw -22725, 29692, 22725, -12299 ; w15 w14 w11 w10
180 :     dw 31521, 26722, 26722, -6270 ; w21 w20 w17 w16
181 :     dw 17855, 6270, -31521, -17855 ; w23 w22 w19 w18
182 :     dw 17855, -31521, 6270, -17855 ; w29 w28 w25 w24
183 :     dw 6270, 26722, 26722, -31521 ; w31 w30 w27 w26
184 : edgomez 851
185 :     ; %3 for rows 2,6 - constants are multiplied by cos_2_16
186 : edgomez 1190 tab_i_26_xmm:
187 : edgomez 1192 dw 21407, 27969, 21407, 11585 ; movq-> w05 w04 w01 w00
188 :     dw 21407, 11585, -21407, -27969 ; w07 w06 w03 w02
189 :     dw 21407, -11585, 21407, -27969 ; w13 w12 w09 w08
190 :     dw -21407, 27969, 21407, -11585 ; w15 w14 w11 w10
191 :     dw 29692, 25172, 25172, -5906 ; w21 w20 w17 w16
192 :     dw 16819, 5906, -29692, -16819 ; w23 w22 w19 w18
193 :     dw 16819, -29692, 5906, -16819 ; w29 w28 w25 w24
194 :     dw 5906, 25172, 25172, -29692 ; w31 w30 w27 w26
195 : edgomez 851
196 :     ; %3 for rows 3,5 - constants are multiplied by cos_3_16
197 : edgomez 1190 tab_i_35_xmm:
198 : edgomez 1192 dw 19266, 25172, 19266, 10426 ; movq-> w05 w04 w01 w00
199 :     dw 19266, 10426, -19266, -25172 ; w07 w06 w03 w02
200 :     dw 19266, -10426, 19266, -25172 ; w13 w12 w09 w08
201 :     dw -19266, 25172, 19266, -10426 ; w15 w14 w11 w10
202 :     dw 26722, 22654, 22654, -5315 ; w21 w20 w17 w16
203 :     dw 15137, 5315, -26722, -15137 ; w23 w22 w19 w18
204 :     dw 15137, -26722, 5315, -15137 ; w29 w28 w25 w24
205 :     dw 5315, 22654, 22654, -26722 ; w31 w30 w27 w26
206 : edgomez 851
207 :     ;=============================================================================
208 : edgomez 1190 ; Code
209 : edgomez 851 ;=============================================================================
210 :    
211 : edgomez 1190 SECTION .text
212 :    
213 : edgomez 1192 cglobal idct_3dne
214 :    
215 : edgomez 851 ;-----------------------------------------------------------------------------
216 : edgomez 1190 ; void idct_3dne(uint16_t block[64]);
217 : edgomez 851 ;-----------------------------------------------------------------------------
218 :    
219 : edgomez 1190 ALIGN 16
220 :     idct_3dne:
221 : edgomez 1192 mov eax, [esp+4]
222 : edgomez 851
223 : edgomez 1192 ; DCT_8_INV_ROW_1_s [eax+64], [eax+64], tab_i_04_sse, rounder_4 ;rounder_4=0
224 :     pshufw mm0, [eax+64],10001000b ; x2 x0 x2 x0
225 :     movq mm3, [tab_i_04_xmm] ; 3 ; w05 w04 w01 w00
226 :     pshufw mm1, [eax+64+8],10001000b ; x6 x4 x6 x4
227 :     movq mm4, [tab_i_04_xmm+8] ; 4 ; w07 w06 w03 w02
228 :     pshufw mm2, [eax+64],11011101b ; x3 x1 x3 x1
229 :     pshufw mm5, [eax+64+8],11011101b ; x7 x5 x7 x5
230 :     movq mm6, [tab_i_04_xmm+32] ; 6 ; w21 w20 w17 w16
231 :     pmaddwd mm3, mm0 ; x2*w05+x0*w04 x2*w01+x0*w00
232 :     movq mm7, [tab_i_04_xmm+40] ; 7 ; w23 w22 w19 w18 ;
233 :     pmaddwd mm0, [tab_i_04_xmm+16]; x2*w13+x0*w12 x2*w09+x0*w08
234 :     pmaddwd mm4, mm1 ; x6*w07+x4*w06 x6*w03+x4*w02
235 :     pmaddwd mm1, [tab_i_04_xmm+24]; x6*w15+x4*w14 x6*w11+x4*w10
236 :     pmaddwd mm6, mm2 ; x3*w21+x1*w20 x3*w17+x1*w16
237 :     pmaddwd mm2, [tab_i_04_xmm+48]; x3*w29+x1*w28 x3*w25+x1*w24
238 :     pmaddwd mm7, mm5 ; 7 ; x7*w23+x5*w22 x7*w19+x5*w18 ; w23 w22 w19 w18
239 :     pmaddwd mm5, [tab_i_04_xmm+56]; x7*w31+x5*w30 x7*w27+x5*w26
240 :     paddd mm3, mm4 ; 4 free ; a1=sum(even1) a0=sum(even0)
241 :     paddd mm0, mm1 ; 1 free ; a3=sum(even3) a2=sum(even2)
242 :     pshufw mm1, [eax+80+8],10001000b ; x6 x4 x6 x4
243 :     movq mm4, mm3 ; 4 ; a1 a0
244 :     paddd mm6, mm7 ; 7 free ; b1=sum(odd1) b0=sum(odd0)
245 :     paddd mm2, mm5 ; 5 free ; b3=sum(odd3) b2=sum(odd2)
246 :     pshufw mm5, [eax+80],10001000b; x2 x0 x2 x0 mm5 & mm0 exchanged for next cycle
247 :     movq mm7, mm0 ; 7 ; a3 a2
248 :     psubd mm4, mm6 ; 6 free ; a1-b1 a0-b0
249 :     paddd mm6, mm3 ; mm6 = mm3+mm6+mm5+mm4; a1+b1 a0+b0
250 :     movq mm3, [tab_i_35_xmm] ; 3 ; w05 w04 w01 w00
251 :     psubd mm7, mm2 ; ; a3-b3 a2-b2
252 :     paddd mm0, mm2 ; 0 free a3+b3 a2+b2
253 :     pshufw mm2, [eax+80],11011101b; x3 x1 x3 x1
254 :     pmaddwd mm3, mm5 ; x2*w05+x0*w04 x2*w01+x0*w00
255 :     pmaddwd mm5, [tab_i_35_xmm+16]; x2*w13+x0*w12 x2*w09+x0*w08
256 :     psrad mm4, SHIFT_INV_ROW ; y6=a1-b1 y7=a0-b0
257 :     psrad mm7, SHIFT_INV_ROW ; y4=a3-b3 y5=a2-b2
258 :     psrad mm6, SHIFT_INV_ROW ; y1=a1+b1 y0=a0+b0
259 :     psrad mm0, SHIFT_INV_ROW ; y3=a3+b3 y2=a2+b2
260 :     packssdw mm7, mm4 ; 4 ; y6 y7 y4 y5
261 :     packssdw mm6, mm0 ; 0 free ; y3 y2 y1 y0
262 :     pshufw mm0, [eax+80+8],11011101b ; x7 x5 x7 x5
263 :     movq [eax+64], mm6 ; 3 ; save y3 y2 y1 y0 stall2
264 : edgomez 851
265 : edgomez 1192 ; DCT_8_INV_ROW_1_s [eax+80], [eax+80], tab_i_35_xmm, rounder_5
266 :     movq mm4, [tab_i_35_xmm+8] ; 4 ; w07 w06 w03 w02
267 :     movq mm6, [tab_i_35_xmm+32] ; 6 ; w21 w20 w17 w16
268 :     pshufw mm7, mm7, 10110001b ; y7 y6 y5 y4
269 :     paddd mm3, [rounder_5] ; +rounder stall 6
270 :     paddd mm5, [rounder_5] ; +rounder
271 :     movq [eax+64+8], mm7 ; 7 ; save y7 y6 y5 y4
272 :     movq mm7, [tab_i_35_xmm+40] ; 7 ; w23 w22 w19 w18
273 :     pmaddwd mm4, mm1 ; x6*w07+x4*w06 x6*w03+x4*w02
274 :     pmaddwd mm1, [tab_i_35_xmm+24]; x6*w15+x4*w14 x6*w11+x4*w10
275 :     pmaddwd mm6, mm2 ; x3*w21+x1*w20 x3*w17+x1*w16
276 :     pmaddwd mm2, [tab_i_35_xmm+48]; x3*w29+x1*w28 x3*w25+x1*w24
277 :     pmaddwd mm7, mm0 ; 7 ; x7*w23+x5*w22 x7*w19+x5*w18 ; w23 w22 w19 w18
278 :     pmaddwd mm0, [tab_i_35_xmm+56]; x7*w31+x5*w30 x7*w27+x5*w26
279 :     paddd mm3, mm4 ; 4 free ; a1=sum(even1) a0=sum(even0)
280 :     paddd mm5, mm1 ; 1 free ; a3=sum(even3) a2=sum(even2)
281 :     pshufw mm1, [eax+96+8],10001000b ; x6 x4 x6 x4
282 :     movq mm4, mm3 ; 4 ; a1 a0
283 :     paddd mm6, mm7 ; 7 free ; b1=sum(odd1) b0=sum(odd0)
284 :     paddd mm2, mm0 ; 5 free ; b3=sum(odd3) b2=sum(odd2)
285 :     pshufw mm0, [eax+96],10001000b ; x2 x0 x2 x0
286 :     movq mm7, mm5 ; 7 ; a3 a2
287 :     psubd mm4, mm6 ; 6 free ; a1-b1 a0-b0 stall 5
288 :     paddd mm6, mm3 ; mm3 = mm3+mm6+mm5+mm4; a1+b1 a0+b0
289 :     movq mm3, [tab_i_26_xmm] ; 3 ; w05 w04 w01 w00
290 :     psubd mm7, mm2 ; ; a3-b3 a2-b2
291 :     paddd mm5, mm2 ; 0 free a3+b3 a2+b2
292 :     pshufw mm2, [eax+96],11011101b; x3 x1 x3 x1
293 :     pmaddwd mm3, mm0 ; x2*w05+x0*w04 x2*w01+x0*w00
294 :     pmaddwd mm0, [tab_i_26_xmm+16]; x2*w13+x0*w12 x2*w09+x0*w08
295 :     psrad mm4, SHIFT_INV_ROW ; y6=a1-b1 y7=a0-b0
296 :     psrad mm7, SHIFT_INV_ROW ; y4=a3-b3 y5=a2-b2
297 :     psrad mm6, SHIFT_INV_ROW ; y1=a1+b1 y0=a0+b0
298 :     psrad mm5, SHIFT_INV_ROW ; y3=a3+b3 y2=a2+b2
299 :     packssdw mm7, mm4 ; 4 ; y6 y7 y4 y5
300 :     packssdw mm6, mm5 ; 0 free ; y3 y2 y1 y0
301 :     pshufw mm5, [eax+96+8],11011101b ; x7 x5 x7 x5
302 :     movq [eax+80], mm6 ; 3 ; save y3 y2 y1 y0
303 : edgomez 851
304 : edgomez 1192 ; DCT_8_INV_ROW_1_s [eax+96], [eax+96], tab_i_26_xmm, rounder_6
305 :     movq mm4, [tab_i_26_xmm+8] ; 4 ; w07 w06 w03 w02
306 :     movq mm6, [tab_i_26_xmm+32] ; 6 ; w21 w20 w17 w16
307 :     pshufw mm7, mm7, 10110001b ; y7 y6 y5 y4 STALL 6
308 :     paddd mm3, [rounder_6] ; +rounder
309 :     paddd mm0, [rounder_6] ; +rounder
310 :     movq [eax+80+8], mm7 ; 7 ; save y7 y6
311 :     movq mm7, [tab_i_26_xmm+40] ; 7 ; w23 w22 w19 w18
312 :     pmaddwd mm4, mm1 ; x6*w07+x4*w06 x6*w03+x4*w02
313 :     pmaddwd mm1, [tab_i_26_xmm+24]; x6*w15+x4*w14 x6*w11+x4*w10
314 :     pmaddwd mm6, mm2 ; x3*w21+x1*w20 x3*w17+x1*w16
315 :     pmaddwd mm2, [tab_i_26_xmm+48]; x3*w29+x1*w28 x3*w25+x1*w24
316 :     pmaddwd mm7, mm5 ; 7 ; x7*w23+x5*w22 x7*w19+x5*w18 ; w23 w22 w19 w18
317 :     pmaddwd mm5, [tab_i_26_xmm+56]; x7*w31+x5*w30 x7*w27+x5*w26
318 :     paddd mm3, mm4 ; 4 free ; a1=sum(even1) a0=sum(even0)
319 :     paddd mm0, mm1 ; 1 free ; a3=sum(even3) a2=sum(even2)
320 :     pshufw mm1, [eax+112+8],10001000b ; x6 x4 x6 x4
321 :     movq mm4, mm3 ; 4 ; a1 a0
322 :     paddd mm6, mm7 ; 7 free ; b1=sum(odd1) b0=sum(odd0)
323 :     paddd mm2, mm5 ; 5 free ; b3=sum(odd3) b2=sum(odd2)
324 :     pshufw mm5, [eax+112],10001000b; x2 x0 x2 x0 mm5 & mm0 exchanged for next cycle
325 :     movq mm7, mm0 ; 7 ; a3 a2
326 :     psubd mm4, mm6 ; 6 free ; a1-b1 a0-b0
327 :     paddd mm6, mm3 ; mm6 = mm3+mm6+mm5+mm4; a1+b1 a0+b0
328 :     movq mm3, [tab_i_17_xmm] ; 3 ; w05 w04 w01 w00
329 :     psubd mm7, mm2 ; ; a3-b3 a2-b2
330 :     paddd mm0, mm2 ; 0 free a3+b3 a2+b2
331 :     pshufw mm2, [eax+112],11011101b; x3 x1 x3 x1
332 :     pmaddwd mm3, mm5 ; x2*w05+x0*w04 x2*w01+x0*w00
333 :     pmaddwd mm5, [tab_i_17_xmm+16]; x2*w13+x0*w12 x2*w09+x0*w08
334 :     psrad mm4, SHIFT_INV_ROW ; y6=a1-b1 y7=a0-b0
335 :     psrad mm7, SHIFT_INV_ROW ; y4=a3-b3 y5=a2-b2
336 :     psrad mm6, SHIFT_INV_ROW ; y1=a1+b1 y0=a0+b0
337 :     psrad mm0, SHIFT_INV_ROW ; y3=a3+b3 y2=a2+b2
338 :     packssdw mm7, mm4 ; 4 ; y6 y7 y4 y5
339 :     packssdw mm6, mm0 ; 0 free ; y3 y2 y1 y0
340 :     pshufw mm0, [eax+112+8],11011101b ; x7 x5 x7 x5
341 :     movq [eax+96], mm6 ; 3 ; save y3 y2 y1 y0 stall2
342 : edgomez 851
343 : edgomez 1192 ; DCT_8_INV_ROW_1_s [eax+112], [eax+112], tab_i_17_xmm, rounder_7
344 :     movq mm4, [tab_i_17_xmm+8] ; 4 ; w07 w06 w03 w02
345 :     movq mm6, [tab_i_17_xmm+32] ; 6 ; w21 w20 w17 w16
346 :     pshufw mm7, mm7, 10110001b ; y7 y6 y5 y4
347 :     paddd mm3, [rounder_7] ; +rounder stall 6
348 :     paddd mm5, [rounder_7] ; +rounder
349 :     movq [eax+96+8], mm7 ; 7 ; save y7 y6 y5 y4
350 :     movq mm7, [tab_i_17_xmm+40] ; 7 ; w23 w22 w19 w18
351 :     pmaddwd mm4, mm1 ; x6*w07+x4*w06 x6*w03+x4*w02
352 :     pmaddwd mm1, [tab_i_17_xmm+24]; x6*w15+x4*w14 x6*w11+x4*w10
353 :     pmaddwd mm6, mm2 ; x3*w21+x1*w20 x3*w17+x1*w16
354 :     pmaddwd mm2, [tab_i_17_xmm+48]; x3*w29+x1*w28 x3*w25+x1*w24
355 :     pmaddwd mm7, mm0 ; 7 ; x7*w23+x5*w22 x7*w19+x5*w18 ; w23 w22 w19 w18
356 :     pmaddwd mm0, [tab_i_17_xmm+56]; x7*w31+x5*w30 x7*w27+x5*w26
357 :     paddd mm3, mm4 ; 4 free ; a1=sum(even1) a0=sum(even0)
358 :     paddd mm5, mm1 ; 1 free ; a3=sum(even3) a2=sum(even2)
359 :     pshufw mm1, [eax+0+8],10001000b; x6 x4 x6 x4
360 :     movq mm4, mm3 ; 4 ; a1 a0
361 :     paddd mm6, mm7 ; 7 free ; b1=sum(odd1) b0=sum(odd0)
362 :     paddd mm2, mm0 ; 5 free ; b3=sum(odd3) b2=sum(odd2)
363 :     pshufw mm0, [eax+0],10001000b ; x2 x0 x2 x0
364 :     movq mm7, mm5 ; 7 ; a3 a2
365 :     psubd mm4, mm6 ; 6 free ; a1-b1 a0-b0 stall 5
366 :     paddd mm6, mm3 ; mm3 = mm3+mm6+mm5+mm4; a1+b1 a0+b0
367 :     movq mm3, [tab_i_04_xmm] ; 3 ; w05 w04 w01 w00
368 :     psubd mm7, mm2 ; ; a3-b3 a2-b2
369 :     paddd mm5, mm2 ; 0 free a3+b3 a2+b2
370 :     pshufw mm2, [eax+0],11011101b ; x3 x1 x3 x1
371 :     pmaddwd mm3, mm0 ; x2*w05+x0*w04 x2*w01+x0*w00
372 :     pmaddwd mm0, [tab_i_04_xmm+16]; x2*w13+x0*w12 x2*w09+x0*w08
373 :     psrad mm4, SHIFT_INV_ROW ; y6=a1-b1 y7=a0-b0
374 :     psrad mm7, SHIFT_INV_ROW ; y4=a3-b3 y5=a2-b2
375 :     psrad mm6, SHIFT_INV_ROW ; y1=a1+b1 y0=a0+b0
376 :     psrad mm5, SHIFT_INV_ROW ; y3=a3+b3 y2=a2+b2
377 :     packssdw mm7, mm4 ; 4 ; y6 y7 y4 y5
378 :     packssdw mm6, mm5 ; 0 free ; y3 y2 y1 y0
379 :     pshufw mm5, [eax+0+8],11011101b; x7 x5 x7 x5
380 :     movq [eax+112], mm6 ; 3 ; save y3 y2 y1 y0
381 : edgomez 851
382 : edgomez 1192 ; DCT_8_INV_ROW_1_s [eax+0], 0, tab_i_04_xmm, rounder_0
383 :     movq mm4, [tab_i_04_xmm+8] ; 4 ; w07 w06 w03 w02
384 :     movq mm6, [tab_i_04_xmm+32] ; 6 ; w21 w20 w17 w16
385 :     pshufw mm7, mm7, 10110001b ; y7 y6 y5 y4 STALL 6
386 :     paddd mm3, [rounder_0] ; +rounder
387 :     paddd mm0, [rounder_0] ; +rounder
388 :     movq [eax+112+8], mm7 ; 7 ; save y7 y6
389 :     movq mm7, [tab_i_04_xmm+40] ; 7 ; w23 w22 w19 w18
390 :     pmaddwd mm4, mm1 ; x6*w07+x4*w06 x6*w03+x4*w02
391 :     pmaddwd mm1, [tab_i_04_xmm+24]; x6*w15+x4*w14 x6*w11+x4*w10
392 :     pmaddwd mm6, mm2 ; x3*w21+x1*w20 x3*w17+x1*w16
393 :     pmaddwd mm2, [tab_i_04_xmm+48]; x3*w29+x1*w28 x3*w25+x1*w24
394 :     pmaddwd mm7, mm5 ; 7 ; x7*w23+x5*w22 x7*w19+x5*w18 ; w23 w22 w19 w18
395 :     pmaddwd mm5, [tab_i_04_xmm+56]; x7*w31+x5*w30 x7*w27+x5*w26
396 :     paddd mm3, mm4 ; 4 free ; a1=sum(even1) a0=sum(even0)
397 :     paddd mm0, mm1 ; 1
398 :     pshufw mm1, [eax+16+8],10001000b ; x6 x4 x6 x4
399 :     movq mm4, mm3 ; 4 ; a1 a0
400 :     paddd mm6, mm7 ; 7 free ; b1=sum(odd1) b0=sum(odd0)
401 :     paddd mm2, mm5 ; 5 free ; b3=sum(odd3) b2=sum(odd2)
402 :     pshufw mm5, [eax+16],10001000b; x2 x0 x2 x0 mm5 & mm0 exchanged for next cycle
403 :     movq mm7, mm0 ; 7 ; a3 a2
404 :     psubd mm4, mm6 ; 6 free ; a1-b1 a0-b0
405 :     paddd mm6, mm3 ; mm6 = mm3+mm6+mm5+mm4; a1+b1 a0+b0
406 :     movq mm3, [tab_i_17_xmm] ; 3 ; w05 w04 w01 w00
407 :     psubd mm7, mm2 ; ; a3-b3 a2-b2
408 :     paddd mm0, mm2 ; 0 free a3+b3 a2+b2
409 :     pshufw mm2, [eax+16],11011101b; x3 x1 x3 x1
410 :     pmaddwd mm3, mm5 ; x2*w05+x0*w04 x2*w01+x0*w00
411 :     pmaddwd mm5, [tab_i_17_xmm+16]; x2*w13+x0*w12 x2*w09+x0*w08
412 :     psrad mm4, SHIFT_INV_ROW ; y6=a1-b1 y7=a0-b0
413 :     psrad mm7, SHIFT_INV_ROW ; y4=a3-b3 y5=a2-b2
414 :     psrad mm6, SHIFT_INV_ROW ; y1=a1+b1 y0=a0+b0
415 :     psrad mm0, SHIFT_INV_ROW ; y3=a3+b3 y2=a2+b2
416 :     packssdw mm7, mm4 ; 4 ; y6 y7 y4 y5
417 :     packssdw mm6, mm0 ; 0 free ; y3 y2 y1 y0
418 :     pshufw mm0, [eax+16+8],11011101b ; x7 x5 x7 x5
419 :     movq [eax+0], mm6 ; 3 ; save y3 y2 y1 y0 stall2
420 : edgomez 851
421 : edgomez 1192 ; DCT_8_INV_ROW_1_s [eax+16], 16, tab_i_17_xmm, rounder_1
422 :     movq mm4, [tab_i_17_xmm+8] ; 4 ; w07 w06 w03 w02
423 :     movq mm6, [tab_i_17_xmm+32] ; 6 ; w21 w20 w17 w16
424 :     pshufw mm7, mm7, 10110001b ; y7 y6 y5 y4
425 :     paddd mm3, [rounder_1] ; +rounder stall 6
426 :     paddd mm5, [rounder_1] ; +rounder
427 :     movq [eax+0+8], mm7 ; 7 ; save y7 y6 y5 y4
428 :     movq mm7, [tab_i_17_xmm+40] ; 7 ; w23 w22 w19 w18
429 :     pmaddwd mm4, mm1 ; x6*w07+x4*w06 x6*w03+x4*w02
430 :     pmaddwd mm1, [tab_i_17_xmm+24]; x6*w15+x4*w14 x6*w11+x4*w10
431 :     pmaddwd mm6, mm2 ; x3*w21+x1*w20 x3*w17+x1*w16
432 :     pmaddwd mm2, [tab_i_17_xmm+48]; x3*w29+x1*w28 x3*w25+x1*w24
433 :     pmaddwd mm7, mm0 ; 7 ; x7*w23+x5*w22 x7*w19+x5*w18 ; w23 w22 w19 w18
434 :     pmaddwd mm0, [tab_i_17_xmm+56]; x7*w31+x5*w30 x7*w27+x5*w26
435 :     paddd mm3, mm4 ; 4 free ; a1=sum(even1) a0=sum(even0)
436 :     paddd mm5, mm1 ; 1 free ; a3=sum(even3) a2=sum(even2)
437 :     pshufw mm1, [eax+32+8],10001000b ; x6 x4 x6 x4
438 :     movq mm4, mm3 ; 4 ; a1 a0
439 :     paddd mm6, mm7 ; 7 free ; b1=sum(odd1) b0=sum(odd0)
440 :     paddd mm2, mm0 ; 5 free ; b3=sum(odd3) b2=sum(odd2)
441 :     pshufw mm0, [eax+32],10001000b; x2 x0 x2 x0
442 :     movq mm7, mm5 ; 7 ; a3 a2
443 :     psubd mm4, mm6 ; 6 free ; a1-b1 a0-b0 stall 5
444 :     paddd mm6, mm3 ; mm3 = mm3+mm6+mm5+mm4; a1+b1 a0+b0
445 :     movq mm3, [tab_i_26_xmm] ; 3 ; w05 w04 w01 w00
446 :     psubd mm7, mm2 ; ; a3-b3 a2-b2
447 :     paddd mm5, mm2 ; 0 free a3+b3 a2+b2
448 :     pshufw mm2, [eax+32],11011101b; x3 x1 x3 x1
449 :     pmaddwd mm3, mm0 ; x2*w05+x0*w04 x2*w01+x0*w00
450 :     pmaddwd mm0, [tab_i_26_xmm+16]; x2*w13+x0*w12 x2*w09+x0*w08
451 :     psrad mm4, SHIFT_INV_ROW ; y6=a1-b1 y7=a0-b0
452 :     psrad mm7, SHIFT_INV_ROW ; y4=a3-b3 y5=a2-b2
453 :     psrad mm6, SHIFT_INV_ROW ; y1=a1+b1 y0=a0+b0
454 :     psrad mm5, SHIFT_INV_ROW ; y3=a3+b3 y2=a2+b2
455 :     packssdw mm7, mm4 ; 4 ; y6 y7 y4 y5
456 :     packssdw mm6, mm5 ; 0 free ; y3 y2 y1 y0
457 :     pshufw mm5, [eax+32+8],11011101b ; x7 x5 x7 x5
458 :     movq [eax+16], mm6 ; 3 ; save y3 y2 y1 y0
459 : edgomez 851
460 : edgomez 1192 ; DCT_8_INV_ROW_1_s [eax+32], 32, tab_i_26_xmm, rounder_2
461 :     movq mm4, [tab_i_26_xmm+8] ; 4 ; w07 w06 w03 w02
462 :     movq mm6, [tab_i_26_xmm+32] ; 6 ; w21 w20 w17 w16
463 :     pshufw mm7, mm7, 10110001b ; y7 y6 y5 y4 STALL 6
464 :     paddd mm3, [rounder_2] ; +rounder
465 :     paddd mm0, [rounder_2] ; +rounder
466 :     movq [eax+16+8], mm7 ; 7 ; save y7 y6
467 :     movq mm7, [tab_i_26_xmm+40] ; 7 ; w23 w22 w19 w18
468 :     pmaddwd mm4, mm1 ; x6*w07+x4*w06 x6*w03+x4*w02
469 :     pmaddwd mm1, [tab_i_26_xmm+24]; x6*w15+x4*w14 x6*w11+x4*w10
470 :     pmaddwd mm6, mm2 ; x3*w21+x1*w20 x3*w17+x1*w16
471 :     pmaddwd mm2, [tab_i_26_xmm+48]; x3*w29+x1*w28 x3*w25+x1*w24
472 :     pmaddwd mm7, mm5 ; 7 ; x7*w23+x5*w22 x7*w19+x5*w18 ; w23 w22 w19 w18
473 :     pmaddwd mm5, [tab_i_26_xmm+56]; x7*w31+x5*w30 x7*w27+x5*w26
474 :     paddd mm3, mm4 ; 4 free ; a1=sum(even1) a0=sum(even0)
475 :     paddd mm0, mm1 ; 1 free ; a3=sum(even3) a2=sum(even2)
476 :     pshufw mm1, [eax+48+8],10001000b ; x6 x4 x6 x4
477 :     movq mm4, mm3 ; 4 ; a1 a0
478 :     paddd mm6, mm7 ; 7 free ; b1=sum(odd1) b0=sum(odd0)
479 :     paddd mm2, mm5 ; 5 free ; b3=sum(odd3) b2=sum(odd2)
480 :     pshufw mm5, [eax+48],10001000b; x2 x0 x2 x0 mm5 & mm0 exchanged for next cycle
481 :     movq mm7, mm0 ; 7 ; a3 a2
482 :     psubd mm4, mm6 ; 6 free ; a1-b1 a0-b0
483 :     paddd mm6, mm3 ; mm6 = mm3+mm6+mm5+mm4; a1+b1 a0+b0
484 :     movq mm3, [tab_i_35_xmm] ; 3 ; w05 w04 w01 w00
485 :     psubd mm7, mm2 ; ; a3-b3 a2-b2
486 :     paddd mm0, mm2 ; 0 free a3+b3 a2+b2
487 :     pshufw mm2, [eax+48],11011101b; x3 x1 x3 x1
488 :     pmaddwd mm3, mm5 ; x2*w05+x0*w04 x2*w01+x0*w00
489 :     pmaddwd mm5, [tab_i_35_xmm+16]; x2*w13+x0*w12 x2*w09+x0*w08
490 :     psrad mm4, SHIFT_INV_ROW ; y6=a1-b1 y7=a0-b0
491 :     psrad mm7, SHIFT_INV_ROW ; y4=a3-b3 y5=a2-b2
492 :     psrad mm6, SHIFT_INV_ROW ; y1=a1+b1 y0=a0+b0
493 :     psrad mm0, SHIFT_INV_ROW ; y3=a3+b3 y2=a2+b2
494 :     packssdw mm7, mm4 ; 4 ; y6 y7 y4 y5
495 :     packssdw mm6, mm0 ; 0 free ; y3 y2 y1 y0
496 :     pshufw mm0, [eax+48+8],11011101b ; x7 x5 x7 x5
497 :     movq [eax+32], mm6 ; 3 ; save y3 y2 y1 y0 stall2
498 : edgomez 851
499 : edgomez 1192 ; DCT_8_INV_ROW_1_s [eax+48], [eax+48], tab_i_35_xmm, rounder_3
500 :     movq mm4, [tab_i_35_xmm+8] ; 4 ; w07 w06 w03 w02
501 :     movq mm6, [tab_i_35_xmm+32] ; 6 ; w21 w20 w17 w16
502 :     pshufw mm7, mm7, 10110001b ; y7 y6 y5 y4
503 :     paddd mm3, [rounder_3] ; +rounder stall 6
504 :     paddd mm5, [rounder_3] ; +rounder
505 :     movq [eax+32+8], mm7 ; 7 ; save y7 y6 y5 y4
506 :     movq mm7, [tab_i_35_xmm+40] ; 7 ; w23 w22 w19 w18
507 :     pmaddwd mm4, mm1 ; x6*w07+x4*w06 x6*w03+x4*w02
508 :     pmaddwd mm1, [tab_i_35_xmm+24]; x6*w15+x4*w14 x6*w11+x4*w10
509 :     pmaddwd mm6, mm2 ; x3*w21+x1*w20 x3*w17+x1*w16
510 :     pmaddwd mm2, [tab_i_35_xmm+48]; x3*w29+x1*w28 x3*w25+x1*w24
511 :     pmaddwd mm7, mm0 ; 7 ; x7*w23+x5*w22 x7*w19+x5*w18 ; w23 w22 w19 w18
512 :     pmaddwd mm0, [tab_i_35_xmm+56]; x7*w31+x5*w30 x7*w27+x5*w26
513 :     paddd mm3, mm4 ; 4 free ; a1=sum(even1) a0=sum(even0)
514 :     paddd mm5, mm1 ; mm1 free ; a3=sum(even3) a2=sum(even2)
515 :     movq mm1, [tg_3_16]
516 :     movq mm4, mm3 ; 4 ; a1 a0
517 :     paddd mm6, mm7 ; 7 free ; b1=sum(odd1) b0=sum(odd0)
518 :     paddd mm2, mm0 ; 5 free ; b3=sum(odd3) b2=sum(odd2)
519 :     movq mm0, [tg_3_16]
520 :     movq mm7, mm5 ; 7 ; a3 a2
521 :     psubd mm4, mm6 ; 6 free ; a1-b1 a0-b0
522 :     paddd mm3, mm6 ; mm3 = mm3+mm6+mm5+mm4; a1+b1 a0+b0
523 :     psubd mm7, mm2 ; ; a3-b3 a2-b2
524 :     paddd mm2, mm5 ; 0 free a3+b3 a2+b2
525 :     movq mm5, [eax+16*5]
526 :     psrad mm4, SHIFT_INV_ROW ; y6=a1-b1 y7=a0-b0
527 :     psrad mm7, SHIFT_INV_ROW ; y4=a3-b3 y5=a2-b2
528 :     psrad mm3, SHIFT_INV_ROW ; y1=a1+b1 y0=a0+b0
529 :     psrad mm2, SHIFT_INV_ROW ; y3=a3+b3 y2=a2+b2
530 :     movq mm6, [eax+16*1]
531 :     packssdw mm7, mm4 ; 4 ; y6 y7 y4 y5
532 :     movq mm4, [tg_1_16]
533 :     packssdw mm3, mm2 ; 0 free ; y3 y2 y1 y0
534 :     pshufw mm2, mm7, 10110001b ; y7 y6 y5 y4
535 : edgomez 851
536 : edgomez 1192 ; DCT_8_INV_COL_4 [eax+0],[eax+0]
537 :     ; movq mm3,mmword ptr [eax+16*3]
538 :     movq mm7, [eax+16*7]
539 :     pmulhw mm0, mm3 ; x3*(tg_3_16-1)
540 :     pmulhw mm1, mm5 ; x5*(tg_3_16-1)
541 :     movq [eax+48+8], mm2 ; 7 ; save y7 y6 y5 y4
542 :     movq mm2, mm4 ; tg_1_16
543 :     pmulhw mm4, mm7 ; x7*tg_1_16
544 :     paddsw mm0, mm3 ; x3*tg_3_16
545 :     pmulhw mm2, mm6 ; x1*tg_1_16
546 :     paddsw mm1, mm3 ; x3+x5*(tg_3_16-1)
547 :     psubsw mm0, mm5 ; x3*tg_3_16-x5 = tm35
548 :     movq [eax+48], mm3 ; 3 ; save y3 y2 y1 y0
549 :     movq mm3, [ocos_4_16]
550 :     paddsw mm1, mm5 ; x3+x5*tg_3_16 = tp35
551 :     paddsw mm4, mm6 ; x1+tg_1_16*x7 = tp17
552 :     psubsw mm2, mm7 ; x1*tg_1_16-x7 = tm17
553 :     movq mm5, mm4 ; tp17
554 :     movq mm6, mm2 ; tm17
555 :     paddsw mm5, mm1 ; tp17+tp35 = b0
556 :     psubsw mm6, mm0 ; tm17-tm35 = b3
557 :     psubsw mm4, mm1 ; tp17-tp35 = t1
558 :     paddsw mm2, mm0 ; tm17+tm35 = t2
559 :     movq mm7, [tg_2_16]
560 :     movq mm1, mm4 ; t1
561 :     movq [eax+3*16], mm5 ; save b0
562 :     paddsw mm1, mm2 ; t1+t2
563 :     movq [eax+5*16], mm6 ; save b3
564 :     psubsw mm4, mm2 ; t1-t2
565 :     movq mm5, [eax+2*16]
566 :     movq mm0, mm7 ; tg_2_16
567 :     movq mm6, [eax+6*16]
568 :     pmulhw mm0, mm5 ; x2*tg_2_16
569 :     pmulhw mm7, mm6 ; x6*tg_2_16
570 : edgomez 851 ; slot
571 : edgomez 1192 pmulhw mm1, mm3 ; ocos_4_16*(t1+t2) = b1/2
572 : edgomez 851 ; slot
573 : edgomez 1192 movq mm2, [eax+0*16]
574 :     pmulhw mm4, mm3 ; ocos_4_16*(t1-t2) = b2/2
575 :     psubsw mm0, mm6 ; t2*tg_2_16-x6 = tm26
576 :     movq mm3, [eax+0*16] ; x0
577 :     movq mm6, [eax+4*16]
578 :     paddsw mm7, mm5 ; x2+x6*tg_2_16 = tp26
579 :     paddsw mm2, mm6 ; x0+x4 = tp04
580 :     psubsw mm3, mm6 ; x0-x4 = tm04
581 :     movq mm5, mm2 ; tp04
582 :     movq mm6, mm3 ; tm04
583 :     psubsw mm2, mm7 ; tp04-tp26 = a3
584 :     paddsw mm3, mm0 ; tm04+tm26 = a1
585 :     paddsw mm1, mm1 ; b1
586 :     paddsw mm4, mm4 ; b2
587 :     paddsw mm5, mm7 ; tp04+tp26 = a0
588 :     psubsw mm6, mm0 ; tm04-tm26 = a2
589 :     movq mm7, mm3 ; a1
590 :     movq mm0, mm6 ; a2
591 :     paddsw mm3, mm1 ; a1+b1
592 :     paddsw mm6, mm4 ; a2+b2
593 :     psraw mm3, SHIFT_INV_COL ; dst1
594 :     psubsw mm7, mm1 ; a1-b1
595 :     psraw mm6, SHIFT_INV_COL ; dst2
596 :     psubsw mm0, mm4 ; a2-b2
597 :     movq mm1, [eax+3*16] ; load b0
598 :     psraw mm7, SHIFT_INV_COL ; dst6
599 :     movq mm4, mm5 ; a0
600 :     psraw mm0, SHIFT_INV_COL ; dst5
601 :     movq [eax+1*16], mm3
602 :     paddsw mm5, mm1 ; a0+b0
603 :     movq [eax+2*16], mm6
604 :     psubsw mm4, mm1 ; a0-b0
605 :     movq mm3, [eax+5*16] ; load b3
606 :     psraw mm5, SHIFT_INV_COL ; dst0
607 :     movq mm6, mm2 ; a3
608 :     psraw mm4, SHIFT_INV_COL ; dst7
609 :     movq [eax+5*16], mm0
610 :     movq mm0, [tg_3_16]
611 :     paddsw mm2, mm3 ; a3+b3
612 :     movq [eax+6*16], mm7
613 :     psubsw mm6, mm3 ; a3-b3
614 :     movq mm3, [eax+8+16*3]
615 :     movq [eax+0*16], mm5
616 :     psraw mm2, SHIFT_INV_COL ; dst3
617 :     movq [eax+7*16], mm4
618 : edgomez 851
619 : edgomez 1192 ; DCT_8_INV_COL_4 [eax+8],[eax+8]
620 :     movq mm1, mm0 ; tg_3_16
621 :     movq mm5, [eax+8+16*5]
622 :     psraw mm6, SHIFT_INV_COL ; dst4
623 :     pmulhw mm0, mm3 ; x3*(tg_3_16-1)
624 :     movq mm4, [tg_1_16]
625 :     pmulhw mm1, mm5 ; x5*(tg_3_16-1)
626 :     movq mm7, [eax+8+16*7]
627 :     movq [eax+3*16], mm2
628 :     movq mm2, mm4 ; tg_1_16
629 :     movq [eax+4*16], mm6
630 :     movq mm6, [eax+8+16*1]
631 :     pmulhw mm4, mm7 ; x7*tg_1_16
632 :     paddsw mm0, mm3 ; x3*tg_3_16
633 :     pmulhw mm2, mm6 ; x1*tg_1_16
634 :     paddsw mm1, mm3 ; x3+x5*(tg_3_16-1)
635 :     psubsw mm0, mm5 ; x3*tg_3_16-x5 = tm35
636 :     movq mm3, [ocos_4_16]
637 :     paddsw mm1, mm5 ; x3+x5*tg_3_16 = tp35
638 :     paddsw mm4, mm6 ; x1+tg_1_16*x7 = tp17
639 :     psubsw mm2, mm7 ; x1*tg_1_16-x7 = tm17
640 :     movq mm5, mm4 ; tp17
641 :     movq mm6, mm2 ; tm17
642 :     paddsw mm5, mm1 ; tp17+tp35 = b0
643 :     psubsw mm4, mm1 ; tp17-tp35 = t1
644 :     paddsw mm2, mm0 ; tm17+tm35 = t2
645 :     movq mm7, [tg_2_16]
646 :     movq mm1, mm4 ; t1
647 :     psubsw mm6, mm0 ; tm17-tm35 = b3
648 :     movq [eax+8+3*16], mm5 ; save b0
649 :     movq [eax+8+5*16], mm6 ; save b3
650 :     psubsw mm4, mm2 ; t1-t2
651 :     movq mm5, [eax+8+2*16]
652 :     movq mm0, mm7 ; tg_2_16
653 :     movq mm6, [eax+8+6*16]
654 :     paddsw mm1, mm2 ; t1+t2
655 :     pmulhw mm0, mm5 ; x2*tg_2_16
656 :     pmulhw mm7, mm6 ; x6*tg_2_16
657 :     movq mm2, [eax+8+0*16]
658 :     pmulhw mm4, mm3 ; ocos_4_16*(t1-t2) = b2/2
659 :     psubsw mm0, mm6 ; t2*tg_2_16-x6 = tm26
660 : edgomez 851 ; slot
661 : edgomez 1192 pmulhw mm1, mm3 ; ocos_4_16*(t1+t2) = b1/2
662 : edgomez 851 ; slot
663 : edgomez 1192 movq mm3, [eax+8+0*16] ; x0
664 :     movq mm6, [eax+8+4*16]
665 :     paddsw mm7, mm5 ; x2+x6*tg_2_16 = tp26
666 :     paddsw mm2, mm6 ; x0+x4 = tp04
667 :     psubsw mm3, mm6 ; x0-x4 = tm04
668 :     movq mm5, mm2 ; tp04
669 :     movq mm6, mm3 ; tm04
670 :     psubsw mm2, mm7 ; tp04-tp26 = a3
671 :     paddsw mm3, mm0 ; tm04+tm26 = a1
672 :     paddsw mm1, mm1 ; b1
673 :     paddsw mm4, mm4 ; b2
674 :     paddsw mm5, mm7 ; tp04+tp26 = a0
675 :     psubsw mm6, mm0 ; tm04-tm26 = a2
676 :     movq mm7, mm3 ; a1
677 :     movq mm0, mm6 ; a2
678 :     paddsw mm3, mm1 ; a1+b1
679 :     paddsw mm6, mm4 ; a2+b2
680 :     psraw mm3, SHIFT_INV_COL ; dst1
681 :     psubsw mm7, mm1 ; a1-b1
682 :     psraw mm6, SHIFT_INV_COL ; dst2
683 :     psubsw mm0, mm4 ; a2-b2
684 :     movq mm1, [eax+8+3*16] ; load b0
685 :     psraw mm7, SHIFT_INV_COL ; dst6
686 :     movq mm4, mm5 ; a0
687 :     psraw mm0, SHIFT_INV_COL ; dst5
688 :     movq [eax+8+1*16], mm3
689 :     paddsw mm5, mm1 ; a0+b0
690 :     movq [eax+8+2*16], mm6
691 :     psubsw mm4, mm1 ; a0-b0
692 :     movq mm3, [eax+8+5*16] ; load b3
693 :     psraw mm5, SHIFT_INV_COL ; dst0
694 :     movq mm6, mm2 ; a3
695 :     psraw mm4, SHIFT_INV_COL ; dst7
696 :     movq [eax+8+5*16], mm0
697 :     paddsw mm2, mm3 ; a3+b3
698 :     movq [eax+8+6*16], mm7
699 :     psubsw mm6, mm3 ; a3-b3
700 :     movq [eax+8+0*16], mm5
701 :     psraw mm2, SHIFT_INV_COL ; dst3
702 :     movq [eax+8+7*16], mm4
703 :     psraw mm6, SHIFT_INV_COL ; dst4
704 :     movq [eax+8+3*16], mm2
705 :     movq [eax+8+4*16], mm6
706 : edgomez 851
707 : edgomez 1192 ret
No admin address has been configured
 ViewVC Help
Powered by ViewVC 1.0.4