From: Denis Davydov Date: Wed, 4 Feb 2015 08:17:56 +0000 (+0100) Subject: fixed KellyErrorEstimator in 1d X-Git-Tag: v8.3.0-rc1~184^2~21 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=592866c5d2566565e86f3c9f7ef293eab84ce491;p=dealii.git fixed KellyErrorEstimator in 1d --- 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