mirror of
https://github.com/KolibriOS/kolibrios.git
synced 2024-12-12 01:53:58 +03:00
402 lines
8.9 KiB
ArmAsm
402 lines
8.9 KiB
ArmAsm
/*
|
|
dct64_sse_float: SSE optimized dct64 (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"
|
|
|
|
#define ARG(n) (8+n*4)(%ebp)
|
|
#define TEMP(n) (4+n*16)(%esp)
|
|
#define TEMP_BYTE(n) (4+n)(%esp)
|
|
|
|
/*
|
|
void dct64_real_sse(real *out0, real *out1, real *samples);
|
|
*/
|
|
|
|
#ifndef __APPLE__
|
|
.section .rodata
|
|
#else
|
|
.data
|
|
#endif
|
|
ALIGN16
|
|
pnpn:
|
|
.long 0
|
|
.long -2147483648
|
|
.long 0
|
|
.long -2147483648
|
|
ALIGN16
|
|
mask:
|
|
.long -1
|
|
.long -1
|
|
.long -1
|
|
.long 0
|
|
|
|
.text
|
|
ALIGN16
|
|
.globl ASM_NAME(dct64_real_sse)
|
|
ASM_NAME(dct64_real_sse):
|
|
pushl %ebp
|
|
movl %esp, %ebp
|
|
|
|
andl $-16, %esp /* align the stack at 16 bytes */
|
|
subl $128, %esp /* reserve space for temporal store */
|
|
pushl %ebx
|
|
|
|
movl ARG(0), %ecx
|
|
movl ARG(1), %ebx
|
|
movl ARG(2), %eax
|
|
|
|
MOVUAPS (%eax), %xmm7
|
|
MOVUAPS 16(%eax), %xmm6
|
|
MOVUAPS 112(%eax), %xmm0
|
|
MOVUAPS 96(%eax), %xmm1
|
|
shufps $0x1b, %xmm0, %xmm0
|
|
shufps $0x1b, %xmm1, %xmm1
|
|
movaps %xmm7, %xmm4
|
|
movaps %xmm6, %xmm5
|
|
addps %xmm0, %xmm4
|
|
addps %xmm1, %xmm5
|
|
subps %xmm0, %xmm7
|
|
subps %xmm1, %xmm6
|
|
movaps %xmm4, TEMP(0)
|
|
movaps %xmm5, TEMP(1)
|
|
|
|
MOVUAPS 32(%eax), %xmm2
|
|
MOVUAPS 48(%eax), %xmm3
|
|
MOVUAPS 80(%eax), %xmm0
|
|
MOVUAPS 64(%eax), %xmm1
|
|
shufps $0x1b, %xmm0, %xmm0
|
|
shufps $0x1b, %xmm1, %xmm1
|
|
movaps %xmm2, %xmm5
|
|
movaps %xmm3, %xmm4
|
|
addps %xmm0, %xmm2
|
|
addps %xmm1, %xmm3
|
|
subps %xmm0, %xmm5
|
|
subps %xmm1, %xmm4
|
|
|
|
mulps ASM_NAME(costab_mmxsse), %xmm7
|
|
mulps ASM_NAME(costab_mmxsse)+16, %xmm6
|
|
mulps ASM_NAME(costab_mmxsse)+32, %xmm5
|
|
mulps ASM_NAME(costab_mmxsse)+48, %xmm4
|
|
|
|
shufps $0x1b, %xmm2, %xmm2
|
|
shufps $0x1b, %xmm3, %xmm3
|
|
shufps $0x1b, %xmm4, %xmm4
|
|
shufps $0x1b, %xmm5, %xmm5
|
|
movaps TEMP(0), %xmm0
|
|
movaps TEMP(1), %xmm1
|
|
subps %xmm3, %xmm0
|
|
subps %xmm2, %xmm1
|
|
addps TEMP(0), %xmm3
|
|
addps TEMP(1), %xmm2
|
|
movaps %xmm3, TEMP(0)
|
|
movaps %xmm2, TEMP(1)
|
|
movaps %xmm6, %xmm2
|
|
movaps %xmm7, %xmm3
|
|
subps %xmm5, %xmm6
|
|
subps %xmm4, %xmm7
|
|
addps %xmm3, %xmm4
|
|
addps %xmm2, %xmm5
|
|
mulps ASM_NAME(costab_mmxsse)+64, %xmm0
|
|
mulps ASM_NAME(costab_mmxsse)+80, %xmm1
|
|
mulps ASM_NAME(costab_mmxsse)+80, %xmm6
|
|
mulps ASM_NAME(costab_mmxsse)+64, %xmm7
|
|
|
|
movaps TEMP(0), %xmm2
|
|
movaps TEMP(1), %xmm3
|
|
shufps $0x1b, %xmm3, %xmm3
|
|
shufps $0x1b, %xmm5, %xmm5
|
|
shufps $0x1b, %xmm1, %xmm1
|
|
shufps $0x1b, %xmm6, %xmm6
|
|
movaps %xmm0, TEMP(1)
|
|
subps %xmm3, %xmm2
|
|
subps %xmm1, %xmm0
|
|
addps TEMP(0), %xmm3
|
|
addps TEMP(1), %xmm1
|
|
movaps %xmm3, TEMP(0)
|
|
movaps %xmm1, TEMP(2)
|
|
movaps %xmm5, %xmm1
|
|
movaps %xmm4, %xmm5
|
|
movaps %xmm7, %xmm3
|
|
subps %xmm1, %xmm5
|
|
subps %xmm6, %xmm7
|
|
addps %xmm1, %xmm4
|
|
addps %xmm3, %xmm6
|
|
mulps ASM_NAME(costab_mmxsse)+96, %xmm2
|
|
mulps ASM_NAME(costab_mmxsse)+96, %xmm0
|
|
mulps ASM_NAME(costab_mmxsse)+96, %xmm5
|
|
mulps ASM_NAME(costab_mmxsse)+96, %xmm7
|
|
movaps %xmm2, TEMP(1)
|
|
movaps %xmm0, TEMP(3)
|
|
|
|
movaps %xmm4, %xmm2
|
|
movaps %xmm5, %xmm3
|
|
shufps $0x44, %xmm6, %xmm2
|
|
shufps $0xbb, %xmm7, %xmm5
|
|
shufps $0xbb, %xmm6, %xmm4
|
|
shufps $0x44, %xmm7, %xmm3
|
|
movaps %xmm2, %xmm6
|
|
movaps %xmm3, %xmm7
|
|
subps %xmm4, %xmm2
|
|
subps %xmm5, %xmm3
|
|
addps %xmm6, %xmm4
|
|
addps %xmm7, %xmm5
|
|
movaps ASM_NAME(costab_mmxsse)+112, %xmm0
|
|
movlhps %xmm0, %xmm0
|
|
mulps %xmm0, %xmm2
|
|
mulps %xmm0, %xmm3
|
|
movaps %xmm0, TEMP(4)
|
|
movaps %xmm4, %xmm6
|
|
movaps %xmm5, %xmm7
|
|
shufps $0x14, %xmm2, %xmm4
|
|
shufps $0xbe, %xmm2, %xmm6
|
|
shufps $0x14, %xmm3, %xmm5
|
|
shufps $0xbe, %xmm3, %xmm7
|
|
movaps %xmm5, TEMP(5)
|
|
movaps %xmm7, TEMP(7)
|
|
|
|
movaps TEMP(0), %xmm0
|
|
movaps TEMP(1), %xmm1
|
|
movaps %xmm0, %xmm2
|
|
movaps %xmm1, %xmm3
|
|
shufps $0x44, TEMP(2), %xmm2
|
|
shufps $0xbb, TEMP(3), %xmm1
|
|
shufps $0xbb, TEMP(2), %xmm0
|
|
shufps $0x44, TEMP(3), %xmm3
|
|
movaps %xmm2, %xmm5
|
|
movaps %xmm3, %xmm7
|
|
subps %xmm0, %xmm2
|
|
subps %xmm1, %xmm3
|
|
addps %xmm5, %xmm0
|
|
addps %xmm7, %xmm1
|
|
mulps TEMP(4), %xmm2
|
|
mulps TEMP(4), %xmm3
|
|
movaps %xmm0, %xmm5
|
|
movaps %xmm1, %xmm7
|
|
shufps $0x14, %xmm2, %xmm0
|
|
shufps $0xbe, %xmm2, %xmm5
|
|
shufps $0x14, %xmm3, %xmm1
|
|
shufps $0xbe, %xmm3, %xmm7
|
|
|
|
movaps %xmm0, TEMP(0)
|
|
movaps %xmm1, TEMP(1)
|
|
movaps %xmm5, TEMP(2)
|
|
movaps %xmm7, TEMP(3)
|
|
|
|
movss ASM_NAME(costab_mmxsse)+120, %xmm5
|
|
shufps $0x00, %xmm5, %xmm5
|
|
xorps pnpn, %xmm5
|
|
|
|
movaps %xmm4, %xmm0
|
|
movaps %xmm6, %xmm1
|
|
unpcklps TEMP(5), %xmm4
|
|
unpckhps TEMP(5), %xmm0
|
|
unpcklps TEMP(7), %xmm6
|
|
unpckhps TEMP(7), %xmm1
|
|
movaps %xmm4, %xmm2
|
|
movaps %xmm6, %xmm3
|
|
unpcklps %xmm0, %xmm4
|
|
unpckhps %xmm0, %xmm2
|
|
unpcklps %xmm1, %xmm6
|
|
unpckhps %xmm1, %xmm3
|
|
movaps %xmm4, %xmm0
|
|
movaps %xmm6, %xmm1
|
|
subps %xmm2, %xmm0
|
|
subps %xmm3, %xmm1
|
|
addps %xmm2, %xmm4
|
|
addps %xmm3, %xmm6
|
|
mulps %xmm5, %xmm0
|
|
mulps %xmm5, %xmm1
|
|
movaps %xmm5, TEMP(5)
|
|
movaps %xmm4, %xmm5
|
|
movaps %xmm6, %xmm7
|
|
unpcklps %xmm0, %xmm4
|
|
unpckhps %xmm0, %xmm5
|
|
unpcklps %xmm1, %xmm6
|
|
unpckhps %xmm1, %xmm7
|
|
|
|
movaps TEMP(0), %xmm0
|
|
movaps TEMP(2), %xmm2
|
|
movaps %xmm4, TEMP(4)
|
|
movaps %xmm6, TEMP(6)
|
|
|
|
movaps %xmm0, %xmm4
|
|
movaps %xmm2, %xmm6
|
|
unpcklps TEMP(1), %xmm0
|
|
unpckhps TEMP(1), %xmm4
|
|
unpcklps TEMP(3), %xmm2
|
|
unpckhps TEMP(3), %xmm6
|
|
movaps %xmm0, %xmm1
|
|
movaps %xmm2, %xmm3
|
|
unpcklps %xmm4, %xmm0
|
|
unpckhps %xmm4, %xmm1
|
|
unpcklps %xmm6, %xmm2
|
|
unpckhps %xmm6, %xmm3
|
|
movaps %xmm0, %xmm4
|
|
movaps %xmm2, %xmm6
|
|
subps %xmm1, %xmm4
|
|
subps %xmm3, %xmm6
|
|
addps %xmm1, %xmm0
|
|
addps %xmm3, %xmm2
|
|
mulps TEMP(5), %xmm4
|
|
mulps TEMP(5), %xmm6
|
|
movaps %xmm0, %xmm1
|
|
movaps %xmm2, %xmm3
|
|
unpcklps %xmm4, %xmm0
|
|
unpckhps %xmm4, %xmm1
|
|
unpcklps %xmm6, %xmm2
|
|
unpckhps %xmm6, %xmm3
|
|
|
|
movaps %xmm0, TEMP(0)
|
|
movaps %xmm1, TEMP(1)
|
|
movaps %xmm2, TEMP(2)
|
|
movaps %xmm3, TEMP(3)
|
|
movaps %xmm5, TEMP(5)
|
|
movaps %xmm7, TEMP(7)
|
|
|
|
movss TEMP_BYTE(12), %xmm0
|
|
movss TEMP_BYTE(28), %xmm1
|
|
movss TEMP_BYTE(44), %xmm2
|
|
movss TEMP_BYTE(60), %xmm3
|
|
addss TEMP_BYTE(8), %xmm0
|
|
addss TEMP_BYTE(24), %xmm1
|
|
addss TEMP_BYTE(40), %xmm2
|
|
addss TEMP_BYTE(56), %xmm3
|
|
movss %xmm0, TEMP_BYTE(8)
|
|
movss %xmm1, TEMP_BYTE(24)
|
|
movss %xmm2, TEMP_BYTE(40)
|
|
movss %xmm3, TEMP_BYTE(56)
|
|
movss TEMP_BYTE(76), %xmm0
|
|
movss TEMP_BYTE(92), %xmm1
|
|
movss TEMP_BYTE(108), %xmm2
|
|
movss TEMP_BYTE(124), %xmm3
|
|
addss TEMP_BYTE(72), %xmm0
|
|
addss TEMP_BYTE(88), %xmm1
|
|
addss TEMP_BYTE(104), %xmm2
|
|
addss TEMP_BYTE(120), %xmm3
|
|
movss %xmm0, TEMP_BYTE(72)
|
|
movss %xmm1, TEMP_BYTE(88)
|
|
movss %xmm2, TEMP_BYTE(104)
|
|
movss %xmm3, TEMP_BYTE(120)
|
|
|
|
movaps TEMP_BYTE(16), %xmm1
|
|
movaps TEMP_BYTE(48), %xmm3
|
|
movaps TEMP_BYTE(80), %xmm5
|
|
movaps TEMP_BYTE(112), %xmm7
|
|
movaps %xmm1, %xmm0
|
|
movaps %xmm3, %xmm2
|
|
movaps %xmm5, %xmm4
|
|
movaps %xmm7, %xmm6
|
|
shufps $0x1e, %xmm0, %xmm0
|
|
shufps $0x1e, %xmm2, %xmm2
|
|
shufps $0x1e, %xmm4, %xmm4
|
|
shufps $0x1e, %xmm6, %xmm6
|
|
andps mask, %xmm0
|
|
andps mask, %xmm2
|
|
andps mask, %xmm4
|
|
andps mask, %xmm6
|
|
addps %xmm0, %xmm1
|
|
addps %xmm2, %xmm3
|
|
addps %xmm4, %xmm5
|
|
addps %xmm6, %xmm7
|
|
|
|
movaps TEMP_BYTE(32), %xmm2
|
|
movaps TEMP_BYTE(96), %xmm6
|
|
movaps %xmm2, %xmm0
|
|
movaps %xmm6, %xmm4
|
|
shufps $0x1e, %xmm0, %xmm0
|
|
shufps $0x1e, %xmm4, %xmm4
|
|
andps mask, %xmm0
|
|
andps mask, %xmm4
|
|
addps %xmm3, %xmm2
|
|
addps %xmm0, %xmm3
|
|
addps %xmm7, %xmm6
|
|
addps %xmm4, %xmm7
|
|
|
|
movaps TEMP_BYTE(0), %xmm0
|
|
movaps TEMP_BYTE(64), %xmm4
|
|
|
|
movss %xmm0, 1024(%ecx)
|
|
movss %xmm2, 896(%ecx)
|
|
movss %xmm1, 768(%ecx)
|
|
movss %xmm3, 640(%ecx)
|
|
|
|
shufps $0xe1, %xmm0, %xmm0
|
|
shufps $0xe1, %xmm2, %xmm2
|
|
shufps $0xe1, %xmm1, %xmm1
|
|
shufps $0xe1, %xmm3, %xmm3
|
|
movss %xmm0, (%ecx)
|
|
movss %xmm0, (%ebx)
|
|
movss %xmm2, 128(%ebx)
|
|
movss %xmm1, 256(%ebx)
|
|
movss %xmm3, 384(%ebx)
|
|
|
|
movhlps %xmm0, %xmm0
|
|
movhlps %xmm2, %xmm2
|
|
movhlps %xmm1, %xmm1
|
|
movhlps %xmm3, %xmm3
|
|
movss %xmm0, 512(%ecx)
|
|
movss %xmm2, 384(%ecx)
|
|
movss %xmm1, 256(%ecx)
|
|
movss %xmm3, 128(%ecx)
|
|
|
|
shufps $0xe1, %xmm0, %xmm0
|
|
shufps $0xe1, %xmm2, %xmm2
|
|
shufps $0xe1, %xmm1, %xmm1
|
|
shufps $0xe1, %xmm3, %xmm3
|
|
movss %xmm0, 512(%ebx)
|
|
movss %xmm2, 640(%ebx)
|
|
movss %xmm1, 768(%ebx)
|
|
movss %xmm3, 896(%ebx)
|
|
|
|
movaps %xmm4, %xmm0
|
|
shufps $0x1e, %xmm0, %xmm0
|
|
movaps %xmm5, %xmm1
|
|
andps mask, %xmm0
|
|
|
|
addps %xmm6, %xmm4
|
|
addps %xmm7, %xmm5
|
|
addps %xmm1, %xmm6
|
|
addps %xmm0, %xmm7
|
|
|
|
movss %xmm4, 960(%ecx)
|
|
movss %xmm6, 832(%ecx)
|
|
movss %xmm5, 704(%ecx)
|
|
movss %xmm7, 576(%ecx)
|
|
movhlps %xmm4, %xmm0
|
|
movhlps %xmm6, %xmm1
|
|
movhlps %xmm5, %xmm2
|
|
movhlps %xmm7, %xmm3
|
|
movss %xmm0, 448(%ecx)
|
|
movss %xmm1, 320(%ecx)
|
|
movss %xmm2, 192(%ecx)
|
|
movss %xmm3, 64(%ecx)
|
|
|
|
shufps $0xe1, %xmm4, %xmm4
|
|
shufps $0xe1, %xmm6, %xmm6
|
|
shufps $0xe1, %xmm5, %xmm5
|
|
shufps $0xe1, %xmm7, %xmm7
|
|
movss %xmm4, 64(%ebx)
|
|
movss %xmm6, 192(%ebx)
|
|
movss %xmm5, 320(%ebx)
|
|
movss %xmm7, 448(%ebx)
|
|
|
|
shufps $0xe1, %xmm0, %xmm0
|
|
shufps $0xe1, %xmm1, %xmm1
|
|
shufps $0xe1, %xmm2, %xmm2
|
|
shufps $0xe1, %xmm3, %xmm3
|
|
movss %xmm0, 576(%ebx)
|
|
movss %xmm1, 704(%ebx)
|
|
movss %xmm2, 832(%ebx)
|
|
movss %xmm3, 960(%ebx)
|
|
|
|
popl %ebx
|
|
movl %ebp, %esp
|
|
popl %ebp
|
|
ret
|
|
|
|
NONEXEC_STACK
|