/* 
 * 
 * Error function for pd
 * Copyright (C) 2006 Tim Blechmann
 * 
 * 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; see the file COPYING.  If not, write to
 * the Free Software Foundation, Inc., 59 Temple Place - Suite 330,
 * Boston, MA 02111-1307, USA.
 */

#include "m_pd.h"
#include "math.h"

static t_class *erf_class;

typedef struct _erf
{
    t_object x_obj;
    float x_f;
} t_erf;

static void *erf_new(void)
{
    t_erf *x = (t_erf *)pd_new(erf_class);
    outlet_new(&x->x_obj, gensym("signal"));
    x->x_f = 0;
    return (x);
}

static t_int *erf_perform(t_int *w)
{
    t_float *in = (t_float *)(w[1]);
    t_float *out = (t_float *)(w[2]);
    int n = (int)(w[3]);

    while (n--)
        *out++ = erff(*in++);

    return w+4;
}

static void erf_dsp(t_erf *x, t_signal **sp)
{
    dsp_add(erf_perform, 3, sp[0]->s_vec, sp[1]->s_vec, sp[0]->s_n);
}

void erf_tilde_setup(void)
{
    erf_class = class_new(gensym("erf~"), (t_newmethod)erf_new, 0,
        sizeof(t_erf), 0, A_DEFFLOAT, 0);
    CLASS_MAINSIGNALIN(erf_class, t_erf, x_f);
    class_addmethod(erf_class, (t_method)erf_dsp, gensym("dsp"), 0);
}