Mercurial > mplayer.hg

changeset 885:cec1562ccf8a

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Improvements
author nickols_k
date Sat, 26 May 2001 10:29:41 +0000
parents 83b5359aee63
children cb432deedb92
files liba52/srfftp_3dnow.h
diffstat 1 files changed, 102 insertions(+), 0 deletions(-) [+]
line wrap: on
line diff
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/liba52/srfftp_3dnow.h	Sat May 26 10:29:41 2001 +0000
@@ -0,0 +1,102 @@
+/* 
+ *  srfftp.h
+ *
+ *  Copyright (C) Yuqing Deng <Yuqing_Deng@brown.edu> - April 2000
+ *
+ *  64 and 128 point split radix fft for ac3dec
+ *
+ *  The algorithm is desribed in the book:
+ *  "Computational Frameworks of the Fast Fourier Transform".
+ *
+ *  The ideas and the the organization of code borrowed from djbfft written by
+ *  D. J. Bernstein <djb@cr.py.to>.  djbff can be found at 
+ *  http://cr.yp.to/djbfft.html.
+ *
+ *  srfftp.h is free software; you can redistribute it and/or modify
+ *  it under the terms of the GNU General Public License as published by
+ *  the Free Software Foundation; either version 2, or (at your option)
+ *  any later version.
+ *
+ *  srfftp.h 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 General Public License for more details.
+ *
+ *  You should have received a copy of the GNU General Public License
+ *  along with GNU Make; see the file COPYING.  If not, write to
+ *  the Free Software Foundation, 675 Mass Ave, Cambridge, MA 02139, USA.
+ *
+ *  Modified for using AMD's 3DNow! - 3DNowEx(DSP)! SIMD operations 
+ *  by Nick Kurshev <nickols_k@mail.ru>
+ */
+
+#ifndef SRFFTP_3DNOW_H__
+#define SRFFTP_3DNOW_H__
+
+#ifdef HAVE_3DNOWEX
+#define TRANS_FILL_MM6_MM7_3DNOW()\
+    asm(\
+	"movl	$-1, %%eax\n\t"\
+	"movd	%%eax, %%mm7\n\t"\
+	"negl	%%eax\n\t"\
+	"movd	%%eax, %%mm6\n\t"\
+	"punpckldq %%mm6, %%mm7\n\t" /* -1.0 | 1.0 */\
+	"pi2fd	%%mm7, %%mm7\n\t"\
+	"pswapd	%%mm7, %%mm6\n\t"/* 1.0 | -1.0 */\
+	:::"eax","memory");
+#else
+#define TRANS_FILL_MM6_MM7_3DNOW()\
+    asm(\
+	"movl	$-1, %%eax\n\t"\
+	"movd	%%eax, %%mm7\n\t"\
+	"negl	%%eax\n\t"\
+	"movd	%%eax, %%mm6\n\t"\
+	"punpckldq %%mm6, %%mm7\n\t" /* -1.0 | 1.0 */\
+	"punpckldq %%mm7, %%mm6\n\t" /* 1.0 | -1.0 */\
+	"pi2fd	%%mm7, %%mm7\n\t"\
+	"pi2fd	%%mm6, %%mm6\n\t"\
+	:::"eax","memory");
+#endif
+
+#ifdef HAVE_3DNOWEX
+#define PSWAP_MM(mm_base,mm_hlp) "pswapd	"##mm_base","##mm_base" \n\t"
+#else
+#define PSWAP_MM(mm_base,mm_hlp)\
+	"movq	"##mm_base","##mm_hlp" \n\t"\
+	"psrlq $32, "##mm_base"\n\t"\
+	"punpckldq "##mm_hlp","##mm_base"\n\t"
+#endif
+
+#define TRANSZERO_3DNOW(A0,A4,A8,A12) \
+{ \
+    asm volatile("femms":::"memory");\
+    TRANS_FILL_MM6_MM7_3DNOW()\
+    asm(\
+	"movq	%4, %%mm0\n\t" /* mm0 = wTB[0]*/\
+	"movq	%5, %%mm1\n\t" /* mm1 = wTB[k*2]*/ \
+	"movq	%%mm0, %%mm5\n\t"/*u.re = wTB[0].re + wTB[k*2].re;*/\
+	"pfadd	%%mm1, %%mm5\n\t"/*u.im = wTB[0].im + wTB[k*2].im; mm5 = u*/\
+	"pfmul  %%mm6, %%mm0\n\t"/*mm0 = wTB[0].re | -wTB[0].im */\
+	"pfmul	%%mm7, %%mm1\n\t"/*mm1 = -wTB[k*2].re | wTB[k*2].im */\
+	"pfadd	%%mm1, %%mm0\n\t"/*v.im = wTB[0].re - wTB[k*2].re;*/\
+	"movq	%%mm0, %%mm4\n\t"/*v.re =-wTB[0].im + wTB[k*2].im;*/\
+	PSWAP_MM("%%mm4","%%mm2")/* mm4 = v*/\
+	"movq	%6, %%mm0\n\t" /* a1 = A0;*/\
+	"movq	%%mm0, %%mm1\n\t"\
+	"pfadd	%%mm5, %%mm0\n\t" /*A0 = a1 + u;*/\
+	"pfsub	%%mm5, %%mm1\n\t" /*A1 = a1 - u;*/\
+	"movq	%%mm0, %0\n\t"\
+	"movq	%%mm1, %1\n\t"\
+	"movq	%7, %%mm2\n\t" /* a1 = A4;*/\
+	"movq	%%mm2, %%mm3\n\t"\
+	"pfadd	%%mm4, %%mm2\n\t" /*A12 = a1 + v;*/\
+	"pfsub	%%mm4, %%mm3\n\t" /*A4  = a1 - v;*/\
+	"movq	%%mm2, %3\n\t"\
+	"movq	%%mm3, %2"\
+	:"=m"(A0), "=m"(A8), "=m"(A4), "=m"(A12)\
+	:"m"(wTB[0]), "m"(wTB[k*2]), "0"(A0), "2"(A4)\
+	:"memory");\
+    asm volatile("femms":::"memory");\
+}
+
+#endif