From da1f0065d07250b86b9a761b012987e4825f491f Mon Sep 17 00:00:00 2001 From: Hans-Christoph Steiner Date: Mon, 28 Nov 2005 16:53:25 +0000 Subject: This commit was generated by cvs2svn to compensate for changes in r4068, which included commits to RCS files with non-trunk default branches. svn path=/trunk/externals/smlib/; revision=4069 --- source/lstd.c | 60 +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 1 file changed, 60 insertions(+) create mode 100644 source/lstd.c (limited to 'source/lstd.c') diff --git a/source/lstd.c b/source/lstd.c new file mode 100644 index 0000000..1873bda --- /dev/null +++ b/source/lstd.c @@ -0,0 +1,60 @@ +#include "defines.h" + +/*--------------- lstd ---------------*/ + +static t_class *lstd_class; + +typedef struct _lstd +{ + t_object x_obj; + float m_avg; + float m_sum_squares; + float m_std; + float m_c_leak; + float m_leak; +} t_lstd; + + +static void lstd_perform(t_lstd *x, t_float in) +{ + float tmp=x->m_avg-in; + x->m_avg= x->m_avg * x->m_c_leak + in * x->m_leak; + x->m_sum_squares=x->m_sum_squares * x->m_c_leak + x->m_leak*tmp*tmp; + x->m_std=(float)sqrtf(x->m_sum_squares); + outlet_float(x->x_obj.ob_outlet, x->m_std); +} + +static void lstd_setHalfDecay(t_lstd *x, t_float halfDecayTime) +{ + x->m_c_leak=(float)powf(.5,(1.0f/halfDecayTime)); + x->m_leak=1.0f-x->m_c_leak; +} + +static void lstd_clear(t_lstd *x) +{ + x->m_sum_squares=0.0f; + x->m_avg=0.0f; +} + +static void *lstd_new( t_float halfDecayTime) +{ + t_lstd *x=(t_lstd *)pd_new(lstd_class); + outlet_new(&x->x_obj, gensym("float")); + + lstd_setHalfDecay(x, halfDecayTime); + return (void *)x; +} + + +void lstd_setup(void) +{ + lstd_class = class_new(gensym("lstd"), + (t_newmethod)lstd_new, 0, + sizeof(t_lstd), + CLASS_DEFAULT, + A_DEFFLOAT, 0); + class_addfloat(lstd_class, (t_method)lstd_perform); + class_addmethod(lstd_class, (t_method)lstd_setHalfDecay, + gensym("decay"), A_DEFFLOAT, NULL); +} + -- cgit v1.2.1