virtualenv/FFmpeg
1
0
Fork 0
mirror of https://github.com/FFmpeg/FFmpeg.git synced 2025-02-04 06:08:26 +02:00
Code Issues Releases Activity

avfilter/asrc_sinc: Use av_bessel_i0()

Browse Source 
The new function is much more precise
For default beta it is slightly slower, but its speed is already at the
worst case in that comparison
while the replaced function becomes much slower for larger beta

Signed-off-by: Michael Niedermayer <michael@niedermayer.cc>
This commit is contained in:
Michael Niedermayer 2023-05-22 23:39:25 +02:00
parent 75918016ab
commit 1e9c337e0f
No known key found for this signature in database
GPG Key ID: B18E8928B3948D64
1 changed files with 2 additions and 17 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

19
libavfilter/asrc_sinc.c
View File

@ -91,27 +91,12 @@ static int query_formats(AVFilterContext *ctx)
return ff_set_common_samplerates_from_list(ctx, sample_rates);
}
static float bessel_I_0(float x)
{
float term = 1, sum = 1, last_sum, x2 = x / 2;
int i = 1;
do {
float y = x2 / i++;
last_sum = sum;
sum += term *= y * y;
} while (sum != last_sum);
return sum;
}
static float *make_lpf(int num_taps, float Fc, float beta, float rho,
float scale, int dc_norm)
{
int i, m = num_taps - 1;
float *h = av_calloc(num_taps, sizeof(*h)), sum = 0;
float mult = scale / bessel_I_0(beta), mult1 = 1.f / (.5f * m + rho);
float mult = scale / av_bessel_i0(beta), mult1 = 1.f / (.5f * m + rho);
if (!h)
return NULL;
@ -121,7 +106,7 @@ static float *make_lpf(int num_taps, float Fc, float beta, float rho,
for (i = 0; i <= m / 2; i++) {
float z = i - .5f * m, x = z * M_PI, y = z * mult1;
h[i] = x ? sinf(Fc * x) / x : Fc;
sum += h[i] *= bessel_I_0(beta * sqrtf(1.f - y * y)) * mult;
sum += h[i] *= av_bessel_i0(beta * sqrtf(1.f - y * y)) * mult;
if (m - i != i) {
h[m - i] = h[i];
sum += h[i];