From: Wolfgang Bangerth Date: Fri, 3 May 2024 17:17:19 +0000 (+0530) Subject: Mark FEValuesExtractors objects as 'const' in the tests/ directory. X-Git-Tag: v9.6.0-rc1~320^2 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=9ec7a5c7ea7e13022e13edc2b3181bb966ed7814;p=dealii.git Mark FEValuesExtractors objects as 'const' in the tests/ directory. --- diff --git a/tests/bits/make_boundary_constraints_03.cc b/tests/bits/make_boundary_constraints_03.cc index a42e9774ad..41f232b276 100644 --- a/tests/bits/make_boundary_constraints_03.cc +++ b/tests/bits/make_boundary_constraints_03.cc @@ -60,9 +60,9 @@ test() // DoFRenumbering::component_wise(dof_handler); deallog << "Number of dofs: " << dof_handler.n_dofs() << std::endl; - AffineConstraints constraints; - FEValuesExtractors::Vector velocities(1); - ComponentMask mask = fe.component_mask(velocities); + AffineConstraints constraints; + const FEValuesExtractors::Vector velocities(1); + ComponentMask mask = fe.component_mask(velocities); deallog << "ComponentMask " << mask[0] << mask[1] << mask[2] << std::endl; DoFTools::make_zero_boundary_constraints(dof_handler, constraints, mask); diff --git a/tests/bits/step-51.cc b/tests/bits/step-51.cc index 4bc4cac7d1..9324e2b960 100644 --- a/tests/bits/step-51.cc +++ b/tests/bits/step-51.cc @@ -891,8 +891,8 @@ namespace Step51 scratch.fe_values_local.reinit(loc_cell); scratch.fe_values.reinit(cell); - FEValuesExtractors::Vector fluxes(0); - FEValuesExtractors::Scalar scalar(dim); + const FEValuesExtractors::Vector fluxes(0); + const FEValuesExtractors::Scalar scalar(dim); const unsigned int n_q_points = scratch.fe_values.get_quadrature().size(); const unsigned int dofs_per_cell = scratch.fe_values.dofs_per_cell; diff --git a/tests/bits/step-51p.cc b/tests/bits/step-51p.cc index 571f78bf11..ceb6e2efea 100644 --- a/tests/bits/step-51p.cc +++ b/tests/bits/step-51p.cc @@ -889,8 +889,8 @@ namespace Step51 scratch.fe_values_local.reinit(loc_cell); scratch.fe_values.reinit(cell); - FEValuesExtractors::Vector fluxes(0); - FEValuesExtractors::Scalar scalar(dim); + const FEValuesExtractors::Vector fluxes(0); + const FEValuesExtractors::Scalar scalar(dim); const unsigned int n_q_points = scratch.fe_values.get_quadrature().size(); const unsigned int dofs_per_cell = scratch.fe_values.dofs_per_cell; diff --git a/tests/dofs/dof_tools_21_b.cc b/tests/dofs/dof_tools_21_b.cc index 5692cf3b68..9dc107e560 100644 --- a/tests/dofs/dof_tools_21_b.cc +++ b/tests/dofs/dof_tools_21_b.cc @@ -254,11 +254,11 @@ print_matching(DoFHandler &dof_handler, dof_handler, support_points); - FEValuesExtractors::Vector v(0); - FEValuesExtractors::Scalar v_1(0); - FEValuesExtractors::Scalar v_2(1); - FEValuesExtractors::Scalar v_3(2); - FEValuesExtractors::Scalar pressure(3); + const FEValuesExtractors::Vector v(0); + const FEValuesExtractors::Scalar v_1(0); + const FEValuesExtractors::Scalar v_2(1); + const FEValuesExtractors::Scalar v_3(2); + const FEValuesExtractors::Scalar pressure(3); ComponentMask velocity_mask; if (constrain_only_velocity) diff --git a/tests/dofs/dof_tools_21_c.cc b/tests/dofs/dof_tools_21_c.cc index 244219b23e..1b7def086a 100644 --- a/tests/dofs/dof_tools_21_c.cc +++ b/tests/dofs/dof_tools_21_c.cc @@ -260,11 +260,11 @@ print_matching(DoFHandler &dof_handler, dof_handler, support_points); - FEValuesExtractors::Vector v(0); - FEValuesExtractors::Scalar v_1(0); - FEValuesExtractors::Scalar v_2(1); - FEValuesExtractors::Scalar v_3(2); - FEValuesExtractors::Scalar pressure(3); + const FEValuesExtractors::Vector v(0); + const FEValuesExtractors::Scalar v_1(0); + const FEValuesExtractors::Scalar v_2(1); + const FEValuesExtractors::Scalar v_3(2); + const FEValuesExtractors::Scalar pressure(3); ComponentMask velocity_mask; if (constrain_only_velocity) diff --git a/tests/fe/fe_abf_gradient_divergence_theorem.cc b/tests/fe/fe_abf_gradient_divergence_theorem.cc index bbc9bfd7d8..2f7525a87e 100644 --- a/tests/fe/fe_abf_gradient_divergence_theorem.cc +++ b/tests/fe/fe_abf_gradient_divergence_theorem.cc @@ -96,7 +96,7 @@ test(const Triangulation &tr, for (unsigned int c = 0; c < fe.n_components(); ++c) { - FEValuesExtractors::Scalar single_component(c); + const FEValuesExtractors::Scalar single_component(c); for (unsigned int i = 0; i < fe_values.dofs_per_cell; ++i) { diff --git a/tests/fe/fe_br.cc b/tests/fe/fe_br.cc index 3e15976cbf..e3bc34075c 100644 --- a/tests/fe/fe_br.cc +++ b/tests/fe/fe_br.cc @@ -85,7 +85,7 @@ test(const Triangulation &tr, for (unsigned int c = 0; c < fe.n_components(); ++c) { - FEValuesExtractors::Scalar single_component(c); + const FEValuesExtractors::Scalar single_component(c); for (unsigned int i = 0; i < fe_values.dofs_per_cell; ++i) { diff --git a/tests/fe/fe_dgp_3rd_derivative_divergence_theorem.cc b/tests/fe/fe_dgp_3rd_derivative_divergence_theorem.cc index 0e6778e399..15a2b29eae 100644 --- a/tests/fe/fe_dgp_3rd_derivative_divergence_theorem.cc +++ b/tests/fe/fe_dgp_3rd_derivative_divergence_theorem.cc @@ -94,7 +94,7 @@ test(const Triangulation &tr, for (unsigned int c = 0; c < fe.n_components(); ++c) { - FEValuesExtractors::Scalar single_component(c); + const FEValuesExtractors::Scalar single_component(c); for (unsigned int i = 0; i < fe_values.dofs_per_cell; ++i) { diff --git a/tests/fe/fe_dgq_3rd_derivative_divergence_theorem.cc b/tests/fe/fe_dgq_3rd_derivative_divergence_theorem.cc index 1d7ea3de62..b9889c6cd6 100644 --- a/tests/fe/fe_dgq_3rd_derivative_divergence_theorem.cc +++ b/tests/fe/fe_dgq_3rd_derivative_divergence_theorem.cc @@ -93,7 +93,7 @@ test(const Triangulation &tr, for (unsigned int c = 0; c < fe.n_components(); ++c) { - FEValuesExtractors::Scalar single_component(c); + const FEValuesExtractors::Scalar single_component(c); for (unsigned int i = 0; i < fe_values.dofs_per_cell; ++i) { diff --git a/tests/fe/fe_dgq_gradient_divergence_theorem.cc b/tests/fe/fe_dgq_gradient_divergence_theorem.cc index f6f66ab0b6..8f5d2b35f6 100644 --- a/tests/fe/fe_dgq_gradient_divergence_theorem.cc +++ b/tests/fe/fe_dgq_gradient_divergence_theorem.cc @@ -95,7 +95,7 @@ test(const Triangulation &tr, for (unsigned int c = 0; c < fe.n_components(); ++c) { - FEValuesExtractors::Scalar single_component(c); + const FEValuesExtractors::Scalar single_component(c); for (unsigned int i = 0; i < fe_values.dofs_per_cell; ++i) { diff --git a/tests/fe/fe_dgq_hessian_divergence_theorem.cc b/tests/fe/fe_dgq_hessian_divergence_theorem.cc index 661bb8df4b..781e1a1eca 100644 --- a/tests/fe/fe_dgq_hessian_divergence_theorem.cc +++ b/tests/fe/fe_dgq_hessian_divergence_theorem.cc @@ -95,7 +95,7 @@ test(const Triangulation &tr, for (unsigned int c = 0; c < fe.n_components(); ++c) { - FEValuesExtractors::Scalar single_component(c); + const FEValuesExtractors::Scalar single_component(c); for (unsigned int i = 0; i < fe_values.dofs_per_cell; ++i) { diff --git a/tests/fe/fe_nedelec_gradient_divergence_theorem.cc b/tests/fe/fe_nedelec_gradient_divergence_theorem.cc index 54b28016f3..3afb3b3a03 100644 --- a/tests/fe/fe_nedelec_gradient_divergence_theorem.cc +++ b/tests/fe/fe_nedelec_gradient_divergence_theorem.cc @@ -95,7 +95,7 @@ test(const Triangulation &tr, for (unsigned int c = 0; c < fe.n_components(); ++c) { - FEValuesExtractors::Scalar single_component(c); + const FEValuesExtractors::Scalar single_component(c); for (unsigned int i = 0; i < fe_values.dofs_per_cell; ++i) { diff --git a/tests/fe/fe_nedelec_hessian_divergence_theorem.cc b/tests/fe/fe_nedelec_hessian_divergence_theorem.cc index 728e6f3bb0..54bfa36fa5 100644 --- a/tests/fe/fe_nedelec_hessian_divergence_theorem.cc +++ b/tests/fe/fe_nedelec_hessian_divergence_theorem.cc @@ -95,7 +95,7 @@ test(const Triangulation &tr, for (unsigned int c = 0; c < fe.n_components(); ++c) { - FEValuesExtractors::Scalar single_component(c); + const FEValuesExtractors::Scalar single_component(c); for (unsigned int i = 0; i < fe_values.dofs_per_cell; ++i) { diff --git a/tests/fe/fe_nedelec_sz_divergence_theorem_hanging_nodes.cc b/tests/fe/fe_nedelec_sz_divergence_theorem_hanging_nodes.cc index a589d43f6d..8303451eac 100644 --- a/tests/fe/fe_nedelec_sz_divergence_theorem_hanging_nodes.cc +++ b/tests/fe/fe_nedelec_sz_divergence_theorem_hanging_nodes.cc @@ -94,7 +94,7 @@ test(const Triangulation &tr, for (unsigned int c = 0; c < fe.n_components(); ++c) { - FEValuesExtractors::Scalar single_component(c); + const FEValuesExtractors::Scalar single_component(c); for (unsigned int i = 0; i < fe_values.dofs_per_cell; ++i) { diff --git a/tests/fe/fe_nedelec_sz_gradient_divergence_theorem.cc b/tests/fe/fe_nedelec_sz_gradient_divergence_theorem.cc index f7b8ddd7a5..b1de9ebaee 100644 --- a/tests/fe/fe_nedelec_sz_gradient_divergence_theorem.cc +++ b/tests/fe/fe_nedelec_sz_gradient_divergence_theorem.cc @@ -94,7 +94,7 @@ test(const Triangulation &tr, for (unsigned int c = 0; c < fe.n_components(); ++c) { - FEValuesExtractors::Scalar single_component(c); + const FEValuesExtractors::Scalar single_component(c); for (unsigned int i = 0; i < fe_values.dofs_per_cell; ++i) { diff --git a/tests/fe/fe_nedelec_sz_hessian_divergence_theorem.cc b/tests/fe/fe_nedelec_sz_hessian_divergence_theorem.cc index e4ad417389..8318d9d6bf 100644 --- a/tests/fe/fe_nedelec_sz_hessian_divergence_theorem.cc +++ b/tests/fe/fe_nedelec_sz_hessian_divergence_theorem.cc @@ -85,7 +85,7 @@ test(const Triangulation &tr, for (unsigned int c = 0; c < fe.n_components(); ++c) { - FEValuesExtractors::Scalar single_component(c); + const FEValuesExtractors::Scalar single_component(c); for (unsigned int i = 0; i < fe_values.dofs_per_cell; ++i) { diff --git a/tests/fe/fe_q_3rd_derivative_divergence_theorem.cc b/tests/fe/fe_q_3rd_derivative_divergence_theorem.cc index ef0cdf3c80..4087ab20ac 100644 --- a/tests/fe/fe_q_3rd_derivative_divergence_theorem.cc +++ b/tests/fe/fe_q_3rd_derivative_divergence_theorem.cc @@ -93,7 +93,7 @@ test(const Triangulation &tr, for (unsigned int c = 0; c < fe.n_components(); ++c) { - FEValuesExtractors::Scalar single_component(c); + const FEValuesExtractors::Scalar single_component(c); for (unsigned int i = 0; i < fe_values.dofs_per_cell; ++i) { diff --git a/tests/fe/fe_q_gradient_divergence_theorem.cc b/tests/fe/fe_q_gradient_divergence_theorem.cc index a5039e9637..60baf89eb6 100644 --- a/tests/fe/fe_q_gradient_divergence_theorem.cc +++ b/tests/fe/fe_q_gradient_divergence_theorem.cc @@ -95,7 +95,7 @@ test(const Triangulation &tr, for (unsigned int c = 0; c < fe.n_components(); ++c) { - FEValuesExtractors::Scalar single_component(c); + const FEValuesExtractors::Scalar single_component(c); for (unsigned int i = 0; i < fe_values.dofs_per_cell; ++i) { diff --git a/tests/fe/fe_q_hessian_divergence_theorem.cc b/tests/fe/fe_q_hessian_divergence_theorem.cc index a3635548ab..c93434fc5c 100644 --- a/tests/fe/fe_q_hessian_divergence_theorem.cc +++ b/tests/fe/fe_q_hessian_divergence_theorem.cc @@ -95,7 +95,7 @@ test(const Triangulation &tr, for (unsigned int c = 0; c < fe.n_components(); ++c) { - FEValuesExtractors::Scalar single_component(c); + const FEValuesExtractors::Scalar single_component(c); for (unsigned int i = 0; i < fe_values.dofs_per_cell; ++i) { diff --git a/tests/fe/fe_q_hierarchical_3rd_derivative_divergence_theorem.cc b/tests/fe/fe_q_hierarchical_3rd_derivative_divergence_theorem.cc index e17fe20228..416c90b25f 100644 --- a/tests/fe/fe_q_hierarchical_3rd_derivative_divergence_theorem.cc +++ b/tests/fe/fe_q_hierarchical_3rd_derivative_divergence_theorem.cc @@ -94,7 +94,7 @@ test(const Triangulation &tr, for (unsigned int c = 0; c < fe.n_components(); ++c) { - FEValuesExtractors::Scalar single_component(c); + const FEValuesExtractors::Scalar single_component(c); for (unsigned int i = 0; i < fe_values.dofs_per_cell; ++i) { diff --git a/tests/fe/fe_rt_gradient_divergence_theorem.cc b/tests/fe/fe_rt_gradient_divergence_theorem.cc index f446c3bf4f..cd37b144af 100644 --- a/tests/fe/fe_rt_gradient_divergence_theorem.cc +++ b/tests/fe/fe_rt_gradient_divergence_theorem.cc @@ -95,7 +95,7 @@ test(const Triangulation &tr, for (unsigned int c = 0; c < fe.n_components(); ++c) { - FEValuesExtractors::Scalar single_component(c); + const FEValuesExtractors::Scalar single_component(c); for (unsigned int i = 0; i < fe_values.dofs_per_cell; ++i) { diff --git a/tests/fe/fe_rt_hessian_divergence_theorem.cc b/tests/fe/fe_rt_hessian_divergence_theorem.cc index 931f8c4477..db2e6de2d3 100644 --- a/tests/fe/fe_rt_hessian_divergence_theorem.cc +++ b/tests/fe/fe_rt_hessian_divergence_theorem.cc @@ -89,7 +89,7 @@ test(const Triangulation &tr, for (unsigned int c = 0; c < fe.n_components(); ++c) { - FEValuesExtractors::Scalar single_component(c); + const FEValuesExtractors::Scalar single_component(c); for (unsigned int i = 0; i < fe_values.dofs_per_cell; ++i) { diff --git a/tests/fe/fe_system_3rd_derivative_divergence_theorem.cc b/tests/fe/fe_system_3rd_derivative_divergence_theorem.cc index 3fb9976933..56892bd989 100644 --- a/tests/fe/fe_system_3rd_derivative_divergence_theorem.cc +++ b/tests/fe/fe_system_3rd_derivative_divergence_theorem.cc @@ -93,7 +93,7 @@ test(const Triangulation &tr, for (unsigned int c = 0; c < fe.n_components(); ++c) { - FEValuesExtractors::Scalar single_component(c); + const FEValuesExtractors::Scalar single_component(c); for (unsigned int i = 0; i < fe_values.dofs_per_cell; ++i) { diff --git a/tests/fe/fe_values_extractor_02.cc b/tests/fe/fe_values_extractor_02.cc index c84a05e9b3..716e4ca128 100644 --- a/tests/fe/fe_values_extractor_02.cc +++ b/tests/fe/fe_values_extractor_02.cc @@ -26,10 +26,10 @@ main() { initlog(); - FEValuesExtractors::Scalar scalar(0); - FEValuesExtractors::Vector vector(1); - FEValuesExtractors::Tensor<2> tensor(3); - FEValuesExtractors::SymmetricTensor<2> symmetric_tensor(4); + const FEValuesExtractors::Scalar scalar(0); + const FEValuesExtractors::Vector vector(1); + const FEValuesExtractors::Tensor<2> tensor(3); + const FEValuesExtractors::SymmetricTensor<2> symmetric_tensor(4); deallog << scalar.get_name() << std::endl diff --git a/tests/fe/fe_values_view_01.cc b/tests/fe/fe_values_view_01.cc index 620199dfe2..7ec93a20c4 100644 --- a/tests/fe/fe_values_view_01.cc +++ b/tests/fe/fe_values_view_01.cc @@ -54,7 +54,7 @@ test(const Triangulation &tr, const FiniteElement &fe) for (unsigned int c = 0; c < fe.n_components(); ++c) { - FEValuesExtractors::Scalar single_component(c); + const FEValuesExtractors::Scalar single_component(c); for (unsigned int i = 0; i < fe_values.dofs_per_cell; ++i) for (const auto q : fe_values.quadrature_point_indices()) diff --git a/tests/fe/fe_values_view_02.cc b/tests/fe/fe_values_view_02.cc index 522d6a38a5..10bde372d2 100644 --- a/tests/fe/fe_values_view_02.cc +++ b/tests/fe/fe_values_view_02.cc @@ -60,7 +60,7 @@ test(const Triangulation &tr, const FiniteElement &fe) // 'c' if (c + dim <= fe.n_components()) { - FEValuesExtractors::Vector vec_components(c); + const FEValuesExtractors::Vector vec_components(c); for (unsigned int i = 0; i < fe_values.dofs_per_cell; ++i) for (const auto q : fe_values.quadrature_point_indices()) diff --git a/tests/fe/fe_values_view_03.cc b/tests/fe/fe_values_view_03.cc index a915e99652..e3d34f3228 100644 --- a/tests/fe/fe_values_view_03.cc +++ b/tests/fe/fe_values_view_03.cc @@ -60,7 +60,7 @@ test(const Triangulation &tr, const FiniteElement &fe) for (unsigned int c = 0; c < fe.n_components(); ++c) { - FEValuesExtractors::Scalar single_component(c); + const FEValuesExtractors::Scalar single_component(c); for (unsigned int i = 0; i < fe_values.dofs_per_cell; ++i) for (const auto q : fe_values.quadrature_point_indices()) diff --git a/tests/fe/fe_values_view_04.cc b/tests/fe/fe_values_view_04.cc index 24a5045877..6dcf482a45 100644 --- a/tests/fe/fe_values_view_04.cc +++ b/tests/fe/fe_values_view_04.cc @@ -65,7 +65,7 @@ test(const Triangulation &tr, const FiniteElement &fe) // 'c' if (c + dim <= fe.n_components()) { - FEValuesExtractors::Vector vec_components(c); + const FEValuesExtractors::Vector vec_components(c); for (unsigned int i = 0; i < fe_values.dofs_per_cell; ++i) for (const auto q : fe_values.quadrature_point_indices()) diff --git a/tests/fe/fe_values_view_05.cc b/tests/fe/fe_values_view_05.cc index 843e7ca09d..3ef48fddc0 100644 --- a/tests/fe/fe_values_view_05.cc +++ b/tests/fe/fe_values_view_05.cc @@ -64,7 +64,7 @@ test(const Triangulation &tr, const FiniteElement &fe) for (unsigned int c = 0; c < fe.n_components(); ++c) { - FEValuesExtractors::Scalar single_component(c); + const FEValuesExtractors::Scalar single_component(c); fe_values[single_component].get_function_values(fe_function, scalar_values); deallog << "component=" << c << std::endl; diff --git a/tests/fe/fe_values_view_06.cc b/tests/fe/fe_values_view_06.cc index 4e8bec3deb..420298d0a9 100644 --- a/tests/fe/fe_values_view_06.cc +++ b/tests/fe/fe_values_view_06.cc @@ -64,7 +64,7 @@ test(const Triangulation &tr, const FiniteElement &fe) for (unsigned int c = 0; c < fe.n_components(); ++c) { - FEValuesExtractors::Scalar single_component(c); + const FEValuesExtractors::Scalar single_component(c); fe_values[single_component].get_function_gradients(fe_function, scalar_values); deallog << "component=" << c << std::endl; diff --git a/tests/fe/fe_values_view_07.cc b/tests/fe/fe_values_view_07.cc index 217d646dc8..57285232a7 100644 --- a/tests/fe/fe_values_view_07.cc +++ b/tests/fe/fe_values_view_07.cc @@ -64,7 +64,7 @@ test(const Triangulation &tr, const FiniteElement &fe) for (unsigned int c = 0; c < fe.n_components(); ++c) { - FEValuesExtractors::Scalar single_component(c); + const FEValuesExtractors::Scalar single_component(c); fe_values[single_component].get_function_hessians(fe_function, scalar_values); deallog << "component=" << c << std::endl; diff --git a/tests/fe/fe_values_view_08.cc b/tests/fe/fe_values_view_08.cc index c40d8113b0..112c0599c2 100644 --- a/tests/fe/fe_values_view_08.cc +++ b/tests/fe/fe_values_view_08.cc @@ -66,7 +66,7 @@ test(const Triangulation &tr, const FiniteElement &fe) for (unsigned int c = 0; c < fe.n_components(); ++c) { - FEValuesExtractors::Scalar single_component(c); + const FEValuesExtractors::Scalar single_component(c); fe_values[single_component].get_function_values(fe_function, scalar_values); deallog << "component=" << c << std::endl; diff --git a/tests/fe/fe_values_view_09.cc b/tests/fe/fe_values_view_09.cc index 2a26a7ab6f..c0b1106439 100644 --- a/tests/fe/fe_values_view_09.cc +++ b/tests/fe/fe_values_view_09.cc @@ -66,7 +66,7 @@ test(const Triangulation &tr, const FiniteElement &fe) for (unsigned int c = 0; c < fe.n_components(); ++c) { - FEValuesExtractors::Scalar single_component(c); + const FEValuesExtractors::Scalar single_component(c); fe_values[single_component].get_function_gradients(fe_function, scalar_values); deallog << "component=" << c << std::endl; diff --git a/tests/fe/fe_values_view_10.cc b/tests/fe/fe_values_view_10.cc index d0592ae015..499eda9b91 100644 --- a/tests/fe/fe_values_view_10.cc +++ b/tests/fe/fe_values_view_10.cc @@ -66,7 +66,7 @@ test(const Triangulation &tr, const FiniteElement &fe) for (unsigned int c = 0; c < fe.n_components(); ++c) { - FEValuesExtractors::Scalar single_component(c); + const FEValuesExtractors::Scalar single_component(c); fe_values[single_component].get_function_hessians(fe_function, scalar_values); deallog << "component=" << c << std::endl; diff --git a/tests/fe/fe_values_view_10_single_01.cc b/tests/fe/fe_values_view_10_single_01.cc index 23ad0b4740..e3b61eb3f4 100644 --- a/tests/fe/fe_values_view_10_single_01.cc +++ b/tests/fe/fe_values_view_10_single_01.cc @@ -69,7 +69,7 @@ test(const Triangulation &tr, const FiniteElement &fe) for (unsigned int c = 0; c < fe.n_components(); ++c) { - FEValuesExtractors::Scalar single_component(c); + const FEValuesExtractors::Scalar single_component(c); fe_values[single_component].get_function_hessians(fe_function, scalar_values); deallog << "component=" << c << std::endl; diff --git a/tests/fe/fe_values_view_10_single_02.cc b/tests/fe/fe_values_view_10_single_02.cc index 420eaf81b3..152e481dcc 100644 --- a/tests/fe/fe_values_view_10_single_02.cc +++ b/tests/fe/fe_values_view_10_single_02.cc @@ -69,7 +69,7 @@ test(const Triangulation &tr, const FiniteElement &fe) for (unsigned int c = 0; c < fe.n_components(); ++c) { - FEValuesExtractors::Scalar single_component(c); + const FEValuesExtractors::Scalar single_component(c); fe_values[single_component].get_function_hessians(fe_function, scalar_values); deallog << "component=" << c << std::endl; diff --git a/tests/fe/fe_values_view_10_single_03.cc b/tests/fe/fe_values_view_10_single_03.cc index e9db1fd77d..f2eb2a755c 100644 --- a/tests/fe/fe_values_view_10_single_03.cc +++ b/tests/fe/fe_values_view_10_single_03.cc @@ -69,7 +69,7 @@ test(const Triangulation &tr, const FiniteElement &fe) for (unsigned int c = 0; c < fe.n_components(); ++c) { - FEValuesExtractors::Scalar single_component(c); + const FEValuesExtractors::Scalar single_component(c); fe_values[single_component].get_function_hessians(fe_function, scalar_values); deallog << "component=" << c << std::endl; diff --git a/tests/fe/fe_values_view_11.cc b/tests/fe/fe_values_view_11.cc index 14d81d2314..faff4c78a3 100644 --- a/tests/fe/fe_values_view_11.cc +++ b/tests/fe/fe_values_view_11.cc @@ -69,7 +69,7 @@ test(const Triangulation &tr, const FiniteElement &fe) // 'c' if (c + dim <= fe.n_components()) { - FEValuesExtractors::Vector vector_components(c); + const FEValuesExtractors::Vector vector_components(c); fe_values[vector_components].get_function_values( fe_function, selected_vector_values); deallog << "component=" << c << std::endl; diff --git a/tests/fe/fe_values_view_12.cc b/tests/fe/fe_values_view_12.cc index a28ea7132b..2f1aab6ba6 100644 --- a/tests/fe/fe_values_view_12.cc +++ b/tests/fe/fe_values_view_12.cc @@ -69,7 +69,7 @@ test(const Triangulation &tr, const FiniteElement &fe) // 'c' if (c + dim <= fe.n_components()) { - FEValuesExtractors::Vector vector_components(c); + const FEValuesExtractors::Vector vector_components(c); fe_values[vector_components].get_function_gradients( fe_function, selected_vector_values); deallog << "component=" << c << std::endl; diff --git a/tests/fe/fe_values_view_13.cc b/tests/fe/fe_values_view_13.cc index 50d9bafc07..cb97651d50 100644 --- a/tests/fe/fe_values_view_13.cc +++ b/tests/fe/fe_values_view_13.cc @@ -69,7 +69,7 @@ test(const Triangulation &tr, const FiniteElement &fe) // 'c' if (c + dim <= fe.n_components()) { - FEValuesExtractors::Vector vector_components(c); + const FEValuesExtractors::Vector vector_components(c); fe_values[vector_components].get_function_hessians( fe_function, selected_vector_values); deallog << "component=" << c << std::endl; diff --git a/tests/fe/fe_values_view_14.cc b/tests/fe/fe_values_view_14.cc index 16f2f69288..f265910c1e 100644 --- a/tests/fe/fe_values_view_14.cc +++ b/tests/fe/fe_values_view_14.cc @@ -71,7 +71,7 @@ test(const Triangulation &tr, const FiniteElement &fe) // 'c' if (c + dim <= fe.n_components()) { - FEValuesExtractors::Vector vector_components(c); + const FEValuesExtractors::Vector vector_components(c); fe_values[vector_components].get_function_values( fe_function, selected_vector_values); deallog << "component=" << c << std::endl; diff --git a/tests/fe/fe_values_view_15.cc b/tests/fe/fe_values_view_15.cc index f71fd85535..2d8ebd5f6a 100644 --- a/tests/fe/fe_values_view_15.cc +++ b/tests/fe/fe_values_view_15.cc @@ -71,7 +71,7 @@ test(const Triangulation &tr, const FiniteElement &fe) // 'c' if (c + dim <= fe.n_components()) { - FEValuesExtractors::Vector vector_components(c); + const FEValuesExtractors::Vector vector_components(c); fe_values[vector_components].get_function_gradients( fe_function, selected_vector_values); deallog << "component=" << c << std::endl; diff --git a/tests/fe/fe_values_view_16.cc b/tests/fe/fe_values_view_16.cc index fe281b40e3..711844f2a5 100644 --- a/tests/fe/fe_values_view_16.cc +++ b/tests/fe/fe_values_view_16.cc @@ -71,7 +71,7 @@ test(const Triangulation &tr, const FiniteElement &fe) // 'c' if (c + dim <= fe.n_components()) { - FEValuesExtractors::Vector vector_components(c); + const FEValuesExtractors::Vector vector_components(c); fe_values[vector_components].get_function_hessians( fe_function, selected_vector_values); deallog << "component=" << c << std::endl; diff --git a/tests/fe/fe_values_view_17.cc b/tests/fe/fe_values_view_17.cc index 4a7d474213..3a67fb73f7 100644 --- a/tests/fe/fe_values_view_17.cc +++ b/tests/fe/fe_values_view_17.cc @@ -71,7 +71,7 @@ test(const Triangulation &tr, const FiniteElement &fe) // 'c' if (c + dim <= fe.n_components()) { - FEValuesExtractors::Vector vector_components(c); + const FEValuesExtractors::Vector vector_components(c); fe_values[vector_components].get_function_symmetric_gradients( fe_function, selected_vector_values); deallog << "component=" << c << std::endl; diff --git a/tests/fe/fe_values_view_18.cc b/tests/fe/fe_values_view_18.cc index d9a26c5dfc..b3824e8726 100644 --- a/tests/fe/fe_values_view_18.cc +++ b/tests/fe/fe_values_view_18.cc @@ -73,7 +73,7 @@ test(const Triangulation &tr, const FiniteElement &fe) // 'c' if (c + dim <= fe.n_components()) { - FEValuesExtractors::Vector vector_components(c); + const FEValuesExtractors::Vector vector_components(c); fe_values[vector_components].get_function_symmetric_gradients( fe_function, selected_vector_values); deallog << "component=" << c << std::endl; diff --git a/tests/fe/fe_values_view_19.cc b/tests/fe/fe_values_view_19.cc index 69a008e61e..f170a3d96a 100644 --- a/tests/fe/fe_values_view_19.cc +++ b/tests/fe/fe_values_view_19.cc @@ -69,7 +69,7 @@ test(const Triangulation &tr, const FiniteElement &fe) // 'c' if (c + dim <= fe.n_components()) { - FEValuesExtractors::Vector vector_components(c); + const FEValuesExtractors::Vector vector_components(c); fe_values[vector_components].get_function_divergences( fe_function, selected_vector_values); deallog << "component=" << c << std::endl; diff --git a/tests/fe/fe_values_view_20.cc b/tests/fe/fe_values_view_20.cc index 90f8b1a553..fe15f17242 100644 --- a/tests/fe/fe_values_view_20.cc +++ b/tests/fe/fe_values_view_20.cc @@ -71,7 +71,7 @@ test(const Triangulation &tr, const FiniteElement &fe) // 'c' if (c + dim <= fe.n_components()) { - FEValuesExtractors::Vector vector_components(c); + const FEValuesExtractors::Vector vector_components(c); fe_values[vector_components].get_function_divergences( fe_function, selected_vector_values); deallog << "component=" << c << std::endl; diff --git a/tests/fe/fe_values_view_23.cc b/tests/fe/fe_values_view_23.cc index f6c327f9b2..c3f3b30612 100644 --- a/tests/fe/fe_values_view_23.cc +++ b/tests/fe/fe_values_view_23.cc @@ -56,8 +56,8 @@ test(const Triangulation &tr, const FiniteElement &fe) // let the FEValues object compute the // divergences at quadrature points - std::vector> divergences(quadrature.size()); - FEValuesExtractors::SymmetricTensor<2> extractor(0); + std::vector> divergences(quadrature.size()); + const FEValuesExtractors::SymmetricTensor<2> extractor(0); fe_values[extractor].get_function_divergences(fe_function, divergences); // now do the same "by hand" diff --git a/tests/fe/fe_values_view_23_nonsymmetric.cc b/tests/fe/fe_values_view_23_nonsymmetric.cc index 6eb7432174..ba574b4d40 100644 --- a/tests/fe/fe_values_view_23_nonsymmetric.cc +++ b/tests/fe/fe_values_view_23_nonsymmetric.cc @@ -55,8 +55,8 @@ test(const Triangulation &tr, const FiniteElement &fe) // let the FEValues object compute the // divergences at quadrature points - std::vector> divergences(quadrature.size()); - FEValuesExtractors::Tensor<2> extractor(0); + std::vector> divergences(quadrature.size()); + const FEValuesExtractors::Tensor<2> extractor(0); fe_values[extractor].get_function_divergences(fe_function, divergences); // now do the same "by hand" diff --git a/tests/fe/fe_values_view_24.cc b/tests/fe/fe_values_view_24.cc index 447f0effd1..f40d320094 100644 --- a/tests/fe/fe_values_view_24.cc +++ b/tests/fe/fe_values_view_24.cc @@ -56,8 +56,8 @@ test(const Triangulation &tr, const FiniteElement &fe) // let the FEValues object compute the // divergences at quadrature points - std::vector> divergences(quadrature.size()); - FEValuesExtractors::SymmetricTensor<2> extractor(0); + std::vector> divergences(quadrature.size()); + const FEValuesExtractors::SymmetricTensor<2> extractor(0); fe_values[extractor].get_function_divergences(fe_function, divergences); // now do the same "by hand" diff --git a/tests/fe/fe_values_view_24_nonsymmetric.cc b/tests/fe/fe_values_view_24_nonsymmetric.cc index 6c1151e08d..14ed8d11c1 100644 --- a/tests/fe/fe_values_view_24_nonsymmetric.cc +++ b/tests/fe/fe_values_view_24_nonsymmetric.cc @@ -55,8 +55,8 @@ test(const Triangulation &tr, const FiniteElement &fe) // let the FEValues object compute the // divergences at quadrature points - std::vector> divergences(quadrature.size()); - FEValuesExtractors::Tensor<2> extractor(0); + std::vector> divergences(quadrature.size()); + const FEValuesExtractors::Tensor<2> extractor(0); fe_values[extractor].get_function_divergences(fe_function, divergences); for (unsigned int i = 0; i < fe.dofs_per_cell; ++i) diff --git a/tests/fe/fe_values_view_25.cc b/tests/fe/fe_values_view_25.cc index 369099a576..e485ae1ee4 100644 --- a/tests/fe/fe_values_view_25.cc +++ b/tests/fe/fe_values_view_25.cc @@ -58,8 +58,8 @@ test(const Triangulation &tr, const FiniteElement &fe) // let the FEValues object compute the // divergences at quadrature points - std::vector> divergences(quadrature.size()); - FEValuesExtractors::SymmetricTensor<2> extractor(0); + std::vector> divergences(quadrature.size()); + const FEValuesExtractors::SymmetricTensor<2> extractor(0); fe_values[extractor].get_function_divergences(fe_function, divergences); // now do the same "by hand" diff --git a/tests/fe/fe_values_view_25_nonsymmetric.cc b/tests/fe/fe_values_view_25_nonsymmetric.cc index d276d8c747..15885e1c4c 100644 --- a/tests/fe/fe_values_view_25_nonsymmetric.cc +++ b/tests/fe/fe_values_view_25_nonsymmetric.cc @@ -58,8 +58,8 @@ test(const Triangulation &tr, const FiniteElement &fe) // let the FEValues object compute the // divergences at quadrature points - std::vector> divergences(quadrature.size()); - FEValuesExtractors::Tensor<2> extractor(0); + std::vector> divergences(quadrature.size()); + const FEValuesExtractors::Tensor<2> extractor(0); fe_values[extractor].get_function_divergences(fe_function, divergences); // now do the same "by hand" diff --git a/tests/fe/fe_values_view_25b_nonsymmetric.cc b/tests/fe/fe_values_view_25b_nonsymmetric.cc index 223aae9af4..c0192856a0 100644 --- a/tests/fe/fe_values_view_25b_nonsymmetric.cc +++ b/tests/fe/fe_values_view_25b_nonsymmetric.cc @@ -98,7 +98,7 @@ test(const Triangulation &tr, const FiniteElement &fe) // let the FEValues object compute the // divergences at quadrature points - FEValuesExtractors::Tensor<2> extractor(0); + const FEValuesExtractors::Tensor<2> extractor(0); // also compare get_function_values/gradients/divergences // the manual evaluation diff --git a/tests/fe/fe_values_view_26.cc b/tests/fe/fe_values_view_26.cc index f66bcdb799..f3ef945e60 100644 --- a/tests/fe/fe_values_view_26.cc +++ b/tests/fe/fe_values_view_26.cc @@ -77,8 +77,8 @@ test(const Triangulation &tr, const FiniteElement &fe) // divergences at quadrature points std::vector::type> curls( quadrature.size()); - std::vector> grads(quadrature.size()); - FEValuesExtractors::Vector extractor(0); + std::vector> grads(quadrature.size()); + const FEValuesExtractors::Vector extractor(0); fe_values[extractor].get_function_curls(fe_function, curls); fe_values[extractor].get_function_gradients(fe_function, grads); diff --git a/tests/fe/fe_values_view_27.cc b/tests/fe/fe_values_view_27.cc index 9147f02580..224db4967a 100644 --- a/tests/fe/fe_values_view_27.cc +++ b/tests/fe/fe_values_view_27.cc @@ -77,8 +77,8 @@ test(const Triangulation &tr, const FiniteElement &fe) // divergences at quadrature points std::vector::type> curls( quadrature.size()); - std::vector> grads(quadrature.size()); - FEValuesExtractors::Vector extractor(0); + std::vector> grads(quadrature.size()); + const FEValuesExtractors::Vector extractor(0); fe_values[extractor].get_function_curls(fe_function, curls); fe_values[extractor].get_function_gradients(fe_function, grads); diff --git a/tests/fe/fe_values_view_28.cc b/tests/fe/fe_values_view_28.cc index 1c42aabb9a..26357c9dd2 100644 --- a/tests/fe/fe_values_view_28.cc +++ b/tests/fe/fe_values_view_28.cc @@ -99,8 +99,8 @@ test(const Triangulation &tr, const FiniteElement &fe) // divergences at quadrature points std::vector::type> curls( quadrature.size()); - std::vector> grads(quadrature.size()); - FEValuesExtractors::Vector extractor(0); + std::vector> grads(quadrature.size()); + const FEValuesExtractors::Vector extractor(0); fe_values[extractor].get_function_curls(fe_function, curls); fe_values[extractor].get_function_gradients(fe_function, grads); diff --git a/tests/fe/fe_values_view_29.cc b/tests/fe/fe_values_view_29.cc index 2998932c42..0a13dc24bc 100644 --- a/tests/fe/fe_values_view_29.cc +++ b/tests/fe/fe_values_view_29.cc @@ -101,8 +101,8 @@ test(const Triangulation &tr, // divergences at quadrature points std::vector::type> curls( quadrature.size()); - std::vector> grads(quadrature.size()); - FEValuesExtractors::Vector extractor(0); + std::vector> grads(quadrature.size()); + const FEValuesExtractors::Vector extractor(0); fe_values[extractor].get_function_curls(fe_function, curls); fe_values[extractor].get_function_gradients(fe_function, grads); diff --git a/tests/fe/fe_values_view_30.cc b/tests/fe/fe_values_view_30.cc index ea2422c5dd..a6e88b3717 100644 --- a/tests/fe/fe_values_view_30.cc +++ b/tests/fe/fe_values_view_30.cc @@ -149,7 +149,7 @@ test(const Triangulation &tr, const FiniteElement &fe) // 'c' if (c + dim <= fe.n_components()) { - FEValuesExtractors::Vector vector_components(c); + const FEValuesExtractors::Vector vector_components(c); fe_values[vector_components].get_function_third_derivatives( function_vals, selected_vector_values); deallog << "component=" << c << std::endl; diff --git a/tests/fe/fe_values_view_get_function_from_local_dof_values_01.cc b/tests/fe/fe_values_view_get_function_from_local_dof_values_01.cc index c811be67a8..1306ed8485 100644 --- a/tests/fe/fe_values_view_get_function_from_local_dof_values_01.cc +++ b/tests/fe/fe_values_view_get_function_from_local_dof_values_01.cc @@ -455,34 +455,34 @@ test() deallog.push("Scalar"); { - FE_Q fe(degree); - FEValuesExtractors::Scalar extractor(0); + FE_Q fe(degree); + const FEValuesExtractors::Scalar extractor(0); test_extractor(fe, extractor); } deallog.pop(); deallog.push("Vector"); { - FESystem fe(FE_Q(degree), dim); - FEValuesExtractors::Vector extractor(0); + FESystem fe(FE_Q(degree), dim); + const FEValuesExtractors::Vector extractor(0); test_extractor(fe, extractor); } deallog.pop(); deallog.push("SymmetricTensor"); { - FESystem fe(FE_Q(degree), + FESystem fe(FE_Q(degree), SymmetricTensor<2, dim>::n_independent_components); - FEValuesExtractors::SymmetricTensor<2> extractor(0); + const FEValuesExtractors::SymmetricTensor<2> extractor(0); test_extractor(fe, extractor); } deallog.pop(); deallog.push("Tensor"); { - FESystem fe(FE_Q(degree), + FESystem fe(FE_Q(degree), Tensor<2, dim>::n_independent_components); - FEValuesExtractors::Tensor<2> extractor(0); + const FEValuesExtractors::Tensor<2> extractor(0); test_extractor(fe, extractor); } deallog.pop(); diff --git a/tests/fe/fe_values_view_get_function_from_local_dof_values_02.cc b/tests/fe/fe_values_view_get_function_from_local_dof_values_02.cc index ae07ac3276..228c505e30 100644 --- a/tests/fe/fe_values_view_get_function_from_local_dof_values_02.cc +++ b/tests/fe/fe_values_view_get_function_from_local_dof_values_02.cc @@ -296,34 +296,34 @@ test() deallog.push("Scalar"); { - FE_Q fe(degree); - FEValuesExtractors::Scalar extractor(0); + FE_Q fe(degree); + const FEValuesExtractors::Scalar extractor(0); test_extractor(fe, extractor); } deallog.pop(); deallog.push("Vector"); { - FESystem fe(FE_Q(degree), dim); - FEValuesExtractors::Vector extractor(0); + FESystem fe(FE_Q(degree), dim); + const FEValuesExtractors::Vector extractor(0); test_extractor(fe, extractor); } deallog.pop(); deallog.push("SymmetricTensor"); { - FESystem fe(FE_Q(degree), + FESystem fe(FE_Q(degree), SymmetricTensor<2, dim>::n_independent_components); - FEValuesExtractors::SymmetricTensor<2> extractor(0); + const FEValuesExtractors::SymmetricTensor<2> extractor(0); test_extractor(fe, extractor); } deallog.pop(); deallog.push("Tensor"); { - FESystem fe(FE_Q(degree), + FESystem fe(FE_Q(degree), Tensor<2, dim>::n_independent_components); - FEValuesExtractors::Tensor<2> extractor(0); + const FEValuesExtractors::Tensor<2> extractor(0); test_extractor(fe, extractor); } deallog.pop(); diff --git a/tests/fe/fe_values_view_invalid_01.cc b/tests/fe/fe_values_view_invalid_01.cc index 38e69cd993..119d86e0d8 100644 --- a/tests/fe/fe_values_view_invalid_01.cc +++ b/tests/fe/fe_values_view_invalid_01.cc @@ -42,7 +42,7 @@ test(const Triangulation &tr, const FiniteElement &fe) FEValues fe_values(fe, quadrature, update_values); fe_values.reinit(dof.begin_active()); - FEValuesExtractors::Scalar extr; // invalid object + const FEValuesExtractors::Scalar extr; // invalid object try { fe_values[extr]; // invalid access diff --git a/tests/fe/fe_values_view_invalid_02.cc b/tests/fe/fe_values_view_invalid_02.cc index 96e38c5300..bdfb7894cf 100644 --- a/tests/fe/fe_values_view_invalid_02.cc +++ b/tests/fe/fe_values_view_invalid_02.cc @@ -42,7 +42,7 @@ test(const Triangulation &tr, const FiniteElement &fe) FEValues fe_values(fe, quadrature, update_values); fe_values.reinit(dof.begin_active()); - FEValuesExtractors::Vector extr; // invalid object + const FEValuesExtractors::Vector extr; // invalid object try { fe_values[extr]; // invalid access diff --git a/tests/fe/fe_values_view_invalid_03.cc b/tests/fe/fe_values_view_invalid_03.cc index cbdb24ee3f..7acdf1886e 100644 --- a/tests/fe/fe_values_view_invalid_03.cc +++ b/tests/fe/fe_values_view_invalid_03.cc @@ -42,7 +42,7 @@ test(const Triangulation &tr, const FiniteElement &fe) FEValues fe_values(fe, quadrature, update_values); fe_values.reinit(dof.begin_active()); - FEValuesExtractors::SymmetricTensor<2> extr; // invalid object + const FEValuesExtractors::SymmetricTensor<2> extr; // invalid object try { fe_values[extr]; // invalid access diff --git a/tests/fe/fe_values_view_invalid_03_nonsymmetric.cc b/tests/fe/fe_values_view_invalid_03_nonsymmetric.cc index 1d9b3c27c7..f3c1deba63 100644 --- a/tests/fe/fe_values_view_invalid_03_nonsymmetric.cc +++ b/tests/fe/fe_values_view_invalid_03_nonsymmetric.cc @@ -42,7 +42,7 @@ test(const Triangulation &tr, const FiniteElement &fe) FEValues fe_values(fe, quadrature, update_values); fe_values.reinit(dof.begin_active()); - FEValuesExtractors::Tensor<2> extr; // invalid object + const FEValuesExtractors::Tensor<2> extr; // invalid object try { fe_values[extr]; // invalid access diff --git a/tests/fe/fe_values_views_types.cc b/tests/fe/fe_values_views_types.cc index b54bac7892..8e77a5fc46 100644 --- a/tests/fe/fe_values_views_types.cc +++ b/tests/fe/fe_values_views_types.cc @@ -42,9 +42,9 @@ main() { initlog(); - FEValuesExtractors::Scalar scalar(0); - FEValuesExtractors::Vector vector(1); - FEValuesExtractors::Tensor<2> tensor(2); + const FEValuesExtractors::Scalar scalar(0); + const FEValuesExtractors::Vector vector(1); + const FEValuesExtractors::Tensor<2> tensor(2); test<1, 1>(scalar); test<1, 1>(vector); diff --git a/tests/fecoupling/fe_coupling_03.cc b/tests/fecoupling/fe_coupling_03.cc index bdc4f2bc8d..131484121d 100644 --- a/tests/fecoupling/fe_coupling_03.cc +++ b/tests/fecoupling/fe_coupling_03.cc @@ -116,7 +116,7 @@ main() << "n_coupling_dofs = " << fcv.n_coupling_dofs() << std::endl; FullMatrix matrix(fcv.n_coupling_dofs(), fcv.n_coupling_dofs()); - FEValuesExtractors::Scalar scalar(0); + const FEValuesExtractors::Scalar scalar(0); const auto first = fcv.get_first_extractor(scalar); const auto second = fcv.get_second_extractor(scalar); diff --git a/tests/fecoupling/fe_coupling_04.cc b/tests/fecoupling/fe_coupling_04.cc index 2342960c9c..bf88ed1223 100644 --- a/tests/fecoupling/fe_coupling_04.cc +++ b/tests/fecoupling/fe_coupling_04.cc @@ -96,7 +96,7 @@ main() QuadratureCouplingType::reorder); FullMatrix matrix(fcv.n_coupling_dofs(), fcv.n_coupling_dofs()); - FEValuesExtractors::Scalar scalar(0); + const FEValuesExtractors::Scalar scalar(0); const auto bulk = fcv.get_first_extractor(scalar); const auto surface = fcv.get_second_extractor(scalar); diff --git a/tests/fecoupling/fe_coupling_05.cc b/tests/fecoupling/fe_coupling_05.cc index b83e4289ef..0b29df7ebc 100644 --- a/tests/fecoupling/fe_coupling_05.cc +++ b/tests/fecoupling/fe_coupling_05.cc @@ -96,7 +96,7 @@ main() QuadratureCouplingType::matching); FullMatrix matrix(fcv.n_coupling_dofs(), fcv.n_coupling_dofs()); - FEValuesExtractors::Scalar scalar(0); + const FEValuesExtractors::Scalar scalar(0); const auto bulk = fcv.get_first_extractor(scalar); const auto surface = fcv.get_second_extractor(scalar); diff --git a/tests/feinterface/fe_interface_values_10.cc b/tests/feinterface/fe_interface_values_10.cc index b3b9bb41f2..374c5a2e19 100644 --- a/tests/feinterface/fe_interface_values_10.cc +++ b/tests/feinterface/fe_interface_values_10.cc @@ -103,8 +103,8 @@ template void test() { - FESystem fe(FE_DGP(4), 1); - FEValuesExtractors::Scalar extractor(0); + FESystem fe(FE_DGP(4), 1); + const FEValuesExtractors::Scalar extractor(0); Triangulation tria; DoFHandler dof_handler(tria); diff --git a/tests/feinterface/fe_interface_view_01.cc b/tests/feinterface/fe_interface_view_01.cc index f952d58d48..67eec23565 100644 --- a/tests/feinterface/fe_interface_view_01.cc +++ b/tests/feinterface/fe_interface_view_01.cc @@ -68,7 +68,7 @@ test(const Triangulation &tr, const FiniteElement &fe) for (unsigned int c = 0; c < fe.n_components(); ++c) { - FEValuesExtractors::Scalar single_component(c); + const FEValuesExtractors::Scalar single_component(c); for (unsigned int i = 0; i < n_dofs_face; ++i) for (unsigned int q = 0; q < n_q_points; ++q) { diff --git a/tests/feinterface/fe_interface_view_02.cc b/tests/feinterface/fe_interface_view_02.cc index c8e6168058..c6964ca860 100644 --- a/tests/feinterface/fe_interface_view_02.cc +++ b/tests/feinterface/fe_interface_view_02.cc @@ -78,7 +78,7 @@ test(const Triangulation &tr, const FiniteElement &fe) // 'c' if (c + dim <= fe.n_components()) { - FEValuesExtractors::Vector vec_components(c); + const FEValuesExtractors::Vector vec_components(c); for (unsigned int i = 0; i < n_dofs_face; ++i) for (unsigned int q = 0; q < n_q_points; ++q) { diff --git a/tests/feinterface/fe_interface_view_get_function_from_local_dof_values_01.cc b/tests/feinterface/fe_interface_view_get_function_from_local_dof_values_01.cc index f215431899..ebb48135ad 100644 --- a/tests/feinterface/fe_interface_view_get_function_from_local_dof_values_01.cc +++ b/tests/feinterface/fe_interface_view_get_function_from_local_dof_values_01.cc @@ -353,8 +353,8 @@ template void test_jump_function() { - FE_DGP fe(0); - FEValuesExtractors::Scalar extractor(0); + FE_DGP fe(0); + const FEValuesExtractors::Scalar extractor(0); Triangulation tria; DoFHandler dof_handler(tria); @@ -421,16 +421,16 @@ test() deallog.push("FEInterfaceViews Scalar"); { - FE_Q fe(degree); - FEValuesExtractors::Scalar extractor(0); + FE_Q fe(degree); + const FEValuesExtractors::Scalar extractor(0); test_feiv_extractor(fe, extractor); } deallog.pop(); deallog.push("FEInterfaceViews Vector"); { - FESystem fe(FE_Q(degree), dim); - FEValuesExtractors::Vector extractor(0); + FESystem fe(FE_Q(degree), dim); + const FEValuesExtractors::Vector extractor(0); test_feiv_extractor(fe, extractor); } deallog.pop(); diff --git a/tests/feinterface/fe_interface_view_get_function_from_local_dof_values_02.cc b/tests/feinterface/fe_interface_view_get_function_from_local_dof_values_02.cc index ed02bf9b5d..1f02218af0 100644 --- a/tests/feinterface/fe_interface_view_get_function_from_local_dof_values_02.cc +++ b/tests/feinterface/fe_interface_view_get_function_from_local_dof_values_02.cc @@ -269,16 +269,16 @@ test() deallog.push("FEInterfaceViews Scalar"); { - FE_Q fe(degree); - FEValuesExtractors::Scalar extractor(0); + FE_Q fe(degree); + const FEValuesExtractors::Scalar extractor(0); test_feiv_extractor(fe, extractor); } deallog.pop(); deallog.push("FEInterfaceViews Vector"); { - FESystem fe(FE_Q(degree), dim); - FEValuesExtractors::Vector extractor(0); + FESystem fe(FE_Q(degree), dim); + const FEValuesExtractors::Vector extractor(0); test_feiv_extractor(fe, extractor); } deallog.pop(); diff --git a/tests/hp/fe_nothing_15.cc b/tests/hp/fe_nothing_15.cc index a9b3d1c5bc..735990e7f2 100644 --- a/tests/hp/fe_nothing_15.cc +++ b/tests/hp/fe_nothing_15.cc @@ -53,8 +53,8 @@ test() DoFHandler dof_handler(triangulation); dof_handler.distribute_dofs(fe); - FEValues fe_values(fe, QGauss(2), update_values); - FEValuesExtractors::Scalar nothing(1); + FEValues fe_values(fe, QGauss(2), update_values); + const FEValuesExtractors::Scalar nothing(1); fe_values.reinit(dof_handler.begin_active()); for (unsigned int i = 0; i < fe.dofs_per_cell; ++i) for (unsigned int q = 0; q < fe_values.n_quadrature_points; ++q) diff --git a/tests/hp/fe_nothing_16.cc b/tests/hp/fe_nothing_16.cc index 2e80bcf33b..19a105b121 100644 --- a/tests/hp/fe_nothing_16.cc +++ b/tests/hp/fe_nothing_16.cc @@ -60,9 +60,9 @@ test() DoFHandler dof_handler(triangulation); dof_handler.distribute_dofs(fe); - QGauss q(2); - FEFaceValues fe_values(fe, q, update_values); - FEValuesExtractors::Scalar nothing(1); + QGauss q(2); + FEFaceValues fe_values(fe, q, update_values); + const FEValuesExtractors::Scalar nothing(1); fe_values.reinit(dof_handler.begin_active(), 0); // the following (second) call to reinit diff --git a/tests/lac/schur_complement_03.cc b/tests/lac/schur_complement_03.cc index 859a5dd2d5..8579c8f9ed 100644 --- a/tests/lac/schur_complement_03.cc +++ b/tests/lac/schur_complement_03.cc @@ -170,7 +170,7 @@ namespace Step22 DoFRenumbering::component_wise(dof_handler, block_component); { constraints.clear(); - FEValuesExtractors::Vector velocities(0); + const FEValuesExtractors::Vector velocities(0); DoFTools::make_hanging_node_constraints(dof_handler, constraints); VectorTools::interpolate_boundary_values(dof_handler, 1, @@ -339,7 +339,7 @@ namespace Step22 StokesProblem::refine_mesh() { Vector estimated_error_per_cell(triangulation.n_active_cells()); - FEValuesExtractors::Scalar pressure(dim); + const FEValuesExtractors::Scalar pressure(dim); KellyErrorEstimator::estimate( dof_handler, QGauss(degree + 1), diff --git a/tests/matrix_free/point_evaluation_03.cc b/tests/matrix_free/point_evaluation_03.cc index 9ff1bb93ad..fed1f11c13 100644 --- a/tests/matrix_free/point_evaluation_03.cc +++ b/tests/matrix_free/point_evaluation_03.cc @@ -100,7 +100,7 @@ test(const unsigned int degree) std::vector> function_values(unit_points.size()); std::vector> function_gradients(unit_points.size()); - FEValuesExtractors::Vector extractor(0); + const FEValuesExtractors::Vector extractor(0); for (const auto &cell : dof_handler.active_cell_iterators()) { diff --git a/tests/matrix_free/point_evaluation_04.cc b/tests/matrix_free/point_evaluation_04.cc index 76ab4f9ab5..bfc0f28eb6 100644 --- a/tests/matrix_free/point_evaluation_04.cc +++ b/tests/matrix_free/point_evaluation_04.cc @@ -110,8 +110,8 @@ test(const unsigned int degree) std::vector function_values_scalar(unit_points.size()); std::vector> function_gradients_scalar(unit_points.size()); - FEValuesExtractors::Vector extractor(0); - FEValuesExtractors::Scalar extractor_s(dim); + const FEValuesExtractors::Vector extractor(0); + const FEValuesExtractors::Scalar extractor_s(dim); for (const auto &cell : dof_handler.active_cell_iterators()) { diff --git a/tests/matrix_free/point_evaluation_07.cc b/tests/matrix_free/point_evaluation_07.cc index a15033002f..b53e8deb35 100644 --- a/tests/matrix_free/point_evaluation_07.cc +++ b/tests/matrix_free/point_evaluation_07.cc @@ -107,7 +107,7 @@ test(const unsigned int degree) std::vector> function_values(unit_points.size()); std::vector> function_gradients(unit_points.size()); - FEValuesExtractors::Vector extractor(0); + const FEValuesExtractors::Vector extractor(0); for (const auto &cell : dof_handler.active_cell_iterators()) { diff --git a/tests/matrix_free/point_evaluation_08.cc b/tests/matrix_free/point_evaluation_08.cc index 39be324fb6..eb32b9496b 100644 --- a/tests/matrix_free/point_evaluation_08.cc +++ b/tests/matrix_free/point_evaluation_08.cc @@ -117,8 +117,8 @@ test(const unsigned int degree) std::vector function_values_scalar(unit_points.size()); std::vector> function_gradients_scalar(unit_points.size()); - FEValuesExtractors::Vector extractor(0); - FEValuesExtractors::Scalar extractor_s(dim); + const FEValuesExtractors::Vector extractor(0); + const FEValuesExtractors::Scalar extractor_s(dim); for (const auto &cell : dof_handler.active_cell_iterators()) { diff --git a/tests/matrix_free/point_evaluation_13.cc b/tests/matrix_free/point_evaluation_13.cc index 56d2ce9585..1396514550 100644 --- a/tests/matrix_free/point_evaluation_13.cc +++ b/tests/matrix_free/point_evaluation_13.cc @@ -98,7 +98,7 @@ test() std::vector> function_values(unit_points.size()); std::vector> function_gradients(unit_points.size()); - FEValuesExtractors::Vector extractor(0); + const FEValuesExtractors::Vector extractor(0); for (const auto &cell : dof_handler.active_cell_iterators()) { diff --git a/tests/matrix_free/point_evaluation_14.cc b/tests/matrix_free/point_evaluation_14.cc index ca0baf459c..bb5d4a3ad0 100644 --- a/tests/matrix_free/point_evaluation_14.cc +++ b/tests/matrix_free/point_evaluation_14.cc @@ -102,7 +102,7 @@ test(const unsigned int degree) std::vector function_values(unit_points.size()); std::vector> function_gradients(unit_points.size()); - FEValuesExtractors::Scalar extractor(1); + const FEValuesExtractors::Scalar extractor(1); for (const auto &cell : dof_handler.active_cell_iterators()) { diff --git a/tests/matrix_free/point_evaluation_15.cc b/tests/matrix_free/point_evaluation_15.cc index 1f384f9f98..214f172ce5 100644 --- a/tests/matrix_free/point_evaluation_15.cc +++ b/tests/matrix_free/point_evaluation_15.cc @@ -97,7 +97,7 @@ test() std::vector solution_values(fe.dofs_per_cell); std::vector solution_values2(fe.dofs_per_cell); - FEValuesExtractors::Vector extractor(0); + const FEValuesExtractors::Vector extractor(0); for (const auto &cell : dof_handler.active_cell_iterators()) { diff --git a/tests/matrix_free/point_evaluation_16.cc b/tests/matrix_free/point_evaluation_16.cc index 102f8a8f1b..bb704d0df9 100644 --- a/tests/matrix_free/point_evaluation_16.cc +++ b/tests/matrix_free/point_evaluation_16.cc @@ -98,7 +98,7 @@ test(const unsigned int degree) std::vector solution_values(fe.dofs_per_cell); std::vector solution_values2(fe.dofs_per_cell); - FEValuesExtractors::Vector extractor(0); + const FEValuesExtractors::Vector extractor(0); for (const auto &cell : dof_handler.active_cell_iterators()) { diff --git a/tests/matrix_free/point_evaluation_17.cc b/tests/matrix_free/point_evaluation_17.cc index ce83f663f5..2fbeb83daa 100644 --- a/tests/matrix_free/point_evaluation_17.cc +++ b/tests/matrix_free/point_evaluation_17.cc @@ -99,7 +99,7 @@ test() std::vector solution_values(fe.dofs_per_cell); std::vector solution_values2(fe.dofs_per_cell); - FEValuesExtractors::Vector extractor(0); + const FEValuesExtractors::Vector extractor(0); for (const auto &cell : dof_handler.active_cell_iterators()) { diff --git a/tests/matrix_free/point_evaluation_31.cc b/tests/matrix_free/point_evaluation_31.cc index 0a97d929b3..2774784641 100644 --- a/tests/matrix_free/point_evaluation_31.cc +++ b/tests/matrix_free/point_evaluation_31.cc @@ -104,7 +104,7 @@ test(const unsigned int degree) std::vector> function_values(unit_points.size()); std::vector> function_gradients(unit_points.size()); - FEValuesExtractors::Vector extractor(0); + const FEValuesExtractors::Vector extractor(0); for (const auto &cell : dof_handler.active_cell_iterators()) { diff --git a/tests/meshworker/scratch_data_04.cc b/tests/meshworker/scratch_data_04.cc index a94a216530..a81d634ca2 100644 --- a/tests/meshworker/scratch_data_04.cc +++ b/tests/meshworker/scratch_data_04.cc @@ -79,8 +79,8 @@ test() {}; - double H1_norm = 0; - FEValuesExtractors::Scalar scalar(0); + double H1_norm = 0; + const FEValuesExtractors::Scalar scalar(0); ScratchData scratch(fe, quad, cell_flags); CopyData copy; diff --git a/tests/meshworker/scratch_data_05.cc b/tests/meshworker/scratch_data_05.cc index 964dfb83a2..041874719d 100644 --- a/tests/meshworker/scratch_data_05.cc +++ b/tests/meshworker/scratch_data_05.cc @@ -113,7 +113,7 @@ test() using Iterator = decltype(cell); - FEValuesExtractors::Scalar scalar(0); + const FEValuesExtractors::Scalar scalar(0); auto cell_worker = [&rhs_function, &solution, &scalar](const Iterator &cell, ScratchData &s, diff --git a/tests/meshworker/scratch_data_06.cc b/tests/meshworker/scratch_data_06.cc index 92c22e74f9..83e26cbbe9 100644 --- a/tests/meshworker/scratch_data_06.cc +++ b/tests/meshworker/scratch_data_06.cc @@ -114,7 +114,7 @@ test() using Iterator = decltype(cell); - FEValuesExtractors::Scalar scalar(0); + const FEValuesExtractors::Scalar scalar(0); auto cell_worker = [&rhs_function, &solution, &scalar](const Iterator &cell, ScratchData &s, diff --git a/tests/meshworker/scratch_data_07.cc b/tests/meshworker/scratch_data_07.cc index 73cc30a305..7c9a7299e7 100644 --- a/tests/meshworker/scratch_data_07.cc +++ b/tests/meshworker/scratch_data_07.cc @@ -115,7 +115,7 @@ test() using Iterator = decltype(cell); - FEValuesExtractors::Scalar scalar(0); + const FEValuesExtractors::Scalar scalar(0); auto cell_worker = [&rhs_function, &solution, &scalar](const Iterator &cell, ScratchData &s, diff --git a/tests/meshworker/scratch_data_09a.cc b/tests/meshworker/scratch_data_09a.cc index c7d75df009..52f5c3f257 100644 --- a/tests/meshworker/scratch_data_09a.cc +++ b/tests/meshworker/scratch_data_09a.cc @@ -92,7 +92,7 @@ main(int argc, char *argv[]) Utilities::MPI::MPI_InitFinalize mpi_initialization( argc, argv, testing_max_num_threads()); - FEValuesExtractors::Scalar extractor(0); + const FEValuesExtractors::Scalar extractor(0); run<2>(extractor); run<3>(extractor); diff --git a/tests/meshworker/scratch_data_09a_hp.cc b/tests/meshworker/scratch_data_09a_hp.cc index 2b29137cb4..a2779ae528 100644 --- a/tests/meshworker/scratch_data_09a_hp.cc +++ b/tests/meshworker/scratch_data_09a_hp.cc @@ -100,7 +100,7 @@ main(int argc, char *argv[]) Utilities::MPI::MPI_InitFinalize mpi_initialization( argc, argv, testing_max_num_threads()); - FEValuesExtractors::Scalar extractor(0); + const FEValuesExtractors::Scalar extractor(0); run<2>(extractor); run<3>(extractor); diff --git a/tests/meshworker/scratch_data_09b.cc b/tests/meshworker/scratch_data_09b.cc index f215a1912d..af3bfcfbb9 100644 --- a/tests/meshworker/scratch_data_09b.cc +++ b/tests/meshworker/scratch_data_09b.cc @@ -101,7 +101,7 @@ main(int argc, char *argv[]) Utilities::MPI::MPI_InitFinalize mpi_initialization( argc, argv, testing_max_num_threads()); - FEValuesExtractors::Vector extractor(0); + const FEValuesExtractors::Vector extractor(0); // Just run in 3d, so we can extract the curl run<3>(extractor); diff --git a/tests/meshworker/scratch_data_09b_hp.cc b/tests/meshworker/scratch_data_09b_hp.cc index 18180f8ee4..2240771968 100644 --- a/tests/meshworker/scratch_data_09b_hp.cc +++ b/tests/meshworker/scratch_data_09b_hp.cc @@ -110,7 +110,7 @@ main(int argc, char *argv[]) Utilities::MPI::MPI_InitFinalize mpi_initialization( argc, argv, testing_max_num_threads()); - FEValuesExtractors::Vector extractor(0); + const FEValuesExtractors::Vector extractor(0); // Just run in 3d, so we can extract the curl run<3>(extractor); diff --git a/tests/meshworker/scratch_data_09c.cc b/tests/meshworker/scratch_data_09c.cc index 6e982b1647..b17c6c9d05 100644 --- a/tests/meshworker/scratch_data_09c.cc +++ b/tests/meshworker/scratch_data_09c.cc @@ -90,7 +90,7 @@ main(int argc, char *argv[]) Utilities::MPI::MPI_InitFinalize mpi_initialization( argc, argv, testing_max_num_threads()); - FEValuesExtractors::Tensor<2> extractor(0); + const FEValuesExtractors::Tensor<2> extractor(0); run<2>(extractor); run<3>(extractor); diff --git a/tests/meshworker/scratch_data_09c_hp.cc b/tests/meshworker/scratch_data_09c_hp.cc index d80eaefb7e..74c4943ce9 100644 --- a/tests/meshworker/scratch_data_09c_hp.cc +++ b/tests/meshworker/scratch_data_09c_hp.cc @@ -100,7 +100,7 @@ main(int argc, char *argv[]) Utilities::MPI::MPI_InitFinalize mpi_initialization( argc, argv, testing_max_num_threads()); - FEValuesExtractors::Tensor<2> extractor(0); + const FEValuesExtractors::Tensor<2> extractor(0); run<2>(extractor); run<3>(extractor); diff --git a/tests/meshworker/scratch_data_09d.cc b/tests/meshworker/scratch_data_09d.cc index ff8e224ef3..a793238495 100644 --- a/tests/meshworker/scratch_data_09d.cc +++ b/tests/meshworker/scratch_data_09d.cc @@ -88,7 +88,7 @@ main(int argc, char *argv[]) Utilities::MPI::MPI_InitFinalize mpi_initialization( argc, argv, testing_max_num_threads()); - FEValuesExtractors::SymmetricTensor<2> extractor(0); + const FEValuesExtractors::SymmetricTensor<2> extractor(0); run<2>(extractor); run<3>(extractor); diff --git a/tests/meshworker/scratch_data_09d_hp.cc b/tests/meshworker/scratch_data_09d_hp.cc index a2f35e9d6e..2ac34b15b7 100644 --- a/tests/meshworker/scratch_data_09d_hp.cc +++ b/tests/meshworker/scratch_data_09d_hp.cc @@ -99,7 +99,7 @@ main(int argc, char *argv[]) Utilities::MPI::MPI_InitFinalize mpi_initialization( argc, argv, testing_max_num_threads()); - FEValuesExtractors::SymmetricTensor<2> extractor(0); + const FEValuesExtractors::SymmetricTensor<2> extractor(0); run<2>(extractor); run<3>(extractor); diff --git a/tests/meshworker/scratch_data_10a.cc b/tests/meshworker/scratch_data_10a.cc index ba022f0635..be0ccbbad7 100644 --- a/tests/meshworker/scratch_data_10a.cc +++ b/tests/meshworker/scratch_data_10a.cc @@ -135,7 +135,7 @@ main(int argc, char *argv[]) Utilities::MPI::MPI_InitFinalize mpi_initialization( argc, argv, testing_max_num_threads()); - FEValuesExtractors::Scalar extractor(0); + const FEValuesExtractors::Scalar extractor(0); run<2>(extractor); run<3>(extractor); diff --git a/tests/meshworker/scratch_data_10a_hp.cc b/tests/meshworker/scratch_data_10a_hp.cc index 9eec1a650e..62a4cbfae2 100644 --- a/tests/meshworker/scratch_data_10a_hp.cc +++ b/tests/meshworker/scratch_data_10a_hp.cc @@ -142,7 +142,7 @@ main(int argc, char *argv[]) Utilities::MPI::MPI_InitFinalize mpi_initialization( argc, argv, testing_max_num_threads()); - FEValuesExtractors::Scalar extractor(0); + const FEValuesExtractors::Scalar extractor(0); run<2>(extractor); run<3>(extractor); diff --git a/tests/meshworker/scratch_data_10b.cc b/tests/meshworker/scratch_data_10b.cc index 6f38390abe..c6fe8b0882 100644 --- a/tests/meshworker/scratch_data_10b.cc +++ b/tests/meshworker/scratch_data_10b.cc @@ -136,7 +136,7 @@ main(int argc, char *argv[]) Utilities::MPI::MPI_InitFinalize mpi_initialization( argc, argv, testing_max_num_threads()); - FEValuesExtractors::Vector extractor(0); + const FEValuesExtractors::Vector extractor(0); run<2>(extractor); run<3>(extractor); diff --git a/tests/meshworker/scratch_data_10b_hp.cc b/tests/meshworker/scratch_data_10b_hp.cc index 7f3668c8f1..212c6957c8 100644 --- a/tests/meshworker/scratch_data_10b_hp.cc +++ b/tests/meshworker/scratch_data_10b_hp.cc @@ -144,7 +144,7 @@ main(int argc, char *argv[]) Utilities::MPI::MPI_InitFinalize mpi_initialization( argc, argv, testing_max_num_threads()); - FEValuesExtractors::Vector extractor(0); + const FEValuesExtractors::Vector extractor(0); run<2>(extractor); run<3>(extractor); diff --git a/tests/mpi/muelu_periodicity.cc b/tests/mpi/muelu_periodicity.cc index 0813130359..fc40d1d080 100644 --- a/tests/mpi/muelu_periodicity.cc +++ b/tests/mpi/muelu_periodicity.cc @@ -339,8 +339,8 @@ namespace Step22 { constraints.clear(); - FEValuesExtractors::Vector velocities(0); - FEValuesExtractors::Scalar pressure(dim); + const FEValuesExtractors::Vector velocities(0); + const FEValuesExtractors::Scalar pressure(dim); DoFTools::make_hanging_node_constraints(dof_handler, constraints); #ifdef PERIODIC @@ -711,7 +711,7 @@ namespace Step22 { Vector estimated_error_per_cell(triangulation.n_active_cells()); - FEValuesExtractors::Scalar pressure(dim); + const FEValuesExtractors::Scalar pressure(dim); KellyErrorEstimator::estimate( dof_handler, QGauss(degree + 1), diff --git a/tests/mpi/periodicity_02.cc b/tests/mpi/periodicity_02.cc index 036962ebd9..c50a6ad665 100644 --- a/tests/mpi/periodicity_02.cc +++ b/tests/mpi/periodicity_02.cc @@ -341,8 +341,8 @@ namespace Step22 constraints.clear(); constraints.reinit(locally_relevant_dofs); - FEValuesExtractors::Vector velocities(0); - FEValuesExtractors::Scalar pressure(dim); + const FEValuesExtractors::Vector velocities(0); + const FEValuesExtractors::Scalar pressure(dim); DoFTools::make_hanging_node_constraints(dof_handler, constraints); #ifdef PERIODIC @@ -730,7 +730,7 @@ namespace Step22 { Vector estimated_error_per_cell(triangulation.n_active_cells()); - FEValuesExtractors::Scalar pressure(dim); + const FEValuesExtractors::Scalar pressure(dim); KellyErrorEstimator::estimate( dof_handler, QGauss(degree + 1), diff --git a/tests/mpi/periodicity_03.cc b/tests/mpi/periodicity_03.cc index 3a7e0b9415..59d2f4b6e7 100644 --- a/tests/mpi/periodicity_03.cc +++ b/tests/mpi/periodicity_03.cc @@ -268,8 +268,8 @@ namespace Step22 constraints.clear(); constraints.reinit(locally_relevant_dofs); - FEValuesExtractors::Vector velocities(0); - FEValuesExtractors::Scalar pressure(dim); + const FEValuesExtractors::Vector velocities(0); + const FEValuesExtractors::Scalar pressure(dim); DoFTools::make_hanging_node_constraints(dof_handler, constraints); @@ -812,7 +812,7 @@ namespace Step22 { Vector estimated_error_per_cell(triangulation.n_active_cells()); - FEValuesExtractors::Scalar pressure(dim); + const FEValuesExtractors::Scalar pressure(dim); KellyErrorEstimator::estimate( dof_handler, QGauss(degree + 1), diff --git a/tests/mpi/project_bv_div_conf.cc b/tests/mpi/project_bv_div_conf.cc index 32e15a81c3..ba4e97f963 100644 --- a/tests/mpi/project_bv_div_conf.cc +++ b/tests/mpi/project_bv_div_conf.cc @@ -199,7 +199,7 @@ namespace ResFlow { constraints.reinit(locally_relevant_dofs); - FEValuesExtractors::Vector velocities(0); + const FEValuesExtractors::Vector velocities(0); DoFTools::make_hanging_node_constraints(dof_handler, constraints); VectorTools::project_boundary_values_div_conforming( diff --git a/tests/physics/elasticity-kinematics_01.cc b/tests/physics/elasticity-kinematics_01.cc index ec5cbbf802..26607e98ad 100644 --- a/tests/physics/elasticity-kinematics_01.cc +++ b/tests/physics/elasticity-kinematics_01.cc @@ -73,11 +73,11 @@ test_kinematic_tensors() dot_soln_t -= soln_t1; dot_soln_t *= (1.0 / delta_t); - FEValuesExtractors::Vector u_fe(0); - std::vector> qp_Grad_u_t; - std::vector> qp_Grad_u_t1; - std::vector> qp_dot_Grad_u_t; - std::vector> qp_dot_grad_u_t; + const FEValuesExtractors::Vector u_fe(0); + std::vector> qp_Grad_u_t; + std::vector> qp_Grad_u_t1; + std::vector> qp_dot_Grad_u_t; + std::vector> qp_dot_grad_u_t; FEValues fe_values(fe, qf, update_gradients); MappingQEulerian q1_mapping(1, dof_handler, soln_t); diff --git a/tests/physics/step-18-rotation_matrix.cc b/tests/physics/step-18-rotation_matrix.cc index df90da7131..9a512d2f6b 100644 --- a/tests/physics/step-18-rotation_matrix.cc +++ b/tests/physics/step-18-rotation_matrix.cc @@ -440,7 +440,7 @@ namespace Step18 } system_matrix.compress(VectorOperation::add); system_rhs.compress(VectorOperation::add); - FEValuesExtractors::Scalar z_component(dim - 1); + const FEValuesExtractors::Scalar z_component(dim - 1); std::map boundary_values; VectorTools::interpolate_boundary_values(dof_handler, 0, diff --git a/tests/physics/step-18.cc b/tests/physics/step-18.cc index 67cae3e4cd..a3b1d9871e 100644 --- a/tests/physics/step-18.cc +++ b/tests/physics/step-18.cc @@ -456,7 +456,7 @@ namespace Step18 } system_matrix.compress(VectorOperation::add); system_rhs.compress(VectorOperation::add); - FEValuesExtractors::Scalar z_component(dim - 1); + const FEValuesExtractors::Scalar z_component(dim - 1); std::map boundary_values; VectorTools::interpolate_boundary_values(dof_handler, 0, diff --git a/tests/simplex/step-18.cc b/tests/simplex/step-18.cc index 1e854cbe2a..0b99edef37 100644 --- a/tests/simplex/step-18.cc +++ b/tests/simplex/step-18.cc @@ -638,7 +638,7 @@ namespace Step18 #endif - FEValuesExtractors::Scalar z_component(dim - 1); + const FEValuesExtractors::Scalar z_component(dim - 1); std::map boundary_values; VectorTools::interpolate_boundary_values(mapping, dof_handler, diff --git a/tests/simplex/step-31.cc b/tests/simplex/step-31.cc index 531eaa59eb..c3dfdd47bc 100644 --- a/tests/simplex/step-31.cc +++ b/tests/simplex/step-31.cc @@ -600,8 +600,8 @@ namespace Step31 deallog << " Rebuilding Stokes preconditioner..." << std::flush; assemble_stokes_preconditioner(); Amg_preconditioner = std::make_shared(); - std::vector> constant_modes; - FEValuesExtractors::Vector velocity_components(0); + std::vector> constant_modes; + const FEValuesExtractors::Vector velocity_components(0); DoFTools::extract_constant_modes(stokes_dof_handler, stokes_fe.component_mask( velocity_components), diff --git a/tests/simplex/step-55.cc b/tests/simplex/step-55.cc index f873f3dae0..26cc9c9c32 100644 --- a/tests/simplex/step-55.cc +++ b/tests/simplex/step-55.cc @@ -467,7 +467,7 @@ namespace Step55 { constraints.reinit(locally_relevant_dofs); - FEValuesExtractors::Vector velocities(0); + const FEValuesExtractors::Vector velocities(0); DoFTools::make_hanging_node_constraints(dof_handler, constraints); VectorTools::interpolate_boundary_values(mapping, dof_handler,