1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
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);
}
|
