Add AssertIndexRange

diff --git a/include/deal.II/matrix_free/operators.h b/include/deal.II/matrix_free/operators.h
index 339e651b18a82a5db83a9d06caced9ef4a02b1e5..771fa19713daf91caa35351979bce20d6a9b6ec0 100644 (file)
--- a/include/deal.II/matrix_free/operators.h
+++ b/include/deal.II/matrix_free/operators.h
@@ -65,6 +65,7 @@ namespace MatrixFreeOperators
                             typename VectorType::BlockType &>::type
     subblock(VectorType &vector, unsigned int block_no)
     {
+      AssertIndexRange(block_no, vector.n_blocks());
       return vector.block(block_no);
     }