From cce2a6d3dde05a4ef5d7fa61ce676382410fdbcd Mon Sep 17 00:00:00 2001 From: bangerth Date: Fri, 15 Dec 2006 14:09:35 +0000 Subject: [PATCH] Check in forgotten hunk git-svn-id: https://svn.dealii.org/trunk@14254 0785d39b-7218-0410-832d-ea1e28bc413d --- .../deal.II/include/numerics/error_estimator.h | 17 ++++++++++++++++- 1 file changed, 16 insertions(+), 1 deletion(-) diff --git a/deal.II/deal.II/include/numerics/error_estimator.h b/deal.II/deal.II/include/numerics/error_estimator.h index 421583b003..0c450fd7ef 100644 --- a/deal.II/deal.II/include/numerics/error_estimator.h +++ b/deal.II/deal.II/include/numerics/error_estimator.h @@ -574,6 +574,15 @@ class KellyErrorEstimator * with and needs not allocate * memory itself, or synchronise * with other threads. + * + * The sizes of the arrays are + * initialized with the maximal number of + * entries necessary for the hp + * case. Within the loop over individual + * cells, we then resize the arrays as + * necessary. Since for std::vector + * resizing to a smaller size doesn't + * imply memory allocation, this is fast. */ struct PerThreadData { @@ -657,9 +666,15 @@ class KellyErrorEstimator */ PerThreadData (const unsigned int n_solution_vectors, const unsigned int n_components, - const unsigned int n_q_points, + const unsigned int max_n_q_points, const unsigned int subdomain_id, const unsigned int material_id); + + /** + * Constructor. + */ + void resize (const unsigned int n_components, + const unsigned int max_n_q_points); }; -- 2.39.5