CarpetLib: add offset3 function identical to index3 without tests

index3 fails with CARPET_DEBUG if it is asked for an index outisde of
the box, which happens in the operators when they try to ask for di
offsets for a box of width 1

1 files changed, 10 insertions, 0 deletions

diff --git a/Carpet/CarpetLib/src/operator_prototypes_3d.hh b/Carpet/CarpetLib/src/operator_prototypes_3d.hh
index b9cb7d343..dd5eb14d9 100644
--- a/Carpet/CarpetLib/src/operator_prototypes_3d.hh
+++ b/Carpet/CarpetLib/src/operator_prototypes_3d.hh
@@ -35,6 +35,16 @@ namespace CarpetLib {
   
   
   
+  static inline
+  size_t
+  offset3 (size_t const i, size_t const j, size_t const k,
+          size_t const exti, size_t const extj, size_t const extk)
+  {
+    return i + exti * (j + extj * k);
+  }
+  
+  
+  
   static int const dim3 = 3;
   
   typedef vect <bool, dim3> bvect3;