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

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

Parent Directory Parent Directory | Revision Log Revision Log


Revision 460 - (view) (download)

1 : chl 460 ;/*****************************************************************************
2 :     ; *
3 :     ; * XVID MPEG-4 VIDEO CODEC
4 :     ; * mmx version - inverse discrete cosine transformation
5 :     ; *
6 :     ; * Initial version provided by Intel at AppNote AP-922
7 :     ; * Copyright (C) 1999 Intel Corporation,
8 :     ; *
9 :     ; * Modifications
10 :     ; * Copyright (c) 2000-2001 Peter Gubanov <peter@elecard.net.ru>
11 :     ; * Copyright (c) 2000 Michel Lespinasse <walken@zoy.org>
12 :     ; *
13 :     ; * ported to NASM and some minor changes
14 :     ; * Copyright (C) 2001 Peter Ross <pross@xvid.org>
15 :     ; *
16 :     ; * This program is an implementation of a part of one or more MPEG-4
17 :     ; * Video tools as specified in ISO/IEC 14496-2 standard. Those intending
18 :     ; * to use this software module in hardware or software products are
19 :     ; * advised that its use may infringe existing patents or copyrights, and
20 :     ; * any such use would be at such party's own risk. The original
21 :     ; * developer of this software module and his/her company, and subsequent
22 :     ; * editors and their companies, will have no liability for use of this
23 :     ; * software or modifications or derivatives thereof.
24 :     ; *
25 :     ; * This program is free software; you can redistribute it and/or modify
26 :     ; * it under the terms of the GNU General Public License as published by
27 :     ; * the Free Software Foundation; either version 2 of the License, or
28 :     ; * (at your option) any later version.
29 :     ; *
30 :     ; * This program is distributed in the hope that it will be useful,
31 :     ; * but WITHOUT ANY WARRANTY; without even the implied warranty of
32 :     ; * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
33 :     ; * GNU General Public License for more details.
34 :     ; *
35 :     ; * You should have received a copy of the GNU General Public License
36 :     ; * along with this program; if not, write to the Free Software
37 :     ; * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
38 :     ; *
39 :     ; *************************************************************************/
40 :    
41 : Isibaar 3 ;=============================================================================
42 :     ;
43 :     ; These examples contain code fragments for first stage iDCT 8x8
44 :     ; (for rows) and first stage DCT 8x8 (for columns)
45 :     ;
46 :     ;=============================================================================
47 :    
48 :     bits 32
49 :    
50 :     %macro cglobal 1
51 :     %ifdef PREFIX
52 :     global _%1
53 :     %define %1 _%1
54 :     %else
55 :     global %1
56 :     %endif
57 :     %endmacro
58 :    
59 :     %define BITS_INV_ACC 5 ; 4 or 5 for IEEE
60 :     %define SHIFT_INV_ROW 16 - BITS_INV_ACC
61 :     %define SHIFT_INV_COL 1 + BITS_INV_ACC
62 :     %define RND_INV_ROW 1024 * (6 - BITS_INV_ACC) ; 1 << (SHIFT_INV_ROW-1)
63 :     %define RND_INV_COL 16 * (BITS_INV_ACC - 3) ; 1 << (SHIFT_INV_COL-1)
64 :     %define RND_INV_CORR RND_INV_COL - 1 ; correction -1.0 and round
65 :    
66 :     %define BITS_FRW_ACC 3 ; 2 or 3 for accuracy
67 :     %define SHIFT_FRW_COL BITS_FRW_ACC
68 :     %define SHIFT_FRW_ROW BITS_FRW_ACC + 17
69 :     %define RND_FRW_ROW 262144 * (BITS_FRW_ACC - 1) ; 1 << (SHIFT_FRW_ROW-1)
70 :    
71 : suxen_drol 207 %ifdef FORMAT_COFF
72 :     section .data data
73 :     %else
74 : Isibaar 154 section .data data align=16
75 : suxen_drol 207 %endif
76 : Isibaar 3
77 :     align 16
78 :    
79 :     one_corr dw 1, 1, 1, 1
80 :     round_inv_row dd RND_INV_ROW, RND_INV_ROW
81 :     round_inv_col dw RND_INV_COL, RND_INV_COL, RND_INV_COL, RND_INV_COL
82 :     round_inv_corr dw RND_INV_CORR, RND_INV_CORR, RND_INV_CORR, RND_INV_CORR
83 :     round_frw_row dd RND_FRW_ROW, RND_FRW_ROW
84 :     tg_1_16 dw 13036, 13036, 13036, 13036 ; tg * (2<<16) + 0.5
85 :     tg_2_16 dw 27146, 27146, 27146, 27146 ; tg * (2<<16) + 0.5
86 :     tg_3_16 dw -21746, -21746, -21746, -21746 ; tg * (2<<16) + 0.5
87 :     cos_4_16 dw -19195, -19195, -19195, -19195 ; cos * (2<<16) + 0.5
88 :     ocos_4_16 dw 23170, 23170, 23170, 23170 ; cos * (2<<15) + 0.5
89 :    
90 :     otg_3_16 dw 21895, 21895, 21895, 21895 ; tg * (2<<16) + 0.5
91 :    
92 :     %if SHIFT_INV_ROW == 12 ; assume SHIFT_INV_ROW == 12
93 :     rounder_0 dd 65536, 65536
94 :     rounder_4 dd 0, 0
95 :     rounder_1 dd 7195, 7195
96 :     rounder_7 dd 1024, 1024
97 :     rounder_2 dd 4520, 4520
98 :     rounder_6 dd 1024, 1024
99 :     rounder_3 dd 2407, 2407
100 :     rounder_5 dd 240, 240
101 :    
102 :     %elif SHIFT_INV_ROW == 11 ; assume SHIFT_INV_ROW == 11
103 :     rounder_0 dd 65536, 65536
104 :     rounder_4 dd 0, 0
105 :     rounder_1 dd 3597, 3597
106 :     rounder_7 dd 512, 512
107 :     rounder_2 dd 2260, 2260
108 :     rounder_6 dd 512, 512
109 :     rounder_3 dd 1203, 1203
110 :     rounder_5 dd 120, 120
111 :     %else
112 :    
113 :     %error invalid _SHIFT_INV_ROW_
114 :    
115 :     %endif
116 :    
117 :     ;=============================================================================
118 :     ;
119 :     ; The first stage iDCT 8x8 - inverse DCTs of rows
120 :     ;
121 :     ;-----------------------------------------------------------------------------
122 :     ; The 8-point inverse DCT direct algorithm
123 :     ;-----------------------------------------------------------------------------
124 :     ;
125 :     ; static const short w[32] = {
126 :     ; FIX(cos_4_16), FIX(cos_2_16), FIX(cos_4_16), FIX(cos_6_16),
127 :     ; FIX(cos_4_16), FIX(cos_6_16), -FIX(cos_4_16), -FIX(cos_2_16),
128 :     ; FIX(cos_4_16), -FIX(cos_6_16), -FIX(cos_4_16), FIX(cos_2_16),
129 :     ; FIX(cos_4_16), -FIX(cos_2_16), FIX(cos_4_16), -FIX(cos_6_16),
130 :     ; FIX(cos_1_16), FIX(cos_3_16), FIX(cos_5_16), FIX(cos_7_16),
131 :     ; FIX(cos_3_16), -FIX(cos_7_16), -FIX(cos_1_16), -FIX(cos_5_16),
132 :     ; FIX(cos_5_16), -FIX(cos_1_16), FIX(cos_7_16), FIX(cos_3_16),
133 :     ; FIX(cos_7_16), -FIX(cos_5_16), FIX(cos_3_16), -FIX(cos_1_16) };
134 :     ;
135 :     ; #define DCT_8_INV_ROW(x, y)
136 :     ; {
137 :     ; int a0, a1, a2, a3, b0, b1, b2, b3;
138 :     ;
139 :     ; a0 =x[0]*w[0]+x[2]*w[1]+x[4]*w[2]+x[6]*w[3];
140 :     ; a1 =x[0]*w[4]+x[2]*w[5]+x[4]*w[6]+x[6]*w[7];
141 :     ; a2 = x[0] * w[ 8] + x[2] * w[ 9] + x[4] * w[10] + x[6] * w[11];
142 :     ; a3 = x[0] * w[12] + x[2] * w[13] + x[4] * w[14] + x[6] * w[15];
143 :     ; b0 = x[1] * w[16] + x[3] * w[17] + x[5] * w[18] + x[7] * w[19];
144 :     ; b1 = x[1] * w[20] + x[3] * w[21] + x[5] * w[22] + x[7] * w[23];
145 :     ; b2 = x[1] * w[24] + x[3] * w[25] + x[5] * w[26] + x[7] * w[27];
146 :     ; b3 = x[1] * w[28] + x[3] * w[29] + x[5] * w[30] + x[7] * w[31];
147 :     ;
148 :     ; y[0] = SHIFT_ROUND ( a0 + b0 );
149 :     ; y[1] = SHIFT_ROUND ( a1 + b1 );
150 :     ; y[2] = SHIFT_ROUND ( a2 + b2 );
151 :     ; y[3] = SHIFT_ROUND ( a3 + b3 );
152 :     ; y[4] = SHIFT_ROUND ( a3 - b3 );
153 :     ; y[5] = SHIFT_ROUND ( a2 - b2 );
154 :     ; y[6] = SHIFT_ROUND ( a1 - b1 );
155 :     ; y[7] = SHIFT_ROUND ( a0 - b0 );
156 :     ; }
157 :     ;
158 :     ;-----------------------------------------------------------------------------
159 :     ;
160 :     ; In this implementation the outputs of the iDCT-1D are multiplied
161 :     ; for rows 0,4 - by cos_4_16,
162 :     ; for rows 1,7 - by cos_1_16,
163 :     ; for rows 2,6 - by cos_2_16,
164 :     ; for rows 3,5 - by cos_3_16
165 :     ; and are shifted to the left for better accuracy
166 :     ;
167 :     ; For the constants used,
168 :     ; FIX(float_const) = (short) (float_const * (1<<15) + 0.5)
169 :     ;
170 :     ;=============================================================================
171 :    
172 :     ;=============================================================================
173 :     ; MMX code
174 :     ;=============================================================================
175 :    
176 :     ; Table for rows 0,4 - constants are multiplied by cos_4_16
177 :    
178 :     tab_i_04 dw 16384, 16384, 16384, -16384 ; movq-> w06 w04 w02 w00
179 :     dw 21407, 8867, 8867, -21407 ; w07 w05 w03 w01
180 :     dw 16384, -16384, 16384, 16384 ; w14 w12 w10 w08
181 :     dw -8867, 21407, -21407, -8867 ; w15 w13 w11 w09
182 :     dw 22725, 12873, 19266, -22725 ; w22 w20 w18 w16
183 :     dw 19266, 4520, -4520, -12873 ; w23 w21 w19 w17
184 :     dw 12873, 4520, 4520, 19266 ; w30 w28 w26 w24
185 :     dw -22725, 19266, -12873, -22725 ; w31 w29 w27 w25
186 :    
187 :     ; Table for rows 1,7 - constants are multiplied by cos_1_16
188 :    
189 :     tab_i_17 dw 22725, 22725, 22725, -22725 ; movq-> w06 w04 w02 w00
190 :     dw 29692, 12299, 12299, -29692 ; w07 w05 w03 w01
191 :     dw 22725, -22725, 22725, 22725 ; w14 w12 w10 w08
192 :     dw -12299, 29692, -29692, -12299 ; w15 w13 w11 w09
193 :     dw 31521, 17855, 26722, -31521 ; w22 w20 w18 w16
194 :     dw 26722, 6270, -6270, -17855 ; w23 w21 w19 w17
195 :     dw 17855, 6270, 6270, 26722 ; w30 w28 w26 w24
196 :     dw -31521, 26722, -17855, -31521 ; w31 w29 w27 w25
197 :    
198 :     ; Table for rows 2,6 - constants are multiplied by cos_2_16
199 :    
200 :     tab_i_26 dw 21407, 21407, 21407, -21407 ; movq-> w06 w04 w02 w00
201 :     dw 27969, 11585, 11585, -27969 ; w07 w05 w03 w01
202 :     dw 21407, -21407, 21407, 21407 ; w14 w12 w10 w08
203 :     dw -11585, 27969, -27969, -11585 ; w15 w13 w11 w09
204 :     dw 29692, 16819, 25172, -29692 ; w22 w20 w18 w16
205 :     dw 25172, 5906, -5906, -16819 ; w23 w21 w19 w17
206 :     dw 16819, 5906, 5906, 25172 ; w30 w28 w26 w24
207 :     dw -29692, 25172, -16819, -29692 ; w31 w29 w27 w25
208 :    
209 :     ; Table for rows 3,5 - constants are multiplied by cos_3_16
210 :    
211 :     tab_i_35 dw 19266, 19266, 19266, -19266 ; movq-> w06 w04 w02 w00
212 :     dw 25172, 10426, 10426, -25172 ; w07 w05 w03 w01
213 :     dw 19266, -19266, 19266, 19266 ; w14 w12 w10 w08
214 :     dw -10426, 25172, -25172, -10426 ; w15 w13 w11 w09
215 :     dw 26722, 15137, 22654, -26722 ; w22 w20 w18 w16
216 :     dw 22654, 5315, -5315, -15137 ; w23 w21 w19 w17
217 :     dw 15137, 5315, 5315, 22654 ; w30 w28 w26 w24
218 :     dw -26722, 22654, -15137, -26722 ; w31 w29 w27 w25
219 :    
220 :     ;-----------------------------------------------------------------------------
221 :    
222 :     ;
223 :     ; DCT_8_INV_ROW_1 INP, OUT, TABLE, ROUNDER
224 :     ;
225 :    
226 :     %macro DCT_8_INV_ROW_1 4
227 :    
228 :     movq mm0, [%1] ; 0 ; x3 x2 x1 x0
229 :    
230 :     movq mm1, [%1+8] ; 1 ; x7 x6 x5 x4
231 :     movq mm2, mm0 ; 2 ; x3 x2 x1 x0
232 :    
233 :     movq mm3, [%3] ; 3 ; w06 w04 w02 w00
234 :     punpcklwd mm0, mm1 ; x5 x1 x4 x0
235 :    
236 :     movq mm5, mm0 ; 5 ; x5 x1 x4 x0
237 :     punpckldq mm0, mm0 ; x4 x0 x4 x0
238 :    
239 :     movq mm4, [%3+8] ; 4 ; w07 w05 w03 w01
240 :     punpckhwd mm2, mm1 ; 1 ; x7 x3 x6 x2
241 :    
242 :     pmaddwd mm3, mm0 ; x4*w06+x0*w04 x4*w02+x0*w00
243 :     movq mm6, mm2 ; 6 ; x7 x3 x6 x2
244 :    
245 :     movq mm1, [%3+32] ; 1 ; w22 w20 w18 w16
246 :     punpckldq mm2, mm2 ; x6 x2 x6 x2
247 :    
248 :     pmaddwd mm4, mm2 ; x6*w07+x2*w05 x6*w03+x2*w01
249 :     punpckhdq mm5, mm5 ; x5 x1 x5 x1
250 :    
251 :     pmaddwd mm0, [%3+16] ; x4*w14+x0*w12 x4*w10+x0*w08
252 :     punpckhdq mm6, mm6 ; x7 x3 x7 x3
253 :    
254 :     movq mm7, [%3+40] ; 7 ; w23 w21 w19 w17
255 :     pmaddwd mm1, mm5 ; x5*w22+x1*w20 x5*w18+x1*w16
256 :    
257 :     paddd mm3, [%4] ; +%4
258 :     pmaddwd mm7, mm6 ; x7*w23+x3*w21 x7*w19+x3*w17
259 :    
260 :     pmaddwd mm2, [%3+24] ; x6*w15+x2*w13 x6*w11+x2*w09
261 :     paddd mm3, mm4 ; 4 ; a1=sum(even1) a0=sum(even0)
262 :    
263 :     pmaddwd mm5, [%3+48] ; x5*w30+x1*w28 x5*w26+x1*w24
264 :     movq mm4, mm3 ; 4 ; a1 a0
265 :    
266 :     pmaddwd mm6, [%3+56] ; x7*w31+x3*w29 x7*w27+x3*w25
267 :     paddd mm1, mm7 ; 7 ; b1=sum(odd1) b0=sum(odd0)
268 :    
269 :     paddd mm0, [%4] ; +%4
270 :     psubd mm3, mm1 ; a1-b1 a0-b0
271 :    
272 :     psrad mm3, SHIFT_INV_ROW ; y6=a1-b1 y7=a0-b0
273 :     paddd mm1, mm4 ; 4 ; a1+b1 a0+b0
274 :    
275 :     paddd mm0, mm2 ; 2 ; a3=sum(even3) a2=sum(even2)
276 :     psrad mm1, SHIFT_INV_ROW ; y1=a1+b1 y0=a0+b0
277 :    
278 :     paddd mm5, mm6 ; 6 ; b3=sum(odd3) b2=sum(odd2)
279 :     movq mm4, mm0 ; 4 ; a3 a2
280 :    
281 :     paddd mm0, mm5 ; a3+b3 a2+b2
282 :     psubd mm4, mm5 ; 5 ; a3-b3 a2-b2
283 :    
284 :     psrad mm0, SHIFT_INV_ROW ; y3=a3+b3 y2=a2+b2
285 :     psrad mm4, SHIFT_INV_ROW ; y4=a3-b3 y5=a2-b2
286 :    
287 :     packssdw mm1, mm0 ; 0 ; y3 y2 y1 y0
288 :     packssdw mm4, mm3 ; 3 ; y6 y7 y4 y5
289 :    
290 :     movq mm7, mm4 ; 7 ; y6 y7 y4 y5
291 :     psrld mm4, 16 ; 0 y6 0 y4
292 :    
293 :     pslld mm7, 16 ; y7 0 y5 0
294 :     movq [%2], mm1 ; 1 ; save y3 y2 y1 y0
295 :    
296 :     por mm7, mm4 ; 4 ; y7 y6 y5 y4
297 :     movq [%2+8], mm7 ; 7 ; save y7 y6 y5 y4
298 :     %endmacro
299 :    
300 :    
301 :    
302 :    
303 :     ;=============================================================================
304 :     ; code for Pentium III
305 :     ;=============================================================================
306 :    
307 :     ; %3 for rows 0,4 - constants are multiplied by cos_4_16
308 :    
309 :     tab_i_04_sse dw 16384, 21407, 16384, 8867 ; movq-> w05 w04 w01 w00
310 :     dw 16384, 8867, -16384, -21407 ; w07 w06 w03 w02
311 :     dw 16384, -8867, 16384, -21407 ; w13 w12 w09 w08
312 :     dw -16384, 21407, 16384, -8867 ; w15 w14 w11 w10
313 :     dw 22725, 19266, 19266, -4520 ; w21 w20 w17 w16
314 :     dw 12873, 4520, -22725, -12873 ; w23 w22 w19 w18
315 :     dw 12873, -22725, 4520, -12873 ; w29 w28 w25 w24
316 :     dw 4520, 19266, 19266, -22725 ; w31 w30 w27 w26
317 :    
318 :     ; %3 for rows 1,7 - constants are multiplied by cos_1_16
319 :    
320 :     tab_i_17_sse dw 22725, 29692, 22725, 12299 ; movq-> w05 w04 w01 w00
321 :     dw 22725, 12299, -22725, -29692 ; w07 w06 w03 w02
322 :     dw 22725, -12299, 22725, -29692 ; w13 w12 w09 w08
323 :     dw -22725, 29692, 22725, -12299 ; w15 w14 w11 w10
324 :     dw 31521, 26722, 26722, -6270 ; w21 w20 w17 w16
325 :     dw 17855, 6270, -31521, -17855 ; w23 w22 w19 w18
326 :     dw 17855, -31521, 6270, -17855 ; w29 w28 w25 w24
327 :     dw 6270, 26722, 26722, -31521 ; w31 w30 w27 w26
328 :    
329 :     ; %3 for rows 2,6 - constants are multiplied by cos_2_16
330 :    
331 :     tab_i_26_sse dw 21407, 27969, 21407, 11585 ; movq-> w05 w04 w01 w00
332 :     dw 21407, 11585, -21407, -27969 ; w07 w06 w03 w02
333 :     dw 21407, -11585, 21407, -27969 ; w13 w12 w09 w08
334 :     dw -21407, 27969, 21407, -11585 ; w15 w14 w11 w10
335 :     dw 29692, 25172, 25172, -5906 ; w21 w20 w17 w16
336 :     dw 16819, 5906, -29692, -16819 ; w23 w22 w19 w18
337 :     dw 16819, -29692, 5906, -16819 ; w29 w28 w25 w24
338 :     dw 5906, 25172, 25172, -29692 ; w31 w30 w27 w26
339 :    
340 :     ; %3 for rows 3,5 - constants are multiplied by cos_3_16
341 :    
342 :     tab_i_35_sse dw 19266, 25172, 19266, 10426 ; movq-> w05 w04 w01 w00
343 :     dw 19266, 10426, -19266, -25172 ; w07 w06 w03 w02
344 :     dw 19266, -10426, 19266, -25172 ; w13 w12 w09 w08
345 :     dw -19266, 25172, 19266, -10426 ; w15 w14 w11 w10
346 :     dw 26722, 22654, 22654, -5315 ; w21 w20 w17 w16
347 :     dw 15137, 5315, -26722, -15137 ; w23 w22 w19 w18
348 :     dw 15137, -26722, 5315, -15137 ; w29 w28 w25 w24
349 :     dw 5315, 22654, 22654, -26722 ; w31 w30 w27 w26
350 :    
351 :     ;-----------------------------------------------------------------------------
352 :    
353 :     ;
354 :     ; DCT_8_INV_ROW_1_sse INP, OUT, TABLE, ROUNDER
355 :     ;
356 :    
357 :     %macro DCT_8_INV_ROW_1_sse 4
358 :    
359 :     movq mm0, [%1] ; 0 ; x3 x2 x1 x0
360 :    
361 :     movq mm1, [%1+8] ; 1 ; x7 x6 x5 x4
362 :     movq mm2, mm0 ; 2 ; x3 x2 x1 x0
363 :    
364 :     movq mm3, [%3] ; 3 ; w05 w04 w01 w00
365 :     pshufw mm0, mm0, 10001000b ; x2 x0 x2 x0
366 :    
367 :     movq mm4, [%3+8] ; 4 ; w07 w06 w03 w02
368 :     movq mm5, mm1 ; 5 ; x7 x6 x5 x4
369 :     pmaddwd mm3, mm0 ; x2*w05+x0*w04 x2*w01+x0*w00
370 :    
371 :     movq mm6, [%3+32] ; 6 ; w21 w20 w17 w16
372 :     pshufw mm1, mm1, 10001000b ; x6 x4 x6 x4
373 :     pmaddwd mm4, mm1 ; x6*w07+x4*w06 x6*w03+x4*w02
374 :    
375 :     movq mm7, [%3+40] ; 7 ; w23 w22 w19 w18
376 :     pshufw mm2, mm2, 11011101b ; x3 x1 x3 x1
377 :     pmaddwd mm6, mm2 ; x3*w21+x1*w20 x3*w17+x1*w16
378 :    
379 :     pshufw mm5, mm5, 11011101b ; x7 x5 x7 x5
380 :     pmaddwd mm7, mm5 ; x7*w23+x5*w22 x7*w19+x5*w18
381 :    
382 :     paddd mm3, [%4] ; +%4
383 :    
384 :     pmaddwd mm0, [%3+16] ; x2*w13+x0*w12 x2*w09+x0*w08
385 :     paddd mm3, mm4 ; 4 ; a1=sum(even1) a0=sum(even0)
386 :    
387 :     pmaddwd mm1, [%3+24] ; x6*w15+x4*w14 x6*w11+x4*w10
388 :     movq mm4, mm3 ; 4 ; a1 a0
389 :    
390 :     pmaddwd mm2, [%3+48] ; x3*w29+x1*w28 x3*w25+x1*w24
391 :     paddd mm6, mm7 ; 7 ; b1=sum(odd1) b0=sum(odd0)
392 :    
393 :     pmaddwd mm5, [%3+56] ; x7*w31+x5*w30 x7*w27+x5*w26
394 :     paddd mm3, mm6 ; a1+b1 a0+b0
395 :    
396 :     paddd mm0, [%4] ; +%4
397 :     psrad mm3, SHIFT_INV_ROW ; y1=a1+b1 y0=a0+b0
398 :    
399 :     paddd mm0, mm1 ; 1 ; a3=sum(even3) a2=sum(even2)
400 :     psubd mm4, mm6 ; 6 ; a1-b1 a0-b0
401 :    
402 :     movq mm7, mm0 ; 7 ; a3 a2
403 :     paddd mm2, mm5 ; 5 ; b3=sum(odd3) b2=sum(odd2)
404 :    
405 :     paddd mm0, mm2 ; a3+b3 a2+b2
406 :     psrad mm4, SHIFT_INV_ROW ; y6=a1-b1 y7=a0-b0
407 :    
408 :     psubd mm7, mm2 ; 2 ; a3-b3 a2-b2
409 :     psrad mm0, SHIFT_INV_ROW ; y3=a3+b3 y2=a2+b2
410 :    
411 :     psrad mm7, SHIFT_INV_ROW ; y4=a3-b3 y5=a2-b2
412 :    
413 :     packssdw mm3, mm0 ; 0 ; y3 y2 y1 y0
414 :    
415 :     packssdw mm7, mm4 ; 4 ; y6 y7 y4 y5
416 :    
417 :     movq [%2], mm3 ; 3 ; save y3 y2 y1 y0
418 :     pshufw mm7, mm7, 10110001b ; y7 y6 y5 y4
419 :    
420 :     movq [%2+8], mm7 ; 7 ; save y7 y6 y5 y4
421 :    
422 :     %endmacro
423 :    
424 :    
425 :     ;=============================================================================
426 :     ;
427 :     ;=============================================================================
428 :    
429 :     ;=============================================================================
430 :     ;
431 :     ; The first stage DCT 8x8 - forward DCTs of columns
432 :     ;
433 :     ; The %2puts are multiplied
434 :     ; for rows 0,4 - on cos_4_16,
435 :     ; for rows 1,7 - on cos_1_16,
436 :     ; for rows 2,6 - on cos_2_16,
437 :     ; for rows 3,5 - on cos_3_16
438 :     ; and are shifted to the left for rise of accuracy
439 :     ;
440 :     ;-----------------------------------------------------------------------------
441 :     ;
442 :     ; The 8-point scaled forward DCT algorithm (26a8m)
443 :     ;
444 :     ;-----------------------------------------------------------------------------
445 :     ;
446 :     ; #define DCT_8_FRW_COL(x, y)
447 :     ;{
448 :     ; short t0, t1, t2, t3, t4, t5, t6, t7;
449 :     ; short tp03, tm03, tp12, tm12, tp65, tm65;
450 :     ; short tp465, tm465, tp765, tm765;
451 :     ;
452 :     ; t0 = LEFT_SHIFT ( x[0] + x[7] );
453 :     ; t1 = LEFT_SHIFT ( x[1] + x[6] );
454 :     ; t2 = LEFT_SHIFT ( x[2] + x[5] );
455 :     ; t3 = LEFT_SHIFT ( x[3] + x[4] );
456 :     ; t4 = LEFT_SHIFT ( x[3] - x[4] );
457 :     ; t5 = LEFT_SHIFT ( x[2] - x[5] );
458 :     ; t6 = LEFT_SHIFT ( x[1] - x[6] );
459 :     ; t7 = LEFT_SHIFT ( x[0] - x[7] );
460 :     ;
461 :     ; tp03 = t0 + t3;
462 :     ; tm03 = t0 - t3;
463 :     ; tp12 = t1 + t2;
464 :     ; tm12 = t1 - t2;
465 :     ;
466 :     ; y[0] = tp03 + tp12;
467 :     ; y[4] = tp03 - tp12;
468 :     ;
469 :     ; y[2] = tm03 + tm12 * tg_2_16;
470 :     ; y[6] = tm03 * tg_2_16 - tm12;
471 :     ;
472 :     ; tp65 =(t6 +t5 )*cos_4_16;
473 :     ; tm65 =(t6 -t5 )*cos_4_16;
474 :     ;
475 :     ; tp765 = t7 + tp65;
476 :     ; tm765 = t7 - tp65;
477 :     ; tp465 = t4 + tm65;
478 :     ; tm465 = t4 - tm65;
479 :     ;
480 :     ; y[1] = tp765 + tp465 * tg_1_16;
481 :     ; y[7] = tp765 * tg_1_16 - tp465;
482 :     ; y[5] = tm765 * tg_3_16 + tm465;
483 :     ; y[3] = tm765 - tm465 * tg_3_16;
484 :     ;}
485 :     ;
486 :     ;=============================================================================
487 :    
488 :    
489 :     ;
490 :     ; DCT_8_FRW_COL_4 INP, OUT
491 :     ;
492 :    
493 :     %macro DCT_8_FRW_COL_4 2
494 :    
495 :     LOCAL x0, x1, x2, x3, x4, x5, x6, x7
496 :     LOCAL y0, y1, y2, y3, y4, y5, y6, y7
497 :     x0 equ [%1 + 0*16]
498 :     x1 equ [%1 + 1*16]
499 :     x2 equ [%1 + 2*16]
500 :     x3 equ [%1 + 3*16]
501 :     x4 equ [%1 + 4*16]
502 :     x5 equ [%1 + 5*16]
503 :     x6 equ [%1 + 6*16]
504 :     x7 equ [%1 + 7*16]
505 :     y0 equ [%2 + 0*16]
506 :     y1 equ [%2 + 1*16]
507 :     y2 equ [%2 + 2*16]
508 :     y3 equ [%2 + 3*16]
509 :     y4 equ [%2 + 4*16]
510 :     y5 equ [%2 + 5*16]
511 :     y6 equ [%2 + 6*16]
512 :     y7 equ [%2 + 7*16]
513 :     movq mm0, x1 ; 0 ; x1
514 :     movq mm1, x6 ; 1 ; x6
515 :     movq mm2, mm0 ; 2 ; x1
516 :     movq mm3, x2 ; 3 ; x2
517 :     paddsw mm0, mm1 ; t1 = x[1] + x[6]
518 :     movq mm4, x5 ; 4 ; x5
519 :     psllw mm0, SHIFT_FRW_COL ; t1
520 :     movq mm5, x0 ; 5 ; x0
521 :     paddsw mm4, mm3 ; t2 = x[2] + x[5]
522 :     paddsw mm5, x7 ; t0 = x[0] + x[7]
523 :     psllw mm4, SHIFT_FRW_COL ; t2
524 :     movq mm6, mm0 ; 6 ; t1
525 :     psubsw mm2, mm1 ; 1 ; t6 = x[1] - x[6]
526 :     movq mm1, [tg_2_16] ; 1 ; tg_2_16
527 :     psubsw mm0, mm4 ; tm12 = t1 - t2
528 :     movq mm7, x3 ; 7 ; x3
529 :     pmulhw mm1, mm0 ; tm12*tg_2_16
530 :     paddsw mm7, x4 ; t3 = x[3] + x[4]
531 :     psllw mm5, SHIFT_FRW_COL ; t0
532 :     paddsw mm6, mm4 ; 4 ; tp12 = t1 + t2
533 :     psllw mm7, SHIFT_FRW_COL ; t3
534 :     movq mm4, mm5 ; 4 ; t0
535 :     psubsw mm5, mm7 ; tm03 = t0 - t3
536 :     paddsw mm1, mm5 ; y2 = tm03 + tm12*tg_2_16
537 :     paddsw mm4, mm7 ; 7 ; tp03 = t0 + t3
538 :     por mm1, [one_corr] ; correction y2 +0.5
539 :     psllw mm2, SHIFT_FRW_COL+1 ; t6
540 :     pmulhw mm5, [tg_2_16] ; tm03*tg_2_16
541 :     movq mm7, mm4 ; 7 ; tp03
542 :     psubsw mm3, x5 ; t5 = x[2] - x[5]
543 :     psubsw mm4, mm6 ; y4 = tp03 - tp12
544 :     movq y2, mm1 ; 1 ; save y2
545 :     paddsw mm7, mm6 ; 6 ; y0 = tp03 + tp12
546 :     movq mm1, x3 ; 1 ; x3
547 :     psllw mm3, SHIFT_FRW_COL+1 ; t5
548 :     psubsw mm1, x4 ; t4 = x[3] - x[4]
549 :     movq mm6, mm2 ; 6 ; t6
550 :     movq y4, mm4 ; 4 ; save y4
551 :     paddsw mm2, mm3 ; t6 + t5
552 :     pmulhw mm2, [ocos_4_16] ; tp65 = (t6 + t5)*cos_4_16
553 :     psubsw mm6, mm3 ; 3 ; t6 - t5
554 :     pmulhw mm6, [ocos_4_16] ; tm65 = (t6 - t5)*cos_4_16
555 :     psubsw mm5, mm0 ; 0 ; y6 = tm03*tg_2_16 - tm12
556 :     por mm5, [one_corr] ; correction y6 +0.5
557 :     psllw mm1, SHIFT_FRW_COL ; t4
558 :     por mm2, [one_corr] ; correction tp65 +0.5
559 :     movq mm4, mm1 ; 4 ; t4
560 :     movq mm3, x0 ; 3 ; x0
561 :     paddsw mm1, mm6 ; tp465 = t4 + tm65
562 :     psubsw mm3, x7 ; t7 = x[0] - x[7]
563 :     psubsw mm4, mm6 ; 6 ; tm465 = t4 - tm65
564 :     movq mm0, [tg_1_16] ; 0 ; tg_1_16
565 :     psllw mm3, SHIFT_FRW_COL ; t7
566 :     movq mm6, [tg_3_16] ; 6 ; tg_3_16
567 :     pmulhw mm0, mm1 ; tp465*tg_1_16
568 :     movq y0, mm7 ; 7 ; save y0
569 :     pmulhw mm6, mm4 ; tm465*tg_3_16
570 :     movq y6, mm5 ; 5 ; save y6
571 :     movq mm7, mm3 ; 7 ; t7
572 :     movq mm5, [tg_3_16] ; 5 ; tg_3_16
573 :     psubsw mm7, mm2 ; tm765 = t7 - tp65
574 :     paddsw mm3, mm2 ; 2 ; tp765 = t7 + tp65
575 :     pmulhw mm5, mm7 ; tm765*tg_3_16
576 :     paddsw mm0, mm3 ; y1 = tp765 + tp465*tg_1_16
577 :     paddsw mm6, mm4 ; tm465*tg_3_16
578 :     pmulhw mm3, [tg_1_16] ; tp765*tg_1_16
579 :     por mm0, [one_corr] ; correction y1 +0.5
580 :     paddsw mm5, mm7 ; tm765*tg_3_16
581 :     psubsw mm7, mm6 ; 6 ; y3 = tm765 - tm465*tg_3_16
582 :     movq y1, mm0 ; 0 ; save y1
583 :     paddsw mm5, mm4 ; 4 ; y5 = tm765*tg_3_16 + tm465
584 :     movq y3, mm7 ; 7 ; save y3
585 :     psubsw mm3, mm1 ; 1 ; y7 = tp765*tg_1_16 - tp465
586 :     movq y5, mm5 ; 5 ; save y5
587 :     movq y7, mm3 ; 3 ; save y7
588 :     %endmacro
589 :    
590 :    
591 :     ;
592 :     ; DCT_8_INV_COL_4 INP,OUT
593 :     ;
594 :    
595 :     %macro DCT_8_INV_COL_4 2
596 :     movq mm0, [tg_3_16]
597 :    
598 :     movq mm3, [%1+16*3]
599 :     movq mm1, mm0 ; tg_3_16
600 :    
601 :     movq mm5, [%1+16*5]
602 :     pmulhw mm0, mm3 ; x3*(tg_3_16-1)
603 :    
604 :     movq mm4, [tg_1_16]
605 :     pmulhw mm1, mm5 ; x5*(tg_3_16-1)
606 :    
607 :     movq mm7, [%1+16*7]
608 :     movq mm2, mm4 ; tg_1_16
609 :    
610 :     movq mm6, [%1+16*1]
611 :     pmulhw mm4, mm7 ; x7*tg_1_16
612 :    
613 :     paddsw mm0, mm3 ; x3*tg_3_16
614 :     pmulhw mm2, mm6 ; x1*tg_1_16
615 :    
616 :     paddsw mm1, mm3 ; x3+x5*(tg_3_16-1)
617 :     psubsw mm0, mm5 ; x3*tg_3_16-x5 = tm35
618 :    
619 :     movq mm3, [ocos_4_16]
620 :     paddsw mm1, mm5 ; x3+x5*tg_3_16 = tp35
621 :    
622 :     paddsw mm4, mm6 ; x1+tg_1_16*x7 = tp17
623 :     psubsw mm2, mm7 ; x1*tg_1_16-x7 = tm17
624 :    
625 :     movq mm5, mm4 ; tp17
626 :     movq mm6, mm2 ; tm17
627 :    
628 :     paddsw mm5, mm1 ; tp17+tp35 = b0
629 :     psubsw mm6, mm0 ; tm17-tm35 = b3
630 :    
631 :     psubsw mm4, mm1 ; tp17-tp35 = t1
632 :     paddsw mm2, mm0 ; tm17+tm35 = t2
633 :    
634 :     movq mm7, [tg_2_16]
635 :     movq mm1, mm4 ; t1
636 :    
637 :     ; movq [SCRATCH+0], mm5 ; save b0
638 :     movq [%2+3*16], mm5 ; save b0
639 :     paddsw mm1, mm2 ; t1+t2
640 :    
641 :     ; movq [SCRATCH+8], mm6 ; save b3
642 :     movq [%2+5*16], mm6 ; save b3
643 :     psubsw mm4, mm2 ; t1-t2
644 :    
645 :     movq mm5, [%1+2*16]
646 :     movq mm0, mm7 ; tg_2_16
647 :    
648 :     movq mm6, [%1+6*16]
649 :     pmulhw mm0, mm5 ; x2*tg_2_16
650 :    
651 :     pmulhw mm7, mm6 ; x6*tg_2_16
652 :     ; slot
653 :     pmulhw mm1, mm3 ; ocos_4_16*(t1+t2) = b1/2
654 :     ; slot
655 :     movq mm2, [%1+0*16]
656 :     pmulhw mm4, mm3 ; ocos_4_16*(t1-t2) = b2/2
657 :    
658 :     psubsw mm0, mm6 ; t2*tg_2_16-x6 = tm26
659 :     movq mm3, mm2 ; x0
660 :    
661 :     movq mm6, [%1+4*16]
662 :     paddsw mm7, mm5 ; x2+x6*tg_2_16 = tp26
663 :    
664 :     paddsw mm2, mm6 ; x0+x4 = tp04
665 :     psubsw mm3, mm6 ; x0-x4 = tm04
666 :    
667 :     movq mm5, mm2 ; tp04
668 :     movq mm6, mm3 ; tm04
669 :    
670 :     psubsw mm2, mm7 ; tp04-tp26 = a3
671 :     paddsw mm3, mm0 ; tm04+tm26 = a1
672 :    
673 :     paddsw mm1, mm1 ; b1
674 :     paddsw mm4, mm4 ; b2
675 :    
676 :     paddsw mm5, mm7 ; tp04+tp26 = a0
677 :     psubsw mm6, mm0 ; tm04-tm26 = a2
678 :    
679 :     movq mm7, mm3 ; a1
680 :     movq mm0, mm6 ; a2
681 :    
682 :     paddsw mm3, mm1 ; a1+b1
683 :     paddsw mm6, mm4 ; a2+b2
684 :    
685 :     psraw mm3, SHIFT_INV_COL ; dst1
686 :     psubsw mm7, mm1 ; a1-b1
687 :    
688 :     psraw mm6, SHIFT_INV_COL ; dst2
689 :     psubsw mm0, mm4 ; a2-b2
690 :    
691 :     ; movq mm1, [SCRATCH+0] ; load b0
692 :     movq mm1, [%2+3*16] ; load b0
693 :     psraw mm7, SHIFT_INV_COL ; dst6
694 :    
695 :     movq mm4, mm5 ; a0
696 :     psraw mm0, SHIFT_INV_COL ; dst5
697 :    
698 :     movq [%2+1*16], mm3
699 :     paddsw mm5, mm1 ; a0+b0
700 :    
701 :     movq [%2+2*16], mm6
702 :     psubsw mm4, mm1 ; a0-b0
703 :    
704 :     ; movq mm3, [SCRATCH+8] ; load b3
705 :     movq mm3, [%2+5*16] ; load b3
706 :     psraw mm5, SHIFT_INV_COL ; dst0
707 :    
708 :     movq mm6, mm2 ; a3
709 :     psraw mm4, SHIFT_INV_COL ; dst7
710 :    
711 :     movq [%2+5*16], mm0
712 :     paddsw mm2, mm3 ; a3+b3
713 :    
714 :     movq [%2+6*16], mm7
715 :     psubsw mm6, mm3 ; a3-b3
716 :    
717 :     movq [%2+0*16], mm5
718 :     psraw mm2, SHIFT_INV_COL ; dst3
719 :    
720 :     movq [%2+7*16], mm4
721 :     psraw mm6, SHIFT_INV_COL ; dst4
722 :    
723 :     movq [%2+3*16], mm2
724 :    
725 :     movq [%2+4*16], mm6
726 :     %endmacro
727 :    
728 :    
729 :    
730 :     section .text
731 :    
732 :     ;=============================================================================
733 :     ;
734 :     ; void idct_mmx (short * const src_result);
735 :     ;
736 :     ;=============================================================================
737 :    
738 :     align 16
739 :     cglobal idct_mmx
740 :     idct_mmx
741 :     mov eax, dword [esp + 4]
742 :    
743 :     DCT_8_INV_ROW_1 eax+0, eax+0, tab_i_04, rounder_0
744 :     DCT_8_INV_ROW_1 eax+16, eax+16, tab_i_17, rounder_1
745 :     DCT_8_INV_ROW_1 eax+32, eax+32, tab_i_26, rounder_2
746 :     DCT_8_INV_ROW_1 eax+48, eax+48, tab_i_35, rounder_3
747 :     DCT_8_INV_ROW_1 eax+64, eax+64, tab_i_04, rounder_4
748 :     DCT_8_INV_ROW_1 eax+80, eax+80, tab_i_35, rounder_5
749 :     DCT_8_INV_ROW_1 eax+96, eax+96, tab_i_26, rounder_6
750 :     DCT_8_INV_ROW_1 eax+112, eax+112, tab_i_17, rounder_7
751 :    
752 :     DCT_8_INV_COL_4 eax+0,eax+0
753 :     DCT_8_INV_COL_4 eax+8,eax+8
754 :    
755 :     ret
756 :    
757 :    
758 :    
759 :     ;=============================================================================
760 :     ;
761 :     ; void idct_sse (short * const src_result);
762 :     ;
763 :     ;=============================================================================
764 :    
765 :     align 16
766 :     cglobal idct_xmm
767 :     idct_xmm
768 :     mov eax, dword [esp + 4]
769 :    
770 :     DCT_8_INV_ROW_1_sse eax+0, eax+0, tab_i_04_sse, rounder_0
771 :     DCT_8_INV_ROW_1_sse eax+16, eax+16, tab_i_17_sse, rounder_1
772 :     DCT_8_INV_ROW_1_sse eax+32, eax+32, tab_i_26_sse, rounder_2
773 :     DCT_8_INV_ROW_1_sse eax+48, eax+48, tab_i_35_sse, rounder_3
774 :     DCT_8_INV_ROW_1_sse eax+64, eax+64, tab_i_04_sse, rounder_4
775 :     DCT_8_INV_ROW_1_sse eax+80, eax+80, tab_i_35_sse, rounder_5
776 :     DCT_8_INV_ROW_1_sse eax+96, eax+96, tab_i_26_sse, rounder_6
777 :     DCT_8_INV_ROW_1_sse eax+112, eax+112, tab_i_17_sse, rounder_7
778 :    
779 :     DCT_8_INV_COL_4 eax+0, eax+0
780 :     DCT_8_INV_COL_4 eax+8, eax+8
781 :    
782 : Isibaar 154 ret
783 :    
784 :     ;=============================================================================
785 :     ; The code below this line is for SSE2-equipped processors
786 :     ; By Dmitry Rozhdestvensky
787 :     ;=============================================================================
788 :    
789 :     section .data
790 :    
791 :     align 16
792 :    
793 :     tab_i_04_s2 dw 16384, 21407, 16384, 8867 ; movq-> w05 w04 w01 w00
794 :     dw 16384, -8867, 16384, -21407 ; w13 w12 w09 w08
795 :     dw 16384, 8867, -16384, -21407 ; w07 w06 w03 w02
796 :     dw -16384, 21407, 16384, -8867 ; w15 w14 w11 w10
797 :     dw 22725, 19266, 19266, -4520 ; w21 w20 w17 w16
798 :     dw 12873, -22725, 4520, -12873 ; w29 w28 w25 w24
799 :     dw 12873, 4520, -22725, -12873 ; w23 w22 w19 w18
800 :     dw 4520, 19266, 19266, -22725 ; w31 w30 w27 w26
801 :    
802 :     ; Table for rows 1,7 - constants are multiplied by cos_1_16
803 :    
804 :     tab_i_17_s2 dw 22725, 29692, 22725, 12299 ; movq-> w05 w04 w01 w00
805 :     dw 22725, -12299, 22725, -29692 ; w13 w12 w09 w08
806 :     dw 22725, 12299, -22725, -29692 ; w07 w06 w03 w02
807 :     dw -22725, 29692, 22725, -12299 ; w15 w14 w11 w10
808 :     dw 31521, 26722, 26722, -6270 ; w21 w20 w17 w16
809 :     dw 17855, -31521, 6270, -17855 ; w29 w28 w25 w24
810 :     dw 17855, 6270, -31521, -17855 ; w23 w22 w19 w18
811 :     dw 6270, 26722, 26722, -31521 ; w31 w30 w27 w26
812 :    
813 :     ; Table for rows 2,6 - constants are multiplied by cos_2_16
814 :    
815 :     tab_i_26_s2 dw 21407, 27969, 21407, 11585 ; movq-> w05 w04 w01 w00
816 :     dw 21407, -11585, 21407, -27969 ; w13 w12 w09 w08
817 :     dw 21407, 11585, -21407, -27969 ; w07 w06 w03 w02
818 :     dw -21407, 27969, 21407, -11585 ; w15 w14 w11 w10
819 :     dw 29692, 25172, 25172, -5906 ; w21 w20 w17 w16
820 :     dw 16819, -29692, 5906, -16819 ; w29 w28 w25 w24
821 :     dw 16819, 5906, -29692, -16819 ; w23 w22 w19 w18
822 :     dw 5906, 25172, 25172, -29692 ; w31 w30 w27 w26
823 :    
824 :     ; Table for rows 3,5 - constants are multiplied by cos_3_16
825 :    
826 :     tab_i_35_s2 dw 19266, 25172, 19266, 10426 ; movq-> w05 w04 w01 w00
827 :     dw 19266, -10426, 19266, -25172 ; w13 w12 w09 w08
828 :     dw 19266, 10426, -19266, -25172 ; w07 w06 w03 w02
829 :     dw -19266, 25172, 19266, -10426 ; w15 w14 w11 w10
830 :     dw 26722, 22654, 22654, -5315 ; w21 w20 w17 w16
831 :     dw 15137, -26722, 5315, -15137 ; w29 w28 w25 w24
832 :     dw 15137, 5315, -26722, -15137 ; w23 w22 w19 w18
833 :     dw 5315, 22654, 22654, -26722 ; w31 w30 w27 w26
834 :    
835 :     %if SHIFT_INV_ROW == 12 ; assume SHIFT_INV_ROW == 12
836 :     rounder_2_0 dd 65536, 65536
837 :     dd 65536, 65536
838 :     rounder_2_4 dd 0, 0
839 :     dd 0, 0
840 :     rounder_2_1 dd 7195, 7195
841 :     dd 7195, 7195
842 :     rounder_2_7 dd 1024, 1024
843 :     dd 1024, 1024
844 :     rounder_2_2 dd 4520, 4520
845 :     dd 4520, 4520
846 :     rounder_2_6 dd 1024, 1024
847 :     dd 1024, 1024
848 :     rounder_2_3 dd 2407, 2407
849 :     dd 2407, 2407
850 :     rounder_2_5 dd 240, 240
851 :     dd 240, 240
852 :    
853 :     %elif SHIFT_INV_ROW == 11 ; assume SHIFT_INV_ROW == 11
854 :     rounder_2_0 dd 65536, 65536
855 :     dd 65536, 65536
856 :     rounder_2_4 dd 0, 0
857 :     dd 0, 0
858 :     rounder_2_1 dd 3597, 3597
859 :     dd 3597, 3597
860 :     rounder_2_7 dd 512, 512
861 :     dd 512, 512
862 :     rounder_2_2 dd 2260, 2260
863 :     dd 2260, 2260
864 :     rounder_2_6 dd 512, 512
865 :     dd 512, 512
866 :     rounder_2_3 dd 1203, 1203
867 :     dd 1203, 1203
868 :     rounder_2_5 dd 120, 120
869 :     dd 120, 120
870 :     %else
871 :    
872 :     %error invalid _SHIFT_INV_ROW_
873 :    
874 :     %endif
875 :    
876 :     tg_1_16_2 dw 13036, 13036, 13036, 13036 ; tg * (2<<16) + 0.5
877 :     dw 13036, 13036, 13036, 13036
878 :     tg_2_16_2 dw 27146, 27146, 27146, 27146 ; tg * (2<<16) + 0.5
879 :     dw 27146, 27146, 27146, 27146
880 :     tg_3_16_2 dw -21746, -21746, -21746, -21746 ; tg * (2<<16) + 0.5
881 :     dw -21746, -21746, -21746, -21746
882 :     ocos_4_16_2 dw 23170, 23170, 23170, 23170 ; cos * (2<<15) + 0.5
883 :     dw 23170, 23170, 23170, 23170
884 :    
885 :     %macro DCT_8_INV_ROW_1_sse2 4
886 :    
887 :     pshufhw xmm1,[%1],11011000b ;x 75643210
888 :     pshuflw xmm1,xmm1,11011000b ;x 75643120
889 :     pshufd xmm0,xmm1,00000000b ;x 20202020
890 :     pmaddwd xmm0,[%3] ;w 13 12 9 8 5410
891 :     ;a 3210 first part
892 :    
893 :     pshufd xmm2,xmm1,10101010b ;x 64646464
894 :     pmaddwd xmm2,[%3+16] ;w 15 14 11 10 7632
895 :     ;a 3210 second part
896 :    
897 :     paddd xmm2,xmm0 ;a 3210 ready
898 :     paddd xmm2,[%4] ;must be 4 dwords long, not 2 as for sse1
899 :     movdqa xmm5,xmm2
900 :    
901 :     pshufd xmm3,xmm1,01010101b ;x 31313131
902 :     pmaddwd xmm3,[%3+32] ;w 29 28 25 24 21 20 17 16
903 :     ;b 3210 first part
904 :    
905 :     pshufd xmm4,xmm1,11111111b ;x 75757575
906 :     pmaddwd xmm4,[%3+48] ;w 31 30 27 26 23 22 19 18
907 :     ;b 3210 second part
908 :     paddd xmm3,xmm4 ;b 3210 ready
909 :    
910 :     paddd xmm2,xmm3 ;will be y 3210
911 :     psubd xmm5,xmm3 ;will be y 4567
912 :     psrad xmm2,SHIFT_INV_ROW
913 :     psrad xmm5,SHIFT_INV_ROW
914 :     packssdw xmm2,xmm5 ;y 45673210
915 :     pshufhw xmm6,xmm2,00011011b ;y 76543210
916 :     movdqa [%2],xmm6
917 :    
918 :     %endmacro
919 :    
920 :     %macro DCT_8_INV_COL_4_sse2 2
921 :    
922 :     movdqa xmm0,[%1+16*0] ;x0 (all columns)
923 :     movdqa xmm2,[%1+16*4] ;x4
924 :     movdqa xmm1,xmm0
925 :    
926 :     movdqa xmm4,[%1+16*2] ;x2
927 :     movdqa xmm5,[%1+16*6] ;x6
928 :     movdqa xmm6,[tg_2_16_2]
929 :     movdqa xmm7,xmm6
930 :    
931 :     paddsw xmm0,xmm2 ;u04=x0+x4
932 :     psubsw xmm1,xmm2 ;v04=x0-x4
933 :     movdqa xmm3,xmm0
934 :     movdqa xmm2,xmm1
935 :    
936 :     pmulhw xmm6,xmm4
937 :     pmulhw xmm7,xmm5
938 :     psubsw xmm6,xmm5 ;v26=x2*T2-x6
939 :     paddsw xmm7,xmm4 ;u26=x6*T2+x2
940 :    
941 :     paddsw xmm1,xmm6 ;a1=v04+v26
942 :     paddsw xmm0,xmm7 ;a0=u04+u26
943 :     psubsw xmm2,xmm6 ;a2=v04-v26
944 :     psubsw xmm3,xmm7 ;a3=u04-u26
945 :    
946 :     movdqa [%2+16*0],xmm0 ;store a3-a0 to
947 :     movdqa [%2+16*6],xmm1 ;free registers
948 :     movdqa [%2+16*2],xmm2
949 :     movdqa [%2+16*4],xmm3
950 :    
951 :     movdqa xmm0,[%1+16*1] ;x1
952 :     movdqa xmm1,[%1+16*7] ;x7
953 :     movdqa xmm2,[tg_1_16_2]
954 :     movdqa xmm3,xmm2
955 :    
956 :     movdqa xmm4,[%1+16*3] ;x3
957 :     movdqa xmm5,[%1+16*5] ;x5
958 :     movdqa xmm6,[tg_3_16_2]
959 :     movdqa xmm7,xmm6
960 :    
961 :     pmulhw xmm2,xmm0
962 :     pmulhw xmm3,xmm1
963 :     psubsw xmm2,xmm1 ;v17=x1*T1-x7
964 :     paddsw xmm3,xmm0 ;u17=x7*T1+x1
965 :     movdqa xmm0,xmm3 ;u17
966 :     movdqa xmm1,xmm2 ;v17
967 :    
968 :     pmulhw xmm6,xmm4 ;x3*(t3-1)
969 :     pmulhw xmm7,xmm5 ;x5*(t3-1)
970 :     paddsw xmm6,xmm4
971 :     paddsw xmm7,xmm5
972 :     psubsw xmm6,xmm5 ;v35=x3*T3-x5
973 :     paddsw xmm7,xmm4 ;u35=x5*T3+x3
974 :    
975 :     movdqa xmm4,[ocos_4_16_2]
976 :    
977 :     paddsw xmm0,xmm7 ;b0=u17+u35
978 :     psubsw xmm1,xmm6 ;b3=v17-v35
979 :     psubsw xmm3,xmm7 ;u12=u17-v35
980 :     paddsw xmm2,xmm6 ;v12=v17+v35
981 :    
982 :     movdqa xmm5,xmm3
983 :     paddsw xmm3,xmm2 ;tb1
984 :     psubsw xmm5,xmm2 ;tb2
985 :     pmulhw xmm5,xmm4
986 :     pmulhw xmm4,xmm3
987 :     paddsw xmm5,xmm5
988 :     paddsw xmm4,xmm4
989 :    
990 :     movdqa xmm6,[%2+16*0] ;a0
991 :     movdqa xmm7,xmm6
992 :     movdqa xmm2,[%2+16*4] ;a3
993 :     movdqa xmm3,xmm2
994 :    
995 :     paddsw xmm6,xmm0
996 :     psubsw xmm7,xmm0
997 :     psraw xmm6,SHIFT_INV_COL ;y0=a0+b0
998 :     psraw xmm7,SHIFT_INV_COL ;y7=a0-b0
999 :     movdqa [%2+16*0],xmm6
1000 :     movdqa [%2+16*7],xmm7
1001 :    
1002 :     paddsw xmm2,xmm1
1003 :     psubsw xmm3,xmm1
1004 :     psraw xmm2,SHIFT_INV_COL ;y3=a3+b3
1005 :     psraw xmm3,SHIFT_INV_COL ;y4=a3-b3
1006 :     movdqa [%2+16*3],xmm2
1007 :     movdqa [%2+16*4],xmm3
1008 :    
1009 :     movdqa xmm0,[%2+16*6] ;a1
1010 :     movdqa xmm1,xmm0
1011 :     movdqa xmm6,[%2+16*2] ;a2
1012 :     movdqa xmm7,xmm6
1013 :    
1014 :    
1015 :     paddsw xmm0,xmm4
1016 :     psubsw xmm1,xmm4
1017 :     psraw xmm0,SHIFT_INV_COL ;y1=a1+b1
1018 :     psraw xmm1,SHIFT_INV_COL ;y6=a1-b1
1019 :     movdqa [%2+16*1],xmm0
1020 :     movdqa [%2+16*6],xmm1
1021 :    
1022 :     paddsw xmm6,xmm5
1023 :     psubsw xmm7,xmm5
1024 :     psraw xmm6,SHIFT_INV_COL ;y2=a2+b2
1025 :     psraw xmm7,SHIFT_INV_COL ;y5=a2-b2
1026 :     movdqa [%2+16*2],xmm6
1027 :     movdqa [%2+16*5],xmm7
1028 :    
1029 :     %endmacro
1030 :    
1031 :     section .text
1032 :    
1033 :     align 16
1034 :     cglobal idct_sse2
1035 :     idct_sse2
1036 :    
1037 :     mov eax, dword [esp + 4]
1038 :    
1039 :     DCT_8_INV_ROW_1_sse2 eax+0, eax+0, tab_i_04_s2, rounder_2_0
1040 :     DCT_8_INV_ROW_1_sse2 eax+16, eax+16, tab_i_17_s2, rounder_2_1
1041 :     DCT_8_INV_ROW_1_sse2 eax+32, eax+32, tab_i_26_s2, rounder_2_2
1042 :     DCT_8_INV_ROW_1_sse2 eax+48, eax+48, tab_i_35_s2, rounder_2_3
1043 :     DCT_8_INV_ROW_1_sse2 eax+64, eax+64, tab_i_04_s2, rounder_2_4
1044 :     DCT_8_INV_ROW_1_sse2 eax+80, eax+80, tab_i_35_s2, rounder_2_5
1045 :     DCT_8_INV_ROW_1_sse2 eax+96, eax+96, tab_i_26_s2, rounder_2_6
1046 :     DCT_8_INV_ROW_1_sse2 eax+112, eax+112, tab_i_17_s2, rounder_2_7
1047 :    
1048 :     DCT_8_INV_COL_4_sse2 eax, eax
1049 :     ;DCT_8_INV_COL_4 eax+8, eax+8
1050 :    
1051 :     ret

No admin address has been configured
ViewVC Help
Powered by ViewVC 1.0.4