FE_Bernstein uses equidistant points

diff --git a/source/fe/fe_bernstein.cc b/source/fe/fe_bernstein.cc
index fdba3c32af1628a5ccf85d2a94079532b4d6d3ee..1e828e207a9769218217ecf16b2bbb58dbe67c28 100644 (file)
--- a/source/fe/fe_bernstein.cc
+++ b/source/fe/fe_bernstein.cc
@@ -86,7 +86,7 @@ get_subface_interpolation_matrix (const FiniteElement<dim,spacedim> &x_source_fe
                ExcInterpolationNotImplemented ()));
 
       const Quadrature<dim-1>
-      quad_face_support(FE_Q<dim,spacedim>(source_fe->degree).get_unit_face_support_points ());
+      quad_face_support(FE_Q<dim,spacedim>(QIterated<1>(QTrapez<1>(),source_fe->degree)).get_unit_face_support_points ());
 
       // Rule of thumb for FP accuracy, that can be expected for a given
       // polynomial degree.  This value is used to cut off values close to