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
|
/*
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 AChaosIkeda : public AChaosBase {
public:
REAL a, b, k, p, nx, ny;
AChaosIkeda(){}
~AChaosIkeda(){}
vector<string> param_labels={"a","b","k","p","",""};
virtual vector<string> &get_param_labels(){return this->param_labels;};
virtual void setup(REAL * params = NULL){
AChaosBase::init(params, 6, 2);
if(params==NULL){
//init
a=0.85;
b=0.9;
k=0.4;
p=7.7;
nx=0;
ny=0;
REAL pa[6] = {nx,ny,a,b,k,p};
set(pa);
} else { set(params); }
}
void reset(){
nx = iv[0];
ny = iv[1];
a = iv[2];
b = iv[3];
k = iv[4];
p = iv[5];
}
void calc(){
REAL x0 = nx, y0 = ny;
REAL temp = (p/(1+x0*x0+y0*y0));
REAL cosx = (COS(b-temp));
REAL sinx = (SIN(b-temp));
/*
xn+1 = a + k*( x0*cos(b-(p/(1+x0*x0+y0*y0))) - y0*sin(b-(p/(1*x0*x0+y0*y0))) )
yn+1 = k*( x0*sin(b-(p/(1+x0*x0+y0*y0))) - y0*cos(b-(p/(1*x0*x0+y0*y0))) )
*/
nx = a + k*(x0*cosx)-y0*sinx;
ny = k*sinx -y0*cosx;
ov[0] = nx;
ov[1] = ny;
}
};
|
