int main() {
  init_grid (N);
  refine (x + y > level/5.); // depth() varies per thread
  scalar s[], w[];
  wavelet (s, w);
}

[Fix a MPI lock bug with `wavelet()`
j.a.v.hooft@gmail.com**20200923124612
 Ignore-this: b1bb5176d1de39ceb4b9238271df49f1
] hunk ./src/grid/multigrid-common.h 89
-  for (int l = depth() - 1; l >= 0; l--) {
+  for (int l = grid->maxdepth - 1; l >= 0; l--) {