+#include <assert.h>
#include <math.h>
#include <stdio.h>
#include <string.h>
-#include <assert.h>
-#include <complex>
#include <algorithm>
+#include <complex>
+
+#include "defs.h"
#ifdef __SSE__
#include <mmintrin.h>
{
omega = M_PI;
resonance = 0.01f;
+ A = 1.0f;
init(FILTER_NONE, 1);
update();
// a2 = 1 - alpha;
break;
+ case FILTER_PEAKING_EQ:
+ b0 = 1 + alpha * A;
+ b1 = -2*cs;
+ b2 = 1 - alpha * A;
+ a0 = 1 + alpha / A;
+ // a1 = -2*cs;
+ a2 = 1 - alpha / A;
+ break;
+
+ case FILTER_LOW_SHELF:
+ b0 = A * ((A + 1) - (A - 1)*cs + 2 * sqrt(A) * alpha);
+ b1 = 2 * A * ((A - 1) - (A + 1)*cs );
+ b2 = A * ((A + 1) - (A - 1)*cs - 2 * sqrt(A) * alpha);
+ a0 = (A + 1) + (A - 1)*cs + 2 * sqrt(A) * alpha ;
+ a1 = -2 * ((A - 1) + (A + 1)*cs );
+ a2 = (A + 1) + (A - 1)*cs - 2 * sqrt(A) * alpha ;
+ break;
+
+ case FILTER_HIGH_SHELF:
+ b0 = A * ((A + 1) + (A - 1)*cs + 2 * sqrt(A) * alpha);
+ b1 = -2 * A * ((A - 1) + (A + 1)*cs );
+ b2 = A * ((A + 1) + (A - 1)*cs - 2 * sqrt(A) * alpha);
+ a0 = (A + 1) - (A - 1)*cs + 2 * sqrt(A) * alpha ;
+ a1 = 2 * ((A - 1) - (A + 1)*cs );
+ a2 = (A + 1) - (A - 1)*cs - 2 * sqrt(A) * alpha ;
+ break;
+
default:
//unknown filter type
assert(false);
memset(feedback, 0, sizeof(feedback));
#else
for (unsigned i = 0; i < 2; ++i) {
- filters[i].init(type, 0, new_order, 0);
+ filters[i].init(type, new_order);
}
#endif
}
-void StereoFilter::render(float *inout_left_ptr, unsigned n_samples, float cutoff, float resonance)
+void StereoFilter::render(float *inout_left_ptr, unsigned n_samples, float cutoff, float resonance, float dbgain_normalized)
{
#ifdef __SSE__
if (parm_filter.filtertype == FILTER_NONE || parm_filter.filter_order == 0)
parm_filter.set_linear_cutoff(cutoff);
parm_filter.set_resonance(resonance);
+ parm_filter.set_dbgain_normalized(dbgain_normalized);
parm_filter.update();
__m128 b0 = _mm_set1_ps(parm_filter.b0);