/*
	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;
 	}
};