FFmpeg
enc_utils.h
Go to the documentation of this file.
1 /*
2  * Opus encoder
3  * Copyright (c) 2017 Rostislav Pehlivanov <atomnuker@gmail.com>
4  *
5  * This file is part of FFmpeg.
6  *
7  * FFmpeg is free software; you can redistribute it and/or
8  * modify it under the terms of the GNU Lesser General Public
9  * License as published by the Free Software Foundation; either
10  * version 2.1 of the License, or (at your option) any later version.
11  *
12  * FFmpeg is distributed in the hope that it will be useful,
13  * but WITHOUT ANY WARRANTY; without even the implied warranty of
14  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
15  * Lesser General Public License for more details.
16  *
17  * You should have received a copy of the GNU Lesser General Public
18  * License along with FFmpeg; if not, write to the Free Software
19  * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
20  */
21 
22 #ifndef AVCODEC_OPUS_ENC_UTILS_H
23 #define AVCODEC_OPUS_ENC_UTILS_H
24 
25 #include <math.h>
26 #include <string.h>
27 
28 #include "opus.h"
29 
30 typedef struct FFBesselFilter {
31  float a[3];
32  float b[2];
33  float x[3];
34  float y[3];
36 
37 /* Fills the coefficients, returns 1 if filter will be unstable */
38 static inline int bessel_reinit(FFBesselFilter *s, float n, float f0, float fs,
39  int highpass)
40 {
41  int unstable;
42  float c, cfreq, w0, k1, k2;
43 
44  if (!highpass) {
45  c = (1.0f/sqrtf(sqrtf(pow(2.0f, 1.0f/n) - 3.0f/4.0f) - 0.5f))/sqrtf(3.0f);
46  cfreq = c*f0/fs;
47  unstable = (cfreq <= 0.0f || cfreq >= 1.0f/4.0f);
48  } else {
49  c = sqrtf(3.0f)*sqrtf(sqrtf(pow(2.0f, 1.0f/n) - 3.0f/4.0f) - 0.5f);
50  cfreq = 0.5f - c*f0/fs;
51  unstable = (cfreq <= 3.0f/8.0f || cfreq >= 1.0f/2.0f);
52  }
53 
54  w0 = tanf(M_PI*cfreq);
55  k1 = 3.0f * w0;
56  k2 = 3.0f * w0;
57 
58  s->a[0] = k2/(1.0f + k1 + k2);
59  s->a[1] = 2.0f * s->a[0];
60  s->a[2] = s->a[0];
61  s->b[0] = 2.0f * s->a[0] * (1.0f/k2 - 1.0f);
62  s->b[1] = 1.0f - (s->a[0] + s->a[1] + s->a[2] + s->b[0]);
63 
64  if (highpass) {
65  s->a[1] *= -1;
66  s->b[0] *= -1;
67  }
68 
69  return unstable;
70 }
71 
72 static inline int bessel_init(FFBesselFilter *s, float n, float f0, float fs,
73  int highpass)
74 {
75  memset(s, 0, sizeof(FFBesselFilter));
76  return bessel_reinit(s, n, f0, fs, highpass);
77 }
78 
79 static inline float bessel_filter(FFBesselFilter *s, float x)
80 {
81  s->x[2] = s->x[1];
82  s->x[1] = s->x[0];
83  s->x[0] = x;
84  s->y[2] = s->y[1];
85  s->y[1] = s->y[0];
86  s->y[0] = s->a[0]*s->x[0] + s->a[1]*s->x[1] + s->a[2]*s->x[2] + s->b[0]*s->y[1] + s->b[1]*s->y[2];
87  return s->y[0];
88 }
89 
90 #endif /* AVCODEC_OPUS_ENC_UTILS_H */
opus.h
FFBesselFilter::a
float a[3]
Definition: enc_utils.h:31
bessel_filter
static float bessel_filter(FFBesselFilter *s, float x)
Definition: enc_utils.h:79
s
#define s(width, name)
Definition: cbs_vp9.c:198
highpass
@ highpass
Definition: af_biquads.c:85
FFBesselFilter::b
float b[2]
Definition: enc_utils.h:32
fs
#define fs(width, name, subs,...)
Definition: cbs_vp9.c:200
FFBesselFilter::y
float y[3]
Definition: enc_utils.h:34
sqrtf
static __device__ float sqrtf(float a)
Definition: cuda_runtime.h:184
c
Undefined Behavior In the C some operations are like signed integer dereferencing freed accessing outside allocated Undefined Behavior must not occur in a C it is not safe even if the output of undefined operations is unused The unsafety may seem nit picking but Optimizing compilers have in fact optimized code on the assumption that no undefined Behavior occurs Optimizing code based on wrong assumptions can and has in some cases lead to effects beyond the output of computations The signed integer overflow problem in speed critical code Code which is highly optimized and works with signed integers sometimes has the problem that often the output of the computation does not c
Definition: undefined.txt:32
f
f
Definition: af_crystalizer.c:122
bessel_init
static int bessel_init(FFBesselFilter *s, float n, float f0, float fs, int highpass)
Definition: enc_utils.h:72
M_PI
#define M_PI
Definition: mathematics.h:67
FFBesselFilter::x
float x[3]
Definition: enc_utils.h:33
FFBesselFilter
Definition: enc_utils.h:30
bessel_reinit
static int bessel_reinit(FFBesselFilter *s, float n, float f0, float fs, int highpass)
Definition: enc_utils.h:38