/* 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);
}