summaryrefslogtreecommitdiffstats
diff options
context:
space:
mode:
authorSebastian Dröge <slomo@circular-chaos.org>2007-08-10 04:44:43 +0000
committerSebastian Dröge <slomo@circular-chaos.org>2007-08-10 04:44:43 +0000
commit11131b7a44ce855d89959c7af1b5773da1630ad7 (patch)
tree7ab8664a9d6a958d9b506a55e7c8aa33640acba0
parent982bf23b353cd825e249d8fc13420f88e9cf4cd4 (diff)
downloadgst-plugins-bad-11131b7a44ce855d89959c7af1b5773da1630ad7.tar.gz
gst-plugins-bad-11131b7a44ce855d89959c7af1b5773da1630ad7.tar.bz2
gst-plugins-bad-11131b7a44ce855d89959c7af1b5773da1630ad7.zip
gst/filter/gstlpwsinc.*: Specify the actual filter length instead of a weird 2N+1. Setting the property will round to...
Original commit message from CVS: * gst/filter/gstlpwsinc.c: (gst_lpwsinc_class_init), (gst_lpwsinc_init), (process_32), (process_64), (lpwsinc_build_kernel), (lpwsinc_set_property), (lpwsinc_get_property): * gst/filter/gstlpwsinc.h: Specify the actual filter length instead of a weird 2N+1. Setting the property will round to the next odd number. Also remove now obsolete FIXMEs.
-rw-r--r--ChangeLog11
-rw-r--r--gst/filter/gstlpwsinc.c50
-rw-r--r--gst/filter/gstlpwsinc.h7
3 files changed, 40 insertions, 28 deletions
diff --git a/ChangeLog b/ChangeLog
index 3707cf27..a62b0116 100644
--- a/ChangeLog
+++ b/ChangeLog
@@ -1,5 +1,16 @@
2007-08-10 Sebastian Dröge <slomo@circular-chaos.org>
+ * gst/filter/gstlpwsinc.c: (gst_lpwsinc_class_init),
+ (gst_lpwsinc_init), (process_32), (process_64),
+ (lpwsinc_build_kernel), (lpwsinc_set_property),
+ (lpwsinc_get_property):
+ * gst/filter/gstlpwsinc.h:
+ Specify the actual filter length instead of a weird
+ 2N+1. Setting the property will round to the next odd number.
+ Also remove now obsolete FIXMEs.
+
+2007-08-10 Sebastian Dröge <slomo@circular-chaos.org>
+
* gst/filter/gstlpwsinc.c: (gst_lpwsinc_window_get_type),
(gst_lpwsinc_class_init), (gst_lpwsinc_init),
(lpwsinc_build_kernel), (lpwsinc_set_property),
diff --git a/gst/filter/gstlpwsinc.c b/gst/filter/gstlpwsinc.c
index fce2b00d..24906d12 100644
--- a/gst/filter/gstlpwsinc.c
+++ b/gst/filter/gstlpwsinc.c
@@ -204,8 +204,8 @@ gst_lpwsinc_class_init (GstLPWSincClass * klass)
g_object_class_install_property (gobject_class, PROP_LENGTH,
g_param_spec_int ("length", "Length",
- "N such that the filter length = 2N + 1",
- 1, G_MAXINT, 1, G_PARAM_READWRITE | GST_PARAM_CONTROLLABLE));
+ "Filter kernel length, will be rounded to the next odd number",
+ 3, G_MAXINT, 101, G_PARAM_READWRITE | GST_PARAM_CONTROLLABLE));
g_object_class_install_property (gobject_class, PROP_MODE,
g_param_spec_enum ("mode", "Mode",
@@ -227,7 +227,7 @@ gst_lpwsinc_init (GstLPWSinc * self, GstLPWSincClass * g_class)
{
self->mode = MODE_LOW_PASS;
self->window = WINDOW_HAMMING;
- self->wing_size = 50;
+ self->kernel_length = 101;
self->frequency = 0.0;
self->kernel = NULL;
self->residue = NULL;
@@ -238,7 +238,7 @@ gst_lpwsinc_init (GstLPWSinc * self, GstLPWSincClass * g_class)
static void
process_32 (GstLPWSinc * self, gfloat * src, gfloat * dst, guint input_samples)
{
- gint kernel_length = self->wing_size * 2 + 1;
+ gint kernel_length = self->kernel_length;
gint i, j, k, l;
gint channels = GST_AUDIO_FILTER (self)->format.channels;
@@ -265,7 +265,7 @@ static void
process_64 (GstLPWSinc * self, gdouble * src, gdouble * dst,
guint input_samples)
{
- gint kernel_length = self->wing_size * 2 + 1;
+ gint kernel_length = self->kernel_length;
gint i, j, k, l;
gint channels = GST_AUDIO_FILTER (self)->format.channels;
@@ -297,8 +297,8 @@ lpwsinc_build_kernel (GstLPWSinc * self)
gdouble w;
/* fill the kernel */
- len = self->wing_size;
- GST_DEBUG ("lpwsinc: initializing filter kernel of length %d", len * 2 + 1);
+ len = self->kernel_length;
+ GST_DEBUG ("lpwsinc: initializing filter kernel of length %d", len);
if (GST_AUDIO_FILTER (self)->format.rate == 0) {
GST_DEBUG ("rate not set yet");
@@ -318,40 +318,40 @@ lpwsinc_build_kernel (GstLPWSinc * self)
if (self->kernel)
g_free (self->kernel);
- self->kernel = g_new (gdouble, 2 * len + 1);
+ self->kernel = g_new (gdouble, len);
- for (i = 0; i <= len * 2; ++i) {
- if (i == len)
+ for (i = 0; i < len; ++i) {
+ if (i == len / 2)
self->kernel[i] = w;
else
- self->kernel[i] = sin (w * (i - len)) / (i - len);
+ self->kernel[i] = sin (w * (i - len / 2)) / (i - len / 2);
/* windowing */
if (self->window == WINDOW_HAMMING)
- self->kernel[i] *= (0.54 - 0.46 * cos (M_PI * i / len));
+ self->kernel[i] *= (0.54 - 0.46 * cos (2 * M_PI * i / len));
else
self->kernel[i] *=
- (0.42 - 0.5 * cos (M_PI * i / len) + 0.08 * cos (2 * M_PI * i / len));
+ (0.42 - 0.5 * cos (2 * M_PI * i / len) +
+ 0.08 * cos (4 * M_PI * i / len));
}
/* normalize for unity gain at DC */
- for (i = 0; i <= len * 2; ++i)
+ for (i = 0; i < len; ++i)
sum += self->kernel[i];
- for (i = 0; i <= len * 2; ++i)
+ for (i = 0; i < len; ++i)
self->kernel[i] /= sum;
/* convert to highpass if specified */
if (self->mode == MODE_HIGH_PASS) {
- for (i = 0; i <= len * 2; ++i)
+ for (i = 0; i < len; ++i)
self->kernel[i] = -self->kernel[i];
- self->kernel[len] += 1.0;
+ self->kernel[len / 2] += 1.0;
}
/* set up the residue memory space */
if (self->residue)
g_free (self->residue);
self->residue =
- g_new0 (gdouble,
- (len * 2 + 1) * GST_AUDIO_FILTER (self)->format.channels);
+ g_new0 (gdouble, len * GST_AUDIO_FILTER (self)->format.channels);
self->have_kernel = TRUE;
}
@@ -434,12 +434,18 @@ lpwsinc_set_property (GObject * object, guint prop_id, const GValue * value,
g_return_if_fail (GST_IS_LPWSINC (self));
switch (prop_id) {
- case PROP_LENGTH:
+ case PROP_LENGTH:{
+ gint val;
+
GST_BASE_TRANSFORM_LOCK (self);
- self->wing_size = g_value_get_int (value);
+ val = g_value_get_int (value);
+ if (val % 2 == 0)
+ val++;
+ self->kernel_length = val;
lpwsinc_build_kernel (self);
GST_BASE_TRANSFORM_UNLOCK (self);
break;
+ }
case PROP_FREQUENCY:
GST_BASE_TRANSFORM_LOCK (self);
self->frequency = g_value_get_double (value);
@@ -472,7 +478,7 @@ lpwsinc_get_property (GObject * object, guint prop_id, GValue * value,
switch (prop_id) {
case PROP_LENGTH:
- g_value_set_int (value, self->wing_size);
+ g_value_set_int (value, self->kernel_length);
break;
case PROP_FREQUENCY:
g_value_set_double (value, self->frequency);
diff --git a/gst/filter/gstlpwsinc.h b/gst/filter/gstlpwsinc.h
index d9838bd2..599bef84 100644
--- a/gst/filter/gstlpwsinc.h
+++ b/gst/filter/gstlpwsinc.h
@@ -25,10 +25,6 @@
* chapter 16
* available at http://www.dspguide.com/
*
- * FIXME:
- * - this filter is totally unoptimized !
- * - we do not destroy the allocated memory for filters and residue
- * - this might be improved upon with bytestream
*/
#ifndef __GST_LPWSINC_H__
@@ -69,8 +65,7 @@ struct _GstLPWSinc {
gint mode;
gint window;
gdouble frequency;
- gint wing_size; /* length of a "wing" of the filter;
- actual length is 2 * wing_size + 1 */
+ gint kernel_length; /* length of the filter kernel */
gdouble *residue; /* buffer for left-over samples from previous buffer */
gdouble *kernel; /* filter kernel */