From 541ab13ab9aaf849014d81d9159bd12a03320c74 Mon Sep 17 00:00:00 2001 From: musil Date: Wed, 8 Nov 2006 13:24:10 +0000 Subject: changed sig* to *_tilde #if MSC_VER is obsolete replaced all float by t_float svn path=/trunk/externals/iemlib/; revision=6231 --- src/iemlib1/prvu~.c | 274 ++++++++++++++++++++++++++++++++++++++++++++++++++++ 1 file changed, 274 insertions(+) create mode 100644 src/iemlib1/prvu~.c (limited to 'src/iemlib1/prvu~.c') diff --git a/src/iemlib1/prvu~.c b/src/iemlib1/prvu~.c new file mode 100644 index 0000000..17fe143 --- /dev/null +++ b/src/iemlib1/prvu~.c @@ -0,0 +1,274 @@ +/* For information on usage and redistribution, and for a DISCLAIMER OF ALL +* WARRANTIES, see the file, "LICENSE.txt," in this distribution. + +iemlib1 written by Thomas Musil, Copyright (c) IEM KUG Graz Austria 2000 - 2005 */ + +#include "m_pd.h" +#include "iemlib.h" +#include + +/* ---------------- prvu~ - simple peak&rms-vu-meter. ----------------- */ + +typedef struct _prvu_tilde +{ + t_object x_obj; + t_atom x_at[3]; + void *x_clock_metro; + t_float x_metro_time; + void *x_clock_hold; + t_float x_hold_time; + t_float x_cur_peak; + t_float x_old_peak; + t_float x_hold_peak; + int x_hold; + t_float x_sum_rms; + t_float x_old_rms; + t_float x_rcp; + t_float x_sr; + t_float x_threshold_over; + int x_overflow_counter; + t_float x_release_time; + t_float x_c1; + int x_started; + t_float x_msi; +} t_prvu_tilde; + +t_class *prvu_tilde_class; +static void prvu_tilde_tick_metro(t_prvu_tilde *x); +static void prvu_tilde_tick_hold(t_prvu_tilde *x); + +static void prvu_tilde_reset(t_prvu_tilde *x) +{ + x->x_at[0].a_w.w_float = -99.9f; + x->x_at[1].a_w.w_float = -99.9f; + x->x_at[2].a_w.w_float = 0.0f; + outlet_list(x->x_obj.ob_outlet, &s_list, 3, x->x_at); + x->x_overflow_counter = 0; + x->x_cur_peak = 0.0f; + x->x_old_peak = 0.0f; + x->x_hold_peak = 0.0f; + x->x_sum_rms = 0.0f; + x->x_old_rms = 0.0f; + x->x_hold = 0; + clock_unset(x->x_clock_hold); + clock_delay(x->x_clock_metro, x->x_metro_time); +} + +static void prvu_tilde_stop(t_prvu_tilde *x) +{ + clock_unset(x->x_clock_metro); + x->x_started = 0; +} + +static void prvu_tilde_start(t_prvu_tilde *x) +{ + clock_delay(x->x_clock_metro, x->x_metro_time); + x->x_started = 1; +} + +static void prvu_tilde_float(t_prvu_tilde *x, t_floatarg f) +{ + if(f == 0.0) + { + clock_unset(x->x_clock_metro); + x->x_started = 0; + } + else + { + clock_delay(x->x_clock_metro, x->x_metro_time); + x->x_started = 1; + } +} + +static void prvu_tilde_t_release(t_prvu_tilde *x, t_floatarg release_time) +{ + if(release_time <= 5.0f) + release_time = 5.0f; + x->x_release_time = release_time; + x->x_c1 = exp(-2.0f*x->x_metro_time/x->x_release_time); +} + +static void prvu_tilde_t_metro(t_prvu_tilde *x, t_floatarg metro_time) +{ + if(metro_time <= 5.0f) + metro_time = 5.0f; + x->x_metro_time = metro_time; + x->x_c1 = exp(-2.0f*x->x_metro_time/x->x_release_time); + x->x_rcp = 1.0f/(x->x_sr*(t_float)x->x_metro_time); +} + +static void prvu_tilde_t_hold(t_prvu_tilde *x, t_floatarg hold_time) +{ + if(hold_time <= 5.0f) + hold_time = 5.0f; + x->x_hold_time = hold_time; +} + +static void prvu_tilde_threshold(t_prvu_tilde *x, t_floatarg thresh) +{ + x->x_threshold_over = thresh; +} + +static t_int *prvu_tilde_perform(t_int *w) +{ + t_float *in = (t_float *)(w[1]); + t_prvu_tilde *x = (t_prvu_tilde *)(w[2]); + int n = (int)(w[3]); + t_float peak = x->x_cur_peak, pow, sum=x->x_sum_rms; + int i; + + if(x->x_started) + { + for(i=0; i peak) + peak = pow; + sum += pow; + } + x->x_cur_peak = peak; + x->x_sum_rms = sum; + } + return(w+4); +} + +static void prvu_tilde_dsp(t_prvu_tilde *x, t_signal **sp) +{ + x->x_sr = 0.001*(t_float)sp[0]->s_sr; + x->x_rcp = 1.0f/(x->x_sr*x->x_metro_time); + dsp_add(prvu_tilde_perform, 3, sp[0]->s_vec, x, sp[0]->s_n); + clock_delay(x->x_clock_metro, x->x_metro_time); +} + +static void prvu_tilde_tick_hold(t_prvu_tilde *x) +{ + x->x_hold = 0; + x->x_hold_peak = x->x_old_peak; +} + +static void prvu_tilde_tick_metro(t_prvu_tilde *x) +{ + t_float dbr, dbp, cur_rms, c1=x->x_c1; + + x->x_old_peak *= c1; + /* NAN protect */ + if(IEM_DENORMAL(x->x_old_peak)) + x->x_old_peak = 0.0f; + + if(x->x_cur_peak > x->x_old_peak) + x->x_old_peak = x->x_cur_peak; + if(x->x_old_peak > x->x_hold_peak) + { + x->x_hold = 1; + x->x_hold_peak = x->x_old_peak; + clock_delay(x->x_clock_hold, x->x_hold_time); + } + if(!x->x_hold) + x->x_hold_peak = x->x_old_peak; + if(x->x_hold_peak <= 0.0000000001f) + dbp = -99.9f; + else if(x->x_hold_peak > 1000000.0f) + { + dbp = 60.0f; + x->x_hold_peak = 1000000.0f; + x->x_old_peak = 1000000.0f; + } + else + dbp = 4.3429448195f*log(x->x_hold_peak); + x->x_cur_peak = 0.0f; + if(dbp >= x->x_threshold_over) + x->x_overflow_counter++; + x->x_at[1].a_w.w_float = dbp; + x->x_at[2].a_w.w_float = (t_float)x->x_overflow_counter; + + cur_rms = (1.0f - c1)*x->x_sum_rms*x->x_rcp + c1*x->x_old_rms; + /* NAN protect */ + if(IEM_DENORMAL(cur_rms)) + cur_rms = 0.0f; + + if(cur_rms <= 0.0000000001f) + dbr = -99.9f; + else if(cur_rms > 1000000.0f) + { + dbr = 60.0f; + x->x_old_rms = 1000000.0f; + } + else + dbr = 4.3429448195f*log(cur_rms); + x->x_sum_rms = 0.0f; + x->x_old_rms = cur_rms; + x->x_at[0].a_w.w_float = dbr; + outlet_list(x->x_obj.ob_outlet, &s_list, 3, x->x_at); + clock_delay(x->x_clock_metro, x->x_metro_time); +} + +static void prvu_tilde_ff(t_prvu_tilde *x) +{ + clock_free(x->x_clock_metro); + clock_free(x->x_clock_hold); +} + +static void *prvu_tilde_new(t_floatarg metro_time, t_floatarg hold_time, + t_floatarg release_time, t_floatarg threshold) +{ + t_prvu_tilde *x; + t_float t; + int i; + + x = (t_prvu_tilde *)pd_new(prvu_tilde_class); + if(metro_time <= 0.0f) + metro_time = 300.0f; + if(metro_time <= 5.0f) + metro_time = 5.0f; + if(release_time <= 0.0f) + release_time = 300.0f; + if(release_time <= 5.0f) + release_time = 5.0f; + if(hold_time <= 0.0f) + hold_time = 1000.0f; + if(hold_time <= 5.0f) + hold_time = 5.0f; + if(threshold == 0.0f) + threshold = -0.01f; + x->x_metro_time = metro_time; + x->x_release_time = release_time; + x->x_hold_time = hold_time; + x->x_threshold_over = threshold; + x->x_c1 = exp(-2.0f*x->x_metro_time/x->x_release_time); + x->x_cur_peak = 0.0f; + x->x_old_peak = 0.0f; + x->x_hold_peak = 0.0f; + x->x_hold = 0; + x->x_sum_rms = 0.0f; + x->x_old_rms = 0.0f; + x->x_sr = 44.1f; + x->x_rcp = 1.0f/(x->x_sr*x->x_metro_time); + x->x_overflow_counter = 0; + x->x_clock_metro = clock_new(x, (t_method)prvu_tilde_tick_metro); + x->x_clock_hold = clock_new(x, (t_method)prvu_tilde_tick_hold); + x->x_started = 1; + outlet_new(&x->x_obj, &s_list); + x->x_at[0].a_type = A_FLOAT; + x->x_at[1].a_type = A_FLOAT; + x->x_at[2].a_type = A_FLOAT; + x->x_msi = 0.0f; + return(x); +} + +void prvu_tilde_setup(void) +{ + prvu_tilde_class = class_new(gensym("prvu~"), (t_newmethod)prvu_tilde_new, + (t_method)prvu_tilde_ff, sizeof(t_prvu_tilde), 0, + A_DEFFLOAT, A_DEFFLOAT, A_DEFFLOAT, A_DEFFLOAT, 0); + CLASS_MAINSIGNALIN(prvu_tilde_class, t_prvu_tilde, x_msi); + class_addmethod(prvu_tilde_class, (t_method)prvu_tilde_dsp, gensym("dsp"), 0); + class_addfloat(prvu_tilde_class, prvu_tilde_float); + class_addmethod(prvu_tilde_class, (t_method)prvu_tilde_reset, gensym("reset"), 0); + class_addmethod(prvu_tilde_class, (t_method)prvu_tilde_start, gensym("start"), 0); + class_addmethod(prvu_tilde_class, (t_method)prvu_tilde_stop, gensym("stop"), 0); + class_addmethod(prvu_tilde_class, (t_method)prvu_tilde_t_release, gensym("t_release"), A_FLOAT, 0); + class_addmethod(prvu_tilde_class, (t_method)prvu_tilde_t_metro, gensym("t_metro"), A_FLOAT, 0); + class_addmethod(prvu_tilde_class, (t_method)prvu_tilde_t_hold, gensym("t_hold"), A_FLOAT, 0); + class_addmethod(prvu_tilde_class, (t_method)prvu_tilde_threshold, gensym("threshold"), A_FLOAT, 0); + class_sethelpsymbol(prvu_tilde_class, gensym("iemhelp/help-prvu~")); +} -- cgit v1.2.1