CarpetLib: Allow empty domains

darcs-hash:20070510222457-dae7b-63f0ce31174188709611001c19feb58f3c516862.gz

1 files changed, 2 insertions, 1 deletions

diff --git a/Carpet/CarpetLib/src/gh.cc b/Carpet/CarpetLib/src/gh.cc
index e85d010d8..66722e11d 100644
--- a/Carpet/CarpetLib/src/gh.cc
+++ b/Carpet/CarpetLib/src/gh.cc
@@ -58,7 +58,8 @@ gh (vector<ivect> const & reffacts_, centering const refcent_,
   for (int ml = 0; ml < (int)baseextents.size(); ++ ml) {
     for (int rl = 0; rl < (int)baseextents.AT(ml).size(); ++ rl) {
       ibbox const & box = baseextents.AT(ml).AT(rl);
-      assert (all (box.shape() / box.stride() >
+      // This condition must hold even for zero-sized grid arrays
+      assert (all (box.shape() / box.stride() >=
                    boundary_width[0] + boundary_width[1]));
     }
   }