mirrors/pipewire
1
0
Fork 0
mirror of https://gitlab.freedesktop.org/pipewire/pipewire.git synced 2025-11-03 09:01:54 -05:00
Code Issues Projects Releases Packages Wiki Activity Actions
pipewire/spa/plugins/audioconvert/resample-native-sse.c

75 lines
2.1 KiB
C
Raw Normal View History

treewide: use SPDX tags to specify copyright information SPDX tags make the licensing information easy to understand and clear, and they are machine parseable. See https://spdx.dev for more information.
2023-02-08 18:12:00 +01:00
/* Spa */
/* SPDX-FileCopyrightText: Copyright © 2019 Wim Taymans */
/* SPDX-License-Identifier: MIT */
audioconvert: add simple and faster resampler Add a sinc based resampler that, unlike speex, avoids memcpy and works directly on the source data. It also allows for ssse3 optimizations and aligned loads. It will later switch to table interpolation when doing variable rate.
2019-03-26 12:58:26 +01:00
audioconvert: handle more optimizations Compile an optimized library for the given CPU with the right flags, then link it with the main library.
2019-03-27 17:58:48 +01:00
#include "resample-native-impl.h"
audioconvert: add simple and faster resampler Add a sinc based resampler that, unlike speex, avoids memcpy and works directly on the source data. It also allows for ssse3 optimizations and aligned loads. It will later switch to table interpolation when doing variable rate.
2019-03-26 12:58:26 +01:00
#include <xmmintrin.h>
audioconvert: refactor peaks resampler Use common code in macro and generate arch specific version. Compile with -Ofast to optimize some fmaxf calls.
2022-09-07 16:00:31 +02:00
static inline void inner_product_sse(float *d, const float * SPA_RESTRICT s,
audioconvert: add simple and faster resampler Add a sinc based resampler that, unlike speex, avoids memcpy and works directly on the source data. It also allows for ssse3 optimizations and aligned loads. It will later switch to table interpolation when doing variable rate.
2019-03-26 12:58:26 +01:00
const float * SPA_RESTRICT taps, uint32_t n_taps)
{
__m128 sum = _mm_setzero_ps();
audioconvert: add avx optimizations
2019-03-28 16:45:57 +01:00
uint32_t i = 0;
#if 0
uint32_t unrolled = n_taps & ~15;
audioconvert: add simple and faster resampler Add a sinc based resampler that, unlike speex, avoids memcpy and works directly on the source data. It also allows for ssse3 optimizations and aligned loads. It will later switch to table interpolation when doing variable rate.
2019-03-26 12:58:26 +01:00
audioconvert: add avx optimizations
2019-03-28 16:45:57 +01:00
for (i = 0; i < unrolled; i += 16) {
sum = _mm_add_ps(sum,
_mm_mul_ps(
_mm_loadu_ps(s + i + 0),
_mm_load_ps(taps + i + 0)));
sum = _mm_add_ps(sum,
_mm_mul_ps(
_mm_loadu_ps(s + i + 4),
_mm_load_ps(taps + i + 4)));
sum = _mm_add_ps(sum,
_mm_mul_ps(
_mm_loadu_ps(s + i + 8),
_mm_load_ps(taps + i + 8)));
sum = _mm_add_ps(sum,
_mm_mul_ps(
_mm_loadu_ps(s + i + 12),
_mm_load_ps(taps + i + 12)));
}
#endif
for (; i < n_taps; i += 8) {
audioconvert: add simple and faster resampler Add a sinc based resampler that, unlike speex, avoids memcpy and works directly on the source data. It also allows for ssse3 optimizations and aligned loads. It will later switch to table interpolation when doing variable rate.
2019-03-26 12:58:26 +01:00
sum = _mm_add_ps(sum,
_mm_mul_ps(
_mm_loadu_ps(s + i + 0),
_mm_load_ps(taps + i + 0)));
sum = _mm_add_ps(sum,
_mm_mul_ps(
_mm_loadu_ps(s + i + 4),
_mm_load_ps(taps + i + 4)));
}
sum = _mm_add_ps(sum, _mm_movehl_ps(sum, sum));
sum = _mm_add_ss(sum, _mm_shuffle_ps(sum, sum, 0x55));
_mm_store_ss(d, sum);
}
audioconvert: refactor peaks resampler Use common code in macro and generate arch specific version. Compile with -Ofast to optimize some fmaxf calls.
2022-09-07 16:00:31 +02:00
static inline void inner_product_ip_sse(float *d, const float * SPA_RESTRICT s,
audioconvert: add simple and faster resampler Add a sinc based resampler that, unlike speex, avoids memcpy and works directly on the source data. It also allows for ssse3 optimizations and aligned loads. It will later switch to table interpolation when doing variable rate.
2019-03-26 12:58:26 +01:00
const float * SPA_RESTRICT t0, const float * SPA_RESTRICT t1, float x,
uint32_t n_taps)
{
__m128 sum[2] = { _mm_setzero_ps (), _mm_setzero_ps () }, t;
uint32_t i;
for (i = 0; i < n_taps; i += 8) {
t = _mm_loadu_ps(s + i + 0);
sum[0] = _mm_add_ps(sum[0], _mm_mul_ps(t, _mm_load_ps(t0 + i + 0)));
sum[1] = _mm_add_ps(sum[1], _mm_mul_ps(t, _mm_load_ps(t1 + i + 0)));
t = _mm_loadu_ps(s + i + 4);
sum[0] = _mm_add_ps(sum[0], _mm_mul_ps(t, _mm_load_ps(t0 + i + 4)));
sum[1] = _mm_add_ps(sum[1], _mm_mul_ps(t, _mm_load_ps(t1 + i + 4)));
}
sum[1] = _mm_mul_ps(_mm_sub_ps(sum[1], sum[0]), _mm_load1_ps(&x));
sum[0] = _mm_add_ps(sum[0], sum[1]);
sum[0] = _mm_add_ps(sum[0], _mm_movehl_ps(sum[0], sum[0]));
sum[0] = _mm_add_ss(sum[0], _mm_shuffle_ps(sum[0], sum[0], 0x55));
_mm_store_ss(d, sum[0]);
}
MAKE_RESAMPLER_FULL(sse);
MAKE_RESAMPLER_INTER(sse);
Reference in a new issue Copy permalink