aboutsummaryrefslogtreecommitdiff
path: root/src/libmpg123/synth_x86_64_float.S
diff options
context:
space:
mode:
authorLucas Eckels <eckels@google.com>2012-08-06 15:17:56 -0700
committerLucas Eckels <eckels@google.com>2012-08-08 09:29:17 -0700
commit3d540f5de5b3a28ce6ad855cef7d9d9a44242c07 (patch)
tree56367e03b9d9fe3cce806508c4d7553647945b91 /src/libmpg123/synth_x86_64_float.S
parent9e51f014ada352e89237a6981093860291c8150d (diff)
downloadmpg123-3d540f5de5b3a28ce6ad855cef7d9d9a44242c07.tar.gz
Add mpg123 1.8.1 source
Change-Id: I602e744b1e24d95e203c05b0e93d15d373a45639
Diffstat (limited to 'src/libmpg123/synth_x86_64_float.S')
-rw-r--r--src/libmpg123/synth_x86_64_float.S262
1 files changed, 262 insertions, 0 deletions
diff --git a/src/libmpg123/synth_x86_64_float.S b/src/libmpg123/synth_x86_64_float.S
new file mode 100644
index 0000000..ef8bb02
--- /dev/null
+++ b/src/libmpg123/synth_x86_64_float.S
@@ -0,0 +1,262 @@
+/*
+ synth_x86_64_float: SSE optimized synth for x86-64 (float output version)
+
+ copyright 1995-2009 by the mpg123 project - free software under the terms of the LGPL 2.1
+ see COPYING and AUTHORS files in distribution or http://mpg123.org
+ initially written by Taihei Monma
+*/
+
+#include "mangle.h"
+
+#ifdef _WIN64
+/* short *window; */
+#define ARG0 %r10
+/* short *b0; */
+#define ARG1 %rdx
+/* short *samples; */
+#define ARG2 %r8
+/* int bo1; */
+#define ARG3 %r9
+#else
+/* real *window; */
+#define ARG0 %rdi
+/* real *b0; */
+#define ARG1 %rsi
+/* real *samples; */
+#define ARG2 %rdx
+/* int bo1; */
+#define ARG3 %rcx
+#endif
+
+#define XMMREG_SCALE %xmm15 /* {1/32768.0, 1/32768.0, 1/32768.0, 1/32768.0} */
+
+/*
+ int synth_1to1_real_x86_64_asm(real *window, real *b0, real *samples, int bo1);
+ return value: number of clipped samples (0)
+*/
+
+#ifndef __APPLE__
+ .section .rodata
+#else
+ .data
+#endif
+ ALIGN32
+ASM_NAME(scale_x86_64):
+ .long 939524096
+ .long 939524096
+ .long 939524096
+ .long 939524096
+ .text
+ ALIGN16,,15
+.globl ASM_NAME(synth_1to1_real_x86_64_asm)
+ASM_NAME(synth_1to1_real_x86_64_asm):
+#ifdef _WIN64 /* should save xmm6-15 */
+ movq %rcx, ARG0
+ subq $120, %rsp /* stack alignment + 7 xmm registers */
+ movaps %xmm6, (%rsp)
+ movaps %xmm7, 16(%rsp)
+ movaps %xmm8, 32(%rsp)
+ movaps %xmm9, 48(%rsp)
+ movaps %xmm10, 64(%rsp)
+ movaps %xmm11, 80(%rsp)
+ movaps %xmm15, 96(%rsp)
+#endif
+
+ leaq ASM_NAME(scale_x86_64)(%rip), %rax
+ movaps (%rax), XMMREG_SCALE
+
+ andq $0xf, ARG3
+ shlq $2, ARG3
+ leaq 64(ARG0), ARG0
+ subq ARG3, ARG0
+
+ movl $4, %ecx
+
+ ALIGN16
+Loop_start_1:
+ movups (ARG0), %xmm8
+ movups 16(ARG0), %xmm1
+ movups 32(ARG0), %xmm2
+ movups 48(ARG0), %xmm3
+ movups 128(ARG0), %xmm9
+ movups 144(ARG0), %xmm5
+ movups 160(ARG0), %xmm6
+ movups 176(ARG0), %xmm7
+ mulps (ARG1), %xmm8
+ mulps 16(ARG1), %xmm1
+ mulps 32(ARG1), %xmm2
+ mulps 48(ARG1), %xmm3
+ mulps 64(ARG1), %xmm9
+ mulps 80(ARG1), %xmm5
+ mulps 96(ARG1), %xmm6
+ mulps 112(ARG1), %xmm7
+
+ addps %xmm1, %xmm8
+ addps %xmm2, %xmm3
+ addps %xmm5, %xmm9
+ addps %xmm7, %xmm6
+ addps %xmm3, %xmm8
+ addps %xmm6, %xmm9
+ leaq 256(ARG0), ARG0
+ leaq 128(ARG1), ARG1
+
+ movups (ARG0), %xmm10
+ movups 16(ARG0), %xmm1
+ movups 32(ARG0), %xmm2
+ movups 48(ARG0), %xmm3
+ movups 128(ARG0), %xmm11
+ movups 144(ARG0), %xmm5
+ movups 160(ARG0), %xmm6
+ movups 176(ARG0), %xmm7
+ mulps (ARG1), %xmm10
+ mulps 16(ARG1), %xmm1
+ mulps 32(ARG1), %xmm2
+ mulps 48(ARG1), %xmm3
+ mulps 64(ARG1), %xmm11
+ mulps 80(ARG1), %xmm5
+ mulps 96(ARG1), %xmm6
+ mulps 112(ARG1), %xmm7
+
+ addps %xmm1, %xmm10
+ addps %xmm2, %xmm3
+ addps %xmm5, %xmm11
+ addps %xmm7, %xmm6
+ addps %xmm3, %xmm10
+ addps %xmm6, %xmm11
+ leaq 256(ARG0), ARG0
+ leaq 128(ARG1), ARG1
+
+ movaps %xmm8, %xmm0
+ movaps %xmm10, %xmm1
+ unpcklps %xmm9, %xmm8
+ unpcklps %xmm11, %xmm10
+ unpckhps %xmm9, %xmm0
+ unpckhps %xmm11, %xmm1
+ movaps %xmm8, %xmm2
+ movaps %xmm0, %xmm3
+ movlhps %xmm10, %xmm8
+ movhlps %xmm2, %xmm10
+ movlhps %xmm1, %xmm0
+ movhlps %xmm3, %xmm1
+ subps %xmm10, %xmm8
+ subps %xmm1, %xmm0
+ addps %xmm8, %xmm0
+
+ movups (ARG2), %xmm1
+ movups 16(ARG2), %xmm2
+ mulps XMMREG_SCALE, %xmm0
+ shufps $0xdd, %xmm2, %xmm1
+ movaps %xmm0, %xmm2
+ unpcklps %xmm1, %xmm0
+ unpckhps %xmm1, %xmm2
+ movups %xmm0, (ARG2)
+ movups %xmm2, 16(ARG2)
+
+ leaq 32(ARG2), ARG2
+ decl %ecx
+ jnz Loop_start_1
+
+ movl $4, %ecx
+
+ ALIGN16
+Loop_start_2:
+ movups (ARG0), %xmm8
+ movups 16(ARG0), %xmm1
+ movups 32(ARG0), %xmm2
+ movups 48(ARG0), %xmm3
+ movups 128(ARG0), %xmm9
+ movups 144(ARG0), %xmm5
+ movups 160(ARG0), %xmm6
+ movups 176(ARG0), %xmm7
+ mulps (ARG1), %xmm8
+ mulps 16(ARG1), %xmm1
+ mulps 32(ARG1), %xmm2
+ mulps 48(ARG1), %xmm3
+ mulps -64(ARG1), %xmm9
+ mulps -48(ARG1), %xmm5
+ mulps -32(ARG1), %xmm6
+ mulps -16(ARG1), %xmm7
+
+ addps %xmm1, %xmm8
+ addps %xmm2, %xmm3
+ addps %xmm5, %xmm9
+ addps %xmm7, %xmm6
+ addps %xmm3, %xmm8
+ addps %xmm6, %xmm9
+ leaq 256(ARG0), ARG0
+ leaq -128(ARG1), ARG1
+
+ movups (ARG0), %xmm10
+ movups 16(ARG0), %xmm1
+ movups 32(ARG0), %xmm2
+ movups 48(ARG0), %xmm3
+ movups 128(ARG0), %xmm11
+ movups 144(ARG0), %xmm5
+ movups 160(ARG0), %xmm6
+ movups 176(ARG0), %xmm7
+ mulps (ARG1), %xmm10
+ mulps 16(ARG1), %xmm1
+ mulps 32(ARG1), %xmm2
+ mulps 48(ARG1), %xmm3
+ mulps -64(ARG1), %xmm11
+ mulps -48(ARG1), %xmm5
+ mulps -32(ARG1), %xmm6
+ mulps -16(ARG1), %xmm7
+
+ addps %xmm1, %xmm10
+ addps %xmm2, %xmm3
+ addps %xmm5, %xmm11
+ addps %xmm7, %xmm6
+ addps %xmm3, %xmm10
+ addps %xmm6, %xmm11
+ leaq 256(ARG0), ARG0
+ leaq -128(ARG1), ARG1
+
+ movaps %xmm8, %xmm0
+ movaps %xmm10, %xmm1
+ unpcklps %xmm9, %xmm8
+ unpcklps %xmm11, %xmm10
+ unpckhps %xmm9, %xmm0
+ unpckhps %xmm11, %xmm1
+ movaps %xmm8, %xmm2
+ movaps %xmm0, %xmm3
+ movlhps %xmm10, %xmm8
+ movhlps %xmm2, %xmm10
+ movlhps %xmm1, %xmm0
+ movhlps %xmm3, %xmm1
+ addps %xmm10, %xmm8
+ addps %xmm1, %xmm0
+ addps %xmm8, %xmm0
+
+ movups (ARG2), %xmm1
+ movups 16(ARG2), %xmm2
+ mulps XMMREG_SCALE, %xmm0
+ shufps $0xdd, %xmm2, %xmm1
+ movaps %xmm0, %xmm2
+ unpcklps %xmm1, %xmm0
+ unpckhps %xmm1, %xmm2
+ movups %xmm0, (ARG2)
+ movups %xmm2, 16(ARG2)
+
+ leaq 32(ARG2), ARG2
+ decl %ecx
+ jnz Loop_start_2
+
+ xorl %eax, %eax
+
+#ifdef _WIN64
+ movaps (%rsp), %xmm6
+ movaps 16(%rsp), %xmm7
+ movaps 32(%rsp), %xmm8
+ movaps 48(%rsp), %xmm9
+ movaps 64(%rsp), %xmm10
+ movaps 80(%rsp), %xmm11
+ movaps 96(%rsp), %xmm15
+ addq $120, %rsp
+#endif
+ ret
+
+/* Mark non-executable stack. */
+#if defined(__linux__) && defined(__ELF__)
+.section .note.GNU-stack,"",%progbits
+#endif