From: kanschat <kanschat@0785d39b-7218-0410-832d-ea1e28bc413d>
Date: Thu, 20 Mar 2008 21:27:04 +0000 (+0000)
Subject: fix problem with zero quedrature weights in FEPolyTensor
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=6c2d2d0738bdb1e866712bb3ead90b5bbe78dc49;p=dealii-svn.git

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 2b7e92b6db..b16f24bddf 100644
--- 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); 
     }