From b1db4325a5cb19d602d8c0e08f9fed19a7ec8b4b Mon Sep 17 00:00:00 2001 From: Wolfgang Bangerth Date: Mon, 30 Oct 2017 16:35:01 -0600 Subject: [PATCH] Rename variables to make their purpose clearer. --- source/numerics/error_estimator_1d.cc | 14 +++++++------- 1 file changed, 7 insertions(+), 7 deletions(-) diff --git a/source/numerics/error_estimator_1d.cc b/source/numerics/error_estimator_1d.cc index d6482821cb..0211ca4c88 100644 --- a/source/numerics/error_estimator_1d.cc +++ b/source/numerics/error_estimator_1d.cc @@ -309,7 +309,7 @@ estimate (const Mapping<1,spacedim> &mapping, std::vector > > > gradients_neighbor (gradients_here); std::vector::type> > - grad_neighbor (n_solution_vectors, Vector::type>(n_components)); + grad_dot_n_neighbor (n_solution_vectors, Vector::type>(n_components)); // reserve some space for coefficient values at one point. if there is no // coefficient, then we fill it by unity once and for all and don't set it @@ -386,7 +386,7 @@ estimate (const Mapping<1,spacedim> &mapping, // extract the gradient in normal direction of all the components. for (unsigned int s=0; s &mapping, v = neumann_bc.find(n)->second->value(cell->vertex(n)); for (unsigned int s=0; s &mapping, neumann_bc.find(n)->second->vector_value(cell->vertex(n), v); for (unsigned int s=0; s &mapping, { // get gradient here const typename ProductType::type - grad_here = gradients_here[s][n][component] * normal; + grad_dot_n_here = gradients_here[s][n][component] * normal; const typename ProductType::type - jump = ((grad_here - grad_neighbor[s](component)) * + jump = ((grad_dot_n_here - grad_dot_n_neighbor[s](component)) * coefficient_values(component)); (*errors[s])(cell->active_cell_index()) += numbers::NumberTraits::type>::abs_square(jump) * cell->diameter(); -- 2.39.5