src/test/mask.c
Checks that “half-mapped” fine/coarse cells do not cause trouble with boundary conditions.
#define TRASH 1
#include "grid/quadtree.h"
int main() {
init_grid (N);
scalar s[];
face vector v[];
// mask (y > L0/4. ? top : none);
mask (sq(x - 0.5) + sq(y - 0.5) > sq(0.3) ? top : none);
output_cells (stdout);
for (int l = depth(); l > 0; l--)
foreach_halo (restriction,l)
fprintf (stderr, "%g %g %d\n", x, y, l);
fflush (stderr);
foreach()
s[] = 1.;
boundary ({s});
foreach_face()
v.x[] = 1.;
boundary ((scalar *){v});
restriction ((scalar *){v});
}