From cdf8ded57019d5c905f16422d40be7b1a18ab3bc Mon Sep 17 00:00:00 2001
From: Tom Schouten <doelie@users.sourceforge.net>
Date: Tue, 21 Jan 2003 10:18:19 +0000
Subject: This commit was generated by cvs2svn to compensate for changes in
 r350, which included commits to RCS files with non-trunk default branches.

svn path=/trunk/externals/creb/; revision=351
---
 modules/permut.c | 177 +++++++++++++++++++++++++++++++++++++++++++++++++++++++
 1 file changed, 177 insertions(+)
 create mode 100644 modules/permut.c

(limited to 'modules/permut.c')

diff --git a/modules/permut.c b/modules/permut.c
new file mode 100644
index 0000000..bc143d2
--- /dev/null
+++ b/modules/permut.c
@@ -0,0 +1,177 @@
+/*
+ *   permut.c  - applies a (random) permutation on a signal block
+ *   intended for spectral processing, dynwav
+ *
+ *   Copyright (c) 2000-2003 by Tom Schouten
+ *
+ *   This program is free software; you can redistribute it and/or modify
+ *   it under the terms of the GNU General Public License as published by
+ *   the Free Software Foundation; either version 2 of the License, or
+ *   (at your option) any later version.
+ *
+ *   This program is distributed in the hope that it will be useful,
+ *   but WITHOUT ANY WARRANTY; without even the implied warranty of
+ *   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
+ *   GNU General Public License for more details.
+ *
+ *   You should have received a copy of the GNU General Public License
+ *   along with this program; if not, write to the Free Software
+ *   Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
+ */
+
+
+#include "m_pd.h"
+#include <math.h>
+#include <stdlib.h>
+                               
+
+
+typedef struct permutctl
+{
+  char c_type;
+  int *c_permutationtable;
+  int c_blocksize;
+} t_permutctl;
+
+typedef struct permut
+{
+  t_object x_obj;
+  t_float x_f;
+  t_permutctl x_ctl;
+} t_permut;
+
+
+static inline void permut_perform_permutation(t_float *S, int n, t_int *f)
+{
+  t_int k,l;
+  t_float swap;
+  for(k=0; k<n; k++)
+    {
+      l = f[k];
+      while (l<k) l = f[l];
+      swap = S[k];
+      S[k] = S[l];
+      S[l] = swap;
+    }
+}
+
+
+static void permut_random(t_permut *x, t_floatarg seed)
+{
+  int i,j;
+  int N = x->x_ctl.c_blocksize;
+  int mask = N-1;
+  int *p = x->x_ctl.c_permutationtable;
+  int r, last = 0;
+
+  srand(* ((unsigned int *)(&seed)));
+
+  if(p)
+    {
+      p[0] = rand() & mask;
+      for (i=1;i<N;i++)
+	{
+	  r = rand() & mask;
+	  j = 0;
+	  while(j<i)
+	    {
+	      if (p[j] == r)
+		{
+		  r = (r + 1) & mask;
+		  j = 0;
+		}
+	      else j++;
+	    }
+	  p[i] = r;
+	  
+	}
+    }
+}
+
+static void permut_bang(t_permut *x)
+{
+    unsigned int r = rand();
+    permut_random(x, *((float *)(&r)));
+}
+
+static void permut_resize_table(t_permut *x, int size)
+{
+  if (x->x_ctl.c_blocksize != size)
+    {
+      if (x->x_ctl.c_permutationtable)
+	free(x->x_ctl.c_permutationtable);
+      x->x_ctl.c_permutationtable = (int *)malloc(sizeof(int)*size);
+      x->x_ctl.c_blocksize = size;
+
+      /* make sure it's initialized */
+      permut_bang(x);
+
+
+    } 
+}
+
+
+
+
+static t_int *permut_perform(t_int *w)
+{
+
+
+  t_float *in    = (float *)(w[3]);
+  t_float *out    = (float *)(w[4]);
+  t_permutctl *ctl  = (t_permutctl *)(w[1]);
+  t_int i;
+  t_int n = (t_int)(w[2]);
+  t_float x,y;
+  t_int *p =  ctl->c_permutationtable;
+
+
+  if (in != out)
+    for (i=0; i<n; i++) out[i] = in[i];
+
+
+  permut_perform_permutation(out, n, p);
+
+  return (w+5);
+}
+
+static void permut_dsp(t_permut *x, t_signal **sp)
+{
+  permut_resize_table(x, sp[0]->s_n);
+  dsp_add(permut_perform, 4, &x->x_ctl, sp[0]->s_n, sp[0]->s_vec, sp[1]->s_vec);
+
+}
+                                  
+static void permut_free(t_permut *x)
+{
+  if (x->x_ctl.c_permutationtable) free(x->x_ctl.c_permutationtable);
+
+}
+
+t_class *permut_class;
+
+static void *permut_new(void)
+{
+    t_permut *x = (t_permut *)pd_new(permut_class);
+    outlet_new(&x->x_obj, gensym("signal")); 
+
+    x->x_ctl.c_permutationtable = 0;
+    x->x_ctl.c_blocksize = 0;
+    permut_resize_table(x, 64);
+    permut_random(x, 0);
+    
+    return (void *)x;
+}
+
+void permut_tilde_setup(void)
+{
+  //post("permut~ v0.1");
+    permut_class = class_new(gensym("permut~"), (t_newmethod)permut_new,
+    	(t_method)permut_free, sizeof(t_permut), 0, A_DEFFLOAT, 0);
+    CLASS_MAINSIGNALIN(permut_class, t_permut, x_f); 
+    class_addmethod(permut_class, (t_method)permut_random, gensym("random"), A_FLOAT, 0);
+    class_addmethod(permut_class, (t_method)permut_bang, gensym("bang"), 0);
+    class_addmethod(permut_class, (t_method)permut_dsp, gensym("dsp"), 0); 
+
+}
+
-- 
cgit v1.2.1