Carpet: Correct error in recomposing

darcs-hash:20070817200016-dae7b-6e8805682638f8255ae53671af861c66be49fc5f.gz

1 files changed, 11 insertions, 1 deletions

diff --git a/Carpet/Carpet/src/Recompose.cc b/Carpet/Carpet/src/Recompose.cc
index b27da0c51..b1301dded 100644
--- a/Carpet/Carpet/src/Recompose.cc
+++ b/Carpet/Carpet/src/Recompose.cc
@@ -983,7 +983,17 @@ namespace Carpet {
             for (int d = 0; d < dim; ++ d) {
               obnds[f][d] = cobnds[f][d].intersects (c);
               if (obnds[f][d]) {
-                assert ((cobnds[f][d] & ibset(c)) == cobnds[f][d]);
+                ivect lo = c.lower();
+                ivect up = c.upper();
+                if (f) lo[d]=up[d]; else up[d]=lo[d];
+                ibbox const cbnds (lo, up, c.stride());
+                if (not ((cobnds[f][d] & ibset(c)) == ibset(cbnds))) {
+                  cout << "cobnds[f][d] = " << cobnds[f][d] << endl
+                       << "ibset(c) = " << ibset(c) << endl
+                       << "(cobnds[f][d] & ibset(c)) = " << (cobnds[f][d] & ibset(c)) << endl
+                       << "ibset(cbnds) = " << ibset(cbnds) << endl;
+                }
+                assert ((cobnds[f][d] & ibset(c)) == ibset(cbnds));
               }
             }
           }