sandbox/bugs/wavelet_bug.c

    wavelet() can lockup with MPI

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

    A fix

    I have send this patch

    [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--) {