1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
| #include "grid/octree.h"
#include "output_vtu_foreach.h"
int main()
{
L0 = 1;
X0 = Y0 = Z0 = -L0 / 2;
N = 1 << 3;
init_grid(N);
scalar s[];
vector u[];
foreach ()
{
s[] = pid();
foreach_dimension()
{
u.x[] = sq(x) + sq(y) + sq(z);
}
}
boundary({s, u});
output_vtu((scalar *){s}, (vector *){u}, "test");
}
|