1cabdff1aSopenharmony_ci;******************************************************************************
2cabdff1aSopenharmony_ci;* Opus SIMD functions
3cabdff1aSopenharmony_ci;*
4cabdff1aSopenharmony_ci;* This file is part of FFmpeg.
5cabdff1aSopenharmony_ci;*
6cabdff1aSopenharmony_ci;* FFmpeg is free software; you can redistribute it and/or
7cabdff1aSopenharmony_ci;* modify it under the terms of the GNU Lesser General Public
8cabdff1aSopenharmony_ci;* License as published by the Free Software Foundation; either
9cabdff1aSopenharmony_ci;* version 2.1 of the License, or (at your option) any later version.
10cabdff1aSopenharmony_ci;*
11cabdff1aSopenharmony_ci;* FFmpeg is distributed in the hope that it will be useful,
12cabdff1aSopenharmony_ci;* but WITHOUT ANY WARRANTY; without even the implied warranty of
13cabdff1aSopenharmony_ci;* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
14cabdff1aSopenharmony_ci;* Lesser General Public License for more details.
15cabdff1aSopenharmony_ci;*
16cabdff1aSopenharmony_ci;* You should have received a copy of the GNU Lesser General Public
17cabdff1aSopenharmony_ci;* License along with FFmpeg; if not, write to the Free Software
18cabdff1aSopenharmony_ci;* Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
19cabdff1aSopenharmony_ci;******************************************************************************
20cabdff1aSopenharmony_ci
21cabdff1aSopenharmony_ci%include "libavutil/x86/x86util.asm"
22cabdff1aSopenharmony_ci
23cabdff1aSopenharmony_ciSECTION_RODATA
24cabdff1aSopenharmony_ci
25cabdff1aSopenharmony_ci         ; 0.85..^1    0.85..^2    0.85..^3    0.85..^4
26cabdff1aSopenharmony_citab_st: dd 0x3f599a00, 0x3f38f671, 0x3f1d382a, 0x3f05a32f
27cabdff1aSopenharmony_ci
28cabdff1aSopenharmony_ciSECTION .text
29cabdff1aSopenharmony_ci
30cabdff1aSopenharmony_ciINIT_XMM fma3
31cabdff1aSopenharmony_ci%if UNIX64
32cabdff1aSopenharmony_cicglobal opus_deemphasis, 3, 3, 8, out, in, len
33cabdff1aSopenharmony_ci%else
34cabdff1aSopenharmony_cicglobal opus_deemphasis, 4, 4, 8, out, in, coeff, len
35cabdff1aSopenharmony_ci%endif
36cabdff1aSopenharmony_ci%if ARCH_X86_32
37cabdff1aSopenharmony_ci    VBROADCASTSS m0, coeffm
38cabdff1aSopenharmony_ci%elif WIN64
39cabdff1aSopenharmony_ci    shufps m0, m2, m2, 0
40cabdff1aSopenharmony_ci%else
41cabdff1aSopenharmony_ci    shufps m0, m0, 0
42cabdff1aSopenharmony_ci%endif
43cabdff1aSopenharmony_ci
44cabdff1aSopenharmony_ci    movaps m4, [tab_st]
45cabdff1aSopenharmony_ci    VBROADCASTSS m5, m4
46cabdff1aSopenharmony_ci    shufps m6, m4, m4, q1111
47cabdff1aSopenharmony_ci    shufps m7, m4, m4, q2222
48cabdff1aSopenharmony_ci
49cabdff1aSopenharmony_ci.loop:
50cabdff1aSopenharmony_ci    movaps  m1, [inq]                ; x0, x1, x2, x3
51cabdff1aSopenharmony_ci
52cabdff1aSopenharmony_ci    pslldq  m2, m1, 4                ;  0, x0, x1, x2
53cabdff1aSopenharmony_ci    pslldq  m3, m1, 8                ;  0,  0, x0, x1
54cabdff1aSopenharmony_ci
55cabdff1aSopenharmony_ci    fmaddps m2, m2, m5, m1           ; x + c1*x[0-2]
56cabdff1aSopenharmony_ci    pslldq  m1, 12                   ;  0,  0,  0, x0
57cabdff1aSopenharmony_ci
58cabdff1aSopenharmony_ci    fmaddps m2, m3, m6, m2           ; x + c1*x[0-2] + c2*x[0-1]
59cabdff1aSopenharmony_ci    fmaddps m1, m1, m7, m2           ; x + c1*x[0-2] + c2*x[0-1] + c3*x[0]
60cabdff1aSopenharmony_ci    fmaddps m0, m0, m4, m1           ; x + c1*x[0-2] + c2*x[0-1] + c3*x[0] + c*s
61cabdff1aSopenharmony_ci
62cabdff1aSopenharmony_ci    movaps [outq], m0
63cabdff1aSopenharmony_ci    shufps m0, m0, q3333             ; new state
64cabdff1aSopenharmony_ci
65cabdff1aSopenharmony_ci    add inq,  mmsize
66cabdff1aSopenharmony_ci    add outq, mmsize
67cabdff1aSopenharmony_ci    sub lend, mmsize >> 2
68cabdff1aSopenharmony_ci    jg .loop
69cabdff1aSopenharmony_ci
70cabdff1aSopenharmony_ci%if ARCH_X86_64 == 0
71cabdff1aSopenharmony_ci    movss r0m, m0
72cabdff1aSopenharmony_ci    fld dword r0m
73cabdff1aSopenharmony_ci%endif
74cabdff1aSopenharmony_ci    RET
75cabdff1aSopenharmony_ci
76cabdff1aSopenharmony_ci
77cabdff1aSopenharmony_ciINIT_XMM fma3
78cabdff1aSopenharmony_cicglobal opus_postfilter, 4, 4, 8, data, period, gains, len
79cabdff1aSopenharmony_ci    VBROADCASTSS m0, [gainsq + 0]
80cabdff1aSopenharmony_ci    VBROADCASTSS m1, [gainsq + 4]
81cabdff1aSopenharmony_ci    VBROADCASTSS m2, [gainsq + 8]
82cabdff1aSopenharmony_ci
83cabdff1aSopenharmony_ci    shl periodd, 2
84cabdff1aSopenharmony_ci    add periodq, 8
85cabdff1aSopenharmony_ci    neg periodq
86cabdff1aSopenharmony_ci
87cabdff1aSopenharmony_ci    movups  m3, [dataq + periodq]
88cabdff1aSopenharmony_ci    mulps   m3, m2
89cabdff1aSopenharmony_ci
90cabdff1aSopenharmony_ci.loop:
91cabdff1aSopenharmony_ci    movups  m4, [dataq + periodq +  4]
92cabdff1aSopenharmony_ci    movups  m5, [dataq + periodq +  8]
93cabdff1aSopenharmony_ci    movups  m6, [dataq + periodq + 12]
94cabdff1aSopenharmony_ci    movups  m7, [dataq + periodq + 16]
95cabdff1aSopenharmony_ci
96cabdff1aSopenharmony_ci    fmaddps m3, m7, m2, m3
97cabdff1aSopenharmony_ci    addps   m6, m4
98cabdff1aSopenharmony_ci
99cabdff1aSopenharmony_ci    fmaddps m5, m5, m0, [dataq]
100cabdff1aSopenharmony_ci    fmaddps m6, m6, m1, m3
101cabdff1aSopenharmony_ci
102cabdff1aSopenharmony_ci    addps   m5, m6
103cabdff1aSopenharmony_ci    mulps   m3, m7, m2
104cabdff1aSopenharmony_ci
105cabdff1aSopenharmony_ci    movaps  [dataq], m5
106cabdff1aSopenharmony_ci
107cabdff1aSopenharmony_ci    add dataq, mmsize
108cabdff1aSopenharmony_ci    sub lend,  mmsize >> 2
109cabdff1aSopenharmony_ci    jg .loop
110cabdff1aSopenharmony_ci
111cabdff1aSopenharmony_ci    RET
112