From: kanschat 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 *element, // gradients, not more for (unsigned int i=0; i plus_quad (diff_points); + const Quadrature plus_quad (diff_points, quadrature.get_weights()); differences[d] = new FEValues (mapping, *element, plus_quad, update_gradients); // now same in minus-direction for (unsigned int i=0; i minus_quad (diff_points); + const Quadrature minus_quad (diff_points, quadrature.get_weights()); differences[d+dim] = new FEValues (mapping, *element, minus_quad, update_gradients); }