/*

dct36_sse: SSE optimized dct36

copyright 1995-2013 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 in %edi

#define out1 %edi

#define out2 %edx

#define w %ecx

#define ts %eax

#define tmp %esi

/*

void dct36_sse(real *inbuf,real *o1,real *o2,real *wintab,real *tsbuf);

*/

#if !defined(__APPLE__) && !defined(__EMX__)

.section .rodata

#else

.data

#endif

ALIGN16

dct36_sse_COS9:

.long 0x3f5db3d7

.long 0x3f5db3d7

.long 0x3f000000

.long 0x3f000000

.long 0x3f7c1c5c

.long 0x3f7c1c5c

.long 0x3f708fb2

.long 0x3f708fb2

.long 0x3f248dbb

.long 0x3f248dbb

.long 0x3e31d0d4

.long 0x3e31d0d4

.long 0x3eaf1d44

.long 0x3eaf1d44

.long 0x3f441b7d

.long 0x3f441b7d

ALIGN16

dct36_sse_tfcos36:

.long 0x3f007d2b

.long 0x3f0483ee

.long 0x3f0d3b7d

.long 0x3f1c4257

.long 0x40b79454

.long 0x3ff746ea

.long 0x3f976fd9

.long 0x3f5f2944

.long 0x3f3504f3

ALIGN16

dct36_sse_mask:

.long 0,0xffffffff,0,0xffffffff

ALIGN16

dct36_sse_sign:

.long 0x80000000,0x80000000,0x80000000,0x80000000

.text

ALIGN16

.globl ASM_NAME(dct36_sse)

ASM_NAME(dct36_sse):

push %ebp

mov %esp, %ebp

and $-16, %esp

sub $80, %esp

push %ebx

push %esi

push %edi

lea 12(%esp), tmp

movl 8(%ebp), in

GET_GOT

lea LOCAL_VAR(dct36_sse_COS9), %eax

lea LOCAL_VAR(dct36_sse_tfcos36), %edx

xorps %xmm0, %xmm0

xorps %xmm5, %xmm5

movlps 64(in), %xmm5

movups 48(in), %xmm4

movups 32(in), %xmm3

movups 16(in), %xmm2

movups (in), %xmm1

movaps %xmm5, %xmm6

shufps $0xe1, %xmm6, %xmm6

movaps %xmm4, %xmm7

shufps $0x93, %xmm7, %xmm7

movss %xmm7, %xmm6

addps %xmm6, %xmm5

movaps %xmm3, %xmm6

shufps $0x93, %xmm6, %xmm6

movss %xmm6, %xmm7

addps %xmm7, %xmm4

movaps %xmm2, %xmm7

shufps $0x93, %xmm7, %xmm7

movss %xmm7, %xmm6

addps %xmm6, %xmm3

movaps %xmm1, %xmm6

shufps $0x93, %xmm6, %xmm6

movss %xmm6, %xmm7

addps %xmm7, %xmm2

movss %xmm0, %xmm6

addps %xmm6, %xmm1

movaps LOCAL_VAR(dct36_sse_mask), %xmm0

movaps %xmm4, %xmm6

shufps $0x4e, %xmm5, %xmm4

movaps %xmm3, %xmm7

shufps $0x4e, %xmm6, %xmm3

andps %xmm0, %xmm6

addps %xmm6, %xmm4

movaps %xmm2, %xmm6

shufps $0x4e, %xmm7, %xmm2

andps %xmm0, %xmm7

addps %xmm7, %xmm3

movaps %xmm1, %xmm7

shufps $0x4e, %xmm6, %xmm1

andps %xmm0, %xmm6

addps %xmm6, %xmm2

movaps %xmm7, %xmm6

andps %xmm0, %xmm7

xorps %xmm0, %xmm0

addps %xmm7, %xmm1

movlhps %xmm6, %xmm0

/*

xmm0 in[-,-,0,1]

xmm1 in[2,3,4,5]

xmm2 in[6,7,8,9]

xmm3 in[10,11,12,13]

xmm4 in[14,15,16,17]

*/

movaps %xmm2, %xmm5

shufps $0xe4, %xmm3, %xmm5

shufps $0xe4, %xmm4, %xmm3

shufps $0xe4, %xmm2, %xmm4

movaps %xmm5, %xmm2

/*

xmm2 in[6,7,12,13]

xmm3 in[10,11,16,17]

xmm4 in[14,15,8,9]

*/

mulps (%eax), %xmm5

addps %xmm0, %xmm5

movaps %xmm0, (tmp)

movaps %xmm2, 16(tmp)

/*

0(tmp) in[-,-,0,1]

xmm5 [ta33,tb33,ta66,tb66]

*/

movaps %xmm1, %xmm6

subps %xmm3, %xmm6

subps %xmm4, %xmm6

xorps %xmm7, %xmm7

shufps $0xe0, %xmm2, %xmm7

mulps (%eax), %xmm6

subps %xmm7, %xmm0

addps %xmm0, %xmm6

movaps %xmm6, 48(tmp)

movaps 16(%eax), %xmm2

movaps %xmm1, %xmm0

movaps %xmm3, %xmm6

movaps %xmm4, %xmm7

mulps %xmm2, %xmm0

mulps 32(%eax), %xmm6

mulps 48(%eax), %xmm7

addps %xmm5, %xmm0

addps %xmm7, %xmm6

addps %xmm6, %xmm0

movaps %xmm0, 32(tmp)

movaps %xmm1, %xmm0

movaps %xmm3, %xmm6

movaps %xmm4, %xmm7

mulps 32(%eax), %xmm0

mulps 48(%eax), %xmm6

mulps %xmm2, %xmm7

subps %xmm5, %xmm0

subps %xmm6, %xmm7

addps %xmm7, %xmm0

movaps %xmm0, 64(tmp)

movaps %xmm1, %xmm6

movaps %xmm4, %xmm7

mulps 48(%eax), %xmm6

mulps %xmm3, %xmm2

mulps 32(%eax), %xmm7

subps %xmm5, %xmm6

subps %xmm7, %xmm2

addps %xmm2, %xmm6

movaps (tmp), %xmm0

movss 32(%edx), %xmm5

subps %xmm1, %xmm0

subps 16(tmp), %xmm4

addps %xmm3, %xmm0

addps %xmm4, %xmm0

shufps $0xaf, %xmm0, %xmm0

mulss %xmm5, %xmm0

movaps %xmm0, (tmp)

movaps 32(tmp), %xmm0

movaps 48(tmp), %xmm1

movaps 64(tmp), %xmm2

/*

xmm0 [1a-0,1b-0, 2a-0, 2b-0]

xmm1 [1a-1,1b-1, 2a-1, 2b-1]

xmm2 [1a-2,1b-2,-2a-2,-2b-2]

xmm6 [1a-3,1b-3,-2a-3,-2b-3]

*/

movaps %xmm0, %xmm3

unpcklps %xmm1, %xmm0

unpckhps %xmm1, %xmm3

movaps %xmm2, %xmm5

unpcklps %xmm6, %xmm2

unpckhps %xmm6, %xmm5

xorps LOCAL_VAR(dct36_sse_sign), %xmm5

/*

xmm0 [1a-0,1a-1,1b-0,1b-1]

xmm3 [2a-0,2a-1,2b-0,2b-1]

xmm2 [1a-2,1a-3,1b-2,1b-3]

xmm5 [2a-2,2a-3,2b-2,2b-3]

*/

movaps %xmm0, %xmm1

movlhps %xmm2, %xmm0

movhlps %xmm1, %xmm2

movaps %xmm3, %xmm4

movlhps %xmm5, %xmm3

movhlps %xmm4, %xmm5

/*

xmm0 tmp1a

xmm3 tmp2a

xmm2 tmp1b

xmm5 tmp2b

*/

movaps (%edx), %xmm6

movaps 16(%edx), %xmm7

movaps %xmm5, %xmm1

addps %xmm2, %xmm5

subps %xmm2, %xmm1

movaps %xmm3, %xmm2

addps %xmm0, %xmm3

subps %xmm0, %xmm2

mulps %xmm6, %xmm5

mulps %xmm1, %xmm7

movaps %xmm2, 16(tmp)

/*

%xmm3 tmp[0,1,2,3]

%xmm5 tmp[17,16,15,14]

16(tmp) tmp[8,7,6,5]

%xmm7 tmp[9,10,11,12]

0(tmp) tmp[13,-,4,-]

*/

movl 12(%ebp), out1

movl 16(%ebp), out2

movl 20(%ebp), w

movl 24(%ebp), ts

movaps %xmm3, %xmm0

movaps %xmm5, %xmm1

movups 108(w), %xmm2

movups 92(w), %xmm3

shufps $0x1b, %xmm3, %xmm3

movups 36(w), %xmm4

movups 20(w), %xmm5

shufps $0x1b, %xmm5, %xmm5

movaps %xmm0, %xmm6

addps %xmm1, %xmm0

subps %xmm1, %xmm6

mulps %xmm0, %xmm2

mulps %xmm3, %xmm0

mulps %xmm6, %xmm4

mulps %xmm5, %xmm6

movups 36(out1), %xmm1

movups 20(out1), %xmm3

shufps $0x1b, %xmm6, %xmm6

addps %xmm4, %xmm1

addps %xmm6, %xmm3

shufps $0x1b, %xmm0, %xmm0

movups %xmm2, 36(out2)

movups %xmm0, 20(out2)

movss %xmm1, 32*36(ts)

movss %xmm3, 32*20(ts)

movhlps %xmm1, %xmm2

movhlps %xmm3, %xmm4

movss %xmm2, 32*44(ts)

movss %xmm4, 32*28(ts)

shufps $0xb1, %xmm1, %xmm1

shufps $0xb1, %xmm3, %xmm3

movss %xmm1, 32*40(ts)

movss %xmm3, 32*24(ts)

movhlps %xmm1, %xmm2

movhlps %xmm3, %xmm4

movss %xmm2, 32*48(ts)

movss %xmm4, 32*32(ts)

movss 8(tmp), %xmm0

movss (tmp), %xmm1

movss 124(w), %xmm2

movss 88(w), %xmm3

movss 52(w), %xmm4

movss 16(w), %xmm5

movss %xmm0, %xmm6

addss %xmm1, %xmm0

subss %xmm1, %xmm6

mulss %xmm0, %xmm2

mulss %xmm3, %xmm0

mulss %xmm6, %xmm4

mulss %xmm5, %xmm6

addss 52(out1), %xmm4

addss 16(out1), %xmm6

movss %xmm2, 52(out2)

movss %xmm0, 16(out2)

movss %xmm4, 32*52(ts)

movss %xmm6, 32*16(ts)

movaps 16(tmp), %xmm0

movaps %xmm7, %xmm1

MOVUAPS 128(w), %xmm2

movups 72(w), %xmm3

shufps $0x1b, %xmm2, %xmm2

movlps 56(w), %xmm4

movhps 64(w), %xmm4

MOVUAPS (w), %xmm5

shufps $0x1b, %xmm4, %xmm4

movaps %xmm0, %xmm6

addps %xmm1, %xmm0

subps %xmm1, %xmm6

mulps %xmm0, %xmm2

mulps %xmm3, %xmm0

mulps %xmm6, %xmm4

mulps %xmm5, %xmm6

movlps 56(out1), %xmm1

movhps 64(out1), %xmm1

movups (out1), %xmm3

shufps $0x1b, %xmm4, %xmm4

addps %xmm6, %xmm3

addps %xmm4, %xmm1

shufps $0x1b, %xmm2, %xmm2

movups %xmm0, (out2)

movlps %xmm2, 56(out2)

movhps %xmm2, 64(out2)

movss %xmm1, 32*56(ts)

movss %xmm3, (ts)

movhlps %xmm1, %xmm2

movhlps %xmm3, %xmm4

movss %xmm2, 32*64(ts)

movss %xmm4, 32*8(ts)

shufps $0xb1, %xmm1, %xmm1

shufps $0xb1, %xmm3, %xmm3

movss %xmm1, 32*60(ts)

movss %xmm3, 32*4(ts)

movhlps %xmm1, %xmm2

movhlps %xmm3, %xmm4

movss %xmm2, 32*68(ts)

movss %xmm4, 32*12(ts)

pop %edi

pop %esi

pop %ebx

mov %ebp, %esp

pop %ebp

ret

NONEXEC_STACK