From 420ada2dde4af8d71bc1d9c6988311961cf47f1b Mon Sep 17 00:00:00 2001
From: daniel aschauer <aschix@users.sourceforge.net>
Date: Fri, 9 Jan 2009 17:23:55 +0000
Subject: svn path=/trunk/externals/algocomp/; revision=10494

---
 fractal.c | 118 ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
 1 file changed, 118 insertions(+)
 create mode 100755 fractal.c

(limited to 'fractal.c')

diff --git a/fractal.c b/fractal.c
new file mode 100755
index 0000000..26c0ccc
--- /dev/null
+++ b/fractal.c
@@ -0,0 +1,118 @@
+#include "m_pd.h"
+#include <math.h>
+#include <stdio.h>
+#include <stdlib.h>
+
+static t_class *oneoverf_class;
+static t_class *chaosgame_class;
+
+
+/* 
+ * Implementation of 1/f noise sequences
+ */
+
+typedef struct _oneoverf {
+  t_object x_obj;
+//  t_int nrp;
+  t_int nrbits;
+  t_float bits;
+  t_float output;
+  t_float rans[32];
+//  t_int count[16];
+  t_int counter;
+//  t_float avg[16];  
+  t_outlet *note_out;//,*r_out,*n_out;
+} t_oneoverf;
+
+void oneoverf_bang(t_oneoverf *x)
+{
+if ((t_int) x->bits != x->nrbits) {
+post("Changed number of bits");
+x->nrbits = (t_int) x->bits;
+x->counter = 0;	
+}
+int i;
+outlet_float(x->note_out, x->output);
+//outlet_float(x->r_out, x->nrp);
+//outlet_float(x->n_out, x->nrp);
+x->output = 0; // ((t_float) rand())/RAND_MAX/x->nrbits;
+x->counter++;
+for (i=x->nrbits-1;i>=0;i--) {
+//  x->avg[i] = ((t_float) (x->avg[i]*(x->count[i]-1)+x->old))/x->count[i];
+  
+  if (x->counter%((int) pow(2,i)) == 0) {
+  if (i == x->nrbits-1) x->counter = 0;
+//  x->count[i] = 1;
+//  x->avg[i] = 0;
+  x->rans[i] = ((t_float) rand())/RAND_MAX/x->nrbits;
+  post("Changing value of %d, counter: %d",i,x->counter);
+  i = -1; // break
+  }
+}
+for (i=0;i<x->nrbits;i++) {
+  x->output += x->rans[i];
+}
+}
+
+void *oneoverf_new(t_floatarg f1)
+{
+  t_oneoverf *x = (t_oneoverf *)pd_new(oneoverf_class);
+  int i;
+  x->bits = f1;
+  x->nrbits = (t_int) f1;
+  if (f1 < 1) x->nrbits = 1;
+  else if (f1 > 32) x->nrbits = 32;
+  x->output = 0;
+  x->counter = 0;
+  for (i=0;i<x->nrbits;i++) {
+  x->rans[i] = (((t_float) rand())/RAND_MAX)/x->nrbits;
+ // x->count[i] = 1;
+//  x->avg[i] = 0;
+  x->output += x->rans[i];
+  post("Random generators: %d", x->nrbits);
+  }
+  x->note_out = outlet_new(&x->x_obj,&s_float);
+    floatinlet_new(&x->x_obj, &x->bits);  
+  //x->r_out = outlet_new(&x->x_obj,&s_float);
+  //x->n_out = outlet_new(&x->x_obj,&s_float);
+  return (void *)x;
+}
+
+/*
+ * Implementation of the chaos game!
+ */
+typedef struct _chaosgame {
+  t_object x_obj;
+  t_float x,y;
+  t_float px[3],py[3];
+  t_outlet *x_out,*y_out;
+} t_chaosgame;
+
+
+void chaosgame_bang(t_chaosgame *x)
+{	
+int i;	
+outlet_float(x->x_out, x->x);
+outlet_float(x->y_out, x->y);
+i = 3*rand()/32765;
+x->x = (x->x + x->px[i])/2;
+x->y = (x->y + x->py[i])/2;
+}
+
+void *chaosgame_new(t_floatarg f1)
+{
+  t_chaosgame *x = (t_chaosgame *)pd_new(chaosgame_class);
+	x->px[0] = -1;
+    x->py[0] = -1;
+    x->px[1] = 0;
+    x->py[1] = 1;
+    x->px[2] = 1;
+    x->py[2] = -1;
+
+	x->x = rand()/RAND_MAX;
+	x->y = rand()/RAND_MAX;
+	
+	x->x_out = outlet_new(&x->x_obj,&s_float);
+  	x->y_out = outlet_new(&x->x_obj,&s_float);
+  	return (void *)x;
+}
-- 
cgit v1.2.1