From: wolf Date: Mon, 8 Mar 1999 12:06:57 +0000 (+0000) Subject: Fix an omission: we need the quadrature not only for computations on the boundary... X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=c0c0b34794eb825b20c68bbc5c2eac34000089cf;p=dealii-svn.git Fix an omission: we need the quadrature not only for computations on the boundary, but also if we have varying coefficients. Do it right this time! git-svn-id: https://svn.dealii.org/trunk@961 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/deal.II/deal.II/source/numerics/error_estimator.cc b/deal.II/deal.II/source/numerics/error_estimator.cc index e374bedd48..cd72620d65 100644 --- a/deal.II/deal.II/source/numerics/error_estimator.cc +++ b/deal.II/deal.II/source/numerics/error_estimator.cc @@ -84,9 +84,9 @@ void KellyErrorEstimator::estimate_error (const DoFHandler &dof, FEFaceValues fe_face_values_cell (fe, quadrature, UpdateFlags(update_gradients | update_JxW_values | - ((neumann_bc.empty() || + ((!neumann_bc.empty() || (coefficient != 0)) ? - 0 : update_q_points) | + update_q_points : 0) | update_normal_vectors)); FEFaceValues fe_face_values_neighbor (fe, quadrature, update_gradients); FESubfaceValues fe_subface_values (fe, quadrature, update_gradients);