From: Wolfgang Bangerth Date: Sat, 14 Aug 2010 03:38:34 +0000 (+0000) Subject: When defining abbreviations for non-scalar data, use references, X-Git-Tag: v8.0.0~5728 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=189e0c257d9cf18568aa86afd4cfd39e410c165e;p=dealii.git When defining abbreviations for non-scalar data, use references, i.e. use const std::vector > & jacobians = fe_values.get_jacobians (); instead of std::vector > jacobians = fe_values.get_jacobians (); git-svn-id: https://svn.dealii.org/trunk@21658 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/deal.II/deal.II/include/numerics/vectors.templates.h b/deal.II/deal.II/include/numerics/vectors.templates.h index 976202d974..1ecd8b47f9 100644 --- a/deal.II/deal.II/include/numerics/vectors.templates.h +++ b/deal.II/deal.II/include/numerics/vectors.templates.h @@ -2947,16 +2947,18 @@ namespace internals // Initialize the required // objects. - std::vector > jacobians = fe_values.get_jacobians (); + const std::vector > & + jacobians = fe_values.get_jacobians (); + const std::vector > & + quadrature_points = quadrature.get_points (); + std::vector > tangentials (fe_values.n_quadrature_points); - std::vector > quadrature_points = fe_values.get_quadrature_points (); std::vector > values (fe_values.n_quadrature_points, Vector (dim)); // Get boundary function values // at quadrature points. boundary_function.vector_value_list (quadrature_points, values); - quadrature_points = quadrature.get_points (); const unsigned int superdegree = cell->get_fe ().degree; const unsigned int degree = superdegree - 1; @@ -3167,8 +3169,11 @@ namespace internals fe_values.reinit (cell); // Initialize the required objects. - std::vector > jacobians = fe_values.get_jacobians (); - std::vector > quadrature_points = fe_values.get_quadrature_points (); + const std::vector > & + jacobians = fe_values.get_jacobians (); + const std::vector > & + quadrature_points = fe_values.get_quadrature_points (); + std::vector > values (fe_values.n_quadrature_points, Vector (dim)); @@ -3416,16 +3421,18 @@ namespace internals fe_values.reinit (cell); // Initialize the required objects. + const std::vector > & + jacobians = fe_values.get_jacobians (); + const std::vector > & + quadrature_points = quadrature.get_points (); + std::vector > tangentials (fe_values.n_quadrature_points); - std::vector > jacobians = fe_values.get_jacobians (); - std::vector > quadrature_points = fe_values.get_quadrature_points (); std::vector > values (fe_values.n_quadrature_points, Vector (dim)); // Get boundary function values // at quadrature points. boundary_function.vector_value_list (quadrature_points, values); - quadrature_points = quadrature.get_points (); const unsigned int degree = cell->get_fe ().degree - 1; Point shifted_reference_point_1;