From: bangerth Date: Fri, 15 Dec 2006 14:09:35 +0000 (+0000) Subject: Check in forgotten hunk X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=cce2a6d3dde05a4ef5d7fa61ce676382410fdbcd;p=dealii-svn.git Check in forgotten hunk git-svn-id: https://svn.dealii.org/trunk@14254 0785d39b-7218-0410-832d-ea1e28bc413d --- 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); };