1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
| #include "utils.h"
int main()
{
int depth = 6;
origin (-0.5, -0.5, -0.5);
init_grid (8);
refine (level < depth - 2 || level <= depth*(1. - sqrt(x*x + y*y + z*z)));
scalar a[];
vector u[];
double k = 2.*pi;
foreach()
u.x[] = u.y[] = a[] = sin(k*x)*cos(k*y);
FILE * fp = fopen ("gfsi.gfs", "w");
output_gfs (fp);
fclose (fp);
}
|