diff mp3lib/decode_i586.c @ 4147:4bbdda22003d

S->C
author nick
date Mon, 14 Jan 2002 09:32:51 +0000
parents
children 3f677202418b
line wrap: on
line diff
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/mp3lib/decode_i586.c	Mon Jan 14 09:32:51 2002 +0000
@@ -0,0 +1,307 @@
+/*
+* mpg123_synth_1to1 works the same way as the c version of this
+* file.  only two types of changes have been made:
+* - reordered floating point instructions to
+*   prevent pipline stalls
+* - made WRITE_SAMPLE use integer instead of
+*   (slower) floating point
+* all kinds of x86 processors should benefit from these
+* modifications.
+*
+* useful sources of information on optimizing x86 code include:
+*
+*     Intel Architecture Optimization Manual
+*     http://www.intel.com/design/pentium/manuals/242816.htm
+*
+*     Cyrix 6x86 Instruction Set Summary
+*     ftp://ftp.cyrix.com/6x86/6x-dbch6.pdf
+*
+*     AMD-K5 Processor Software Development
+*     http://www.amd.com/products/cpg/techdocs/appnotes/20007e.pdf
+*
+* Stefan Bieschewski <stb@acm.org>
+*
+* $Id$
+*/
+#define real float /* ugly - but only way */
+
+static long buffs[1088];
+static long bo=1;
+
+int synth_1to1_pent(real *bandPtr, int channel, short *samples)
+{
+  real tmp[3];
+  register int retval;
+    __asm __volatile(
+"        movl %1,%%eax\n\t"/*bandPtr*/
+"        movl %3,%%esi\n\t"
+"        xorl %%edi,%%edi\n\t"
+"        movl bo,%%ebp\n\t"
+"        cmpl %%edi,%2\n\t"
+"        jne .L48\n\t"
+"        decl %%ebp\n\t"
+"        andl $15,%%ebp\n\t"
+"        movl %%ebp,bo\n\t"
+"        movl $buffs,%%ecx\n\t"
+"        jmp .L49\n\t"
+".L48:\n\t"
+"        addl $2,%%esi\n\t"
+"        movl $buffs+2176,%%ecx\n\t"
+".L49:\n\t"
+"        testl $1,%%ebp\n\t"
+"        je .L50\n\t"
+"        movl %%ecx,%%ebx\n\t"
+"        movl %%ebp,%4\n\t"
+"        pushl %%eax\n\t"
+"        movl 4+%4,%%edx\n\t"
+"        leal (%%ebx,%%edx,4),%%eax\n\t"
+"        pushl %%eax\n\t"
+"        movl 8+%4,%%eax\n\t"
+"        incl %%eax\n\t"
+"        andl $15,%%eax\n\t"
+"        leal 1088(,%%eax,4),%%eax\n\t"
+"        addl %%ebx,%%eax\n\t"
+"        jmp .L74\n\t"
+".L50:\n\t"
+"        leal 1088(%%ecx),%%ebx\n\t"
+"        leal 1(%%ebp),%%edx\n\t"
+"        movl %%edx,%4\n\t"
+"        pushl %%eax\n\t"
+"        leal 1092(%%ecx,%%ebp,4),%%eax\n\t"
+"        pushl %%eax\n\t"
+"        leal (%%ecx,%%ebp,4),%%eax\n\t"
+".L74:\n\t"
+"        pushl %%eax\n\t"
+"        call dct64\n\t"
+"        addl $12,%%esp\n\t"
+"        movl %4,%%edx\n\t"
+"        leal 0(,%%edx,4),%%edx\n\t"
+"        movl $decwin+64,%%eax\n\t"
+"        movl %%eax,%%ecx\n\t"
+"        subl %%edx,%%ecx\n\t"
+"        movl $16,%%ebp\n\t"
+".L55:\n\t"
+"        flds (%%ecx)\n\t"
+"        fmuls (%%ebx)\n\t"
+"        flds 4(%%ecx)\n\t"
+"        fmuls 4(%%ebx)\n\t"
+"        fxch %%st(1)\n\t"
+"        flds 8(%%ecx)\n\t"
+"        fmuls 8(%%ebx)\n\t"
+"        fxch %%st(2)\n\t"
+"        fsubrp %%st,%%st(1)\n\t"
+"        flds 12(%%ecx)\n\t"
+"        fmuls 12(%%ebx)\n\t"
+"        fxch %%st(2)\n\t"
+"        faddp %%st,%%st(1)\n\t"
+"        flds 16(%%ecx)\n\t"
+"        fmuls 16(%%ebx)\n\t"
+"        fxch %%st(2)\n\t"
+"        fsubrp %%st,%%st(1)\n\t"
+"        flds 20(%%ecx)\n\t"
+"        fmuls 20(%%ebx)\n\t"
+"        fxch %%st(2)\n\t"
+"        faddp %%st,%%st(1)\n\t"
+"        flds 24(%%ecx)\n\t"
+"        fmuls 24(%%ebx)\n\t"
+"        fxch %%st(2)\n\t"
+"        fsubrp %%st,%%st(1)\n\t"
+"        flds 28(%%ecx)\n\t"
+"        fmuls 28(%%ebx)\n\t"
+"        fxch %%st(2)\n\t"
+"        faddp %%st,%%st(1)\n\t"
+"        flds 32(%%ecx)\n\t"
+"        fmuls 32(%%ebx)\n\t"
+"        fxch %%st(2)\n\t"
+"        fsubrp %%st,%%st(1)\n\t"
+"        flds 36(%%ecx)\n\t"
+"        fmuls 36(%%ebx)\n\t"
+"        fxch %%st(2)\n\t"
+"        faddp %%st,%%st(1)\n\t"
+"        flds 40(%%ecx)\n\t"
+"        fmuls 40(%%ebx)\n\t"
+"        fxch %%st(2)\n\t"
+"        fsubrp %%st,%%st(1)\n\t"
+"        flds 44(%%ecx)\n\t"
+"        fmuls 44(%%ebx)\n\t"
+"        fxch %%st(2)\n\t"
+"        faddp %%st,%%st(1)\n\t"
+"        flds 48(%%ecx)\n\t"
+"        fmuls 48(%%ebx)\n\t"
+"        fxch %%st(2)\n\t"
+"        fsubrp %%st,%%st(1)\n\t"
+"        flds 52(%%ecx)\n\t"
+"        fmuls 52(%%ebx)\n\t"
+"        fxch %%st(2)\n\t"
+"        faddp %%st,%%st(1)\n\t"
+"        flds 56(%%ecx)\n\t"
+"        fmuls 56(%%ebx)\n\t"
+"        fxch %%st(2)\n\t"
+"        fsubrp %%st,%%st(1)\n\t"
+"        flds 60(%%ecx)\n\t"
+"        fmuls 60(%%ebx)\n\t"
+"        fxch %%st(2)\n\t"
+"        subl $4,%%esp\n\t"
+"        faddp %%st,%%st(1)\n\t"
+"        fxch %%st(1)\n\t"
+"        fsubrp %%st,%%st(1)\n\t"
+"        fistpl (%%esp)\n\t"
+"        popl %%eax\n\t"
+"        cmpl $32767,%%eax\n\t"
+"        jg 1f\n\t"
+"        cmpl $-32768,%%eax\n\t"
+"        jl 2f\n\t"
+"        movw %%ax,(%%esi)\n\t"
+"        jmp 4f\n\t"
+"1:      movw $32767,(%%esi)\n\t"
+"        jmp 3f\n\t"
+"2:      movw $-32768,(%%esi)\n\t"
+"3:      incl %%edi\n\t"
+"4:\n\t"
+".L54:\n\t"
+"        addl $64,%%ebx\n\t"
+"        subl $-128,%%ecx\n\t"
+"        addl $4,%%esi\n\t"
+"        decl %%ebp\n\t"
+"        jnz .L55\n\t"
+"        flds (%%ecx)\n\t"
+"        fmuls (%%ebx)\n\t"
+"        flds 8(%%ecx)\n\t"
+"        fmuls 8(%%ebx)\n\t"
+"        flds 16(%%ecx)\n\t"
+"        fmuls 16(%%ebx)\n\t"
+"        fxch %%st(2)\n\t"
+"        faddp %%st,%%st(1)\n\t"
+"        flds 24(%%ecx)\n\t"
+"        fmuls 24(%%ebx)\n\t"
+"        fxch %%st(2)\n\t"
+"        faddp %%st,%%st(1)\n\t"
+"        flds 32(%%ecx)\n\t"
+"        fmuls 32(%%ebx)\n\t"
+"        fxch %%st(2)\n\t"
+"        faddp %%st,%%st(1)\n\t"
+"        flds 40(%%ecx)\n\t"
+"        fmuls 40(%%ebx)\n\t"
+"        fxch %%st(2)\n\t"
+"        faddp %%st,%%st(1)\n\t"
+"        flds 48(%%ecx)\n\t"
+"        fmuls 48(%%ebx)\n\t"
+"        fxch %%st(2)\n\t"
+"        faddp %%st,%%st(1)\n\t"
+"        flds 56(%%ecx)\n\t"
+"        fmuls 56(%%ebx)\n\t"
+"        fxch %%st(2)\n\t"
+"        subl $4,%%esp\n\t"
+"        faddp %%st,%%st(1)\n\t"
+"        fxch %%st(1)\n\t"
+"        faddp %%st,%%st(1)\n\t"
+"        fistpl (%%esp)\n\t"
+"        popl %%eax\n\t"
+"        cmpl $32767,%%eax\n\t"
+"        jg 1f\n\t"
+"        cmpl $-32768,%%eax\n\t"
+"        jl 2f\n\t"
+"        movw %%ax,(%%esi)\n\t"
+"        jmp 4f\n\t"
+"1:      movw $32767,(%%esi)\n\t"
+"        jmp 3f\n\t"
+"2:      movw $-32768,(%%esi)\n\t"
+"3:      incl %%edi\n\t"
+"4:\n\t"
+".L62:\n\t"
+"        addl $-64,%%ebx\n\t"
+"        addl $4,%%esi\n\t"
+"        movl %4,%%edx\n\t"
+"        leal -128(%%ecx,%%edx,8),%%ecx\n\t"
+"        movl $15,%%ebp\n\t"
+".L68:\n\t"
+"        flds -4(%%ecx)\n\t"
+"        fchs\n\t"
+"        fmuls (%%ebx)\n\t"
+"        flds -8(%%ecx)\n\t"
+"        fmuls 4(%%ebx)\n\t"
+"        fxch %%st(1)\n\t"
+"        flds -12(%%ecx)\n\t"
+"        fmuls 8(%%ebx)\n\t"
+"        fxch %%st(2)\n\t"
+"        fsubrp %%st,%%st(1)\n\t"
+"        flds -16(%%ecx)\n\t"
+"        fmuls 12(%%ebx)\n\t"
+"        fxch %%st(2)\n\t"
+"        fsubrp %%st,%%st(1)\n\t"
+"        flds -20(%%ecx)\n\t"
+"        fmuls 16(%%ebx)\n\t"
+"        fxch %%st(2)\n\t"
+"        fsubrp %%st,%%st(1)\n\t"
+"        flds -24(%%ecx)\n\t"
+"        fmuls 20(%%ebx)\n\t"
+"        fxch %%st(2)\n\t"
+"        fsubrp %%st,%%st(1)\n\t"
+"        flds -28(%%ecx)\n\t"
+"        fmuls 24(%%ebx)\n\t"
+"        fxch %%st(2)\n\t"
+"        fsubrp %%st,%%st(1)\n\t"
+"        flds -32(%%ecx)\n\t"
+"        fmuls 28(%%ebx)\n\t"
+"        fxch %%st(2)\n\t"
+"        fsubrp %%st,%%st(1)\n\t"
+"        flds -36(%%ecx)\n\t"
+"        fmuls 32(%%ebx)\n\t"
+"        fxch %%st(2)\n\t"
+"        fsubrp %%st,%%st(1)\n\t"
+"        flds -40(%%ecx)\n\t"
+"        fmuls 36(%%ebx)\n\t"
+"        fxch %%st(2)\n\t"
+"        fsubrp %%st,%%st(1)\n\t"
+"        flds -44(%%ecx)\n\t"
+"        fmuls 40(%%ebx)\n\t"
+"        fxch %%st(2)\n\t"
+"        fsubrp %%st,%%st(1)\n\t"
+"        flds -48(%%ecx)\n\t"
+"        fmuls 44(%%ebx)\n\t"
+"        fxch %%st(2)\n\t"
+"        fsubrp %%st,%%st(1)\n\t"
+"        flds -52(%%ecx)\n\t"
+"        fmuls 48(%%ebx)\n\t"
+"        fxch %%st(2)\n\t"
+"        fsubrp %%st,%%st(1)\n\t"
+"        flds -56(%%ecx)\n\t"
+"        fmuls 52(%%ebx)\n\t"
+"        fxch %%st(2)\n\t"
+"        fsubrp %%st,%%st(1)\n\t"
+"        flds -60(%%ecx)\n\t"
+"        fmuls 56(%%ebx)\n\t"
+"        fxch %%st(2)\n\t"
+"        fsubrp %%st,%%st(1)\n\t"
+"        flds (%%ecx)\n\t"
+"        fmuls 60(%%ebx)\n\t"
+"        fxch %%st(2)\n\t"
+"        subl $4,%%esp\n\t"
+"        fsubrp %%st,%%st(1)\n\t"
+"        fxch %%st(1)\n\t"
+"        fsubrp %%st,%%st(1)\n\t"
+"        fistpl (%%esp)\n\t"
+"        popl %%eax\n\t"
+"        cmpl $32767,%%eax\n\t"
+"        jg 1f\n\t"
+"        cmpl $-32768,%%eax\n\t"
+"        jl 2f\n\t"
+"        movw %%ax,(%%esi)\n\t"
+"        jmp 4f\n\t"
+"1:      movw $32767,(%%esi)\n\t"
+"        jmp 3f\n\t"
+"2:      movw $-32768,(%%esi)\n\t"
+"3:      incl %%edi\n\t"
+"4:\n\t"
+".L67:\n\t"
+"        addl $-64,%%ebx\n\t"
+"        addl $-128,%%ecx\n\t"
+"        addl $4,%%esi\n\t"
+"        decl %%ebp\n\t"
+"        jnz .L68\n\t"
+"        movl %%edi,%%eax\n\t"
+	:"=a"(retval)
+	:"m"(bandPtr),"m"(channel),"m"(samples),"m"(tmp[0])
+	:"memory","%ebp","%edi","%esi","%ebx");
+  return retval;
+}