diff options
author | Tom Schouten <doelie@users.sourceforge.net> | 2003-09-12 22:26:57 +0000 |
---|---|---|
committer | Tom Schouten <doelie@users.sourceforge.net> | 2003-09-12 22:26:57 +0000 |
commit | 494a07a361fe4ee0e54f77468a976b1a77818770 (patch) | |
tree | b784b32b2a29e379c6ae5947cdcdcdaadde72b5e /modules/blocknorm.c | |
parent | 9f8757ccf742d984bd1c1a5d1e5906e9fd50784e (diff) |
creb 0.9.0
svn path=/trunk/externals/creb/; revision=956
Diffstat (limited to 'modules/blocknorm.c')
-rw-r--r-- | modules/blocknorm.c | 129 |
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); +} + |