sandbox/bugs/QuadtreeLevelZeroGhostcells.c

Testing the Quadtree Implementation at level Zero for implementation of boundary conditions at Ghost cells.

#define BGHOSTS 2

scalar A[];

int main() {

  init_grid(2);

  foreach_level(0)
    A[] = 1.;

  foreach_dimension() {
    A[left] = neumann (0);
    A[right] = neumann (0);
  }
  boundary_level({A}, 0);

  foreach_level(0)
    foreach_neighbor()
      fprintf (stderr, "%g %g %g\n", x, y, A[]);
}

There are NaNs in some locations.