/*
	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 AChaosNavierStokes : public AChaosBase {
public:

	REAL a, b, c, d, e, r, dt;

	vector<string> param_labels={"nv","nw","nx","ny","nz","r","dt"};
	
	AChaosNavierStokes(){}
	~AChaosNavierStokes(){}	

	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];
	}
	
#define nava(a, b, c, d, e)		((-2*(a)) + (4*(b)*(c)) + (4*(d)*(e)))
#define navb(a, b, c)			((-9*(b)) + (3*(a)*(c)))
#define navc(a, b, c, r)		((-5*(c)) + (-7*(a)*(b)) + (r))
#define navd(a, d, e)			((-5*(d)) - ((a)*(e)))
#define nave(a, d, e)			((-3*(a)*(d)) - (e))

	void calc(){ 	

		REAL dt2=dt/2., afa, bfb, cfc, dfd, efe;
		REAL fa, fb, fc, fd, fe;
		REAL ffa, ffb, ffc, ffd, ffe;

		fa=nava(a, b, c, d, e);
		fb=navb(a, b, c);
		fc=navc(a, b, c, r);
		fd=navd(a, d, e);
		fe=nave(a, d, e);
		
		afa = a +(dt*fa);
		bfb = b +(dt*fb);
		cfc = c +(dt*fc);
		dfd = d +(dt*fd);
		efe = e +(dt*fe);
		
		ffa = nava(afa,bfb,cfc,dfd,efe);
		ffb = navb(afa,bfb,cfc);
		ffc = navc(afa,bfb,cfc,r);
		ffd = navd(afa,dfd,efe);
		ffe = nave(afa,dfd,efe);
		
		a += (dt2*(fa+ffa));
		b += (dt2*(fb+ffb));
		c += (dt2*(fc+ffc));
		d += (dt2*(fd+ffd));
		e += (dt2*(fe+ffe));

		ov[0] = a;
		ov[1] = b;
		ov[2] = c;
		ov[3] = d;
		ov[4] = e;

 	}
};