diff options
author | Erik Schnetter <schnetter@cct.lsu.edu> | 2008-10-01 17:24:56 -0500 |
---|---|---|
committer | Erik Schnetter <schnetter@cct.lsu.edu> | 2008-10-01 17:24:56 -0500 |
commit | 28bdc89b87b8e6f64c747034c46b1414dcaba224 (patch) | |
tree | 71f93c729b99696ae31db8bbf2031f4bb11809dd /Carpet/CarpetLib/src/bbox.cc | |
parent | f4b007242eaa3a0701957535e6bfedabf99931c1 (diff) |
CarpetLib: Allow expanding empty bboxes
Diffstat (limited to 'Carpet/CarpetLib/src/bbox.cc')
-rw-r--r-- | Carpet/CarpetLib/src/bbox.cc | 3 |
1 files changed, 2 insertions, 1 deletions
diff --git a/Carpet/CarpetLib/src/bbox.cc b/Carpet/CarpetLib/src/bbox.cc index 6183f17e5..a14e1e9a3 100644 --- a/Carpet/CarpetLib/src/bbox.cc +++ b/Carpet/CarpetLib/src/bbox.cc @@ -160,7 +160,8 @@ bool bbox<T,D>::operator> (const bbox& b) const { template<class T, int D> bbox<T,D> bbox<T,D>::expand (const vect<T,D>& lo, const vect<T,D>& hi) const { // Allow expansion only into directions where the extent is not negative - assert (all(lower()<=upper() or (lo==T(0) and hi==T(0)))); + // assert (all(lower()<=upper() or (lo==T(0) and hi==T(0)))); + assert (all(shape()>=vect<T,D>(0) or (lo==T(0) and hi==T(0)))); const vect<T,D> str = stride(); const vect<T,D> lb = lower() - lo * str; const vect<T,D> ub = upper() + hi * str; |