summaryrefslogtreecommitdiff
path: root/gui/libs/ofxAChaosLib/src/AChaosRossler.h
diff options
context:
space:
mode:
authorTim Redfern <tim@getdrop.com>2018-07-08 12:22:06 +0100
committerTim Redfern <tim@getdrop.com>2018-07-08 12:22:06 +0100
commitfb68eed64f548d090eb550047fd0d898e4e033fc (patch)
tree890fa917ca6bd88505b45e5c261d05e0c5c9b6ec /gui/libs/ofxAChaosLib/src/AChaosRossler.h
parentc10665e02e49da2d8f7e1a1eada151f17f507a21 (diff)
add altered AChaoslib
Diffstat (limited to 'gui/libs/ofxAChaosLib/src/AChaosRossler.h')
-rw-r--r--gui/libs/ofxAChaosLib/src/AChaosRossler.h74
1 files changed, 74 insertions, 0 deletions
diff --git a/gui/libs/ofxAChaosLib/src/AChaosRossler.h b/gui/libs/ofxAChaosLib/src/AChaosRossler.h
new file mode 100644
index 0000000..3bbbc42
--- /dev/null
+++ b/gui/libs/ofxAChaosLib/src/AChaosRossler.h
@@ -0,0 +1,74 @@
+/*
+ 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"
+
+#define rosx(y,z) -((y)+(z))
+#define rosy(x,y,a) (x)+((a)*(y))
+#define rosz(x,z,b,c) (b)+((x)*(z))-((c)*(z))
+
+class AChaosRossler : public AChaosBase {
+public:
+
+ REAL a, b, c, nx, ny, nz,dt;
+
+ vector<string> param_labels={"a","b","c","","","","dt"};
+
+ AChaosRossler(){}
+ ~AChaosRossler(){}
+
+ virtual void setup(REAL * params = NULL){
+
+ AChaosBase::init(params, 7, 3);
+ if(params==NULL){
+ //classic rossler
+ a = 0.2;
+ b = 0.2;
+ c = 5.7;
+ nx = 1;
+ ny = 1;
+ nz = 1;
+ dt = 0.01;
+ REAL p[7] = {a,b,c,nx,ny,nz,dt};
+ set(p);
+ } else { set(params); }
+
+ }
+
+ void reset(){
+ a = iv[0];
+ b = iv[1];
+ c = iv[2];
+ nx = iv[3];
+ ny = iv[4];
+ nz = iv[5];
+ dt = iv[6];
+ }
+
+
+ void calc(){
+ REAL fx, fy, fz, ffx, ffy, ffz, dt2 = dt/2.;
+ REAL dtfx, dtfy, dtfz;
+
+ fx = rosx(ny,nz);
+ fy = rosy(nx,ny,a);
+ fz = rosz(nx,nz,b,c);
+
+ ffx = rosx((ny+dt*fy), (nz+(dt*fz)));
+ ffy = rosy((nx+dt*nx), (ny+dt*ny), a);
+ ffz = rosz((nx+dt*fx), (nz+dt*nz), b, c);
+
+ nx += (dt2*(fx+ffx));
+ ny += (dt2*(fy+ffy));
+ nz += (dt2*(fz+ffz));
+
+ ov[0] = nx;
+ ov[1] = ny;
+ ov[2] = nz;
+ }
+}; \ No newline at end of file