From: Peter Munch <peterrmuench@gmail.com>
Date: Wed, 12 May 2021 21:41:06 +0000 (+0200)
Subject: FEPointEvaluation for degree = 0
X-Git-Tag: v9.3.0-rc1~77^2
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=dc7dfc18afe814837af30e3a298a8af37adbe772;p=dealii.git

FEPointEvaluation for degree = 0
---

diff --git a/include/deal.II/fe/fe_point_evaluation.h b/include/deal.II/fe/fe_point_evaluation.h
index 33bc273475..00dbb6ddda 100644
--- a/include/deal.II/fe/fe_point_evaluation.h
+++ b/include/deal.II/fe/fe_point_evaluation.h
@@ -586,8 +586,11 @@ FEPointEvaluation<n_components, dim, spacedim, Number>::FEPointEvaluation(
       shape_info.reinit(QMidpoint<1>(), fe, base_element_number);
       renumber           = shape_info.lexicographic_numbering;
       dofs_per_component = shape_info.dofs_per_component_on_cell;
-      poly               = Polynomials::generate_complete_Lagrange_basis(
-        QGaussLobatto<1>(shape_info.data[0].fe_degree + 1).get_points());
+      poly =
+        shape_info.data[0].fe_degree == 0 ?
+          Polynomials::LagrangeEquidistant::generate_complete_basis(0) :
+          Polynomials::generate_complete_Lagrange_basis(
+            QGaussLobatto<1>(shape_info.data[0].fe_degree + 1).get_points());
     }
   if (true /*TODO: as long as the fast path of integrate() is not working*/)
     {