diff options
Diffstat (limited to 'Carpet/CarpetInterp')
| -rw-r--r-- | Carpet/CarpetInterp/src/interp.cc | 14 |
1 files changed, 12 insertions, 2 deletions
diff --git a/Carpet/CarpetInterp/src/interp.cc b/Carpet/CarpetInterp/src/interp.cc index ae9461698..9ebb438f2 100644 --- a/Carpet/CarpetInterp/src/interp.cc +++ b/Carpet/CarpetInterp/src/interp.cc @@ -502,12 +502,22 @@ namespace CarpetInterp { { vector<bool> filled(tmp.size(), false); for (size_t n=0; n<dist::size(); ++n) { -#pragma omp parallel for + //#pragma omp parallel for for (int i=0; i<recvcnt.AT(n); ++i) { assert (not filled.AT(recvdispl.AT(n)+i)); filled.AT(recvdispl.AT(n)+i) = true; } } + bool error = false; + for (int i=0; i<filled.size(); ++i) { + error = error or not (filled.AT(i)); + } + if (error) { + cerr << "error" << endl; + cerr << "recvdispl: " << recvdispl << endl; + cerr << "recvcnt: " << recvcnt << endl; + cerr << "filled: " << filled << endl; + } #pragma omp parallel for for (int i=0; i<filled.size(); ++i) { assert (filled.AT(i)); @@ -603,7 +613,7 @@ namespace CarpetInterp { { vector<bool> filled(source_map.size(), false); for (size_t n=0; n<dist::size(); ++n) { -#pragma omp parallel for + //#pragma omp parallel for for (int i=0; i<recvcnt.AT(n); ++i) { assert (not filled.AT(recvdispl.AT(n)+i)); filled.AT(recvdispl.AT(n)+i) = true; |