From: Peter Munch 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::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*/) {