1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
|
/* 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 <string.h>
#include "m_pd.h"
#include "shared.h"
#include "sickle/sic.h"
#include "sickle/arsic.h"
/* CHECKME (the refman): if the buffer~ has more channels, channels are mixed */
typedef t_arsic t_play;
static t_class *play_class;
static void play_set(t_play *x, t_symbol *s)
{
arsic_setarray((t_arsic *)x, s, 1);
}
/* LATER optimize */
static t_int *play_perform(t_int *w)
{
t_arsic *sic = (t_arsic *)(w[1]);
int nblock = (int)(w[2]);
int nch = sic->s_nchannels;
t_int *outp = w + 4;
if (sic->s_playable)
{
t_play *x = (t_play *)sic;
t_float *xin = (t_float *)(w[3]);
int vecsize = sic->s_vecsize;
t_word **vectable = sic->s_vectors;
float ksr = sic->s_ksr;
int nointerp = 0;
int maxindex = (nointerp ? vecsize - 1 : vecsize - 3);
int iblock;
for (iblock = 0; iblock < nblock; iblock++)
{
float phase = *xin++ * ksr;
int ndx;
int ch = nch;
float frac, a, b, c, d, cminusb;
if (phase < 0 || phase > maxindex)
phase = 0; /* CHECKED: a value 0, not ndx 0 */
ndx = (int)phase;
/* CHECKME: what kind of interpolation? (CHECKED: multi-point) */
if (ndx < 1)
ndx = 1, frac = 0;
else if (ndx > maxindex)
ndx = maxindex, frac = 1;
else frac = phase - ndx;
while (ch--)
{
t_word *vp = vectable[ch];
t_float *out = (t_float *)(outp[ch]);
if (vp)
{
vp += ndx;
a = vp[-1].w_float;
b = vp[0].w_float;
c = vp[1].w_float;
d = vp[2].w_float;
cminusb = c-b;
out[iblock] = b + frac * (
cminusb - 0.1666667f * (1. - frac) * (
(d - a - 3.0f * cminusb) * frac
+ (d + 2.0f * a - 3.0f * b)
)
);
}
else out[iblock] = 0;
}
}
}
else
{
int ch = nch;
while (ch--)
{
t_float *out = (t_float *)outp[ch];
int n = nblock;
while (n--) *out++ = 0;
}
}
return (w + sic->s_nperfargs + 1);
}
static void play_dsp(t_play *x, t_signal **sp)
{
arsic_dsp((t_arsic *)x, sp, play_perform, 1);
}
static void play_free(t_play *x)
{
arsic_free((t_arsic *)x);
}
static void *play_new(t_symbol *s, t_floatarg f)
{
/* one auxiliary signal: position input */
t_play *x = (t_play *)arsic_new(play_class, s, (int)f, 0, 1);
if (x)
{
int nch = arsic_getnchannels((t_arsic *)x);
while (nch--)
outlet_new((t_object *)x, &s_signal);
}
return (x);
}
void play_tilde_setup(void)
{
play_class = class_new(gensym("play~"),
(t_newmethod)play_new,
(t_method)play_free,
sizeof(t_play), 0,
A_DEFSYM, A_DEFFLOAT, 0);
arsic_setup(play_class, play_dsp, SIC_FLOATTOSIGNAL);
class_addmethod(play_class, (t_method)play_set,
gensym("set"), A_SYMBOL, 0);
logpost(NULL, 4, "this is cyclone/play~ %s, %dth %s build",
CYCLONE_VERSION, CYCLONE_BUILD, CYCLONE_RELEASE);
}
|
