From c9fd663a781feb45be8f226d153bab1b4123f3f8 Mon Sep 17 00:00:00 2001
From: Timo Heister <timo.heister@gmail.com>
Date: Fri, 20 Sep 2019 15:07:39 -0400
Subject: [PATCH] make compile

---
 source/base/tensor_product_polynomials.cc | 24 ++++++++++++++++++++++-
 1 file changed, 23 insertions(+), 1 deletion(-)

diff --git a/source/base/tensor_product_polynomials.cc b/source/base/tensor_product_polynomials.cc
index 7db3763f10..d10d94b28c 100644
--- a/source/base/tensor_product_polynomials.cc
+++ b/source/base/tensor_product_polynomials.cc
@@ -164,7 +164,7 @@ TensorProductPolynomials<dim, PolynomialType>::compute_grad(
   const unsigned int i,
   const Point<dim> & p) const
 {
-  unsigned int indices[(dim > 0) ? dim : 1];
+  unsigned int indices[dim];
   compute_index(i, indices);
 
   // compute values and
@@ -195,6 +195,17 @@ TensorProductPolynomials<dim, PolynomialType>::compute_grad(
 
 
 
+template <>
+Tensor<1, 0>
+TensorProductPolynomials<0, Polynomials::Polynomial<double>>::compute_grad(
+  const unsigned int,
+  const Point<0> &) const
+{
+  return Tensor<1, 0>();
+}
+
+
+
 template <int dim, typename PolynomialType>
 Tensor<2, dim>
 TensorProductPolynomials<dim, PolynomialType>::compute_grad_grad(
@@ -240,6 +251,17 @@ TensorProductPolynomials<dim, PolynomialType>::compute_grad_grad(
 
 
 
+template <>
+Tensor<2, 0>
+TensorProductPolynomials<0, Polynomials::Polynomial<double>>::compute_grad_grad(
+  const unsigned int,
+  const Point<0> &) const
+{
+  return Tensor<2, 0>();
+}
+
+
+
 template <int dim, typename PolynomialType>
 void
 TensorProductPolynomials<dim, PolynomialType>::evaluate(
-- 
2.39.5