2008-08-12 02:26:58 +02:00
|
|
|
;******************************************************************************
|
|
|
|
;* FFT transform with SSE/3DNow optimizations
|
|
|
|
;* Copyright (c) 2008 Loren Merritt
|
|
|
|
;*
|
2010-07-18 22:06:42 +02:00
|
|
|
;* This algorithm (though not any of the implementation details) is
|
|
|
|
;* based on libdjbfft by D. J. Bernstein.
|
|
|
|
;*
|
2008-08-12 02:26:58 +02:00
|
|
|
;* This file is part of FFmpeg.
|
|
|
|
;*
|
|
|
|
;* FFmpeg is free software; you can redistribute it and/or
|
|
|
|
;* modify it under the terms of the GNU Lesser General Public
|
|
|
|
;* License as published by the Free Software Foundation; either
|
|
|
|
;* version 2.1 of the License, or (at your option) any later version.
|
|
|
|
;*
|
|
|
|
;* FFmpeg is distributed in the hope that it will be useful,
|
|
|
|
;* but WITHOUT ANY WARRANTY; without even the implied warranty of
|
|
|
|
;* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
|
|
|
|
;* Lesser General Public License for more details.
|
|
|
|
;*
|
|
|
|
;* You should have received a copy of the GNU Lesser General Public
|
|
|
|
;* License along with FFmpeg; if not, write to the Free Software
|
|
|
|
;* 51, Inc., Foundation Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
|
|
|
|
;******************************************************************************
|
|
|
|
|
|
|
|
; These functions are not individually interchangeable with the C versions.
|
|
|
|
; While C takes arrays of FFTComplex, SSE/3DNow leave intermediate results
|
|
|
|
; in blocks as conventient to the vector size.
|
|
|
|
; i.e. {4x real, 4x imaginary, 4x real, ...} (or 2x respectively)
|
|
|
|
|
|
|
|
%include "x86inc.asm"
|
|
|
|
|
2010-08-22 16:39:58 +02:00
|
|
|
%ifdef ARCH_X86_64
|
|
|
|
%define pointer resq
|
|
|
|
%else
|
|
|
|
%define pointer resd
|
|
|
|
%endif
|
|
|
|
|
|
|
|
struc FFTContext
|
|
|
|
.nbits: resd 1
|
|
|
|
.reverse: resd 1
|
|
|
|
.revtab: pointer 1
|
|
|
|
.tmpbuf: pointer 1
|
|
|
|
.mdctsize: resd 1
|
|
|
|
.mdctbits: resd 1
|
|
|
|
.tcos: pointer 1
|
|
|
|
.tsin: pointer 1
|
|
|
|
endstruc
|
|
|
|
|
2008-08-12 02:26:58 +02:00
|
|
|
SECTION_RODATA
|
|
|
|
|
|
|
|
%define M_SQRT1_2 0.70710678118654752440
|
|
|
|
ps_root2: times 4 dd M_SQRT1_2
|
|
|
|
ps_root2mppm: dd -M_SQRT1_2, M_SQRT1_2, M_SQRT1_2, -M_SQRT1_2
|
|
|
|
ps_m1p1: dd 1<<31, 0
|
|
|
|
|
|
|
|
%assign i 16
|
|
|
|
%rep 13
|
2010-06-23 21:20:46 +02:00
|
|
|
cextern cos_ %+ i
|
2008-08-12 02:26:58 +02:00
|
|
|
%assign i i<<1
|
|
|
|
%endrep
|
|
|
|
|
|
|
|
%ifdef ARCH_X86_64
|
|
|
|
%define pointer dq
|
|
|
|
%else
|
|
|
|
%define pointer dd
|
|
|
|
%endif
|
|
|
|
|
|
|
|
%macro IF0 1+
|
|
|
|
%endmacro
|
|
|
|
%macro IF1 1+
|
|
|
|
%1
|
|
|
|
%endmacro
|
|
|
|
|
|
|
|
section .text align=16
|
|
|
|
|
|
|
|
%macro T2_3DN 4 ; z0, z1, mem0, mem1
|
|
|
|
mova %1, %3
|
|
|
|
mova %2, %1
|
|
|
|
pfadd %1, %4
|
|
|
|
pfsub %2, %4
|
|
|
|
%endmacro
|
|
|
|
|
|
|
|
%macro T4_3DN 6 ; z0, z1, z2, z3, tmp0, tmp1
|
|
|
|
mova %5, %3
|
|
|
|
pfsub %3, %4
|
|
|
|
pfadd %5, %4 ; {t6,t5}
|
2010-06-23 21:20:46 +02:00
|
|
|
pxor %3, [ps_m1p1] ; {t8,t7}
|
2008-08-12 02:26:58 +02:00
|
|
|
mova %6, %1
|
|
|
|
pswapd %3, %3
|
|
|
|
pfadd %1, %5 ; {r0,i0}
|
|
|
|
pfsub %6, %5 ; {r2,i2}
|
|
|
|
mova %4, %2
|
|
|
|
pfadd %2, %3 ; {r1,i1}
|
|
|
|
pfsub %4, %3 ; {r3,i3}
|
|
|
|
SWAP %3, %6
|
|
|
|
%endmacro
|
|
|
|
|
|
|
|
; in: %1={r0,i0,r1,i1} %2={r2,i2,r3,i3}
|
|
|
|
; out: %1={r0,r1,r2,r3} %2={i0,i1,i2,i3}
|
|
|
|
%macro T4_SSE 3
|
|
|
|
mova %3, %1
|
|
|
|
shufps %1, %2, 0x64 ; {r0,i0,r3,i2}
|
|
|
|
shufps %3, %2, 0xce ; {r1,i1,r2,i3}
|
|
|
|
mova %2, %1
|
|
|
|
addps %1, %3 ; {t1,t2,t6,t5}
|
|
|
|
subps %2, %3 ; {t3,t4,t8,t7}
|
|
|
|
mova %3, %1
|
|
|
|
shufps %1, %2, 0x44 ; {t1,t2,t3,t4}
|
|
|
|
shufps %3, %2, 0xbe ; {t6,t5,t7,t8}
|
|
|
|
mova %2, %1
|
|
|
|
addps %1, %3 ; {r0,i0,r1,i1}
|
|
|
|
subps %2, %3 ; {r2,i2,r3,i3}
|
|
|
|
mova %3, %1
|
|
|
|
shufps %1, %2, 0x88 ; {r0,r1,r2,r3}
|
|
|
|
shufps %3, %2, 0xdd ; {i0,i1,i2,i3}
|
|
|
|
SWAP %2, %3
|
|
|
|
%endmacro
|
|
|
|
|
|
|
|
%macro T8_SSE 6 ; r0,i0,r1,i1,t0,t1
|
|
|
|
mova %5, %3
|
|
|
|
shufps %3, %4, 0x44 ; {r4,i4,r6,i6}
|
|
|
|
shufps %5, %4, 0xee ; {r5,i5,r7,i7}
|
|
|
|
mova %6, %3
|
|
|
|
subps %3, %5 ; {r5,i5,r7,i7}
|
|
|
|
addps %6, %5 ; {t1,t2,t3,t4}
|
|
|
|
mova %5, %3
|
|
|
|
shufps %5, %5, 0xb1 ; {i5,r5,i7,r7}
|
2010-06-23 21:20:46 +02:00
|
|
|
mulps %3, [ps_root2mppm] ; {-r5,i5,r7,-i7}
|
|
|
|
mulps %5, [ps_root2]
|
2008-08-12 02:26:58 +02:00
|
|
|
addps %3, %5 ; {t8,t7,ta,t9}
|
|
|
|
mova %5, %6
|
|
|
|
shufps %6, %3, 0x36 ; {t3,t2,t9,t8}
|
|
|
|
shufps %5, %3, 0x9c ; {t1,t4,t7,ta}
|
|
|
|
mova %3, %6
|
|
|
|
addps %6, %5 ; {t1,t2,t9,ta}
|
|
|
|
subps %3, %5 ; {t6,t5,tc,tb}
|
|
|
|
mova %5, %6
|
|
|
|
shufps %6, %3, 0xd8 ; {t1,t9,t5,tb}
|
|
|
|
shufps %5, %3, 0x8d ; {t2,ta,t6,tc}
|
|
|
|
mova %3, %1
|
|
|
|
mova %4, %2
|
|
|
|
addps %1, %6 ; {r0,r1,r2,r3}
|
|
|
|
addps %2, %5 ; {i0,i1,i2,i3}
|
|
|
|
subps %3, %6 ; {r4,r5,r6,r7}
|
|
|
|
subps %4, %5 ; {i4,i5,i6,i7}
|
|
|
|
%endmacro
|
|
|
|
|
|
|
|
; scheduled for cpu-bound sizes
|
|
|
|
%macro PASS_SMALL 3 ; (to load m4-m7), wre, wim
|
|
|
|
IF%1 mova m4, Z(4)
|
|
|
|
IF%1 mova m5, Z(5)
|
|
|
|
mova m0, %2 ; wre
|
|
|
|
mova m2, m4
|
|
|
|
mova m1, %3 ; wim
|
|
|
|
mova m3, m5
|
|
|
|
mulps m2, m0 ; r2*wre
|
2010-09-26 11:08:11 +02:00
|
|
|
IF%1 mova m6, Z2(6)
|
2008-08-12 02:26:58 +02:00
|
|
|
mulps m3, m1 ; i2*wim
|
2010-09-26 11:08:11 +02:00
|
|
|
IF%1 mova m7, Z2(7)
|
2008-08-12 02:26:58 +02:00
|
|
|
mulps m4, m1 ; r2*wim
|
|
|
|
mulps m5, m0 ; i2*wre
|
|
|
|
addps m2, m3 ; r2*wre + i2*wim
|
|
|
|
mova m3, m1
|
|
|
|
mulps m1, m6 ; r3*wim
|
|
|
|
subps m5, m4 ; i2*wre - r2*wim
|
|
|
|
mova m4, m0
|
|
|
|
mulps m3, m7 ; i3*wim
|
|
|
|
mulps m4, m6 ; r3*wre
|
|
|
|
mulps m0, m7 ; i3*wre
|
|
|
|
subps m4, m3 ; r3*wre - i3*wim
|
|
|
|
mova m3, Z(0)
|
|
|
|
addps m0, m1 ; i3*wre + r3*wim
|
|
|
|
mova m1, m4
|
|
|
|
addps m4, m2 ; t5
|
|
|
|
subps m1, m2 ; t3
|
|
|
|
subps m3, m4 ; r2
|
|
|
|
addps m4, Z(0) ; r0
|
|
|
|
mova m6, Z(2)
|
|
|
|
mova Z(4), m3
|
|
|
|
mova Z(0), m4
|
|
|
|
mova m3, m5
|
|
|
|
subps m5, m0 ; t4
|
|
|
|
mova m4, m6
|
|
|
|
subps m6, m5 ; r3
|
|
|
|
addps m5, m4 ; r1
|
2010-09-26 11:08:11 +02:00
|
|
|
mova Z2(6), m6
|
2008-08-12 02:26:58 +02:00
|
|
|
mova Z(2), m5
|
|
|
|
mova m2, Z(3)
|
|
|
|
addps m3, m0 ; t6
|
|
|
|
subps m2, m1 ; i3
|
|
|
|
mova m7, Z(1)
|
|
|
|
addps m1, Z(3) ; i1
|
2010-09-26 11:08:11 +02:00
|
|
|
mova Z2(7), m2
|
2008-08-12 02:26:58 +02:00
|
|
|
mova Z(3), m1
|
|
|
|
mova m4, m7
|
|
|
|
subps m7, m3 ; i2
|
|
|
|
addps m3, m4 ; i0
|
|
|
|
mova Z(5), m7
|
|
|
|
mova Z(1), m3
|
|
|
|
%endmacro
|
|
|
|
|
|
|
|
; scheduled to avoid store->load aliasing
|
|
|
|
%macro PASS_BIG 1 ; (!interleave)
|
|
|
|
mova m4, Z(4) ; r2
|
|
|
|
mova m5, Z(5) ; i2
|
|
|
|
mova m2, m4
|
|
|
|
mova m0, [wq] ; wre
|
|
|
|
mova m3, m5
|
|
|
|
mova m1, [wq+o1q] ; wim
|
|
|
|
mulps m2, m0 ; r2*wre
|
2010-09-26 11:08:11 +02:00
|
|
|
mova m6, Z2(6) ; r3
|
2008-08-12 02:26:58 +02:00
|
|
|
mulps m3, m1 ; i2*wim
|
2010-09-26 11:08:11 +02:00
|
|
|
mova m7, Z2(7) ; i3
|
2008-08-12 02:26:58 +02:00
|
|
|
mulps m4, m1 ; r2*wim
|
|
|
|
mulps m5, m0 ; i2*wre
|
|
|
|
addps m2, m3 ; r2*wre + i2*wim
|
|
|
|
mova m3, m1
|
|
|
|
mulps m1, m6 ; r3*wim
|
|
|
|
subps m5, m4 ; i2*wre - r2*wim
|
|
|
|
mova m4, m0
|
|
|
|
mulps m3, m7 ; i3*wim
|
|
|
|
mulps m4, m6 ; r3*wre
|
|
|
|
mulps m0, m7 ; i3*wre
|
|
|
|
subps m4, m3 ; r3*wre - i3*wim
|
|
|
|
mova m3, Z(0)
|
|
|
|
addps m0, m1 ; i3*wre + r3*wim
|
|
|
|
mova m1, m4
|
|
|
|
addps m4, m2 ; t5
|
|
|
|
subps m1, m2 ; t3
|
|
|
|
subps m3, m4 ; r2
|
|
|
|
addps m4, Z(0) ; r0
|
|
|
|
mova m6, Z(2)
|
|
|
|
mova Z(4), m3
|
|
|
|
mova Z(0), m4
|
|
|
|
mova m3, m5
|
|
|
|
subps m5, m0 ; t4
|
|
|
|
mova m4, m6
|
|
|
|
subps m6, m5 ; r3
|
|
|
|
addps m5, m4 ; r1
|
2010-09-26 11:08:11 +02:00
|
|
|
IF%1 mova Z2(6), m6
|
2008-08-12 02:26:58 +02:00
|
|
|
IF%1 mova Z(2), m5
|
|
|
|
mova m2, Z(3)
|
|
|
|
addps m3, m0 ; t6
|
|
|
|
subps m2, m1 ; i3
|
|
|
|
mova m7, Z(1)
|
|
|
|
addps m1, Z(3) ; i1
|
2010-09-26 11:08:11 +02:00
|
|
|
IF%1 mova Z2(7), m2
|
2008-08-12 02:26:58 +02:00
|
|
|
IF%1 mova Z(3), m1
|
|
|
|
mova m4, m7
|
|
|
|
subps m7, m3 ; i2
|
|
|
|
addps m3, m4 ; i0
|
|
|
|
IF%1 mova Z(5), m7
|
|
|
|
IF%1 mova Z(1), m3
|
|
|
|
%if %1==0
|
|
|
|
mova m4, m5 ; r1
|
|
|
|
mova m0, m6 ; r3
|
|
|
|
unpcklps m5, m1
|
|
|
|
unpckhps m4, m1
|
|
|
|
unpcklps m6, m2
|
|
|
|
unpckhps m0, m2
|
|
|
|
mova m1, Z(0)
|
|
|
|
mova m2, Z(4)
|
|
|
|
mova Z(2), m5
|
|
|
|
mova Z(3), m4
|
2010-09-26 11:08:11 +02:00
|
|
|
mova Z2(6), m6
|
|
|
|
mova Z2(7), m0
|
2008-08-12 02:26:58 +02:00
|
|
|
mova m5, m1 ; r0
|
|
|
|
mova m4, m2 ; r2
|
|
|
|
unpcklps m1, m3
|
|
|
|
unpckhps m5, m3
|
|
|
|
unpcklps m2, m7
|
|
|
|
unpckhps m4, m7
|
|
|
|
mova Z(0), m1
|
|
|
|
mova Z(1), m5
|
|
|
|
mova Z(4), m2
|
|
|
|
mova Z(5), m4
|
|
|
|
%endif
|
|
|
|
%endmacro
|
|
|
|
|
|
|
|
%macro PUNPCK 3
|
|
|
|
mova %3, %1
|
|
|
|
punpckldq %1, %2
|
|
|
|
punpckhdq %3, %2
|
|
|
|
%endmacro
|
|
|
|
|
|
|
|
INIT_XMM
|
2009-10-25 04:09:53 +01:00
|
|
|
%define mova movaps
|
2008-08-12 02:26:58 +02:00
|
|
|
|
|
|
|
%define Z(x) [r0+mmsize*x]
|
2010-09-26 11:08:11 +02:00
|
|
|
%define Z2(x) [r0+mmsize*x]
|
2008-08-12 02:26:58 +02:00
|
|
|
|
|
|
|
align 16
|
|
|
|
fft4_sse:
|
|
|
|
mova m0, Z(0)
|
|
|
|
mova m1, Z(1)
|
|
|
|
T4_SSE m0, m1, m2
|
|
|
|
mova Z(0), m0
|
|
|
|
mova Z(1), m1
|
|
|
|
ret
|
|
|
|
|
|
|
|
align 16
|
|
|
|
fft8_sse:
|
|
|
|
mova m0, Z(0)
|
|
|
|
mova m1, Z(1)
|
|
|
|
T4_SSE m0, m1, m2
|
|
|
|
mova m2, Z(2)
|
|
|
|
mova m3, Z(3)
|
|
|
|
T8_SSE m0, m1, m2, m3, m4, m5
|
|
|
|
mova Z(0), m0
|
|
|
|
mova Z(1), m1
|
|
|
|
mova Z(2), m2
|
|
|
|
mova Z(3), m3
|
|
|
|
ret
|
|
|
|
|
|
|
|
align 16
|
|
|
|
fft16_sse:
|
|
|
|
mova m0, Z(0)
|
|
|
|
mova m1, Z(1)
|
|
|
|
T4_SSE m0, m1, m2
|
|
|
|
mova m2, Z(2)
|
|
|
|
mova m3, Z(3)
|
|
|
|
T8_SSE m0, m1, m2, m3, m4, m5
|
|
|
|
mova m4, Z(4)
|
|
|
|
mova m5, Z(5)
|
|
|
|
mova Z(0), m0
|
|
|
|
mova Z(1), m1
|
|
|
|
mova Z(2), m2
|
|
|
|
mova Z(3), m3
|
|
|
|
T4_SSE m4, m5, m6
|
2010-09-26 11:08:11 +02:00
|
|
|
mova m6, Z2(6)
|
|
|
|
mova m7, Z2(7)
|
2008-08-12 02:26:58 +02:00
|
|
|
T4_SSE m6, m7, m0
|
2010-06-23 21:20:46 +02:00
|
|
|
PASS_SMALL 0, [cos_16], [cos_16+16]
|
2008-08-12 02:26:58 +02:00
|
|
|
ret
|
|
|
|
|
|
|
|
|
|
|
|
INIT_MMX
|
|
|
|
|
|
|
|
%macro FFT48_3DN 1
|
|
|
|
align 16
|
|
|
|
fft4%1:
|
|
|
|
T2_3DN m0, m1, Z(0), Z(1)
|
|
|
|
mova m2, Z(2)
|
|
|
|
mova m3, Z(3)
|
|
|
|
T4_3DN m0, m1, m2, m3, m4, m5
|
|
|
|
PUNPCK m0, m1, m4
|
|
|
|
PUNPCK m2, m3, m5
|
|
|
|
mova Z(0), m0
|
|
|
|
mova Z(1), m4
|
|
|
|
mova Z(2), m2
|
|
|
|
mova Z(3), m5
|
|
|
|
ret
|
|
|
|
|
|
|
|
align 16
|
|
|
|
fft8%1:
|
|
|
|
T2_3DN m0, m1, Z(0), Z(1)
|
|
|
|
mova m2, Z(2)
|
|
|
|
mova m3, Z(3)
|
|
|
|
T4_3DN m0, m1, m2, m3, m4, m5
|
|
|
|
mova Z(0), m0
|
|
|
|
mova Z(2), m2
|
2010-09-26 11:08:11 +02:00
|
|
|
T2_3DN m4, m5, Z(4), Z(5)
|
|
|
|
T2_3DN m6, m7, Z2(6), Z2(7)
|
2008-08-12 02:26:58 +02:00
|
|
|
pswapd m0, m5
|
|
|
|
pswapd m2, m7
|
2010-06-23 21:20:46 +02:00
|
|
|
pxor m0, [ps_m1p1]
|
|
|
|
pxor m2, [ps_m1p1]
|
2008-08-12 02:26:58 +02:00
|
|
|
pfsub m5, m0
|
|
|
|
pfadd m7, m2
|
2010-06-23 21:20:46 +02:00
|
|
|
pfmul m5, [ps_root2]
|
|
|
|
pfmul m7, [ps_root2]
|
2008-08-12 02:26:58 +02:00
|
|
|
T4_3DN m1, m3, m5, m7, m0, m2
|
|
|
|
mova Z(5), m5
|
2010-09-26 11:08:11 +02:00
|
|
|
mova Z2(7), m7
|
2008-08-12 02:26:58 +02:00
|
|
|
mova m0, Z(0)
|
|
|
|
mova m2, Z(2)
|
|
|
|
T4_3DN m0, m2, m4, m6, m5, m7
|
|
|
|
PUNPCK m0, m1, m5
|
|
|
|
PUNPCK m2, m3, m7
|
|
|
|
mova Z(0), m0
|
|
|
|
mova Z(1), m5
|
|
|
|
mova Z(2), m2
|
|
|
|
mova Z(3), m7
|
2010-09-26 11:08:11 +02:00
|
|
|
PUNPCK m4, Z(5), m5
|
|
|
|
PUNPCK m6, Z2(7), m7
|
2008-08-12 02:26:58 +02:00
|
|
|
mova Z(4), m4
|
|
|
|
mova Z(5), m5
|
2010-09-26 11:08:11 +02:00
|
|
|
mova Z2(6), m6
|
|
|
|
mova Z2(7), m7
|
2008-08-12 02:26:58 +02:00
|
|
|
ret
|
|
|
|
%endmacro
|
|
|
|
|
|
|
|
FFT48_3DN _3dn2
|
|
|
|
|
|
|
|
%macro pswapd 2
|
|
|
|
%ifidn %1, %2
|
|
|
|
movd [r0+12], %1
|
|
|
|
punpckhdq %1, [r0+8]
|
|
|
|
%else
|
|
|
|
movq %1, %2
|
|
|
|
psrlq %1, 32
|
|
|
|
punpckldq %1, %2
|
|
|
|
%endif
|
|
|
|
%endmacro
|
|
|
|
|
|
|
|
FFT48_3DN _3dn
|
|
|
|
|
|
|
|
|
2010-09-26 11:08:11 +02:00
|
|
|
%define Z(x) [zq + o1q*(x&6) + mmsize*(x&1)]
|
|
|
|
%define Z2(x) [zq + o3q + mmsize*(x&1)]
|
2008-08-12 02:26:58 +02:00
|
|
|
|
|
|
|
%macro DECL_PASS 2+ ; name, payload
|
|
|
|
align 16
|
|
|
|
%1:
|
|
|
|
DEFINE_ARGS z, w, n, o1, o3
|
|
|
|
lea o3q, [nq*3]
|
|
|
|
lea o1q, [nq*8]
|
|
|
|
shl o3q, 4
|
|
|
|
.loop:
|
|
|
|
%2
|
|
|
|
add zq, mmsize*2
|
|
|
|
add wq, mmsize
|
|
|
|
sub nd, mmsize/8
|
|
|
|
jg .loop
|
|
|
|
rep ret
|
|
|
|
%endmacro
|
|
|
|
|
|
|
|
INIT_XMM
|
2009-10-25 04:09:53 +01:00
|
|
|
%define mova movaps
|
2008-08-12 02:26:58 +02:00
|
|
|
DECL_PASS pass_sse, PASS_BIG 1
|
|
|
|
DECL_PASS pass_interleave_sse, PASS_BIG 0
|
|
|
|
|
|
|
|
INIT_MMX
|
|
|
|
%define mulps pfmul
|
|
|
|
%define addps pfadd
|
|
|
|
%define subps pfsub
|
|
|
|
%define unpcklps punpckldq
|
|
|
|
%define unpckhps punpckhdq
|
|
|
|
DECL_PASS pass_3dn, PASS_SMALL 1, [wq], [wq+o1q]
|
|
|
|
DECL_PASS pass_interleave_3dn, PASS_BIG 0
|
|
|
|
%define pass_3dn2 pass_3dn
|
|
|
|
%define pass_interleave_3dn2 pass_interleave_3dn
|
|
|
|
|
2010-01-30 20:26:47 +01:00
|
|
|
%ifdef PIC
|
|
|
|
%define SECTION_REL - $$
|
|
|
|
%else
|
|
|
|
%define SECTION_REL
|
|
|
|
%endif
|
2008-08-12 02:26:58 +02:00
|
|
|
|
2010-08-22 16:39:58 +02:00
|
|
|
%macro FFT_DISPATCH 2; clobbers 5 GPRs, 8 XMMs
|
|
|
|
lea r2, [dispatch_tab%1]
|
|
|
|
mov r2, [r2 + (%2q-2)*gprsize]
|
|
|
|
%ifdef PIC
|
|
|
|
lea r3, [$$]
|
|
|
|
add r2, r3
|
|
|
|
%endif
|
|
|
|
call r2
|
|
|
|
%endmacro ; FFT_DISPATCH
|
|
|
|
|
2008-08-12 02:26:58 +02:00
|
|
|
%macro DECL_FFT 2-3 ; nbits, cpu, suffix
|
2010-01-30 20:26:47 +01:00
|
|
|
%xdefine list_of_fft fft4%2 SECTION_REL, fft8%2 SECTION_REL
|
2008-08-12 02:26:58 +02:00
|
|
|
%if %1==5
|
2010-01-30 20:26:47 +01:00
|
|
|
%xdefine list_of_fft list_of_fft, fft16%2 SECTION_REL
|
2008-08-12 02:26:58 +02:00
|
|
|
%endif
|
|
|
|
|
|
|
|
%assign n 1<<%1
|
|
|
|
%rep 17-%1
|
|
|
|
%assign n2 n/2
|
|
|
|
%assign n4 n/4
|
2010-01-30 20:26:47 +01:00
|
|
|
%xdefine list_of_fft list_of_fft, fft %+ n %+ %3%2 SECTION_REL
|
2008-08-12 02:26:58 +02:00
|
|
|
|
|
|
|
align 16
|
|
|
|
fft %+ n %+ %3%2:
|
|
|
|
call fft %+ n2 %+ %2
|
|
|
|
add r0, n*4 - (n&(-2<<%1))
|
|
|
|
call fft %+ n4 %+ %2
|
|
|
|
add r0, n*2 - (n2&(-2<<%1))
|
|
|
|
call fft %+ n4 %+ %2
|
|
|
|
sub r0, n*6 + (n2&(-2<<%1))
|
2010-06-23 21:20:46 +02:00
|
|
|
lea r1, [cos_ %+ n]
|
2008-08-12 02:26:58 +02:00
|
|
|
mov r2d, n4/2
|
|
|
|
jmp pass%3%2
|
|
|
|
|
|
|
|
%assign n n*2
|
|
|
|
%endrep
|
|
|
|
%undef n
|
|
|
|
|
|
|
|
align 8
|
|
|
|
dispatch_tab%3%2: pointer list_of_fft
|
|
|
|
|
2009-02-13 01:57:22 +01:00
|
|
|
section .text
|
|
|
|
|
2008-08-12 02:26:58 +02:00
|
|
|
; On x86_32, this function does the register saving and restoring for all of fft.
|
|
|
|
; The others pass args in registers and don't spill anything.
|
2009-08-04 09:42:55 +02:00
|
|
|
cglobal fft_dispatch%3%2, 2,5,8, z, nbits
|
2010-08-22 16:39:58 +02:00
|
|
|
FFT_DISPATCH %3%2, nbits
|
2008-08-12 02:26:58 +02:00
|
|
|
RET
|
|
|
|
%endmacro ; DECL_FFT
|
|
|
|
|
|
|
|
DECL_FFT 5, _sse
|
|
|
|
DECL_FFT 5, _sse, _interleave
|
|
|
|
DECL_FFT 4, _3dn
|
|
|
|
DECL_FFT 4, _3dn, _interleave
|
|
|
|
DECL_FFT 4, _3dn2
|
|
|
|
DECL_FFT 4, _3dn2, _interleave
|
|
|
|
|
2010-08-22 16:39:58 +02:00
|
|
|
INIT_XMM
|
|
|
|
%undef mulps
|
|
|
|
%undef addps
|
|
|
|
%undef subps
|
|
|
|
%undef unpcklps
|
|
|
|
%undef unpckhps
|
|
|
|
|
|
|
|
%macro PREROTATER 5 ;-2*k, 2*k, input+n4, tcos+n8, tsin+n8
|
|
|
|
movaps xmm0, [%3+%2*4]
|
|
|
|
movaps xmm1, [%3+%1*4-0x10]
|
|
|
|
movaps xmm2, xmm0
|
|
|
|
shufps xmm0, xmm1, 0x88
|
|
|
|
shufps xmm1, xmm2, 0x77
|
|
|
|
movlps xmm4, [%4+%2*2]
|
|
|
|
movlps xmm5, [%5+%2*2+0x0]
|
|
|
|
movhps xmm4, [%4+%1*2-0x8]
|
|
|
|
movhps xmm5, [%5+%1*2-0x8]
|
|
|
|
movaps xmm2, xmm0
|
|
|
|
movaps xmm3, xmm1
|
|
|
|
mulps xmm0, xmm5
|
|
|
|
mulps xmm1, xmm4
|
|
|
|
mulps xmm2, xmm4
|
|
|
|
mulps xmm3, xmm5
|
|
|
|
subps xmm1, xmm0
|
|
|
|
addps xmm2, xmm3
|
|
|
|
movaps xmm0, xmm1
|
|
|
|
unpcklps xmm1, xmm2
|
|
|
|
unpckhps xmm0, xmm2
|
|
|
|
%endmacro
|
|
|
|
|
|
|
|
%macro CMUL 6 ;j, xmm0, xmm1, 3, 4, 5
|
|
|
|
movaps xmm6, [%4+%1*2]
|
|
|
|
movaps %2, [%4+%1*2+0x10]
|
|
|
|
movaps %3, xmm6
|
|
|
|
movaps xmm7, %2
|
2010-08-28 23:03:13 +02:00
|
|
|
mulps xmm6, [%5+%1]
|
|
|
|
mulps %2, [%6+%1]
|
|
|
|
mulps %3, [%6+%1]
|
|
|
|
mulps xmm7, [%5+%1]
|
2010-08-22 16:39:58 +02:00
|
|
|
subps %2, xmm6
|
|
|
|
addps %3, xmm7
|
|
|
|
%endmacro
|
|
|
|
|
|
|
|
%macro POSROTATESHUF 5 ;j, k, z+n8, tcos+n8, tsin+n8
|
|
|
|
.post:
|
|
|
|
CMUL %1, xmm0, xmm1, %3, %4, %5
|
|
|
|
CMUL %2, xmm4, xmm5, %3, %4, %5
|
|
|
|
shufps xmm1, xmm1, 0x1b
|
|
|
|
shufps xmm5, xmm5, 0x1b
|
|
|
|
movaps xmm6, xmm4
|
|
|
|
unpckhps xmm4, xmm1
|
|
|
|
unpcklps xmm6, xmm1
|
|
|
|
movaps xmm2, xmm0
|
|
|
|
unpcklps xmm0, xmm5
|
|
|
|
unpckhps xmm2, xmm5
|
|
|
|
movaps [%3+%2*2], xmm6
|
|
|
|
movaps [%3+%2*2+0x10], xmm4
|
|
|
|
movaps [%3+%1*2], xmm0
|
|
|
|
movaps [%3+%1*2+0x10], xmm2
|
|
|
|
sub %2, 0x10
|
|
|
|
add %1, 0x10
|
|
|
|
jl .post
|
|
|
|
%endmacro
|
|
|
|
|
|
|
|
cglobal imdct_half_sse, 3,7,8; FFTContext *s, FFTSample *output, const FFTSample *input
|
|
|
|
%ifdef ARCH_X86_64
|
|
|
|
%define rrevtab r10
|
|
|
|
%define rtcos r11
|
|
|
|
%define rtsin r12
|
|
|
|
push r12
|
|
|
|
push r13
|
|
|
|
push r14
|
|
|
|
%else
|
|
|
|
%define rrevtab r6
|
|
|
|
%define rtsin r6
|
|
|
|
%define rtcos r5
|
|
|
|
%endif
|
|
|
|
mov r3d, [r0+FFTContext.mdctsize]
|
|
|
|
add r2, r3
|
|
|
|
shr r3, 1
|
|
|
|
mov rtcos, [r0+FFTContext.tcos]
|
|
|
|
mov rtsin, [r0+FFTContext.tsin]
|
|
|
|
add rtcos, r3
|
|
|
|
add rtsin, r3
|
|
|
|
%ifndef ARCH_X86_64
|
|
|
|
push rtcos
|
|
|
|
push rtsin
|
|
|
|
%endif
|
|
|
|
shr r3, 1
|
|
|
|
mov rrevtab, [r0+FFTContext.revtab]
|
|
|
|
add rrevtab, r3
|
|
|
|
%ifndef ARCH_X86_64
|
|
|
|
push rrevtab
|
|
|
|
%endif
|
|
|
|
|
|
|
|
sub r3, 4
|
|
|
|
%ifdef ARCH_X86_64
|
|
|
|
xor r4, r4
|
|
|
|
sub r4, r3
|
|
|
|
%endif
|
|
|
|
.pre:
|
|
|
|
%ifndef ARCH_X86_64
|
|
|
|
;unspill
|
|
|
|
xor r4, r4
|
|
|
|
sub r4, r3
|
|
|
|
mov rtsin, [esp+4]
|
|
|
|
mov rtcos, [esp+8]
|
|
|
|
%endif
|
|
|
|
|
|
|
|
PREROTATER r4, r3, r2, rtcos, rtsin
|
|
|
|
%ifdef ARCH_X86_64
|
2010-08-28 23:03:13 +02:00
|
|
|
movzx r5, word [rrevtab+r4-4]
|
|
|
|
movzx r6, word [rrevtab+r4-2]
|
|
|
|
movzx r13, word [rrevtab+r3]
|
|
|
|
movzx r14, word [rrevtab+r3+2]
|
|
|
|
movlps [r1+r5 *8], xmm0
|
|
|
|
movhps [r1+r6 *8], xmm0
|
|
|
|
movlps [r1+r13*8], xmm1
|
|
|
|
movhps [r1+r14*8], xmm1
|
2010-08-22 16:39:58 +02:00
|
|
|
add r4, 4
|
|
|
|
%else
|
|
|
|
mov r6, [esp]
|
2010-08-28 23:03:13 +02:00
|
|
|
movzx r5, word [r6+r4-4]
|
|
|
|
movzx r4, word [r6+r4-2]
|
|
|
|
movlps [r1+r5*8], xmm0
|
|
|
|
movhps [r1+r4*8], xmm0
|
|
|
|
movzx r5, word [r6+r3]
|
|
|
|
movzx r4, word [r6+r3+2]
|
|
|
|
movlps [r1+r5*8], xmm1
|
|
|
|
movhps [r1+r4*8], xmm1
|
2010-08-22 16:39:58 +02:00
|
|
|
%endif
|
|
|
|
sub r3, 4
|
|
|
|
jns .pre
|
|
|
|
|
|
|
|
mov r5, r0
|
|
|
|
mov r6, r1
|
|
|
|
mov r0, r1
|
|
|
|
mov r1d, [r5+FFTContext.nbits]
|
|
|
|
|
|
|
|
FFT_DISPATCH _sse, r1
|
|
|
|
|
|
|
|
mov r0d, [r5+FFTContext.mdctsize]
|
|
|
|
add r6, r0
|
|
|
|
shr r0, 1
|
|
|
|
%ifndef ARCH_X86_64
|
|
|
|
%define rtcos r2
|
|
|
|
%define rtsin r3
|
|
|
|
mov rtcos, [esp+8]
|
|
|
|
mov rtsin, [esp+4]
|
|
|
|
%endif
|
|
|
|
neg r0
|
|
|
|
mov r1, -16
|
|
|
|
sub r1, r0
|
|
|
|
POSROTATESHUF r0, r1, r6, rtcos, rtsin
|
|
|
|
%ifdef ARCH_X86_64
|
|
|
|
pop r14
|
|
|
|
pop r13
|
|
|
|
pop r12
|
|
|
|
%else
|
|
|
|
add esp, 12
|
|
|
|
%endif
|
|
|
|
RET
|