fix problem with zero quedrature weights in FEPolyTensor

git-svn-id: https://svn.dealii.org/trunk@15922 0785d39b-7218-0410-832d-ea1e28bc413d

diff --git a/deal.II/deal.II/source/fe/fe.cc b/deal.II/deal.II/source/fe/fe.cc
index 2b7e92b6db97a3c209dce97c1a7c4672507b59c2..b16f24bddf76665346da8bc2176076d03028707c 100644 (file)
--- a/deal.II/deal.II/source/fe/fe.cc
+++ b/deal.II/deal.II/source/fe/fe.cc
@@ -73,14 +73,14 @@ InternalDataBase::initialize_2nd (const FiniteElement<dim> *element,
                                       // gradients, not more
       for (unsigned int i=0; i<diff_points.size(); ++i)
        diff_points[i] = quadrature.point(i) + shift;
-      const Quadrature<dim> plus_quad (diff_points);
+      const Quadrature<dim> plus_quad (diff_points, quadrature.get_weights());
       differences[d] = new FEValues<dim> (mapping, *element,
                                          plus_quad, update_gradients);
 
                                       // now same in minus-direction
       for (unsigned int i=0; i<diff_points.size(); ++i)
        diff_points[i] = quadrature.point(i) - shift;
-      const Quadrature<dim> minus_quad (diff_points);
+      const Quadrature<dim> minus_quad (diff_points, quadrature.get_weights());
       differences[d+dim] = new FEValues<dim> (mapping, *element,
                                              minus_quad, update_gradients); 
     }