/* Copyright (c) 2003 krzYszcz, and others.
 * For information on usage and redistribution, and for a DISCLAIMER OF ALL
 * WARRANTIES, see the file, "LICENSE.txt," in this distribution.  */

/* Filter: Paul Kellet's from music-dsp.  This is the most popular one,
   used in jMax, Csound, etc.  LATER compare to McCartney's (sc).
   Rng: noise~ code from d_osc.c. */

#include "m_pd.h"
#include "common/loud.h"
#include "sickle/sic.h"

/* more like 0.085 in c74's */
#define PINK_GAIN  .105

typedef struct _pink
{
    t_sic  x_sic;
    int    x_state;
    float  x_b0;
    float  x_b1;
    float  x_b2;
    float  x_b3;
    float  x_b4;
    float  x_b5;
    float  x_b6;
} t_pink;

static t_class *pink_class;

static void pink_float(t_pink *x, t_float f)
{
    loud_nomethod((t_pd *)x, &s_float);  /* LATER rethink */
}

static t_int *pink_perform(t_int *w)
{
    t_pink *x = (t_pink *)(w[1]);
    int nblock = (int)(w[2]);
    t_float *out = (t_float *)(w[3]);
    int state = x->x_state;
    float b0 = x->x_b0;
    float b1 = x->x_b1;
    float b2 = x->x_b2;
    float b3 = x->x_b3;
    float b4 = x->x_b4;
    float b5 = x->x_b5;
    float b6 = x->x_b6;
    while (nblock--)
    {
    	float white = ((float)((state & 0x7fffffff) - 0x40000000)) *
    	    (float)(1.0 / 0x40000000);
    	state = state * 435898247 + 382842987;
	b0 = 0.99886 * b0 + white * 0.0555179;
	b1 = 0.99332 * b1 + white * 0.0750759;
	b2 = 0.96900 * b2 + white * 0.1538520;
	b3 = 0.86650 * b3 + white * 0.3104856;
	b4 = 0.55000 * b4 + white * 0.5329522;
	b5 = -0.7616 * b5 - white * 0.0168980;
	*out++ = (b0 + b1 + b2 + b3 + b4 + b5 + b6 + white * 0.5362)
	    * PINK_GAIN;
	b6 = white * 0.115926;
    }
    x->x_state = state;
    x->x_b0 = b0;
    x->x_b1 = b1;
    x->x_b2 = b2;
    x->x_b3 = b3;
    x->x_b4 = b4;
    x->x_b5 = b5;
    x->x_b6 = b6;
    return (w + 4);
}

static void pink_dsp(t_pink *x, t_signal **sp)
{
    dsp_add(pink_perform, 3, x, sp[1]->s_n, sp[1]->s_vec);
}

static void *pink_new(void)
{
    t_pink *x = (t_pink *)pd_new(pink_class);
    /* borrowed from d_osc.c, LATER rethink */
    static int init = 307;
    x->x_state = (init *= 1319);
    /* all coefs set to zero */
    outlet_new((t_object *)x, &s_signal);
    return (x);
}

void pink_tilde_setup(void)
{
    pink_class = class_new(gensym("pink~"),
			   (t_newmethod)pink_new, 0,
			   sizeof(t_pink), 0, 0);
    /* dummy float method: we need signal in (for begin~), but neither
       float-to-signal conversion, nor a float method (the only control
       input is 'enable').  LATER rethink. */
    sic_setup(pink_class, pink_dsp, pink_float);
}