From 086e6b358ad45bda0e6a948cfb7ea2d33ab7da5a Mon Sep 17 00:00:00 2001
From: Hans-Christoph Steiner <eighthave@users.sourceforge.net>
Date: Wed, 10 Feb 2010 00:34:05 +0000
Subject: re-arranged into libdir and setup with the template Makefile+debian
 stuff. renabled vv+ and vv- as vvplus and vvminus

svn path=/trunk/externals/smlib/; revision=13158
---
 lhist.c | 174 ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
 1 file changed, 174 insertions(+)
 create mode 100644 lhist.c

(limited to 'lhist.c')

diff --git a/lhist.c b/lhist.c
new file mode 100644
index 0000000..ea20edc
--- /dev/null
+++ b/lhist.c
@@ -0,0 +1,174 @@
+#include "defines.h"
+
+/*--------------- lhist ---------------*/
+
+static t_class *lhist_class;
+
+typedef struct _lhist
+{
+    t_object x_obj;
+	float m_lo;
+	float m_hi;
+	float m_scale;
+	float m_c_leak;
+	float m_leak;
+	int m_nbins;
+//	int m_n_observations;
+	float *m_lhist;
+} t_lhist;
+
+static void lhist_setHalfDecay(t_lhist *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 lhist_perform_float(t_lhist *x, t_float f)
+{
+	int j;
+	j=(int)(.5+(f-x->m_lo)*x->m_scale);
+	j=(j>0)?(j<x->m_nbins?j:x->m_nbins-1):0; // limit without IF
+	x->m_lhist[j]++;
+//	x->m_n_observations++;
+}
+
+static void lhist_perform_list(t_lhist *x, t_symbol *s, int argc, t_atom *argv)
+{
+	int i,j;
+	for (i = 0; i < argc; i++)
+	{
+		j=(int)(.5f+(atom_getfloat(&argv[i])-x->m_lo)*x->m_scale);
+		j=(j>0)?(j<x->m_nbins?j:x->m_nbins-1):0; // limit without IF
+		x->m_lhist[j]++;
+	}
+//	x->m_n_observations+=argc;
+}
+
+static void lhist_leak(t_lhist *x)
+{
+	int i;
+	float *f;
+	float sc;
+	f=x->m_lhist;
+    sc=x->m_c_leak;
+	i=x->m_nbins;
+	while(i--)
+      *f++*=sc;
+}
+
+static void lhist_bang(t_lhist *x)
+{
+	int i,n;
+	float *f;
+	t_atom *ap,*app;
+	n=x->m_nbins;
+    ap = (t_atom *)getbytes(sizeof(t_atom)*n);
+	app=ap;
+	i=n;
+	f=x->m_lhist;
+
+    while(i--){
+      SETFLOAT(app, *f);
+	  f++;
+      app++;
+    }
+	outlet_list(x->x_obj.ob_outlet,gensym("list"),n,ap);
+	freebytes(ap,n);
+}
+
+static void lhist_relative(t_lhist *x)
+{
+	int i,n;
+	float *f;
+	float invn,sum;
+	t_atom *ap,*app;
+
+	n=x->m_nbins;
+    ap = (t_atom *)getbytes(sizeof(t_atom)*n);
+	app=ap;
+	i=x->m_nbins;
+	f=x->m_lhist;
+	sum=0.0f;
+    while(i--) sum+=*f++;
+	invn=1.0f/(1e-10f+sum);
+
+	i=x->m_nbins;
+	f=x->m_lhist;
+
+    while(i--){
+      SETFLOAT(app, (*f*invn));
+	  f++;
+      app++;
+    }
+	outlet_list(x->x_obj.ob_outlet,gensym("list"),n,ap);
+	freebytes(ap,n);
+}
+
+static void lhist_clear(t_lhist *x)
+{
+	int i;
+	float *f;
+	f=x->m_lhist;
+	for (i=0;i<x->m_nbins;i++)
+		*f++=0.0f;
+//	x->m_n_observations=0;
+}
+
+static void lhist_set(t_lhist *x, t_float lo, t_float hi, t_float nbins)
+{
+	if (nbins<1)
+	{
+		nbins=1;
+		post("lhist: number of bins is minimum 1...");
+	}
+	if (hi<=lo)
+	{
+		post("lhist: higher bound must be higher than lower bound...");	
+		hi=lo+1.0f;
+	}
+	freebytes(x->m_lhist, x->m_nbins);
+	x->m_hi=hi;
+	x->m_lo=lo;
+	x->m_nbins=(int)nbins;
+	x->m_scale=(float)x->m_nbins/(hi-lo);
+    x->m_lhist = (float*)getbytes(sizeof(float)*x->m_nbins);
+
+	lhist_clear(x);
+}
+
+static void *lhist_new(t_float lo, t_float hi, t_float nbins, t_float decay)
+{
+	t_lhist *x=(t_lhist *)pd_new(lhist_class);
+	outlet_new(&x->x_obj, gensym("list"));
+	lhist_setHalfDecay(x,decay);
+	x->m_nbins=0;
+	x->m_lhist=0;
+	lhist_set(x, lo, hi, nbins);
+	return (void *)x;
+}
+
+static void lhist_free(t_lhist *x)
+{
+	freebytes(x->m_lhist, x->m_nbins);
+}
+
+void lhist_setup(void)
+{
+    lhist_class = class_new(gensym("lhist"),
+    	(t_newmethod)lhist_new, (t_method)lhist_free,
+		sizeof(t_lhist), 
+		CLASS_DEFAULT,
+	    A_DEFFLOAT, A_DEFFLOAT, A_DEFFLOAT,A_DEFFLOAT,0);
+
+	class_addmethod(lhist_class, (t_method)lhist_clear, gensym("clear"),0);
+	class_addmethod(lhist_class, (t_method)lhist_bang, gensym("absolute"),0);
+	class_addmethod(lhist_class, (t_method)lhist_relative, gensym("relative"),0);
+	class_addmethod(lhist_class, (t_method)lhist_leak, gensym("leak"),0);
+	class_addmethod(lhist_class, (t_method)lhist_setHalfDecay,
+    	gensym("decay"), A_DEFFLOAT, NULL);
+
+    class_addlist(lhist_class, (t_method)lhist_perform_list);
+    class_addfloat(lhist_class, (t_method)lhist_perform_float);
+    class_addbang(lhist_class, (t_method)lhist_bang);
+}
+
-- 
cgit v1.2.1