1 files changed, 117 insertions, 0 deletions
diff --git a/cyclone/sickle/train.c b/cyclone/sickle/train.c
new file mode 100644
index 0000000..b07ba9a
--- /dev/null
+++ b/cyclone/sickle/train.c
@@ -0,0 +1,117 @@
+/* Copyright (c) 2002-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.  */
+
+#include "m_pd.h"
+#include "shared.h"
+#include "sickle/sic.h"
+
+#define TRAIN_DEFPERIOD  1000
+#define TRAIN_DEFWIDTH      0.5
+#define TRAIN_DEFOFFSET     0
+
+typedef struct _train
+{
+    t_sic      x_sic;
+    int        x_on;
+    double     x_phase;
+    float      x_rcpksr;
+    t_outlet  *x_bangout;
+    t_clock   *x_clock;
+} t_train;
+
+static t_class *train_class;
+
+static void train_tick(t_train *x)
+{
+    outlet_bang(x->x_bangout);
+}
+
+static t_int *train_perform(t_int *w)
+{
+    t_train *x = (t_train *)(w[1]);
+    int nblock = (int)(w[2]);
+    t_float *in1 = (t_float *)(w[3]);
+    t_float *in2 = (t_float *)(w[4]);
+    t_float *in3 = (t_float *)(w[5]);
+    t_float *out = (t_float *)(w[6]);
+    float rcpksr = x->x_rcpksr;
+    double ph = x->x_phase;
+    double tfph = ph + SHARED_UNITBIT32;
+    t_shared_wrappy wrappy;
+    int32 normhipart;
+    int on = x->x_on;
+    int edge = 0;
+
+    wrappy.w_d = SHARED_UNITBIT32;
+    normhipart = wrappy.w_i[SHARED_HIOFFSET];
+
+    while (nblock--)
+    {
+	double onph, offph;
+	float period = *in1++;
+
+	wrappy.w_d = *in3++ + SHARED_UNITBIT32;
+    	wrappy.w_i[SHARED_HIOFFSET] = normhipart;
+	onph = wrappy.w_d - SHARED_UNITBIT32;
+
+	wrappy.w_d = onph + *in2++ + SHARED_UNITBIT32;
+    	wrappy.w_i[SHARED_HIOFFSET] = normhipart;
+	offph = wrappy.w_d - SHARED_UNITBIT32;
+
+	if (offph > onph ? ph < offph && ph >= onph : ph < offph || ph >= onph)
+	{
+	    if (!on) on = edge = 1;
+	    *out++ = 1.;
+	}
+	else
+	{
+	    on = 0;
+	    *out++ = 0.;
+	}
+	if (period > rcpksr)  /* LATER rethink */
+	    tfph += rcpksr / period;  /* LATER revisit (profiling?) */
+	wrappy.w_d = tfph;
+    	wrappy.w_i[SHARED_HIOFFSET] = normhipart;
+	ph = wrappy.w_d - SHARED_UNITBIT32;
+    }
+    x->x_phase = ph;
+    x->x_on = on;
+    if (edge) clock_delay(x->x_clock, 0);
+    return (w + 7);
+}
+
+static void train_dsp(t_train *x, t_signal **sp)
+{
+    x->x_rcpksr = 1000. / sp[0]->s_sr;
+    dsp_add(train_perform, 6, x, sp[0]->s_n,
+	    sp[0]->s_vec, sp[1]->s_vec, sp[2]->s_vec, sp[3]->s_vec);
+}
+
+static void train_free(t_train *x)
+{
+    if (x->x_clock) clock_free(x->x_clock);
+}
+
+static void *train_new(t_symbol *s, int ac, t_atom *av)
+{
+    t_train *x = (t_train *)pd_new(train_class);
+    x->x_on = 0;
+    x->x_phase = 0;
+    sic_inlet((t_sic *)x, 0, TRAIN_DEFPERIOD, 0, ac, av);
+    sic_inlet((t_sic *)x, 1, TRAIN_DEFWIDTH, 1, ac, av);
+    sic_inlet((t_sic *)x, 2, TRAIN_DEFOFFSET, 2, ac, av);
+    outlet_new((t_object *)x, &s_signal);
+    x->x_bangout = outlet_new((t_object *)x, &s_bang);
+    x->x_clock = clock_new(x, (t_method)train_tick);
+    return (x);
+}
+
+void train_tilde_setup(void)
+{
+    train_class = class_new(gensym("train~"),
+			    (t_newmethod)train_new,
+			    (t_method)train_free,
+			    sizeof(t_train), 0, A_GIMME, 0);
+    sic_setup(train_class, train_dsp, SIC_FLOATTOSIGNAL);
+}