From 5542f02701603b5efdd7538a9b88d3ffc9610978 Mon Sep 17 00:00:00 2001
From: Erik Schnetter <schnetter@cct.lsu.edu>
Date: Fri, 25 Jul 2008 13:21:56 -0500
Subject: CarpetLib: Allow negative array indices in vect<T,D>::index

---
 Carpet/CarpetLib/src/vect.hh | 4 +++-
 1 file changed, 3 insertions(+), 1 deletion(-)

(limited to 'Carpet')

diff --git a/Carpet/CarpetLib/src/vect.hh b/Carpet/CarpetLib/src/vect.hh
index 3c4232525..56898d970 100644
--- a/Carpet/CarpetLib/src/vect.hh
+++ b/Carpet/CarpetLib/src/vect.hh
@@ -442,7 +442,9 @@ inline T index (const vect<T,D>& lsh, const vect<T,D>& ind) {
   T r(0);
   for (int d=D-1; d>=0; --d) {
     assert (lsh[d]>=0);
-    assert (ind[d]>=0 and ind[d]<lsh[d]);
+    // Be generous, and allow relative indices which may be negtive
+    // assert (ind[d]>=0 and ind[d]<lsh[d]);
+    assert (abs(ind[d])<=lsh[d]);
     r = r * lsh[d] + ind[d];
   }
   return r;
-- 
cgit v1.2.3