1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
| #define BGHOSTS 2
scalar s[];
int main (int argc, char * argv[])
{
init_grid (64);
unrefine ((x < 0.75 || y < 0.75) && level > 2);
output_cells (stdout);
foreach()
s[] = 1.;
foreach()
foreach_neighbor()
assert ((s[] == 1.));
// rebalancing
int nf = 0;
foreach(reduction(+:nf))
nf++;
int npe;
MPI_Comm_size (MPI_COMM_WORLD, &npe);
nf = max(1, nf/npe);
scalar index[];
z_indexing (index, true);
foreach()
fprintf (stderr, "%g %g %d\n", x, y, min(npe - 1, (int)(index[]/nf)));
}
|