From ee64c0be3e982f3dee65a80cb758d40f485f8340 Mon Sep 17 00:00:00 2001 From: wolf Date: Mon, 8 Mar 1999 11:59:37 +0000 Subject: [PATCH] Fix an omission: we need the quadrature not only for computations on the boundary, but also if we have varying coefficients. git-svn-id: https://svn.dealii.org/trunk@960 0785d39b-7218-0410-832d-ea1e28bc413d --- deal.II/deal.II/source/numerics/error_estimator.cc | 3 ++- 1 file changed, 2 insertions(+), 1 deletion(-) diff --git a/deal.II/deal.II/source/numerics/error_estimator.cc b/deal.II/deal.II/source/numerics/error_estimator.cc index 4c4a508280..e374bedd48 100644 --- a/deal.II/deal.II/source/numerics/error_estimator.cc +++ b/deal.II/deal.II/source/numerics/error_estimator.cc @@ -84,7 +84,8 @@ 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_normal_vectors)); FEFaceValues fe_face_values_neighbor (fe, quadrature, update_gradients); -- 2.39.5