From 592866c5d2566565e86f3c9f7ef293eab84ce491 Mon Sep 17 00:00:00 2001 From: Denis Davydov Date: Wed, 4 Feb 2015 09:17:56 +0100 Subject: [PATCH] fixed KellyErrorEstimator in 1d --- source/numerics/error_estimator_1d.cc | 10 +++++----- 1 file changed, 5 insertions(+), 5 deletions(-) diff --git a/source/numerics/error_estimator_1d.cc b/source/numerics/error_estimator_1d.cc index 74e8aae9a0..1cc92477f9 100644 --- a/source/numerics/error_estimator_1d.cc +++ b/source/numerics/error_estimator_1d.cc @@ -300,13 +300,13 @@ estimate (const Mapping<1,spacedim> &mapping, // // for the neighbor gradient, we need several auxiliary fields, depending on // the way we get it (see below) - std::vector > > > + std::vector > > > gradients_here (n_solution_vectors, - std::vector > >(2, std::vector >(n_components))); - std::vector > > > + std::vector > >(2, std::vector >(n_components))); + std::vector > > > gradients_neighbor (gradients_here); - std::vector > - grad_neighbor (n_solution_vectors, Vector(n_components)); + std::vector > + grad_neighbor (n_solution_vectors, Vector(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 -- 2.39.5