#include "config.h"
#endif
+#include <string.h>
+
#include <gst/gst.h>
#include <gst/base/gstbasetransform.h>
#include <gst/audio/audio.h>
#include <gst/audio/gstaudiofilter.h>
-#include <gst/controller/gstcontroller.h>
#include <math.h>
#include "audiochebband.h"
+#include "gst/glib-compat-private.h"
+
#define GST_CAT_DEFAULT gst_audio_cheb_band_debug
GST_DEBUG_CATEGORY_STATIC (GST_CAT_DEFAULT);
filter->poles = 4;
filter->ripple = 0.25;
- filter->lock = g_mutex_new ();
+ g_mutex_init (&filter->lock);
}
static void
generate_biquad_coefficients (GstAudioChebBand * filter,
- gint p, gdouble * a0, gdouble * a1, gdouble * a2, gdouble * a3,
- gdouble * a4, gdouble * b1, gdouble * b2, gdouble * b3, gdouble * b4)
+ gint p, gdouble * b0, gdouble * b1, gdouble * b2, gdouble * b3,
+ gdouble * b4, gdouble * a1, gdouble * a2, gdouble * a3, gdouble * a4)
{
gint np = filter->poles / 2;
gdouble ripple = filter->ripple;
d = 1.0 + beta * (y1 - beta * y2);
- *a0 = (x0 + beta * (-x1 + beta * x2)) / d;
- *a1 = (alpha * (-2.0 * x0 + x1 + beta * x1 - 2.0 * beta * x2)) / d;
- *a2 =
+ *b0 = (x0 + beta * (-x1 + beta * x2)) / d;
+ *b1 = (alpha * (-2.0 * x0 + x1 + beta * x1 - 2.0 * beta * x2)) / d;
+ *b2 =
(-x1 - beta * beta * x1 + 2.0 * beta * (x0 + x2) +
alpha * alpha * (x0 - x1 + x2)) / d;
- *a3 = (alpha * (x1 + beta * (-2.0 * x0 + x1) - 2.0 * x2)) / d;
- *a4 = (beta * (beta * x0 - x1) + x2) / d;
- *b1 = (alpha * (2.0 + y1 + beta * y1 - 2.0 * beta * y2)) / d;
- *b2 =
+ *b3 = (alpha * (x1 + beta * (-2.0 * x0 + x1) - 2.0 * x2)) / d;
+ *b4 = (beta * (beta * x0 - x1) + x2) / d;
+ *a1 = (alpha * (2.0 + y1 + beta * y1 - 2.0 * beta * y2)) / d;
+ *a2 =
(-y1 - beta * beta * y1 - alpha * alpha * (1.0 + y1 - y2) +
2.0 * beta * (-1.0 + y2)) / d;
- *b3 = (alpha * (y1 + beta * (2.0 + y1) - 2.0 * y2)) / d;
- *b4 = (-beta * beta - beta * y1 + y2) / d;
+ *a3 = (alpha * (y1 + beta * (2.0 + y1) - 2.0 * y2)) / d;
+ *a4 = (-beta * beta - beta * y1 + y2) / d;
} else {
a = cos ((w1 + w0) / 2.0) / cos ((w1 - w0) / 2.0);
b = tan (1.0 / 2.0) * tan ((w1 - w0) / 2.0);
d = -1.0 + beta * (beta * y2 + y1);
- *a0 = (-x0 - beta * x1 - beta * beta * x2) / d;
- *a1 = (alpha * (2.0 * x0 + x1 + beta * x1 + 2.0 * beta * x2)) / d;
- *a2 =
+ *b0 = (-x0 - beta * x1 - beta * beta * x2) / d;
+ *b1 = (alpha * (2.0 * x0 + x1 + beta * x1 + 2.0 * beta * x2)) / d;
+ *b2 =
(-x1 - beta * beta * x1 - 2.0 * beta * (x0 + x2) -
alpha * alpha * (x0 + x1 + x2)) / d;
- *a3 = (alpha * (x1 + beta * (2.0 * x0 + x1) + 2.0 * x2)) / d;
- *a4 = (-beta * beta * x0 - beta * x1 - x2) / d;
- *b1 = (alpha * (-2.0 + y1 + beta * y1 + 2.0 * beta * y2)) / d;
- *b2 =
+ *b3 = (alpha * (x1 + beta * (2.0 * x0 + x1) + 2.0 * x2)) / d;
+ *b4 = (-beta * beta * x0 - beta * x1 - x2) / d;
+ *a1 = (alpha * (-2.0 + y1 + beta * y1 + 2.0 * beta * y2)) / d;
+ *a2 =
-(y1 + beta * beta * y1 + 2.0 * beta * (-1.0 + y2) +
alpha * alpha * (-1.0 + y1 + y2)) / d;
- *b3 = (alpha * (beta * (-2.0 + y1) + y1 + 2.0 * y2)) / d;
- *b4 = -(-beta * beta + beta * y1 + y2) / d;
+ *a3 = (alpha * (beta * (-2.0 + y1) + y1 + 2.0 * y2)) / d;
+ *a4 = -(-beta * beta + beta * y1 + y2) / d;
}
}
}
if (rate == 0) {
gdouble *a = g_new0 (gdouble, 1);
+ gdouble *b = g_new0 (gdouble, 1);
a[0] = 1.0;
+ b[0] = 1.0;
gst_audio_fx_base_iir_filter_set_coefficients (GST_AUDIO_FX_BASE_IIR_FILTER
- (filter), a, 1, NULL, 0);
+ (filter), a, 1, b, 1);
GST_LOG_OBJECT (filter, "rate was not set yet");
return;
}
if (filter->upper_frequency <= filter->lower_frequency) {
gdouble *a = g_new0 (gdouble, 1);
+ gdouble *b = g_new0 (gdouble, 1);
- a[0] = (filter->mode == MODE_BAND_PASS) ? 0.0 : 1.0;
+ a[0] = 1.0;
+ b[0] = (filter->mode == MODE_BAND_PASS) ? 0.0 : 1.0;
gst_audio_fx_base_iir_filter_set_coefficients (GST_AUDIO_FX_BASE_IIR_FILTER
- (filter), a, 1, NULL, 0);
+ (filter), a, 1, b, 1);
GST_LOG_OBJECT (filter, "frequency band had no or negative dimension");
return;
b[4] = 1.0;
for (p = 1; p <= np / 4; p++) {
- gdouble a0, a1, a2, a3, a4, b1, b2, b3, b4;
+ gdouble b0, b1, b2, b3, b4, a1, a2, a3, a4;
gdouble *ta = g_new0 (gdouble, np + 5);
gdouble *tb = g_new0 (gdouble, np + 5);
- generate_biquad_coefficients (filter, p, &a0, &a1, &a2, &a3, &a4, &b1,
- &b2, &b3, &b4);
+ generate_biquad_coefficients (filter, p, &b0, &b1, &b2, &b3, &b4, &a1,
+ &a2, &a3, &a4);
memcpy (ta, a, sizeof (gdouble) * (np + 5));
memcpy (tb, b, sizeof (gdouble) * (np + 5));
* to the cascade by multiplication of the transfer
* functions */
for (i = 4; i < np + 5; i++) {
- a[i] =
- a0 * ta[i] + a1 * ta[i - 1] + a2 * ta[i - 2] + a3 * ta[i - 3] +
- a4 * ta[i - 4];
b[i] =
- tb[i] - b1 * tb[i - 1] - b2 * tb[i - 2] - b3 * tb[i - 3] -
+ b0 * tb[i] + b1 * tb[i - 1] + b2 * tb[i - 2] + b3 * tb[i - 3] +
b4 * tb[i - 4];
+ a[i] =
+ ta[i] - a1 * ta[i - 1] - a2 * ta[i - 2] - a3 * ta[i - 3] -
+ a4 * ta[i - 4];
}
g_free (ta);
g_free (tb);
}
- /* Move coefficients to the beginning of the array
- * and multiply the b coefficients with -1 to move from
+ /* Move coefficients to the beginning of the array to move from
* the transfer function's coefficients to the difference
* equation's coefficients */
- b[4] = 0.0;
for (i = 0; i <= np; i++) {
a[i] = a[i + 4];
- b[i] = -b[i + 4];
+ b[i] = b[i + 4];
}
/* Normalize to unity gain at frequency 0 and frequency
gain1 = sqrt (gain1 * gain2);
for (i = 0; i <= np; i++) {
- a[i] /= gain1;
+ b[i] /= gain1;
}
} else {
/* gain is H(wc), wc = center frequency */
zi);
for (i = 0; i <= np; i++) {
- a[i] /= gain;
+ b[i] /= gain;
}
}
{
GstAudioChebBand *filter = GST_AUDIO_CHEB_BAND (object);
- g_mutex_free (filter->lock);
- filter->lock = NULL;
+ g_mutex_clear (&filter->lock);
G_OBJECT_CLASS (parent_class)->finalize (object);
}
switch (prop_id) {
case PROP_MODE:
- g_mutex_lock (filter->lock);
+ g_mutex_lock (&filter->lock);
filter->mode = g_value_get_enum (value);
generate_coefficients (filter);
- g_mutex_unlock (filter->lock);
+ g_mutex_unlock (&filter->lock);
break;
case PROP_TYPE:
- g_mutex_lock (filter->lock);
+ g_mutex_lock (&filter->lock);
filter->type = g_value_get_int (value);
generate_coefficients (filter);
- g_mutex_unlock (filter->lock);
+ g_mutex_unlock (&filter->lock);
break;
case PROP_LOWER_FREQUENCY:
- g_mutex_lock (filter->lock);
+ g_mutex_lock (&filter->lock);
filter->lower_frequency = g_value_get_float (value);
generate_coefficients (filter);
- g_mutex_unlock (filter->lock);
+ g_mutex_unlock (&filter->lock);
break;
case PROP_UPPER_FREQUENCY:
- g_mutex_lock (filter->lock);
+ g_mutex_lock (&filter->lock);
filter->upper_frequency = g_value_get_float (value);
generate_coefficients (filter);
- g_mutex_unlock (filter->lock);
+ g_mutex_unlock (&filter->lock);
break;
case PROP_RIPPLE:
- g_mutex_lock (filter->lock);
+ g_mutex_lock (&filter->lock);
filter->ripple = g_value_get_float (value);
generate_coefficients (filter);
- g_mutex_unlock (filter->lock);
+ g_mutex_unlock (&filter->lock);
break;
case PROP_POLES:
- g_mutex_lock (filter->lock);
+ g_mutex_lock (&filter->lock);
filter->poles = GST_ROUND_UP_4 (g_value_get_int (value));
generate_coefficients (filter);
- g_mutex_unlock (filter->lock);
+ g_mutex_unlock (&filter->lock);
break;
default:
G_OBJECT_WARN_INVALID_PROPERTY_ID (object, prop_id, pspec);