Add generalized support points to support_point tests

diff --git a/tests/fe/fe_support_points_common.h b/tests/fe/fe_support_points_common.h
index 0d7a63f756abed889e49150687204274ac36728f..dff05f6714617a0451fd68452d01a8cca6f21a4b 100644 (file)
--- a/tests/fe/fe_support_points_common.h
+++ b/tests/fe/fe_support_points_common.h
@@ -74,6 +74,31 @@ check_support (const FiniteElement<dim> &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<Point<dim> > &cell_g_points = finel.get_generalized_support_points ();
+
+  for (unsigned int k=0; k<cell_g_points.size(); ++k)
+    deallog << std::setprecision(3) << cell_g_points[k] << std::endl;
+
+  const std::vector<Point<dim-1> > &face_g_points = finel.get_generalized_face_support_points ();
+  const std::vector<double> dummy_g_weights (face_g_points.size());
+
+  Quadrature<dim-1> qg(face_g_points, dummy_g_weights);
+
+  for (unsigned int i=0; i<GeometryInfo<dim>::faces_per_cell; ++i)
+    {
+      std::vector<Point<dim> > q_points (qg.get_points().size());
+      QProjector<dim>::project_to_face (qg, i, q_points);
+      Quadrature<dim> qp(q_points);
+      deallog << name << '<' << dim << '>' << " face " << i
+              << " generalized support points" << std::endl;
+
+      for (unsigned int k=0; k<face_points.size(); ++k)
+        deallog << std::setprecision(3) << qp.point(k)
+                << std::endl;
+    }
 }