1 files changed, 69 insertions, 0 deletions
diff --git a/gui/libs/ofxAChaosLib/src/AChaosLorenzEuler.h b/gui/libs/ofxAChaosLib/src/AChaosLorenzEuler.h
new file mode 100644
index 0000000..d23889d
--- /dev/null
+++ b/gui/libs/ofxAChaosLib/src/AChaosLorenzEuler.h
@@ -0,0 +1,69 @@
+/*
+	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 lorx(x, y, a) 		(a) * ((y) - (x));		
+#define lory(x, y, z, a) 	((x) * ((r) - (z))) - (y);
+#define lorz(x, y, z, c) 	((x) * (y)) - ((z) * (c));
+
+class AChaosLorenzEuler : public AChaosBase {
+public:
+
+	REAL	a, r, c, nx, ny, nz,dt; //eq variables
+	
+	AChaosLorenzEuler(){}
+	~AChaosLorenzEuler(){}	
+
+	vector<string> param_labels={"a","r","c","nx","ny","nz","dt"};
+	
+	virtual void setup(REAL * params = NULL){
+	
+		AChaosBase::init(params, 7, 3);	
+		if(params==NULL){		
+			//classic lorenz
+			 a=10.0;
+			 r=28.0; 
+			 c= 2.67; 
+			 nx=0.1; 
+			 ny=0.1; 
+			 nz=0.1;
+			 dt=0.01;
+			REAL p[7] = {nx,ny,nz,a,r,c,dt};
+			set(p);
+		} else { set(params); }
+
+	}
+
+
+	void reset(){
+		nx = iv[0];
+		ny = iv[1];
+		nz = iv[2];
+		a = iv[3];
+		r = iv[4];
+		c = iv[5];
+		dt = iv[6];
+	}
+	
+ 	void calc(){ 		
+		REAL dx, dy, dz;
+		
+		dx =  ((ny*a)-(nx*a))*dt;
+		dy = ((nx*r)-ny-(nx*nz))*dt;
+		dz = ((nx*ny)-(nz*c))*dt;
+		
+		nx += dx;
+		ny += dy;
+		nz += dz;
+
+		ov[0] = nx;
+		ov[1] = ny;
+		ov[2] = nz;
+ 	}
+};
+\ No newline at end of file