1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
|
/*
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;
}
};
|
