From f4c1f1f44756e5daa030ff65ce14fdda124d2f8e Mon Sep 17 00:00:00 2001
From: Timo Heister <timo.heister@gmail.com>
Date: Mon, 27 Jul 2015 17:10:15 -0400
Subject: [PATCH] fix out of bounds access

---
 source/base/tensor_product_polynomials.cc | 12 +++++++++---
 1 file changed, 9 insertions(+), 3 deletions(-)

diff --git a/source/base/tensor_product_polynomials.cc b/source/base/tensor_product_polynomials.cc
index 61423c94bb..aa8ce7230c 100644
--- a/source/base/tensor_product_polynomials.cc
+++ b/source/base/tensor_product_polynomials.cc
@@ -357,13 +357,19 @@ AnisotropicPolynomials<dim>::
 compute_index (const unsigned int i,
                unsigned int       (&indices)[dim]) const
 {
+#ifdef DEBUG
   unsigned int n_poly = 1;
   for (unsigned int d=0; d<dim; ++d)
     n_poly *= polynomials[d].size();
   Assert (i < n_poly, ExcInternalError());
-
-  internal::compute_tensor_index(i, polynomials[0].size(),
-                                 polynomials[1].size(), indices);
+#endif
+
+  if (dim==1)
+    internal::compute_tensor_index(i, polynomials[0].size(),
+                                   0 /*not used*/, indices);
+  else
+    internal::compute_tensor_index(i, polynomials[0].size(),
+                                   polynomials[1].size(), indices);
 }
 
 
-- 
2.39.5