aboutsummaryrefslogtreecommitdiff
path: root/modules/blocknorm~.c
diff options
context:
space:
mode:
authorHans-Christoph Steiner <eighthave@users.sourceforge.net>2006-05-25 16:40:19 +0000
committerHans-Christoph Steiner <eighthave@users.sourceforge.net>2006-05-25 16:40:19 +0000
commitd962ed33f6b5fe3040cfbd798f63ab9aa36d1f9e (patch)
treef700234efa8d57c28ce9be56e0db5f6526921fbe /modules/blocknorm~.c
parent2994949cca515f63fb82fa1d114ce07eb213220d (diff)
renamed files to match their class names
svn path=/trunk/externals/creb/; revision=5127
Diffstat (limited to 'modules/blocknorm~.c')
-rw-r--r--modules/blocknorm~.c129
1 files changed, 129 insertions, 0 deletions
diff --git a/modules/blocknorm~.c b/modules/blocknorm~.c
new file mode 100644
index 0000000..bcf695a
--- /dev/null
+++ b/modules/blocknorm~.c
@@ -0,0 +1,129 @@
+/*
+ * blocknorm.c - normalize an array of dsp blocks (for spectral processing)
+ * Copyright (c) 2000-2003 by Tom Schouten
+ *
+ * This program is free software; you can redistribute it and/or modify
+ * it under the terms of the GNU General Public License as published by
+ * the Free Software Foundation; either version 2 of the License, or
+ * (at your option) any later version.
+ *
+ * This program is distributed in the hope that it will be useful,
+ * but WITHOUT ANY WARRANTY; without even the implied warranty of
+ * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
+ * GNU General Public License for more details.
+ *
+ * You should have received a copy of the GNU General Public License
+ * along with this program; if not, write to the Free Software
+ * Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
+ */
+
+#include "m_pd.h"
+#include <math.h>
+#include <stdlib.h>
+#include <stdio.h>
+#include <math.h>
+
+#define MAXCHANNELS 32
+
+typedef struct blocknormctl
+{
+ t_int c_channels;
+ t_float **c_input;
+ t_float **c_output;
+} t_blocknormctl;
+
+typedef struct blocknorm
+{
+ t_object x_obj;
+ t_float x_f;
+ t_blocknormctl x_ctl;
+} t_blocknorm;
+
+
+static t_int *blocknorm_perform(t_int *word)
+{
+
+ t_blocknormctl *ctl = (t_blocknormctl *)(word[1]);
+ t_int n = (t_int)(word[2]);
+ t_float **in = ctl->c_input;
+ t_float **out = ctl->c_output;
+ t_int c = ctl->c_channels;
+ t_int i,j;
+
+ t_float p = 0.0f;
+ t_float x, s;
+
+ /* get power */
+ for (j=0;j<c;j++){
+ for (i=0;i<n;i++){
+ x = in[j][i];
+ p += x*x;
+ }
+ }
+
+ /* compute normalization */
+ if (p == 0.0f) s = 1.0f;
+ else s =sqrt(((float)(c * n)) / p);
+
+ /* normalize */
+ for (j=0;j<c;j++){
+ for (i=0;i<n;i++){
+ out[j][i] *= s;
+ }
+ }
+
+
+
+ return (word+3);
+}
+
+
+
+static void blocknorm_dsp(t_blocknorm *x, t_signal **sp)
+{
+
+ int i;
+ int c = x->x_ctl.c_channels;
+ for (i=0;i<c;i++){
+ x->x_ctl.c_input[i] = sp[i]->s_vec;
+ x->x_ctl.c_output[i] = sp[c+i]->s_vec;
+ }
+ dsp_add(blocknorm_perform, 2, &x->x_ctl, sp[0]->s_n);
+}
+
+
+static void blocknorm_free(t_blocknorm *x)
+{
+ free (x->x_ctl.c_output);
+ free (x->x_ctl.c_input);
+}
+
+t_class *blocknorm_class;
+
+static void *blocknorm_new(t_floatarg channels)
+{
+ int i = (int)channels;
+ int j;
+ t_blocknorm *x = (t_blocknorm *)pd_new(blocknorm_class);
+
+ if (i<1) i = 1;
+ if (i>MAXCHANNELS) i = MAXCHANNELS;
+ x->x_ctl.c_channels = i;
+ x->x_ctl.c_input = malloc(sizeof(float)*i);
+ x->x_ctl.c_output = malloc(sizeof(float)*i);
+
+ j = i;
+ while (--j) inlet_new(&x->x_obj, &x->x_obj.ob_pd, gensym("signal"), gensym("signal"));
+ while (i--) outlet_new(&x->x_obj, gensym("signal"));
+
+ return (void *)x;
+}
+
+void blocknorm_tilde_setup(void)
+{
+ blocknorm_class = class_new(gensym("blocknorm~"), (t_newmethod)blocknorm_new,
+ (t_method)blocknorm_free, sizeof(t_blocknorm), 0, A_DEFFLOAT, 0);
+ CLASS_MAINSIGNALIN(blocknorm_class, t_blocknorm, x_f);
+ class_addmethod(blocknorm_class, (t_method)blocknorm_dsp, gensym("dsp"), 0);
+}
+