[svn] / trunk / xvidcore / src / dct / x86_asm / idct_3dne.asm Repository:
ViewVC logotype

Annotation of /trunk/xvidcore/src/dct/x86_asm/idct_3dne.asm

Parent Directory Parent Directory | Revision Log Revision Log

Revision 1382 - (view) (download)

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