1 files changed, 75 insertions, 0 deletions
diff --git a/gui/libs/ofxAChaosLib/src/AChaosHenonF.h b/gui/libs/ofxAChaosLib/src/AChaosHenonF.h
new file mode 100644
index 0000000..384a6ab
--- /dev/null
+++ b/gui/libs/ofxAChaosLib/src/AChaosHenonF.h
@@ -0,0 +1,75 @@
+/*
+	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 AChaosHenonF : public AChaosBase {
+public:
+
+	REAL a, b, nx, ny, dt, t;
+	
+	AChaosHenonF(){}
+	~AChaosHenonF(){}	
+
+		vector<string> param_labels={"a","b","","","dt","t"};
+
+
+	virtual void setup(REAL * params = NULL){
+		AChaosBase::init(params, 5,2);
+		if(params==NULL){
+			//init	
+			a = 1.4f;
+			b = 0.3f;
+			nx = 1.0f;
+			ny = 1.0f;
+			dt = 0.01f;
+			t=0.0f;
+			REAL p[5] = {nx,ny,a,b,dt};
+			set(p);
+		} else {
+			set(params);
+		}
+	}
+
+
+	void reset(){
+		nx = iv[0];
+		ny = iv[1];
+		a = iv[2];
+		b = iv[3];
+		dt = iv[4];
+	}
+	
+ 	void calc(){ 		
+		// Henon Phase Diagrams
+ 		REAL cosa=COS(a),sina=SIN(a), lx1, ly1, nx2=nx*nx;
+		
+		//henon calc
+		lx1 = (ny + 1) - (a * nx2);
+		ly1 = b * nx;
+
+		//advance 'a' by dt _AFTER_ calculating henon
+		t += dt;
+		a *= t;
+
+		//henon phase
+		nx = nx*cosa - (ny -nx2)*sina;
+		ny = nx*sina - (ny -nx2)*cosa;
+
+
+/*
+	// henonf attractor
+
+	lx1 = nx*cos (a) - (ny - (nx*nx))*sin(a);
+	ly1 = nx*sin(a) - (ny - (nx*nx))*cos(a);
+*/
+
+		ov[0] = nx;
+		ov[1] = ny;
+ 	}
+};
+\ No newline at end of file