From: Jonathan Robey <class4kayaker@gmail.com>
Date: Sun, 18 Sep 2016 03:28:00 +0000 (-0700)
Subject: Add generalized support points to support_point tests
X-Git-Tag: v8.5.0-rc1~658^2~1
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=4982c4517d2db6dd78938e92eb8a33a257f594c0;p=dealii.git

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 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<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;
+    }
 }