1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
|
--- libaf/af_equalizer.c.orig 2009-05-13 04:58:57.000000000 +0200
+++ libaf/af_equalizer.c 2010-03-09 22:03:10.000000000 +0100
@@ -51,7 +51,7 @@
8 8.000 kHz
9 16.00 kHz
*/
-#define CF {31.25,62.5,125,250,500,1000,2000,4000,8000,16000}
+#define CF {23.4375,46.875,93.75,187.5,375,750,1500,3000,6000,12000}
// Maximum and minimum gain for the bands
#define G_MAX +12.0
@@ -60,9 +60,9 @@
// Data for specific instances of this filter
typedef struct af_equalizer_s
{
- float a[KM][L]; // A weights
- float b[KM][L]; // B weights
- float wq[AF_NCH][KM][L]; // Circular buffer for W data
+ double a[KM][L]; // A weights
+ double b[KM][L]; // B weights
+ double wq[AF_NCH][KM][L]; // Circular buffer for W data
float g[AF_NCH][KM]; // Gain factor for each channel and band
int K; // Number of used eq bands
int channels; // Number of channels
@@ -70,9 +70,9 @@
} af_equalizer_t;
// 2nd order Band-pass Filter design
-static void bp2(float* a, float* b, float fc, float q){
+static void bp2(double* a, double* b, double fc, double q){
double th= 2.0 * M_PI * fc;
- double C = (1.0 - tan(th*q/2.0))/(1.0 + tan(th*q/2.0));
+ double C = (1.0 - tan(th*q/4.0))/(1.0 + tan(th*q/4.0));
a[0] = (1.0 + C) * cos(th);
a[1] = -1 * C;
@@ -89,7 +89,7 @@
switch(cmd){
case AF_CONTROL_REINIT:{
int k =0, i =0;
- float F[KM] = CF;
+ double F[KM] = CF;
s->gain_factor=0.0;
@@ -103,7 +103,7 @@
// Calculate number of active filters
s->K=KM;
- while(F[s->K-1] > (float)af->data->rate/2.2)
+ while(F[s->K-1] > (double)af->data->rate/2.2)
s->K--;
if(s->K != KM)
@@ -112,7 +112,7 @@
// Generate filter taps
for(k=0;k<s->K;k++)
- bp2(s->a[k],s->b[k],F[k]/((float)af->data->rate),Q);
+ bp2(s->a[k],s->b[k],F[k]/((double)af->data->rate),Q);
// Calculate how much this plugin adds to the overall time delay
af->delay = 2 * af->data->nch * af->data->bps;
@@ -209,9 +209,9 @@
// Run the filters
for(;k<s->K;k++){
// Pointer to circular buffer wq
- register float* wq = s->wq[ci][k];
+ double* wq = s->wq[ci][k];
// Calculate output from AR part of current filter
- register float w=yt*s->b[k][0] + wq[0]*s->a[k][0] + wq[1]*s->a[k][1];
+ double w=yt*s->b[k][0] + wq[0]*s->a[k][0] + wq[1]*s->a[k][1];
// Calculate output form MA part of current filter
yt+=(w + wq[1]*s->b[k][1])*g[k];
// Update circular buffer
|