From 880150c4b4853788c00cffe3e4223cc59cd73d86 Mon Sep 17 00:00:00 2001 From: wolf Date: Mon, 9 Apr 2001 14:30:17 +0000 Subject: [PATCH] Take over changes 1.58->1.60. git-svn-id: https://svn.dealii.org/branches/Branch-3-1@4411 0785d39b-7218-0410-832d-ea1e28bc413d --- .../source/numerics/error_estimator.cc | 53 ++++++++++++++----- 1 file changed, 40 insertions(+), 13 deletions(-) diff --git a/deal.II/deal.II/source/numerics/error_estimator.cc b/deal.II/deal.II/source/numerics/error_estimator.cc index a026e12f07..d635e60787 100644 --- a/deal.II/deal.II/source/numerics/error_estimator.cc +++ b/deal.II/deal.II/source/numerics/error_estimator.cc @@ -302,9 +302,24 @@ void KellyErrorEstimator<1>::estimate (const DoFHandler<1> &dof_ // if Neumann b.c., then fill // the gradients field which // will be used later on. - for (unsigned int s=0; ssecond->vector_value(cell->vertex(0), - grad_neighbor[s]); + { + if (n_components==1) + { + double v; + neumann_bc.find(n)->second->value(cell->vertex(0), v); + + for (unsigned int s=0; s v(n_components); + neumann_bc.find(n)->second->vector_value(cell->vertex(0), v); + + for (unsigned int s=0; s > g(n_q_points, Vector(n_components)); - data.neumann_bc.find(boundary_indicator)->second - ->vector_value_list (fe_face_values_cell.get_quadrature_points(), - g); - - for (unsigned int n=0; n g(n_q_points); + data.neumann_bc.find(boundary_indicator)->second + ->value_list (fe_face_values_cell.get_quadrature_points(), g); + + for (unsigned int n=0; n > g(n_q_points, Vector(n_components)); + data.neumann_bc.find(boundary_indicator)->second + ->vector_value_list (fe_face_values_cell.get_quadrature_points(), + g); + + for (unsigned int n=0; n