1 ;*****************************************************************************
3 ;*****************************************************************************
4 ;* Copyright (C) 2008 x264 project
6 ;* Authors: Holger Lubitz <holger@lubitz.org>
7 ;* Loren Merritt <lorenm@u.washington.edu>
9 ;* This program is free software; you can redistribute it and/or modify
10 ;* it 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.
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.
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., 51 Franklin Street, Fifth Floor, Boston, MA 02111, USA.
22 ;*****************************************************************************
24 %assign FENC_STRIDE 16
25 %assign FDEC_STRIDE 32
41 %macro TRANSPOSE4x4W 5
42 SBUTTERFLY wd, %1, %2, %5
43 SBUTTERFLY wd, %3, %4, %5
44 SBUTTERFLY dq, %1, %3, %5
45 SBUTTERFLY dq, %2, %4, %5
49 %macro TRANSPOSE2x4x4W 5
50 SBUTTERFLY wd, %1, %2, %5
51 SBUTTERFLY wd, %3, %4, %5
52 SBUTTERFLY dq, %1, %3, %5
53 SBUTTERFLY dq, %2, %4, %5
54 SBUTTERFLY qdq, %1, %2, %5
55 SBUTTERFLY qdq, %3, %4, %5
58 %macro TRANSPOSE4x4D 5
59 SBUTTERFLY dq, %1, %2, %5
60 SBUTTERFLY dq, %3, %4, %5
61 SBUTTERFLY qdq, %1, %3, %5
62 SBUTTERFLY qdq, %2, %4, %5
66 %macro TRANSPOSE8x8W 9-11
68 SBUTTERFLY wd, %1, %2, %9
69 SBUTTERFLY wd, %3, %4, %9
70 SBUTTERFLY wd, %5, %6, %9
71 SBUTTERFLY wd, %7, %8, %9
72 SBUTTERFLY dq, %1, %3, %9
73 SBUTTERFLY dq, %2, %4, %9
74 SBUTTERFLY dq, %5, %7, %9
75 SBUTTERFLY dq, %6, %8, %9
76 SBUTTERFLY qdq, %1, %5, %9
77 SBUTTERFLY qdq, %2, %6, %9
78 SBUTTERFLY qdq, %3, %7, %9
79 SBUTTERFLY qdq, %4, %8, %9
83 ; in: m0..m7, unless %11 in which case m6 is in %9
84 ; out: m0..m7, unless %11 in which case m4 is in %10
85 ; spills into %9 and %10
89 SBUTTERFLY wd, %1, %2, %7
92 SBUTTERFLY wd, %3, %4, %2
93 SBUTTERFLY wd, %5, %6, %2
94 SBUTTERFLY wd, %7, %8, %2
95 SBUTTERFLY dq, %1, %3, %2
98 SBUTTERFLY dq, %2, %4, %3
99 SBUTTERFLY dq, %5, %7, %3
100 SBUTTERFLY dq, %6, %8, %3
101 SBUTTERFLY qdq, %1, %5, %3
102 SBUTTERFLY qdq, %2, %6, %3
105 SBUTTERFLY qdq, %3, %7, %2
106 SBUTTERFLY qdq, %4, %8, %2
115 %macro ABS1_MMX 2 ; a, tmp
121 %macro ABS2_MMX 4 ; a, b, tmp0, tmp1
139 %define ABS1 ABS1_MMX
140 %define ABS2 ABS2_MMX
148 movd %1, [%2-3] ;to avoid crossing a cacheline
158 %macro SPLATB_SSSE3 3
177 %macro PALIGNR_SSSE3 4
181 %macro DEINTB 5 ; mask, reg1, mask, reg2, optional src to fill masks from
189 pand m%1, m%2 ; dst .. y6 .. y4
190 pand m%3, m%4 ; src .. y6 .. y4
191 psrlw m%2, 8 ; dst .. y7 .. y5
192 psrlw m%4, 8 ; src .. y7 .. y5
207 %macro SUMSUB_BADC 4-5
221 %macro HADAMARD4_V 4+
222 SUMSUB_BADC %1, %2, %3, %4
223 SUMSUB_BADC %1, %3, %2, %4
226 %macro HADAMARD8_V 8+
227 SUMSUB_BADC %1, %2, %3, %4
228 SUMSUB_BADC %5, %6, %7, %8
229 SUMSUB_BADC %1, %3, %2, %4
230 SUMSUB_BADC %5, %7, %6, %8
231 SUMSUB_BADC %1, %5, %2, %6
232 SUMSUB_BADC %3, %7, %4, %8
235 %macro TRANS_SSE2 5-6
237 ; %1: transpose width (d/q) - use SBUTTERFLY qdq for dq
238 ; %2: ord/unord (for compat with sse4, unused)
242 %define mask [mask_10 GLOBAL]
245 %define mask [mask_1100 GLOBAL]
248 %if %0==6 ; less dependency if we have two tmp
249 mova m%5, mask ; ff00
251 psll%1 m%4, shift ; x4..
253 pandn m%5, m%3 ; ..x0
254 psrl%1 m%3, shift ; ..x1
257 %else ; more dependency, one insn less. sometimes faster, sometimes not
259 psll%1 m%4, shift ; x4..
260 pxor m%4, m%3 ; (x4^x1)x0
261 pand m%4, mask ; (x4^x1)..
263 psrl%1 m%4, shift ; ..(x1^x4)
269 %macro TRANS_SSE4 5-6 ; see above
275 pblendw m%3, m%4, 10101010b
278 pblendw m%4, m%5, 01010101b
285 shufps m%3, m%4, 10001000b
286 shufps m%5, m%4, 11011101b
292 ; %1=distance in words (0 for vertical pass, 1/2/4 for horizontal passes)
293 ; %2=sumsub/max/amax (sum and diff / maximum / maximum of absolutes)
296 %if %1!=0 ; have to reorder stuff for horizontal op
299 ; sumsub needs order because a-b != b-a unless a=b
302 ; if we just max, order doesn't matter (allows pblendw+or in sse4)
305 TRANS d, ORDER, %3, %4, %5, %6
308 SBUTTERFLY dq, %3, %4, %5
310 TRANS q, ORDER, %3, %4, %5, %6
313 SBUTTERFLY qdq, %3, %4, %5
317 SUMSUB_BA m%3, m%4, m%5
321 ABS2 m%3, m%4, m%5, m%6
332 %macro HADAMARD2_2D 6-7 sumsub
333 HADAMARD 0, sumsub, %1, %2, %5
334 HADAMARD 0, sumsub, %3, %4, %5
335 SBUTTERFLY %6, %1, %2, %5
337 HADAMARD 0, amax, %1, %2, %5, %7
339 HADAMARD 0, %7, %1, %2, %5
341 SBUTTERFLY %6, %3, %4, %5
343 HADAMARD 0, amax, %3, %4, %5, %7
345 HADAMARD 0, %7, %3, %4, %5
349 %macro HADAMARD4_2D 5-6 sumsub
350 HADAMARD2_2D %1, %2, %3, %4, %5, wd
351 HADAMARD2_2D %1, %3, %2, %4, %5, dq, %6
355 %macro HADAMARD4_2D_SSE 5-6 sumsub
356 HADAMARD 0, sumsub, %1, %2, %5 ; 1st V row 0 + 1
357 HADAMARD 0, sumsub, %3, %4, %5 ; 1st V row 2 + 3
358 SBUTTERFLY wd, %1, %2, %5 ; %1: m0 1+0 %2: m1 1+0
359 SBUTTERFLY wd, %3, %4, %5 ; %3: m0 3+2 %4: m1 3+2
360 HADAMARD2_2D %1, %3, %2, %4, %5, dq
361 SBUTTERFLY qdq, %1, %2, %5
362 HADAMARD 0, %6, %1, %2, %5 ; 2nd H m1/m0 row 0+1
363 SBUTTERFLY qdq, %3, %4, %5
364 HADAMARD 0, %6, %3, %4, %5 ; 2nd H m1/m0 row 2+3
367 %macro HADAMARD8_2D 9-10 sumsub
368 HADAMARD2_2D %1, %2, %3, %4, %9, wd
369 HADAMARD2_2D %5, %6, %7, %8, %9, wd
370 HADAMARD2_2D %1, %3, %2, %4, %9, dq
371 HADAMARD2_2D %5, %7, %6, %8, %9, dq
372 HADAMARD2_2D %1, %5, %3, %7, %9, qdq, %10
373 HADAMARD2_2D %2, %6, %4, %8, %9, qdq, %10
399 psraw %2, 1 ; %2: %2>>1
400 psraw %1, 1 ; %1: %1>>1
401 paddw %2, %4 ; %2: %2>>1+%1
402 psubw %1, %3 ; %1: %1>>1-%2
407 SUMSUB_BADC m%4, m%1, m%3, m%2; m%5
408 SUMSUB_BA m%3, m%4, m%5
409 SUMSUB2_AB m%1, m%2, m%5
410 SWAP %1, %3, %4, %5, %2
412 SUMSUB_BADC m%4, m%1, m%3, m%2
415 SUMSUB2_AB m%1, [%5], m%2
422 SUMSUBD2_AB m%2, m%4, m%6, m%5
423 ; %2: %2>>1-%4 %4: %2+%4>>1
424 SUMSUB_BA m%3, m%1, m%6
425 ; %3: %1+%3 %1: %1-%3
426 SUMSUB_BADC m%4, m%3, m%2, m%1, m%6
427 ; %4: %1+%3 + (%2+%4>>1)
428 ; %3: %1+%3 - (%2+%4>>1)
429 ; %2: %1-%3 + (%2>>1-%4)
430 ; %1: %1-%3 - (%2>>1-%4)
432 SUMSUBD2_AB m%2, m%4, [%5], [%5+16]
434 SUMSUB_BADC m%4, m%3, m%2, m%1
437 ; %1: %1+%3 + (%2+%4>>1) row0
438 ; %2: %1-%3 + (%2>>1-%4) row1
439 ; %3: %1-%3 - (%2>>1-%4) row2
440 ; %4: %1+%3 - (%2+%4>>1) row3
460 %macro LOAD_DIFF8x4_SSE2 8
461 LOAD_DIFF m%1, m%5, m%6, [%7+%1*FENC_STRIDE], [%8+%1*FDEC_STRIDE]
462 LOAD_DIFF m%2, m%5, m%6, [%7+%2*FENC_STRIDE], [%8+%2*FDEC_STRIDE]
463 LOAD_DIFF m%3, m%5, m%6, [%7+%3*FENC_STRIDE], [%8+%3*FDEC_STRIDE]
464 LOAD_DIFF m%4, m%5, m%6, [%7+%4*FENC_STRIDE], [%8+%4*FDEC_STRIDE]
467 %macro LOAD_DIFF8x4_SSSE3 8 ; 4x dst, 1x tmp, 1x mul, 2x ptr
468 movh m%2, [%8+%1*FDEC_STRIDE]
469 movh m%1, [%7+%1*FENC_STRIDE]
471 movh m%3, [%8+%2*FDEC_STRIDE]
472 movh m%2, [%7+%2*FENC_STRIDE]
474 movh m%4, [%8+%3*FDEC_STRIDE]
475 movh m%3, [%7+%3*FENC_STRIDE]
477 movh m%5, [%8+%4*FDEC_STRIDE]
478 movh m%4, [%7+%4*FENC_STRIDE]
491 movhps [%5+%6+32], m%1
492 movhps [%5+%6+40], m%2
493 movhps [%5+%6+48], m%3
494 movhps [%5+%6+56], m%4
498 movhps [r0-4*FDEC_STRIDE], %1
499 movh [r0-3*FDEC_STRIDE], %1
500 movhps [r0-2*FDEC_STRIDE], %2
501 movh [r0-1*FDEC_STRIDE], %2
502 movhps [r0+0*FDEC_STRIDE], %3
503 movh [r0+1*FDEC_STRIDE], %3
504 movhps [r0+2*FDEC_STRIDE], %4
505 movh [r0+3*FDEC_STRIDE], %4
508 %macro LOAD_DIFF_8x4P 7-10 r0,r2,0 ; 4x dest, 2x temp, 2x pointer, increment?
509 LOAD_DIFF m%1, m%5, m%7, [%8], [%9]
510 LOAD_DIFF m%2, m%6, m%7, [%8+r1], [%9+r3]
511 LOAD_DIFF m%3, m%5, m%7, [%8+2*r1], [%9+2*r3]
512 LOAD_DIFF m%4, m%6, m%7, [%8+r4], [%9+r5]