sandbox/bugs/wavelet_bug.c
wavelet()
can lockup with MPI
int main() {
init_grid (N);
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--) {