From: Wolfgang Bangerth Date: Wed, 5 Aug 1998 15:24:54 +0000 (+0000) Subject: Move the multiplication with the coefficient to the right place. X-Git-Tag: v8.0.0~22777 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=dfe648449bf8c1456fbdeaedffe56ec5147398f6;p=dealii.git Move the multiplication with the coefficient to the right place. git-svn-id: https://svn.dealii.org/trunk@477 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 a03e614147..647482a501 100644 --- a/deal.II/deal.II/source/numerics/error_estimator.cc +++ b/deal.II/deal.II/source/numerics/error_estimator.cc @@ -319,6 +319,17 @@ integrate_over_regular_face (const active_cell_iterator &cell, for (unsigned int point=0; point coefficient_values (n_q_points); + coefficient->value_list (fe_face_values_cell.get_quadrature_points(), + coefficient_values); + for (unsigned int point=0; pointat_boundary() == true) // neumann boundary face. compute @@ -341,17 +352,6 @@ integrate_over_regular_face (const active_cell_iterator &cell, phi[point] -= g[point]; }; - // if a coefficient was given: use that - // to scale the jump in the gradient - if (coefficient != 0) - { - vector coefficient_values (n_q_points); - coefficient->value_list (fe_face_values_cell.get_quadrature_points(), - coefficient_values); - for (unsigned int point=0; point