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
| #define BGHOSTS 2
#include "grid/multigrid.h"
int main()
{
init_grid (4);
output_cells (stdout);
vector v[];
scalar s[];
face vector f[];
s[left] = dirichlet(x + y);
s[top] = dirichlet(x + y);
s[right] = dirichlet(x + y);
s[bottom] = dirichlet(x + y);
v.t[bottom] = dirichlet(0);
f.n[bottom] = 0.;
foreach() {
s[] = x + y;
v.x[] = v.y[] = 1.;
}
foreach_face()
f.x[] = 1.;
foreach()
foreach_neighbor()
fprintf (stderr, "%g %g %g %g %g\n", x, y, s[], v.x[], v.y[]);
foreach_face(x)
for (int i = -2; i <= 2; i++)
fprintf (stderr, "%g %g %g\n", x, y + i*Delta, f.x[0,i]);
foreach_face(y)
for (int i = -2; i <= 2; i++)
fprintf (stderr, "%g %g %g\n", x + i*Delta, y, f.y[i]);
}
|