mirror of
https://github.com/FFmpeg/FFmpeg.git
synced 2024-11-26 19:01:44 +02:00
3ade6a8644
Old one was written with the assumption only even inputs would be given. This very messy replacement supports even and odd inputs, and supports AVX2 for extra speed. The buffers given are usually quite big (4k samples), so the speedup is worth it. The new SSE version is still faster than the old inline asm version by 33%. Also checkasm is provided to make sure this monstrosity works. This fixes some FATE tests.
121 lines
4.7 KiB
C
121 lines
4.7 KiB
C
/*
|
|
* SIMD-optimized LPC functions
|
|
* Copyright (c) 2007 Loren Merritt
|
|
*
|
|
* This file is part of FFmpeg.
|
|
*
|
|
* FFmpeg is free software; you can redistribute it and/or
|
|
* modify it under the terms of the GNU Lesser General Public
|
|
* License as published by the Free Software Foundation; either
|
|
* version 2.1 of the License, or (at your option) any later version.
|
|
*
|
|
* FFmpeg is distributed in the hope that it will be useful,
|
|
* but WITHOUT ANY WARRANTY; without even the implied warranty of
|
|
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
|
|
* Lesser General Public License for more details.
|
|
*
|
|
* You should have received a copy of the GNU Lesser General Public
|
|
* License along with FFmpeg; if not, write to the Free Software
|
|
* Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
|
|
*/
|
|
|
|
#include "libavutil/attributes.h"
|
|
#include "libavutil/x86/asm.h"
|
|
#include "libavutil/x86/cpu.h"
|
|
#include "libavcodec/lpc.h"
|
|
|
|
void ff_lpc_apply_welch_window_sse2(const int32_t *data, int len,
|
|
double *w_data);
|
|
void ff_lpc_apply_welch_window_avx2(const int32_t *data, int len,
|
|
double *w_data);
|
|
|
|
DECLARE_ASM_CONST(16, double, pd_1)[2] = { 1.0, 1.0 };
|
|
|
|
#if HAVE_SSE2_INLINE
|
|
|
|
static void lpc_compute_autocorr_sse2(const double *data, int len, int lag,
|
|
double *autoc)
|
|
{
|
|
int j;
|
|
|
|
if((x86_reg)data & 15)
|
|
data++;
|
|
|
|
for(j=0; j<lag; j+=2){
|
|
x86_reg i = -len*sizeof(double);
|
|
if(j == lag-2) {
|
|
__asm__ volatile(
|
|
"movsd "MANGLE(pd_1)", %%xmm0 \n\t"
|
|
"movsd "MANGLE(pd_1)", %%xmm1 \n\t"
|
|
"movsd "MANGLE(pd_1)", %%xmm2 \n\t"
|
|
"1: \n\t"
|
|
"movapd (%2,%0), %%xmm3 \n\t"
|
|
"movupd -8(%3,%0), %%xmm4 \n\t"
|
|
"movapd (%3,%0), %%xmm5 \n\t"
|
|
"mulpd %%xmm3, %%xmm4 \n\t"
|
|
"mulpd %%xmm3, %%xmm5 \n\t"
|
|
"mulpd -16(%3,%0), %%xmm3 \n\t"
|
|
"addpd %%xmm4, %%xmm1 \n\t"
|
|
"addpd %%xmm5, %%xmm0 \n\t"
|
|
"addpd %%xmm3, %%xmm2 \n\t"
|
|
"add $16, %0 \n\t"
|
|
"jl 1b \n\t"
|
|
"movhlps %%xmm0, %%xmm3 \n\t"
|
|
"movhlps %%xmm1, %%xmm4 \n\t"
|
|
"movhlps %%xmm2, %%xmm5 \n\t"
|
|
"addsd %%xmm3, %%xmm0 \n\t"
|
|
"addsd %%xmm4, %%xmm1 \n\t"
|
|
"addsd %%xmm5, %%xmm2 \n\t"
|
|
"movsd %%xmm0, (%1) \n\t"
|
|
"movsd %%xmm1, 8(%1) \n\t"
|
|
"movsd %%xmm2, 16(%1) \n\t"
|
|
:"+&r"(i)
|
|
:"r"(autoc+j), "r"(data+len), "r"(data+len-j)
|
|
NAMED_CONSTRAINTS_ARRAY_ADD(pd_1)
|
|
:"memory"
|
|
);
|
|
} else {
|
|
__asm__ volatile(
|
|
"movsd "MANGLE(pd_1)", %%xmm0 \n\t"
|
|
"movsd "MANGLE(pd_1)", %%xmm1 \n\t"
|
|
"1: \n\t"
|
|
"movapd (%3,%0), %%xmm3 \n\t"
|
|
"movupd -8(%4,%0), %%xmm4 \n\t"
|
|
"mulpd %%xmm3, %%xmm4 \n\t"
|
|
"mulpd (%4,%0), %%xmm3 \n\t"
|
|
"addpd %%xmm4, %%xmm1 \n\t"
|
|
"addpd %%xmm3, %%xmm0 \n\t"
|
|
"add $16, %0 \n\t"
|
|
"jl 1b \n\t"
|
|
"movhlps %%xmm0, %%xmm3 \n\t"
|
|
"movhlps %%xmm1, %%xmm4 \n\t"
|
|
"addsd %%xmm3, %%xmm0 \n\t"
|
|
"addsd %%xmm4, %%xmm1 \n\t"
|
|
"movsd %%xmm0, %1 \n\t"
|
|
"movsd %%xmm1, %2 \n\t"
|
|
:"+&r"(i), "=m"(autoc[j]), "=m"(autoc[j+1])
|
|
:"r"(data+len), "r"(data+len-j)
|
|
NAMED_CONSTRAINTS_ARRAY_ADD(pd_1)
|
|
);
|
|
}
|
|
}
|
|
}
|
|
|
|
#endif /* HAVE_SSE2_INLINE */
|
|
|
|
av_cold void ff_lpc_init_x86(LPCContext *c)
|
|
{
|
|
int cpu_flags = av_get_cpu_flags();
|
|
|
|
#if HAVE_SSE2_INLINE
|
|
if (INLINE_SSE2_SLOW(cpu_flags))
|
|
c->lpc_compute_autocorr = lpc_compute_autocorr_sse2;
|
|
#endif
|
|
|
|
if (EXTERNAL_SSE2(cpu_flags))
|
|
c->lpc_apply_welch_window = ff_lpc_apply_welch_window_sse2;
|
|
|
|
if (EXTERNAL_AVX2(cpu_flags))
|
|
c->lpc_apply_welch_window = ff_lpc_apply_welch_window_avx2;
|
|
}
|