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
|
/*
32/64 bits A-Chaos Lib in openFrameworks
(c) s373.net/x 2004, 2012, 2015
http://s373.net/code/A-Chaos-Lib/A-Chaos.html
programmed by Andre Sier, revised 2015
License: MIT
*/
#pragma once
#include "AChaosBase.h"
class AChaosHenonHeilles : public AChaosBase {
public:
REAL nx, ny, nxdot, nydot, e, dt;
AChaosHenonHeilles(){}
~AChaosHenonHeilles(){}
vector<string> param_labels={"","","","","e","dt"};
virtual void setup(REAL * params = NULL){
AChaosBase::init(params, 5, 4);
if(params==NULL){
//init
nx = 1.0f;
ny = 1.0f;
nydot = 1.4f;
e = 0.125f;
dt = 0.02;
REAL p[5] = {nx,ny, nydot,e,dt};
set(p);
} else { set(params); }
}
void reset(){
nx = iv[0];
ny = iv[1];
nydot = iv[2];
e = iv[3];
dt = iv[4];
}
void calc(){
// henon_heilles_calc
REAL xdd,ydd;
xdd = -nx-(2*nx*ny);
ydd = ny*ny -ny -(nx*nx);
nxdot += xdd*dt;
nydot += ydd*dt;
nx += nxdot*dt;
ny += nydot*dt;
ov[0] = nx;
ov[1] = ny;
ov[2] = nxdot;
ov[3] = nydot;
}
};
|
