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
94
95
96
97
98
99
100
101
102
103
104
105
106
|
function buildsource(in_N)
Globals2D;
NODETOL = 1e-8;
N = in_N;
Np = (N+1)*(N+2)/2; Nfp = N+1; Nfaces = 3;
[x,y] = Nodes2D(N);
[r,s] = xytors(x,y);
% find all the nodes that lie on each edge
fmask1 = find( abs(s+1) < NODETOL)';
fmask2 = find( abs(r+s) < NODETOL)';
fmask3 = find( abs(r+1) < NODETOL)';
Fmask = [fmask1;fmask2;fmask3]';
% vandermonde matrix
V = Vandermonde2D(N, r, s);
% derivative matrices
[Dr, Ds] = Dmatrices2D(N, r, s, V);
% lift matrix
LIFT = Lift2D();
% correct to 0-index
Fmask = Fmask-1;
fid = fopen(sprintf('dataN%02d.h', N), 'w');
fprintf(fid, sprintf('#ifndef DATAN%02d \n\n', N));
fprintf(fid, sprintf('#define DATAN%02d 1\n\n', N));
fprintf(fid, 'double p_r[%d] = { ', Np);
fprintf(fid, '%17.15g ', r(1));
for n=2:Np
fprintf(fid, ', %17.15g ', r(n));
end
fprintf(fid, '};\n')
fprintf(fid, 'double p_s[%d] = {', Np);
fprintf(fid, '%17.15g ', s(1));
for n=2:Np
fprintf(fid, ', %17.15g ', s(n));
end
fprintf(fid, '};\n')
fprintf(fid, 'double p_Dr[%d][%d] = {', Np, Np);
for m=1:Np
fprintf(fid, '{%17.15g ', Dr(m,1));
for n=2:Np
fprintf(fid, ', %17.15g ', Dr(m,n));
end
if(m<Np)
fprintf(fid, '},\n')
else
fprintf(fid, '}};\n')
end
end
fprintf(fid, 'double p_Ds[%d][%d] = {', Np, Np);
for m=1:Np
fprintf(fid, '{%17.15g ', Ds(m,1));
for n=2:Np
fprintf(fid, ', %17.15g ', Ds(m,n));
end
if(m<Np)
fprintf(fid, '},\n')
else
fprintf(fid, '}};\n')
end
end
fprintf(fid, 'double p_LIFT[%d][%d] = {', Np, Nfp*Nfaces);
for m=1:Np
fprintf(fid, '{%17.15g ', LIFT(m,1));
for n=2:Nfp*Nfaces
fprintf(fid, ', %17.15g ', LIFT(m,n));
end
if(m<Np)
fprintf(fid, '},\n')
else
fprintf(fid, '}};\n')
end
end
fprintf(fid, 'int p_Fmask[%d][%d] = {', Nfaces, Nfp);
for m=1:Nfaces
fprintf(fid, '{%d ', Fmask(1,m));
for n=2:Nfp
fprintf(fid, ', %d ', Fmask(n,m));
end
if(m<Nfaces)
fprintf(fid, '},\n')
else
fprintf(fid, '}};\n')
end
end
fprintf(fid, sprintf('\n#endif\n'));
fclose(fid)
|
