From dc7dfc18afe814837af30e3a298a8af37adbe772 Mon Sep 17 00:00:00 2001 From: Peter Munch Date: Wed, 12 May 2021 23:41:06 +0200 Subject: [PATCH] FEPointEvaluation for degree = 0 --- include/deal.II/fe/fe_point_evaluation.h | 7 +++++-- 1 file changed, 5 insertions(+), 2 deletions(-) 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*/) { -- 2.39.5