/* See ../figures/mpi-refine1.svg */
#include "refine_unbalanced.h"
int main (int argc, char * argv[])
{
#if 1
init_grid (16);
#else
init_grid(1);
foreach_cell() {
cell.pid = pid();
cell.flags |= active;
}
tree->dirty = true;
refine_unbalanced (level < 4, NULL);
mpi_partitioning();
#endif
refine_unbalanced (level < 5 && y < 0.315 && (1. - x) < 0.438 && y > 0.25 &&
(1. - x) > 0.375, NULL);
unrefine (y < 0.25 && (1. - x) > 0.5);
refine_unbalanced (level < 5 && y < 0.315 && (1. - x) < 0.5 && y > 0.25 &&
(1. - x) > 0.438, NULL);
scalar s[];
foreach()
s[] = 0.;
boundary ({s});
}