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
|
///////////////////////////////////////////////////////////////////////////////////
/* Lorenz's Attractor PD External */
/* Copyright Ben Bogart 2002 */
/* This program is distributed under the terms of the GNU General Public License */
///////////////////////////////////////////////////////////////////////////////////
///////////////////////////////////////////////////////////////////////////////////
/* This file is part of Chaos PD Externals. */
/* */
/* Chaos PD Externals are free software; you can redistribute them and/or modify */
/* them 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. */
/* */
/* Chaos PD Externals are distributed in the hope that they 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 the Chaos PD Externals; 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>
t_class *lorenz_class;
typedef struct lorenz_struct
{
t_object lorenz_obj;
double h, a, b, c, lx0, ly0, lz0;
t_outlet *y_outlet;
t_outlet *z_outlet;
} lorenz_struct;
static void calculate(lorenz_struct *x)
{
double lx0, ly0, lz0, lx1, ly1, lz1;
double h, a, b, c;
h = x->h;
a = x->a;
b = x->b;
c = x->c;
lx0 = x->lx0;
ly0 = x->ly0;
lz0 = x->lz0;
lx1 = lx0 + h * a * (ly0 - lx0);
ly1 = ly0 + h * (lx0 * (b - lz0) - ly0);
lz1 = lz0 + h * (lx0 * ly0 - c * lz0);
x->lx0 = lx1;
x->ly0 = ly1;
x->lz0 = lz1;
outlet_float(x->lorenz_obj.ob_outlet, (t_float)lx1);
outlet_float(x->y_outlet, (t_float)ly1);
outlet_float(x->z_outlet, (t_float)lz1);
}
static void reset(lorenz_struct *x, t_floatarg lx0, t_floatarg ly0, t_floatarg lz0)
{
x->lx0 = lx0;
x->ly0 = ly0;
x->lz0 = lz0;
}
static void param(lorenz_struct *x, t_floatarg h, t_floatarg a, t_floatarg b, t_floatarg c)
{
x->h = (double)h;
x->a = (double)a;
x->b = (double)b;
x->c = (double)c;
}
void *lorenz_new(void)
{
lorenz_struct *x = (lorenz_struct *)pd_new(lorenz_class);
x->h = 0.01;
x->a = 10.0;
x->b = 28.0;
x->c = 8.0/3.0;
x->lx0 = 0.1;
x->ly0 = 0;
x->lz0 = 0;
outlet_new(&x->lorenz_obj, &s_float); /* Default float outlet */
x->y_outlet = outlet_new(&x->lorenz_obj, &s_float); /* Two New Outlets */
x->z_outlet = outlet_new(&x->lorenz_obj, &s_float);
return (void *)x;
}
void lorenz_setup(void)
{
post("lorenz");
lorenz_class = class_new(gensym("lorenz"), /* symname is the symbolic name */
(t_newmethod)lorenz_new, /* Constructor Function */
0, /* Destructor Function */
sizeof(lorenz_struct), /* Size of the structure */
CLASS_DEFAULT, /* Graphical Representation */
0); /* 0 Terminates Argument List */
class_addbang(lorenz_class, (t_method)calculate);
class_addmethod(lorenz_class,
(t_method)reset,
gensym("reset"),
A_DEFFLOAT,
A_DEFFLOAT,
A_DEFFLOAT,
0);
class_addmethod(lorenz_class,
(t_method)param,
gensym("param"),
A_DEFFLOAT,
A_DEFFLOAT,
A_DEFFLOAT,
A_DEFFLOAT,
0);
}
|
