add altered AChaoslib

1 files changed, 73 insertions, 0 deletions

diff --git a/gui/libs/ofxAChaosLib/src/AChaosNavierStokesEuler.h b/gui/libs/ofxAChaosLib/src/AChaosNavierStokesEuler.h
new file mode 100644
index 0000000..9018ce4
--- /dev/null
+++ b/gui/libs/ofxAChaosLib/src/AChaosNavierStokesEuler.h
@@ -0,0 +1,73 @@
+/*
+	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 AChaosNavierStokesEuler : public AChaosBase {
+public:
+
+	REAL a, b, c, d, e, r, dt;
+	
+	AChaosNavierStokesEuler(){}
+	~AChaosNavierStokesEuler(){}	
+
+	vector<string> param_labels={"nv","nw","nx","ny","nz","r","dt"};
+
+	virtual void setup(REAL * params = NULL){
+	
+		AChaosBase::init(params, 7, 5);	
+		if(params==NULL){		
+			//classic navierstokes
+			a = 1.0;
+			b = 1.0;
+			c = 1.0;
+			d = 1.0;
+			e = 1.0;
+			r = 28.0;
+			dt = 0.01;
+			REAL p[7] = {a,b,c,d,e,r,dt};
+			set(p);
+		} else { set(params); }
+
+	}
+
+	void reset(){
+		a = iv[0];
+		b = iv[1];
+		c = iv[2];
+		d = iv[3];
+		e = iv[4];
+		r = iv[5];
+		dt = iv[6];
+	}
+	
+
+	void calc(){ 	
+
+		REAL da, db, dc, dd, de;
+	
+		da=((a * -2) + (4 * b * c) + (4 * d * e))*dt;
+		db=((b * -9) + (3 * a * c))*dt;
+		dc=((c * -5) + (-7 * a * b) + r)*dt;
+		dd=((d * -5) - (a * e))*dt;
+		de=((a * d * -3) - e)*dt;
+
+		a += da;
+		b += db;
+		c += dc;
+		d += dd;
+		e += de;
+
+		ov[0] = a;
+		ov[1] = b;
+		ov[2] = c;
+		ov[3] = d;
+		ov[4] = e;
+
+ 	}
+};
\ No newline at end of file