Added ihisto Integer histogram external

svn path=/trunk/externals/ekext/; revision=14775

1 files changed, 147 insertions, 0 deletions

diff --git a/ihisto.c b/ihisto.c
new file mode 100644
index 0000000..07cf95a
--- /dev/null
+++ b/ihisto.c
@@ -0,0 +1,147 @@
+/* Histogram for integers - copyright 2011 Edward Kelly */
+/* released under the GNU GPL license */
+#include "m_pd.h"
+
+#define MAXENTRIES 13823
+#define LASTENTRY 13822
+
+static t_class *ihisto_class;
+
+/* mode = 0: most recent when equal, 1: highest, 2: lowest */
+
+typedef struct _hist
+{
+  t_atom histogram[MAXENTRIES];
+} t_hist;
+
+typedef struct _ihisto
+{
+  t_object x_obj;
+  t_hist histog;
+  int intinput, size;
+  t_float max, maxvalue, count, exclude, mode;
+  t_outlet *peak, *counted, *current, *counting, *histo;
+} t_ihisto;
+
+void ihisto_float(t_ihisto *x, t_floatarg fin)
+{
+  x->intinput = fin > 0 ? (int)fin < x->size ? (int)fin : x->size - 1 : 0;
+  x->count = atom_getfloatarg(x->intinput, MAXENTRIES, x->histog.histogram);
+  x->count++;
+  SETFLOAT(&x->histog.histogram[x->intinput],x->count);
+
+  if(x->mode == 0)
+    {
+      if(x->count >= x->maxvalue && x->intinput != (int)x->exclude)
+	{
+	  x->max = (float)x->intinput;
+	  x->maxvalue = x->count;
+	}
+    }
+  else if (x->mode == 1)
+    {
+      if(x->count >= x->maxvalue && x->intinput >= (int)x->max && x->intinput != (int)x->exclude || x->count > x->maxvalue && x->intinput != (int)x->exclude)
+	{
+	  x->max = (float)x->intinput;
+	  x->maxvalue = x->count;
+	}
+    }
+  else if (x->mode == 2)
+    {
+      if(x->count >= x->maxvalue && x->intinput <= (int)x->max && x->intinput != (int)x->exclude || x->count > x->maxvalue && x->intinput != (int)x->exclude)
+	{
+	  x->max = (float)x->intinput;
+	  x->maxvalue = x->count;
+	}
+    }
+
+  outlet_float(x->current, x->max);
+  outlet_float(x->counting, x->maxvalue);
+}
+
+void ihisto_bang(t_ihisto *x)
+{
+  if(x->max)
+    {
+      outlet_float(x->counted, x->maxvalue);
+      outlet_float(x->peak, x->max);
+    }
+
+  int i;
+  for(i=0;i<x->size;i++)
+    {
+      SETFLOAT(x->histog.histogram+i, 0);
+    }
+
+  x->max = -1;
+  x->maxvalue = 0;
+  x->count = 0;
+  x->intinput = -1;
+  x->exclude = -1;
+}
+
+void ihisto_poll(t_ihisto *x)
+{
+  if(x->max)
+    {
+      outlet_float(x->counted, x->maxvalue);
+      outlet_float(x->peak, x->max);
+    }
+}
+
+void ihisto_mode(t_ihisto *x, t_floatarg f)
+{
+  x->mode = f < 0 ? 0 : f > 2 ? 2 : (int)f;
+}
+
+void ihisto_get(t_ihisto *x)
+{
+  outlet_list(x->histo, gensym("list"), x->size, x->histog.histogram);
+}
+
+void ihisto_exclude(t_ihisto *x, t_floatarg f)
+{
+  x->exclude = f;
+}
+
+void *ihisto_new(t_floatarg f)
+{
+  t_ihisto *x = (t_ihisto *)pd_new(ihisto_class);
+  int init = (int)f;
+  x->size = init > MAXENTRIES ? MAXENTRIES : init <= 0 ? 128 : init;
+
+  x->max = -1;
+  x->maxvalue = 0;
+  x->count = 0;
+  x->intinput = 0;
+  x->exclude = -1;
+
+  int i;
+  for(i=0;i<MAXENTRIES;i++) 
+    {
+      SETFLOAT(x->histog.histogram+i, 0);
+    }
+  x->peak = outlet_new(&x->x_obj, &s_float);
+  x->counted = outlet_new(&x->x_obj, &s_float);
+  x->current = outlet_new(&x->x_obj, &s_float);
+  x->counting = outlet_new(&x->x_obj, &s_float);
+  x->histo = outlet_new(&x->x_obj, &s_list);
+
+  return (void *)x;
+}
+
+void ihisto_setup(void) 
+{
+  ihisto_class = class_new(gensym("ihisto"),
+  (t_newmethod)ihisto_new,
+  0, sizeof(t_ihisto),
+  0, A_DEFFLOAT, 0);
+  post("Integer Histogram by Ed Kelly 2011");
+
+  class_addfloat(ihisto_class, ihisto_float);
+  class_addbang(ihisto_class, ihisto_bang);
+  class_addmethod(ihisto_class, (t_method)ihisto_poll, gensym("poll"), A_DEFFLOAT, 0);
+  class_addmethod(ihisto_class, (t_method)ihisto_mode, gensym("mode"), A_DEFFLOAT, 0);
+  class_addmethod(ihisto_class, (t_method)ihisto_get, gensym("get"), A_GIMME, 0);
+  class_addmethod(ihisto_class, (t_method)ihisto_exclude, gensym("exclude"), A_DEFFLOAT, 0);
+}