/** # Inconsistent 5x5 stencil on root level of trees */ #define BGHOSTS 2 int main() { init_grid(2); scalar A[]; 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. ~~~gnuplot set size ratio -1 plot 'log' u 1:2:3 w labels ~~~ */