array: more checks

src/array.hpp

@@ -89,6 +89,11 @@ public:
             std::array<size_t, NDim> stride;
             size_t offset;
 
+#ifndef NO_BOUNDS_CHECK
+            if (end < begin || begin >= shape[axis] || end > shape[axis])
+                throw std::out_of_range("begin/end indices out of bounds");
+#endif
+
             offset = offset_ + stride_[axis] * begin;
 
             for (size_t i = 0; i < NDim; ++i) {