From 4982c4517d2db6dd78938e92eb8a33a257f594c0 Mon Sep 17 00:00:00 2001 From: Jonathan Robey Date: Sat, 17 Sep 2016 20:28:00 -0700 Subject: [PATCH] Add generalized support points to support_point tests --- tests/fe/fe_support_points_common.h | 25 +++++++++++++++++++++++++ 1 file changed, 25 insertions(+) diff --git a/tests/fe/fe_support_points_common.h b/tests/fe/fe_support_points_common.h index 0d7a63f756..dff05f6714 100644 --- a/tests/fe/fe_support_points_common.h +++ b/tests/fe/fe_support_points_common.h @@ -74,6 +74,31 @@ check_support (const FiniteElement &finel, const char *name) deallog << std::setprecision(3) << qp.point(k) << std::endl; } + + deallog << name << '<' << dim << '>' << " cell generalized support points" << std::endl; + + const std::vector > &cell_g_points = finel.get_generalized_support_points (); + + for (unsigned int k=0; k > &face_g_points = finel.get_generalized_face_support_points (); + const std::vector dummy_g_weights (face_g_points.size()); + + Quadrature qg(face_g_points, dummy_g_weights); + + for (unsigned int i=0; i::faces_per_cell; ++i) + { + std::vector > q_points (qg.get_points().size()); + QProjector::project_to_face (qg, i, q_points); + Quadrature qp(q_points); + deallog << name << '<' << dim << '>' << " face " << i + << " generalized support points" << std::endl; + + for (unsigned int k=0; k