From da1f0065d07250b86b9a761b012987e4825f491f Mon Sep 17 00:00:00 2001
From: Hans-Christoph Steiner <eighthave@users.sourceforge.net>
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/lhisti.c | 174 ++++++++++++++++++++++++++++++++++++++++++++++++++++++++
 1 file changed, 174 insertions(+)
 create mode 100644 source/lhisti.c

(limited to 'source/lhisti.c')

diff --git a/source/lhisti.c b/source/lhisti.c
new file mode 100644
index 0000000..826cce4
--- /dev/null
+++ b/source/lhisti.c
@@ -0,0 +1,174 @@
+#include "defines.h"
+
+/*--------------- lhisti ---------------*/
+// ignores samples outside bins
+
+static t_class *lhisti_class;
+
+typedef struct _lhisti
+{
+    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_lhisti;
+} t_lhisti;
+
+static void lhisti_setHalfDecay(t_lhisti *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 lhisti_perform_float(t_lhisti *x, t_float f)
+{
+	int j;
+	j=(int)(.5+(f-x->m_lo)*x->m_scale);
+	if ((j>=0)&&(j<x->m_nbins))
+		x->m_lhisti[j]++;
+}
+
+static void lhisti_perform_list(t_lhisti *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);
+		if ((j>=0)&&(j<x->m_nbins))
+			x->m_lhisti[j]++;
+	}
+//	x->m_n_observations+=argc;
+}
+
+static void lhisti_leak(t_lhisti *x)
+{
+	int i;
+	float *f;
+	float sc;
+	f=x->m_lhisti;
+    sc=x->m_c_leak;
+	i=x->m_nbins;
+	while(i--)
+      *f++*=sc;
+}
+
+static void lhisti_bang(t_lhisti *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_lhisti;
+
+    while(i--){
+      SETFLOAT(app, *f);
+	  f++;
+      app++;
+    }
+	outlet_list(x->x_obj.ob_outlet,gensym("list"),n,ap);
+	freebytes(ap,n);
+}
+
+static void lhisti_relative(t_lhisti *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_lhisti;
+	sum=0.0f;
+    while(i--) sum+=*f++;
+	invn=1.0f/(1e-10f+sum);
+
+	i=x->m_nbins;
+	f=x->m_lhisti;
+
+    while(i--){
+      SETFLOAT(app, (*f*invn));
+	  f++;
+      app++;
+    }
+	outlet_list(x->x_obj.ob_outlet,gensym("list"),n,ap);
+	freebytes(ap,n);
+}
+
+static void lhisti_clear(t_lhisti *x)
+{
+	int i;
+	float *f;
+	f=x->m_lhisti;
+	for (i=0;i<x->m_nbins;i++)
+		*f++=0.0f;
+//	x->m_n_observations=0;
+}
+
+static void lhisti_set(t_lhisti *x, t_float lo, t_float hi, t_float nbins)
+{
+	if (nbins<1)
+	{
+		nbins=1;
+		post("lhisti: number of bins is minimum 1...");
+	}
+	if (hi<=lo)
+	{
+		post("lhisti: higher bound must be higher than lower bound...");	
+		hi=lo+1.0f;
+	}
+	freebytes(x->m_lhisti, 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_lhisti = (float*)getbytes(sizeof(float)*x->m_nbins);
+
+	lhisti_clear(x);
+}
+
+static void *lhisti_new(t_float lo, t_float hi, t_float nbins, t_float decay)
+{
+	t_lhisti *x=(t_lhisti *)pd_new(lhisti_class);
+	outlet_new(&x->x_obj, gensym("list"));
+	lhisti_setHalfDecay(x,decay);
+	x->m_nbins=0;
+	x->m_lhisti=0;
+	lhisti_set(x, lo, hi, nbins);
+	return (void *)x;
+}
+
+static void lhisti_free(t_lhisti *x)
+{
+	freebytes(x->m_lhisti, x->m_nbins);
+}
+
+void lhisti_setup(void)
+{
+    lhisti_class = class_new(gensym("lhisti"),
+    	(t_newmethod)lhisti_new, (t_method)lhisti_free,
+		sizeof(t_lhisti), 
+		CLASS_DEFAULT,
+	    A_DEFFLOAT, A_DEFFLOAT, A_DEFFLOAT,A_DEFFLOAT,0);
+
+	class_addmethod(lhisti_class, (t_method)lhisti_clear, gensym("clear"),0);
+	class_addmethod(lhisti_class, (t_method)lhisti_bang, gensym("absolute"),0);
+	class_addmethod(lhisti_class, (t_method)lhisti_relative, gensym("relative"),0);
+	class_addmethod(lhisti_class, (t_method)lhisti_leak, gensym("leak"),0);
+	class_addmethod(lhisti_class, (t_method)lhisti_setHalfDecay,
+    	gensym("decay"), A_DEFFLOAT, NULL);
+
+    class_addlist(lhisti_class, (t_method)lhisti_perform_list);
+    class_addfloat(lhisti_class, (t_method)lhisti_perform_float);
+    class_addbang(lhisti_class, (t_method)lhisti_bang);
+}
+
-- 
cgit v1.2.1