From a7e6303e3e7e87bf90b21821d2e419ed4799783c Mon Sep 17 00:00:00 2001
From: musil <tmusil@users.sourceforge.net>
Date: Fri, 28 Jun 2013 17:25:22 +0000
Subject: initial check in of double precision library of iem

svn path=/trunk/externals/iem/iem_dp/; revision=17167
---
 src/tabread~~.c | 117 ++++++++++++++++++++++++++++++++++++++++++++++++++++++++
 1 file changed, 117 insertions(+)
 create mode 100755 src/tabread~~.c

(limited to 'src/tabread~~.c')

diff --git a/src/tabread~~.c b/src/tabread~~.c
new file mode 100755
index 0000000..6f5a3ce
--- /dev/null
+++ b/src/tabread~~.c
@@ -0,0 +1,117 @@
+/* For information on usage and redistribution, and for a DISCLAIMER OF ALL
+* WARRANTIES, see the file, "LICENSE.txt," in this distribution.
+
+iem_dp written by IOhannes m zmoelnig, Thomas Musil, Copyright (c) IEM KUG Graz Austria 1999 - 2007 */
+/* double precision library */
+
+
+#include "m_pd.h"
+#include "iemlib.h"
+#include "iem_dp.h"
+
+/* -------------------------- tabread~~ ------------------------------ */
+/* based on miller's tabread~ which is part of pd */
+
+/******************** tabread~ ***********************/
+
+static t_class *tabread_tilde_tilde_class;
+
+typedef struct _tabread_tilde_tilde
+  {
+    t_object x_obj;
+    int x_npoints;
+    iemarray_t  *x_vec;
+    t_symbol *x_arrayname;
+    t_float x_f;
+  } t_tabread_tilde_tilde;
+
+static void *tabread_tilde_tilde_new(t_symbol *s)
+{
+  t_tabread_tilde_tilde *x = (t_tabread_tilde_tilde *)pd_new(tabread_tilde_tilde_class);
+  x->x_arrayname = s;
+  x->x_npoints=0;
+  x->x_vec = 0;
+  inlet_new(&x->x_obj, &x->x_obj.ob_pd, &s_signal, &s_signal);
+  outlet_new(&x->x_obj, gensym("signal"));
+  x->x_f = 0;
+  return (x);
+}
+
+static t_int *tabread_tilde_tilde_perform(t_int *w)
+{
+  t_tabread_tilde_tilde *x = (t_tabread_tilde_tilde *)(w[1]);
+  t_sample *in0 = (t_sample *)(w[2]);
+  t_sample *in1 = (t_sample *)(w[3]);
+  t_sample *out = (t_sample *)(w[4]);
+  int n = (int)(w[5]);    
+  int maxindexminusone;
+  iemarray_t *buf = x->x_vec;
+  int i;
+  
+  maxindexminusone = x->x_npoints - 1;
+  if(!buf)
+  {
+    while(n--)
+      *out++ = 0;
+    return(w+6);
+  }
+  
+  for (i = 0; i < n; i++)
+  {
+    t_sample in0_s=*in0++;
+    t_sample in1_s=*in1++;
+    double findex = iem_dp_calc_sum(in0_s, in1_s);
+    int index = (int)findex;
+    
+    if (index < 0)
+      index = 0;
+    else if (index > maxindexminusone)
+      index = maxindexminusone;
+    *out++ = iemarray_getfloat(buf, index);
+  }
+  return (w+6);
+}
+
+void tabread_tilde_tilde_set(t_tabread_tilde_tilde *x, t_symbol *s)
+{
+  t_garray *a;
+  
+  x->x_arrayname = s;
+  if (!(a = (t_garray *)pd_findbyclass(x->x_arrayname, garray_class)))
+  {
+    if (*s->s_name)
+      pd_error(x, "tabread~: %s: no such array", x->x_arrayname->s_name);
+    x->x_vec = 0;
+  }
+  else if (!iemarray_getarray(a, &x->x_npoints, &x->x_vec))
+  {
+    pd_error(x, "%s: bad template for tabread~", x->x_arrayname->s_name);
+    x->x_vec = 0;
+  }
+  else garray_usedindsp(a);
+}
+
+static void tabread_tilde_tilde_dsp(t_tabread_tilde_tilde *x, t_signal **sp)
+{
+  tabread_tilde_tilde_set(x, x->x_arrayname);
+  
+  dsp_add(tabread_tilde_tilde_perform, 5, x,
+          sp[0]->s_vec, sp[1]->s_vec, sp[2]->s_vec, sp[0]->s_n);
+  
+}
+
+static void tabread_tilde_tilde_free(t_tabread_tilde_tilde *x)
+{
+}
+
+void tabread_tilde_tilde_setup(void)
+{
+  tabread_tilde_tilde_class = class_new(gensym("tabread~~"),
+                                  (t_newmethod)tabread_tilde_tilde_new, (t_method)tabread_tilde_tilde_free,
+                                  sizeof(t_tabread_tilde_tilde), 0, A_DEFSYM, 0);
+  CLASS_MAINSIGNALIN(tabread_tilde_tilde_class, t_tabread_tilde_tilde, x_f);
+  class_addmethod(tabread_tilde_tilde_class, (t_method)tabread_tilde_tilde_dsp,
+                  gensym("dsp"), 0);
+  class_addmethod(tabread_tilde_tilde_class, (t_method)tabread_tilde_tilde_set,
+        gensym("set"), A_SYMBOL, 0);
+}
-- 
cgit v1.2.1