-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathhenon~.c
executable file
·165 lines (125 loc) · 2.99 KB
/
henon~.c
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
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
/* code by Chikashi Miyama 2009 */
#include "m_pd.h"
#include <math.h>
#include <stdlib.h>
#ifdef NT
#pragma warning( disable : 4244 )
#pragma warning( disable : 4305 )
#endif
/* ----------------------------- henon----------------------------- */
static t_class *henon_class;
typedef struct _henon
{
t_object x_obj;
t_float x_nextx;
t_float x_nexty;
t_float x_previousx;
t_float x_previousy;
t_float x_gapx;
t_float x_gapy;
t_float x_hold;
t_float x_head;
t_float x_sr;
t_float x_f;
} t_henon;
static void *henon_new(void)
{
t_henon*x = (t_henon*)pd_new(henon_class);
inlet_new(&x->x_obj, &x->x_obj.ob_pd, &s_signal, &s_signal);
inlet_new(&x->x_obj, &x->x_obj.ob_pd, &s_signal, &s_signal);
inlet_new(&x->x_obj, &x->x_obj.ob_pd, &s_signal, &s_signal);
outlet_new(&x->x_obj, &s_signal);
outlet_new(&x->x_obj, &s_signal);
x->x_previousx = 0.6;
x->x_previousy = 0.6;
x->x_nextx = 0.0;
x->x_nexty = 0.0;
x->x_head = 1.1;
return (x);
}
t_int *henon_perform(t_int *w)
{
t_henon *x = (t_henon*)(w[1]);
t_sample *freq = (t_sample *)(w[2]);
t_sample *coef1 = (t_sample *)(w[3]);
t_sample *coef2 = (t_sample *)(w[4]);
t_sample *curve = (t_sample *)(w[5]);
t_sample *out1 = (t_sample *)(w[6]);
t_sample *out2 = (t_sample *)(w[7]);
t_float nextx = x->x_nextx;
t_float nexty = x->x_nexty;
t_float previousx = x->x_previousx;
t_float previousy = x->x_previousy;
t_float gapx = x->x_gapx;
t_float gapy = x->x_gapy;
t_float head = x->x_head;
t_float hold = x->x_hold;
t_float sr = x->x_sr;
t_float co1,co2;
int n = (int)(w[8]);
float step,cfrq, distorted;
while (n--){
if(head > 1.0)
{
head = 0.0;
previousx = nextx;
previousy = nexty;
co1 = *coef1;
co2 = *coef2;
nextx = previousy + 1 - co1 * previousx * previousx;
nexty = co2 * previousx;
gapx = nextx - previousx;
gapy = nexty - previousy;
hold = *curve;
}
cfrq = *freq++;
if(cfrq == 0)
{
step = 0;
}
else
{
step = 1.0/(sr / cfrq);
}
head += step;
//power
if(hold == 0.0)
{
distorted = 0.0;
}
else if(hold == 1.0)
{
distorted = head;
}
else
{
distorted = pow(head, hold);
}
*curve++;
*coef1++;
*coef2++;
*out1++ = distorted * gapx + previousx;
*out2++ = distorted * gapy + previousy;
}
x->x_nextx = nextx;
x->x_nexty = nexty;
x->x_previousx = previousx;
x->x_previousy = previousy;
x->x_head = head;
x->x_gapx = gapx;
x->x_gapy = gapy;
x->x_hold = hold;
return (w+9);
}
static void henon_dsp(t_henon*x, t_signal **sp)
{
x->x_sr = sp[0]->s_sr;
dsp_add(henon_perform, 8, x, sp[0]->s_vec, sp[1]->s_vec, sp[2]->s_vec, sp[3]->s_vec, sp[4]->s_vec, sp[5]->s_vec, sp[0]->s_n);
}
void henon_tilde_setup(void)
{
henon_class = class_new(gensym("henon~"), (t_newmethod)henon_new, 0, sizeof(t_henon), 0, 0);
class_addmethod(henon_class, (t_method)henon_dsp, gensym("dsp"), 0);
CLASS_MAINSIGNALIN(henon_class, t_henon, x_f);
class_sethelpsymbol(henon_class, gensym("sigbinops"));
}