1 files changed, 61 insertions, 0 deletions
diff --git a/gui/libs/ofxAChaosLib/src/AChaosHenonHeilles.h b/gui/libs/ofxAChaosLib/src/AChaosHenonHeilles.h
new file mode 100644
index 0000000..c1aeb98
--- /dev/null
+++ b/gui/libs/ofxAChaosLib/src/AChaosHenonHeilles.h
@@ -0,0 +1,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;
+ 	}
+};
+\ No newline at end of file