]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Mark FEValuesExtractors objects as 'const' in the tests/ directory. 16961/head
authorWolfgang Bangerth <bangerth@colostate.edu>
Fri, 3 May 2024 17:17:19 +0000 (22:47 +0530)
committerWolfgang Bangerth <bangerth@colostate.edu>
Fri, 3 May 2024 17:17:19 +0000 (22:47 +0530)
113 files changed:
tests/bits/make_boundary_constraints_03.cc
tests/bits/step-51.cc
tests/bits/step-51p.cc
tests/dofs/dof_tools_21_b.cc
tests/dofs/dof_tools_21_c.cc
tests/fe/fe_abf_gradient_divergence_theorem.cc
tests/fe/fe_br.cc
tests/fe/fe_dgp_3rd_derivative_divergence_theorem.cc
tests/fe/fe_dgq_3rd_derivative_divergence_theorem.cc
tests/fe/fe_dgq_gradient_divergence_theorem.cc
tests/fe/fe_dgq_hessian_divergence_theorem.cc
tests/fe/fe_nedelec_gradient_divergence_theorem.cc
tests/fe/fe_nedelec_hessian_divergence_theorem.cc
tests/fe/fe_nedelec_sz_divergence_theorem_hanging_nodes.cc
tests/fe/fe_nedelec_sz_gradient_divergence_theorem.cc
tests/fe/fe_nedelec_sz_hessian_divergence_theorem.cc
tests/fe/fe_q_3rd_derivative_divergence_theorem.cc
tests/fe/fe_q_gradient_divergence_theorem.cc
tests/fe/fe_q_hessian_divergence_theorem.cc
tests/fe/fe_q_hierarchical_3rd_derivative_divergence_theorem.cc
tests/fe/fe_rt_gradient_divergence_theorem.cc
tests/fe/fe_rt_hessian_divergence_theorem.cc
tests/fe/fe_system_3rd_derivative_divergence_theorem.cc
tests/fe/fe_values_extractor_02.cc
tests/fe/fe_values_view_01.cc
tests/fe/fe_values_view_02.cc
tests/fe/fe_values_view_03.cc
tests/fe/fe_values_view_04.cc
tests/fe/fe_values_view_05.cc
tests/fe/fe_values_view_06.cc
tests/fe/fe_values_view_07.cc
tests/fe/fe_values_view_08.cc
tests/fe/fe_values_view_09.cc
tests/fe/fe_values_view_10.cc
tests/fe/fe_values_view_10_single_01.cc
tests/fe/fe_values_view_10_single_02.cc
tests/fe/fe_values_view_10_single_03.cc
tests/fe/fe_values_view_11.cc
tests/fe/fe_values_view_12.cc
tests/fe/fe_values_view_13.cc
tests/fe/fe_values_view_14.cc
tests/fe/fe_values_view_15.cc
tests/fe/fe_values_view_16.cc
tests/fe/fe_values_view_17.cc
tests/fe/fe_values_view_18.cc
tests/fe/fe_values_view_19.cc
tests/fe/fe_values_view_20.cc
tests/fe/fe_values_view_23.cc
tests/fe/fe_values_view_23_nonsymmetric.cc
tests/fe/fe_values_view_24.cc
tests/fe/fe_values_view_24_nonsymmetric.cc
tests/fe/fe_values_view_25.cc
tests/fe/fe_values_view_25_nonsymmetric.cc
tests/fe/fe_values_view_25b_nonsymmetric.cc
tests/fe/fe_values_view_26.cc
tests/fe/fe_values_view_27.cc
tests/fe/fe_values_view_28.cc
tests/fe/fe_values_view_29.cc
tests/fe/fe_values_view_30.cc
tests/fe/fe_values_view_get_function_from_local_dof_values_01.cc
tests/fe/fe_values_view_get_function_from_local_dof_values_02.cc
tests/fe/fe_values_view_invalid_01.cc
tests/fe/fe_values_view_invalid_02.cc
tests/fe/fe_values_view_invalid_03.cc
tests/fe/fe_values_view_invalid_03_nonsymmetric.cc
tests/fe/fe_values_views_types.cc
tests/fecoupling/fe_coupling_03.cc
tests/fecoupling/fe_coupling_04.cc
tests/fecoupling/fe_coupling_05.cc
tests/feinterface/fe_interface_values_10.cc
tests/feinterface/fe_interface_view_01.cc
tests/feinterface/fe_interface_view_02.cc
tests/feinterface/fe_interface_view_get_function_from_local_dof_values_01.cc
tests/feinterface/fe_interface_view_get_function_from_local_dof_values_02.cc
tests/hp/fe_nothing_15.cc
tests/hp/fe_nothing_16.cc
tests/lac/schur_complement_03.cc
tests/matrix_free/point_evaluation_03.cc
tests/matrix_free/point_evaluation_04.cc
tests/matrix_free/point_evaluation_07.cc
tests/matrix_free/point_evaluation_08.cc
tests/matrix_free/point_evaluation_13.cc
tests/matrix_free/point_evaluation_14.cc
tests/matrix_free/point_evaluation_15.cc
tests/matrix_free/point_evaluation_16.cc
tests/matrix_free/point_evaluation_17.cc
tests/matrix_free/point_evaluation_31.cc
tests/meshworker/scratch_data_04.cc
tests/meshworker/scratch_data_05.cc
tests/meshworker/scratch_data_06.cc
tests/meshworker/scratch_data_07.cc
tests/meshworker/scratch_data_09a.cc
tests/meshworker/scratch_data_09a_hp.cc
tests/meshworker/scratch_data_09b.cc
tests/meshworker/scratch_data_09b_hp.cc
tests/meshworker/scratch_data_09c.cc
tests/meshworker/scratch_data_09c_hp.cc
tests/meshworker/scratch_data_09d.cc
tests/meshworker/scratch_data_09d_hp.cc
tests/meshworker/scratch_data_10a.cc
tests/meshworker/scratch_data_10a_hp.cc
tests/meshworker/scratch_data_10b.cc
tests/meshworker/scratch_data_10b_hp.cc
tests/mpi/muelu_periodicity.cc
tests/mpi/periodicity_02.cc
tests/mpi/periodicity_03.cc
tests/mpi/project_bv_div_conf.cc
tests/physics/elasticity-kinematics_01.cc
tests/physics/step-18-rotation_matrix.cc
tests/physics/step-18.cc
tests/simplex/step-18.cc
tests/simplex/step-31.cc
tests/simplex/step-55.cc

index a42e9774ad48287833a387f9c05381a705675d28..41f232b27655eea0d2a6a04a02b627d93976f3d9 100644 (file)
@@ -60,9 +60,9 @@ test()
   // DoFRenumbering::component_wise(dof_handler);
   deallog << "Number of dofs: " << dof_handler.n_dofs() << std::endl;
 
-  AffineConstraints<double>  constraints;
-  FEValuesExtractors::Vector velocities(1);
-  ComponentMask              mask = fe.component_mask(velocities);
+  AffineConstraints<double>        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);
index 4bc4cac7d148e3b93ef986db124618b06a903a23..9324e2b9605ce5b334dcb1afeb5dec6e19f9abae 100644 (file)
@@ -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;
index 571f78bf114505a91c673f9f40d4492f6a1e075d..ceb6e2efea880f6e853fdd8e1f93d8c697cbe346 100644 (file)
@@ -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;
index 5692cf3b6886bd9e6fc5832bf9fa327fc581379d..9dc107e5604afa1fffc62a3c4ab9756a5e33ef18 100644 (file)
@@ -254,11 +254,11 @@ print_matching(DoFHandler<dim, spacedim> &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)
index 244219b23e29ffcf3ceba2b1c79d058f9f26a7de..1b7def086a1f8ff45a5a4e8635811482d28c40f8 100644 (file)
@@ -260,11 +260,11 @@ print_matching(DoFHandler<dim, spacedim> &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)
index bbc9bfd7d811f0e6543a2ee14cacdf146343f9af..2f7525a87e0a467e0af3ef8e906a5019454d3a8c 100644 (file)
@@ -96,7 +96,7 @@ test(const Triangulation<dim> &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)
             {
index 3e15976cbfca855c98b583bb9528b2006a6afd54..e3bc34075cea7ffde05349b5d5b91e1b0fe19f1f 100644 (file)
@@ -85,7 +85,7 @@ test(const Triangulation<dim> &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)
             {
index 0e6778e39954963da85f947544bc850ac421c76f..15a2b29eae6831085bcba04327cbf42c3923c65b 100644 (file)
@@ -94,7 +94,7 @@ test(const Triangulation<dim> &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)
             {
index 1d7ea3de62b92afab3144ccff2667ee1d5c69f70..b9889c6cd6695591fe918f9d96986fc7e7919363 100644 (file)
@@ -93,7 +93,7 @@ test(const Triangulation<dim> &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)
             {
index f6f66ab0b6a58dd7a55e69d08268d3d42820ac0b..8f5d2b35f618da17f400b55a462414fa3378c514 100644 (file)
@@ -95,7 +95,7 @@ test(const Triangulation<dim> &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)
             {
index 661bb8df4bfd15c75e7e91624f06a4ae124479c4..781e1a1ecad70e2c9c21500bc3f86b46bde8b3cb 100644 (file)
@@ -95,7 +95,7 @@ test(const Triangulation<dim> &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)
             {
index 54b28016f3bd756c7065033946576efb5c11713e..3afb3b3a033cdd9f02cd6d57112b09954ccc9801 100644 (file)
@@ -95,7 +95,7 @@ test(const Triangulation<dim> &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)
             {
index 728e6f3bb0089524f2f65b3b1a504376ece07e66..54bfa36fa57a5f020e089722c86b4e512b1f80aa 100644 (file)
@@ -95,7 +95,7 @@ test(const Triangulation<dim> &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)
             {
index a589d43f6dc23b581975bd3f1d2fd99930a5e68c..8303451eac31837e0bc10204e15710673be16ffb 100644 (file)
@@ -94,7 +94,7 @@ test(const Triangulation<dim> &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)
             {
index f7b8ddd7a5f1ef45e8dcaccab1889299331035bc..b1de9ebaee3d75f04159fb6fccb4fe1ee2f891f2 100644 (file)
@@ -94,7 +94,7 @@ test(const Triangulation<dim> &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)
             {
index e4ad41738944a112aa585b7774fe0933f1323a39..8318d9d6bff40b1c995d14dfab6204948d287ab3 100644 (file)
@@ -85,7 +85,7 @@ test(const Triangulation<dim> &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)
             {
index ef0cdf3c80a09b82ad7c4ac5a5461b5e1fa710f5..4087ab20acd8ff339a84e4d00dcc254fe761fe43 100644 (file)
@@ -93,7 +93,7 @@ test(const Triangulation<dim> &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)
             {
index a5039e963728785ed1af0df79293f2b5fbb94647..60baf89eb6e5242c82250ec8dae05d487660904b 100644 (file)
@@ -95,7 +95,7 @@ test(const Triangulation<dim> &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)
             {
index a3635548abca851bfda383dc4d1a9a83dfabe38c..c93434fc5ce283153ba17a7ecf9068096bfa5640 100644 (file)
@@ -95,7 +95,7 @@ test(const Triangulation<dim> &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)
             {
index e17fe2022850b441a5fe32b7ea0e09b605b44dcd..416c90b25fd5872f65bc28e48e50f4b1498405ca 100644 (file)
@@ -94,7 +94,7 @@ test(const Triangulation<dim> &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)
             {
index f446c3bf4f3a9df4dd866b2087a9181422aeb0fd..cd37b144af1dbf3ab97505f6f0b044a1fabf452d 100644 (file)
@@ -95,7 +95,7 @@ test(const Triangulation<dim> &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)
             {
index 931f8c44776e60dce1e4e88b8d6abcdc20d16cc9..db2e6de2d3494a32833ab81e17ea7e47112d03ef 100644 (file)
@@ -89,7 +89,7 @@ test(const Triangulation<dim> &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)
             {
index 3fb9976933beccacaf1f4b41831f11511bc65dea..56892bd98973cfbab1cdedb9e52c64b57e067d9c 100644 (file)
@@ -93,7 +93,7 @@ test(const Triangulation<dim> &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)
             {
index c84a05e9b38df21847bb75cea8ab039585fa9db5..716e4ca1283bedfb83f23b3e22f21da4681bc21d 100644 (file)
@@ -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
index 620199dfe2a5fdc5528803c95e95af6835d735cd..7ec93a20c46735083aa658a04073047995aa5396 100644 (file)
@@ -54,7 +54,7 @@ test(const Triangulation<dim> &tr, const FiniteElement<dim> &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())
index 522d6a38a570ff488bb68418bf7bea5f0d4b327d..10bde372d28c845c9c92e839202844990b59ef48 100644 (file)
@@ -60,7 +60,7 @@ test(const Triangulation<dim> &tr, const FiniteElement<dim> &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())
index a915e996525b530558bf0ec9e365b1ab1df61b70..e3d34f3228223ea6be0f18b3e5c11319c744194c 100644 (file)
@@ -60,7 +60,7 @@ test(const Triangulation<dim> &tr, const FiniteElement<dim> &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())
index 24a5045877ac3cbf4053f78cbf0d76844cd8f258..6dcf482a45295988a46260d0843368482a6640e4 100644 (file)
@@ -65,7 +65,7 @@ test(const Triangulation<dim> &tr, const FiniteElement<dim> &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())
index 843e7ca09df2c478edd4eaa73276e69d66028acd..3ef48fddc0fc1bfb07932f2ff3fc368033819334 100644 (file)
@@ -64,7 +64,7 @@ test(const Triangulation<dim> &tr, const FiniteElement<dim> &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;
index 4e8bec3deba41a1a49ad0fac7a3add78041bc906..420298d0a906dfafc15656228a7ee16f3f058201 100644 (file)
@@ -64,7 +64,7 @@ test(const Triangulation<dim> &tr, const FiniteElement<dim> &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;
index 217d646dc8ba7d06fa4c27a4e6001ff4219f8a6c..57285232a74b187bcd55714518b812f8083b8c31 100644 (file)
@@ -64,7 +64,7 @@ test(const Triangulation<dim> &tr, const FiniteElement<dim> &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;
index c40d8113b020d8b4671612c283004d5be565b9ca..112c0599c28025ea7bffeca5c859108151981133 100644 (file)
@@ -66,7 +66,7 @@ test(const Triangulation<dim> &tr, const FiniteElement<dim> &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;
index 2a26a7ab6fe6ebe156358fd45235bfd69f5b8e9f..c0b11064395539ba7352050652c7dd6028100fc0 100644 (file)
@@ -66,7 +66,7 @@ test(const Triangulation<dim> &tr, const FiniteElement<dim> &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;
index d0592ae0155905add8d95580f8ee4ce4a5c9d1b8..499eda9b9173d1b91e3c9261d62fe8f2cdb7f38a 100644 (file)
@@ -66,7 +66,7 @@ test(const Triangulation<dim> &tr, const FiniteElement<dim> &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;
index 23ad0b4740dee74185beeaf3e2758ac08d98133c..e3b61eb3f4cb8a667345a485a2a723eff855464a 100644 (file)
@@ -69,7 +69,7 @@ test(const Triangulation<dim> &tr, const FiniteElement<dim> &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;
index 420eaf81b39d6ca50acae5ce13b71bd31cf0fa6c..152e481dccfa831c98157b2f0d2fa6306dc26678 100644 (file)
@@ -69,7 +69,7 @@ test(const Triangulation<dim> &tr, const FiniteElement<dim> &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;
index e9db1fd77db28536e7542476cecb2338e7f57338..f2eb2a755c5901734b73e5e41253892e4e0dfbdd 100644 (file)
@@ -69,7 +69,7 @@ test(const Triangulation<dim> &tr, const FiniteElement<dim> &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;
index 14d81d2314e7874067fe062b8d306c903e99d4ef..faff4c78a38cd2e6f19a490ef5798f69b6a4555e 100644 (file)
@@ -69,7 +69,7 @@ test(const Triangulation<dim> &tr, const FiniteElement<dim> &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;
index a28ea7132b7df3ef6f06f7f1b76fa7c83952bdf5..2f1aab6ba61c5454e24402ddfd4dd4810e5e270e 100644 (file)
@@ -69,7 +69,7 @@ test(const Triangulation<dim> &tr, const FiniteElement<dim> &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;
index 50d9bafc07034b79e331cc8de41ad75a90320ceb..cb97651d502eff5a9fdee1da630d930cd8bd33e2 100644 (file)
@@ -69,7 +69,7 @@ test(const Triangulation<dim> &tr, const FiniteElement<dim> &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;
index 16f2f6928823e2aa759b726dbc26b44a4733ef97..f265910c1e2e99882f01370f54d329e9dace8e20 100644 (file)
@@ -71,7 +71,7 @@ test(const Triangulation<dim> &tr, const FiniteElement<dim> &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;
index f71fd85535c6a7ebb32ed67aae1b27d98363c95b..2d8ebd5f6a3a033ad65ef1fc1c8be7ce499a7ff3 100644 (file)
@@ -71,7 +71,7 @@ test(const Triangulation<dim> &tr, const FiniteElement<dim> &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;
index fe281b40e3c5d101b38033e79b8c75d1a215154b..711844f2a5b3d50ffd3134b18c3720c129f2816b 100644 (file)
@@ -71,7 +71,7 @@ test(const Triangulation<dim> &tr, const FiniteElement<dim> &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;
index 4a7d47421384d9dba810b37ae59136448284320f..3a67fb73f71a0659be0e7b93178192f0b1a0e564 100644 (file)
@@ -71,7 +71,7 @@ test(const Triangulation<dim> &tr, const FiniteElement<dim> &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;
index d9a26c5dfc712142a200a11f05cd5b5abbf04ae1..b3824e8726c240969e6e248502c981a7f0cba22a 100644 (file)
@@ -73,7 +73,7 @@ test(const Triangulation<dim> &tr, const FiniteElement<dim> &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;
index 69a008e61e3d29adfd16df653b4c1b604b1bd072..f170a3d96a4b126ba4d24ca7720bd18688fe74fa 100644 (file)
@@ -69,7 +69,7 @@ test(const Triangulation<dim> &tr, const FiniteElement<dim> &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;
index 90f8b1a55366a17f61851beb9614f254ed9b9354..fe15f1724233f1fdc83d6115206ac68382a58b79 100644 (file)
@@ -71,7 +71,7 @@ test(const Triangulation<dim> &tr, const FiniteElement<dim> &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;
index f6c327f9b2535350f62b52fec6d6893bf31ecc74..c3f3b306121b6ec0c955f248cd3eb542aeec77fe 100644 (file)
@@ -56,8 +56,8 @@ test(const Triangulation<dim> &tr, const FiniteElement<dim> &fe)
 
   // let the FEValues object compute the
   // divergences at quadrature points
-  std::vector<Tensor<1, dim>>            divergences(quadrature.size());
-  FEValuesExtractors::SymmetricTensor<2> extractor(0);
+  std::vector<Tensor<1, dim>>                  divergences(quadrature.size());
+  const FEValuesExtractors::SymmetricTensor<2> extractor(0);
   fe_values[extractor].get_function_divergences(fe_function, divergences);
 
   // now do the same "by hand"
index 6eb7432174db1af8be5c56da63de1104964cac98..ba574b4d401367a37a53406009524b996267bd28 100644 (file)
@@ -55,8 +55,8 @@ test(const Triangulation<dim> &tr, const FiniteElement<dim> &fe)
 
   // let the FEValues object compute the
   // divergences at quadrature points
-  std::vector<Tensor<1, dim>>   divergences(quadrature.size());
-  FEValuesExtractors::Tensor<2> extractor(0);
+  std::vector<Tensor<1, dim>>         divergences(quadrature.size());
+  const FEValuesExtractors::Tensor<2> extractor(0);
   fe_values[extractor].get_function_divergences(fe_function, divergences);
 
   // now do the same "by hand"
index 447f0effd192c8dc6cd297d72ce9da43a3d61190..f40d3200941086d1b8493374c5023560f1fe9f87 100644 (file)
@@ -56,8 +56,8 @@ test(const Triangulation<dim> &tr, const FiniteElement<dim> &fe)
 
   // let the FEValues object compute the
   // divergences at quadrature points
-  std::vector<Tensor<1, dim>>            divergences(quadrature.size());
-  FEValuesExtractors::SymmetricTensor<2> extractor(0);
+  std::vector<Tensor<1, dim>>                  divergences(quadrature.size());
+  const FEValuesExtractors::SymmetricTensor<2> extractor(0);
   fe_values[extractor].get_function_divergences(fe_function, divergences);
 
   // now do the same "by hand"
index 6c1151e08d7a53d98383ebf9d203fcd06792263d..14ed8d11c1c29d69613ad268e8586bb6c09237ea 100644 (file)
@@ -55,8 +55,8 @@ test(const Triangulation<dim> &tr, const FiniteElement<dim> &fe)
 
   // let the FEValues object compute the
   // divergences at quadrature points
-  std::vector<Tensor<1, dim>>   divergences(quadrature.size());
-  FEValuesExtractors::Tensor<2> extractor(0);
+  std::vector<Tensor<1, dim>>         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)
index 369099a5765cf81d926ba7426bab6e0ee3799a97..e485ae1ee48a6817f330136bc2a89c3c704da1bb 100644 (file)
@@ -58,8 +58,8 @@ test(const Triangulation<dim> &tr, const FiniteElement<dim> &fe)
 
   // let the FEValues object compute the
   // divergences at quadrature points
-  std::vector<Tensor<1, dim>>            divergences(quadrature.size());
-  FEValuesExtractors::SymmetricTensor<2> extractor(0);
+  std::vector<Tensor<1, dim>>                  divergences(quadrature.size());
+  const FEValuesExtractors::SymmetricTensor<2> extractor(0);
   fe_values[extractor].get_function_divergences(fe_function, divergences);
 
   // now do the same "by hand"
index d276d8c7472217d60c58b2f43181c772f05bb9d9..15885e1c4ca6a4a8bcd8c7e1c841a29a0de08fdd 100644 (file)
@@ -58,8 +58,8 @@ test(const Triangulation<dim> &tr, const FiniteElement<dim> &fe)
 
   // let the FEValues object compute the
   // divergences at quadrature points
-  std::vector<Tensor<1, dim>>   divergences(quadrature.size());
-  FEValuesExtractors::Tensor<2> extractor(0);
+  std::vector<Tensor<1, dim>>         divergences(quadrature.size());
+  const FEValuesExtractors::Tensor<2> extractor(0);
   fe_values[extractor].get_function_divergences(fe_function, divergences);
 
   // now do the same "by hand"
index 223aae9af4399093d8a57ebff9daa7e1b622d7ef..c0192856a076e28958a06f43811866b1fbe52611 100644 (file)
@@ -98,7 +98,7 @@ test(const Triangulation<dim> &tr, const FiniteElement<dim> &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
index f66bcdb7996693ff8fc11fd87f2045928b90e84e..f3ef945e60939f80a6d545302230de66eb31f6e5 100644 (file)
@@ -77,8 +77,8 @@ test(const Triangulation<dim> &tr, const FiniteElement<dim> &fe)
   // divergences at quadrature points
   std::vector<typename dealii::internal::CurlType<dim>::type> curls(
     quadrature.size());
-  std::vector<Tensor<2, dim>> grads(quadrature.size());
-  FEValuesExtractors::Vector  extractor(0);
+  std::vector<Tensor<2, dim>>      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);
 
index 9147f02580a0f1ec9b920d7f882d09262dcbbd14..224db4967a8d3acbea362610bdc7120a88abc6c5 100644 (file)
@@ -77,8 +77,8 @@ test(const Triangulation<dim> &tr, const FiniteElement<dim> &fe)
   // divergences at quadrature points
   std::vector<typename dealii::internal::CurlType<dim>::type> curls(
     quadrature.size());
-  std::vector<Tensor<2, dim>> grads(quadrature.size());
-  FEValuesExtractors::Vector  extractor(0);
+  std::vector<Tensor<2, dim>>      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);
 
index 1c42aabb9ad4ff3d27778af863e40d9aa3328cf0..26357c9dd26599899fdfef5f77b48d0997eca11b 100644 (file)
@@ -99,8 +99,8 @@ test(const Triangulation<dim> &tr, const FiniteElement<dim> &fe)
   // divergences at quadrature points
   std::vector<typename dealii::internal::CurlType<dim>::type> curls(
     quadrature.size());
-  std::vector<Tensor<2, dim>> grads(quadrature.size());
-  FEValuesExtractors::Vector  extractor(0);
+  std::vector<Tensor<2, dim>>      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);
 
index 2998932c421861fea65272fe44a4092db249ccc7..0a13dc24bc4c1ac4064c9d9ccc3b9b1db005ac79 100644 (file)
@@ -101,8 +101,8 @@ test(const Triangulation<dim> &tr,
   // divergences at quadrature points
   std::vector<typename dealii::internal::CurlType<dim>::type> curls(
     quadrature.size());
-  std::vector<Tensor<2, dim>> grads(quadrature.size());
-  FEValuesExtractors::Vector  extractor(0);
+  std::vector<Tensor<2, dim>>      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);
 
index ea2422c5dd2190757d5bd48d0efb2d565511dbda..a6e88b37176d6eea2c8bdf96f913363fb4685779 100644 (file)
@@ -149,7 +149,7 @@ test(const Triangulation<dim> &tr, const FiniteElement<dim> &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;
index c811be67a80b4656e9b87c816a0d78bb97f15a70..1306ed848587ef22165936c72bf239ee3ec20990 100644 (file)
@@ -455,34 +455,34 @@ test()
 
   deallog.push("Scalar");
   {
-    FE_Q<dim>                  fe(degree);
-    FEValuesExtractors::Scalar extractor(0);
+    FE_Q<dim>                        fe(degree);
+    const FEValuesExtractors::Scalar extractor(0);
     test_extractor<NumberType, dim>(fe, extractor);
   }
   deallog.pop();
 
   deallog.push("Vector");
   {
-    FESystem<dim>              fe(FE_Q<dim>(degree), dim);
-    FEValuesExtractors::Vector extractor(0);
+    FESystem<dim>                    fe(FE_Q<dim>(degree), dim);
+    const FEValuesExtractors::Vector extractor(0);
     test_extractor<NumberType, dim>(fe, extractor);
   }
   deallog.pop();
 
   deallog.push("SymmetricTensor");
   {
-    FESystem<dim>                          fe(FE_Q<dim>(degree),
+    FESystem<dim>                                fe(FE_Q<dim>(degree),
                      SymmetricTensor<2, dim>::n_independent_components);
-    FEValuesExtractors::SymmetricTensor<2> extractor(0);
+    const FEValuesExtractors::SymmetricTensor<2> extractor(0);
     test_extractor<NumberType, dim>(fe, extractor);
   }
   deallog.pop();
 
   deallog.push("Tensor");
   {
-    FESystem<dim>                 fe(FE_Q<dim>(degree),
+    FESystem<dim>                       fe(FE_Q<dim>(degree),
                      Tensor<2, dim>::n_independent_components);
-    FEValuesExtractors::Tensor<2> extractor(0);
+    const FEValuesExtractors::Tensor<2> extractor(0);
     test_extractor<NumberType, dim>(fe, extractor);
   }
   deallog.pop();
index ae07ac32766eb9df2ce3cb738d1f24dd23607e94..228c505e3011ae63f28a24ba9e6f2f1128cd2208 100644 (file)
@@ -296,34 +296,34 @@ test()
 
   deallog.push("Scalar");
   {
-    FE_Q<dim>                  fe(degree);
-    FEValuesExtractors::Scalar extractor(0);
+    FE_Q<dim>                        fe(degree);
+    const FEValuesExtractors::Scalar extractor(0);
     test_extractor<NumberType, dim>(fe, extractor);
   }
   deallog.pop();
 
   deallog.push("Vector");
   {
-    FESystem<dim>              fe(FE_Q<dim>(degree), dim);
-    FEValuesExtractors::Vector extractor(0);
+    FESystem<dim>                    fe(FE_Q<dim>(degree), dim);
+    const FEValuesExtractors::Vector extractor(0);
     test_extractor<NumberType, dim>(fe, extractor);
   }
   deallog.pop();
 
   deallog.push("SymmetricTensor");
   {
-    FESystem<dim>                          fe(FE_Q<dim>(degree),
+    FESystem<dim>                                fe(FE_Q<dim>(degree),
                      SymmetricTensor<2, dim>::n_independent_components);
-    FEValuesExtractors::SymmetricTensor<2> extractor(0);
+    const FEValuesExtractors::SymmetricTensor<2> extractor(0);
     test_extractor<NumberType, dim>(fe, extractor);
   }
   deallog.pop();
 
   deallog.push("Tensor");
   {
-    FESystem<dim>                 fe(FE_Q<dim>(degree),
+    FESystem<dim>                       fe(FE_Q<dim>(degree),
                      Tensor<2, dim>::n_independent_components);
-    FEValuesExtractors::Tensor<2> extractor(0);
+    const FEValuesExtractors::Tensor<2> extractor(0);
     test_extractor<NumberType, dim>(fe, extractor);
   }
   deallog.pop();
index 38e69cd993ddceaed9a06dc135a194b72b16d3f4..119d86e0d85e1e18ef440ac4a7c8ab647fc4d32f 100644 (file)
@@ -42,7 +42,7 @@ test(const Triangulation<dim> &tr, const FiniteElement<dim> &fe)
   FEValues<dim>     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
index 96e38c5300152546d1187415bc09facc292d6883..bdfb7894cf1bdeca75e3cfb96838c220beb4b8f6 100644 (file)
@@ -42,7 +42,7 @@ test(const Triangulation<dim> &tr, const FiniteElement<dim> &fe)
   FEValues<dim>     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
index cbdb24ee3fb9d60be150b0d924e3b721b531f50c..7acdf1886e84d86a687ff9786f1c6768b0b89a09 100644 (file)
@@ -42,7 +42,7 @@ test(const Triangulation<dim> &tr, const FiniteElement<dim> &fe)
   FEValues<dim>     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
index 1d9b3c27c7673a7f669909352e9a9c5fc3f3ac11..f3c1deba630481db0be759ab8df34635df26c3d9 100644 (file)
@@ -42,7 +42,7 @@ test(const Triangulation<dim> &tr, const FiniteElement<dim> &fe)
   FEValues<dim>     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
index b54bac78928e45acb6cafa893b56ab15d593aa71..8e77a5fc4681985ca15c4253269df7b96a2d04cf 100644 (file)
@@ -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);
index bdc4f2bc8de148385e439d7e5bf46207b8d7c72b..131484121d090dbc3484254f9809d7e62181ab4e 100644 (file)
@@ -116,7 +116,7 @@ main()
           << "n_coupling_dofs = " << fcv.n_coupling_dofs() << std::endl;
 
   FullMatrix<double> 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);
index 2342960c9c4bd09fa66cc66839eaa941dd33ad0f..bf88ed12234af1af846a17077b9816760dc639e3 100644 (file)
@@ -96,7 +96,7 @@ main()
                                         QuadratureCouplingType::reorder);
 
   FullMatrix<double> 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);
index b83e4289efb86892e05fb3c7e701ea80ddbc1cd1..0b29df7ebcecae392e49e18c6a475df04c6dd4f1 100644 (file)
@@ -96,7 +96,7 @@ main()
                                         QuadratureCouplingType::matching);
 
   FullMatrix<double> 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);
index b3b9bb41f26254f718e9f8bc345906a3ba484b82..374c5a2e190af698f5687b035f87eb6686255dd3 100644 (file)
@@ -103,8 +103,8 @@ template <int dim = 2>
 void
 test()
 {
-  FESystem<dim>              fe(FE_DGP<dim>(4), 1);
-  FEValuesExtractors::Scalar extractor(0);
+  FESystem<dim>                    fe(FE_DGP<dim>(4), 1);
+  const FEValuesExtractors::Scalar extractor(0);
 
   Triangulation<dim> tria;
   DoFHandler<dim>    dof_handler(tria);
index f952d58d48c00253d30deb8d6553de167b6f21f7..67eec235658be390c3c5b77618975fb27061a8fe 100644 (file)
@@ -68,7 +68,7 @@ test(const Triangulation<dim> &tr, const FiniteElement<dim> &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)
                 {
index c8e6168058a522b7a42bb2b0c50dde87caaf5dd7..c6964ca860556095fb908c2b34c1a7733971f09d 100644 (file)
@@ -78,7 +78,7 @@ test(const Triangulation<dim> &tr, const FiniteElement<dim> &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)
                     {
index f21543189902cd37a6df8a50357249fe1afff57b..ebb48135adc37a282c1ccf13bd6de3639f6edf3d 100644 (file)
@@ -353,8 +353,8 @@ template <int dim = 2>
 void
 test_jump_function()
 {
-  FE_DGP<dim>                fe(0);
-  FEValuesExtractors::Scalar extractor(0);
+  FE_DGP<dim>                      fe(0);
+  const FEValuesExtractors::Scalar extractor(0);
 
   Triangulation<dim> tria;
   DoFHandler<dim>    dof_handler(tria);
@@ -421,16 +421,16 @@ test()
 
   deallog.push("FEInterfaceViews Scalar");
   {
-    FE_Q<dim>                  fe(degree);
-    FEValuesExtractors::Scalar extractor(0);
+    FE_Q<dim>                        fe(degree);
+    const FEValuesExtractors::Scalar extractor(0);
     test_feiv_extractor<NumberType, dim>(fe, extractor);
   }
   deallog.pop();
 
   deallog.push("FEInterfaceViews Vector");
   {
-    FESystem<dim>              fe(FE_Q<dim>(degree), dim);
-    FEValuesExtractors::Vector extractor(0);
+    FESystem<dim>                    fe(FE_Q<dim>(degree), dim);
+    const FEValuesExtractors::Vector extractor(0);
     test_feiv_extractor<NumberType, dim>(fe, extractor);
   }
   deallog.pop();
index ed02bf9b5d16dd29db6e1ad9578da63168425dd8..1f02218af0bd6afb0906dec8f8ff385b3f46bad2 100644 (file)
@@ -269,16 +269,16 @@ test()
 
   deallog.push("FEInterfaceViews Scalar");
   {
-    FE_Q<dim>                  fe(degree);
-    FEValuesExtractors::Scalar extractor(0);
+    FE_Q<dim>                        fe(degree);
+    const FEValuesExtractors::Scalar extractor(0);
     test_feiv_extractor<NumberType, dim>(fe, extractor);
   }
   deallog.pop();
 
   deallog.push("FEInterfaceViews Vector");
   {
-    FESystem<dim>              fe(FE_Q<dim>(degree), dim);
-    FEValuesExtractors::Vector extractor(0);
+    FESystem<dim>                    fe(FE_Q<dim>(degree), dim);
+    const FEValuesExtractors::Vector extractor(0);
     test_feiv_extractor<NumberType, dim>(fe, extractor);
   }
   deallog.pop();
index a9b3d1c5bc8abf221eb8b82a51fb580a5d694f85..735990e7f235616c0c4c330a988cf5f2dc0917d4 100644 (file)
@@ -53,8 +53,8 @@ test()
   DoFHandler<dim> dof_handler(triangulation);
   dof_handler.distribute_dofs(fe);
 
-  FEValues<dim>              fe_values(fe, QGauss<dim>(2), update_values);
-  FEValuesExtractors::Scalar nothing(1);
+  FEValues<dim>                    fe_values(fe, QGauss<dim>(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)
index 2e80bcf33bbcc2856d6e338d3ef3126e2c10e43e..19a105b121df35a8e21d7fbf3ed8308188c68ead 100644 (file)
@@ -60,9 +60,9 @@ test()
   DoFHandler<dim> dof_handler(triangulation);
   dof_handler.distribute_dofs(fe);
 
-  QGauss<dim - 1>            q(2);
-  FEFaceValues<dim>          fe_values(fe, q, update_values);
-  FEValuesExtractors::Scalar nothing(1);
+  QGauss<dim - 1>                  q(2);
+  FEFaceValues<dim>                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
index 859a5dd2d58576afddd015cb9f296a975fa6a02b..8579c8f9ed30895fb376f4db8d3ac823dc0f3767 100644 (file)
@@ -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<dim>::refine_mesh()
   {
     Vector<float> estimated_error_per_cell(triangulation.n_active_cells());
-    FEValuesExtractors::Scalar pressure(dim);
+    const FEValuesExtractors::Scalar pressure(dim);
     KellyErrorEstimator<dim>::estimate(
       dof_handler,
       QGauss<dim - 1>(degree + 1),
index 9ff1bb93ad3975e90a5c122bfe64154bdb9b8503..fed1f11c134cbd8991e08f1f5c8cec1001d391af 100644 (file)
@@ -100,7 +100,7 @@ test(const unsigned int degree)
   std::vector<Tensor<1, dim>> function_values(unit_points.size());
   std::vector<Tensor<2, dim>> function_gradients(unit_points.size());
 
-  FEValuesExtractors::Vector extractor(0);
+  const FEValuesExtractors::Vector extractor(0);
 
   for (const auto &cell : dof_handler.active_cell_iterators())
     {
index 76ab4f9ab506211375058d6b93a75c38810db45e..bfc0f28eb6910fa294541a247a6d09d2cea80fa7 100644 (file)
@@ -110,8 +110,8 @@ test(const unsigned int degree)
   std::vector<double>         function_values_scalar(unit_points.size());
   std::vector<Tensor<1, dim>> 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())
     {
index a15033002fede270dd718561d97b2f95ad6ce0ce..b53e8deb350e34ae60c8169ce2d79b361f484ba8 100644 (file)
@@ -107,7 +107,7 @@ test(const unsigned int degree)
   std::vector<Tensor<1, dim>> function_values(unit_points.size());
   std::vector<Tensor<2, dim>> function_gradients(unit_points.size());
 
-  FEValuesExtractors::Vector extractor(0);
+  const FEValuesExtractors::Vector extractor(0);
 
   for (const auto &cell : dof_handler.active_cell_iterators())
     {
index 39be324fb6c0d7fa9fb0eb88d520f1971259f951..eb32b9496b6fd310c4d244cc2bf34815ca91bd6e 100644 (file)
@@ -117,8 +117,8 @@ test(const unsigned int degree)
   std::vector<double>         function_values_scalar(unit_points.size());
   std::vector<Tensor<1, dim>> 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())
     {
index 56d2ce958544c327772943708ed83cac5cfe45bc..13965145504ebf80f3ede5adb91e29bee5960a1e 100644 (file)
@@ -98,7 +98,7 @@ test()
   std::vector<Tensor<1, dim>> function_values(unit_points.size());
   std::vector<Tensor<2, dim>> function_gradients(unit_points.size());
 
-  FEValuesExtractors::Vector extractor(0);
+  const FEValuesExtractors::Vector extractor(0);
 
   for (const auto &cell : dof_handler.active_cell_iterators())
     {
index ca0baf459c116c67c39431b1a75266e131ed690e..bb5d4a3ad05f39d2f7fb341e5822a1e99e6514c7 100644 (file)
@@ -102,7 +102,7 @@ test(const unsigned int degree)
   std::vector<double>         function_values(unit_points.size());
   std::vector<Tensor<1, dim>> function_gradients(unit_points.size());
 
-  FEValuesExtractors::Scalar extractor(1);
+  const FEValuesExtractors::Scalar extractor(1);
 
   for (const auto &cell : dof_handler.active_cell_iterators())
     {
index 1f384f9f989457f0ac46823b390354574ea25fa5..214f172ce57eec9cae7fc17ececa721fee067749 100644 (file)
@@ -97,7 +97,7 @@ test()
   std::vector<double> solution_values(fe.dofs_per_cell);
   std::vector<double> solution_values2(fe.dofs_per_cell);
 
-  FEValuesExtractors::Vector extractor(0);
+  const FEValuesExtractors::Vector extractor(0);
 
   for (const auto &cell : dof_handler.active_cell_iterators())
     {
index 102f8a8f1b098031424db7100809d6df4c618a86..bb704d0df97262af403adc850a5a3b60481f6789 100644 (file)
@@ -98,7 +98,7 @@ test(const unsigned int degree)
   std::vector<double> solution_values(fe.dofs_per_cell);
   std::vector<double> solution_values2(fe.dofs_per_cell);
 
-  FEValuesExtractors::Vector extractor(0);
+  const FEValuesExtractors::Vector extractor(0);
 
   for (const auto &cell : dof_handler.active_cell_iterators())
     {
index ce83f663f56cd44933ab387b8a0a46ef2d8e16da..2fbeb83daadd7a70b73983c3d8d6fb00c3ac9b85 100644 (file)
@@ -99,7 +99,7 @@ test()
   std::vector<double> solution_values(fe.dofs_per_cell);
   std::vector<double> solution_values2(fe.dofs_per_cell);
 
-  FEValuesExtractors::Vector extractor(0);
+  const FEValuesExtractors::Vector extractor(0);
 
   for (const auto &cell : dof_handler.active_cell_iterators())
     {
index 0a97d929b310510a45f1b941fe360eb8f30d1d10..2774784641f1fe640473aacc9f1e8b9209bdc649 100644 (file)
@@ -104,7 +104,7 @@ test(const unsigned int degree)
   std::vector<Tensor<1, spacedim>> function_values(unit_points.size());
   std::vector<Tensor<2, spacedim>> function_gradients(unit_points.size());
 
-  FEValuesExtractors::Vector extractor(0);
+  const FEValuesExtractors::Vector extractor(0);
 
   for (const auto &cell : dof_handler.active_cell_iterators())
     {
index a94a216530dcdc20f2ca157affe1697f761e2671..a81d634ca2ddf8adf7099002d844edd2d53d60a1 100644 (file)
@@ -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;
index 964dfb83a245e1d1e02a30311f16aad3fe5154a3..041874719d5a0895a3cf5f353c84bb4ca083673f 100644 (file)
@@ -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,
index 92c22e74f91567c3c7d2241ea7ce9a925aa50703..83e26cbbe9bd5c829c32d96fcae58a215080db9a 100644 (file)
@@ -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,
index 73cc30a3054c34d25fe465fa51f146623c0b32ad..7c9a7299e76a281743eb3668d5bebc8b94c494f5 100644 (file)
@@ -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,
index c7d75df009562e3d1e6716f46e05ebd67ca9dff5..52f5c3f257601af4a08245d9d790b2ed71fdc808 100644 (file)
@@ -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);
index 2b29137cb4b5c1b80770729459f18b5814151ad0..a2779ae528c778892cc799be2452e1251880fedc 100644 (file)
@@ -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);
index f215a1912de181a00561c4bf1187b7b2222d2451..af3bfcfbb9e37fe69b604d15a68eefe0cadecfec 100644 (file)
@@ -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);
index 18180f8ee4311d1cbb26755a9daaa3a33511b84c..22407719684c2609e363b0afb0f6d41b0b48fea1 100644 (file)
@@ -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);
index 6e982b16472c87e7e996aa46e1bb09ddbcad5b9f..b17c6c9d05ff68b933b0d80f66fc6bccfc30d7f8 100644 (file)
@@ -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);
index d80eaefb7e7a27b18a6e603ad85111d0a811967a..74c4943ce94348a4457d12d99e4f7dc831c67dd1 100644 (file)
@@ -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);
index ff8e224ef3b0e24f6146a610f9127ffae5d51f36..a793238495b41ecc574e323d11da609abd9c0882 100644 (file)
@@ -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);
index a2f35e9d6e5d2f9e2c7d7ac177129080bc79c7be..2ac34b15b7755e211bb01a6081f0b269b3e06eb9 100644 (file)
@@ -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);
index ba022f06356f6a66a7b96cbbfb54631c2e901369..be0ccbbad7a552e9c9f50204f41d5e450535c956 100644 (file)
@@ -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);
index 9eec1a650ef4cc0257f39d17db068e100a5ce6ce..62a4cbfae2c56931b1b23aac654f809646a04742 100644 (file)
@@ -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);
index 6f38390abec8532eeb51e114e91b797c27a93327..c6fe8b0882c99aa6f4b726547c7c655c7259a491 100644 (file)
@@ -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);
index 7f3668c8f13c6502add5660c78b5f9231e22c5bd..212c6957c82bb59e9a57395877d24f7c7ed30fc1 100644 (file)
@@ -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);
index 0813130359b4b4b24ede65d75027e6f7f9662f01..fc40d1d080b103eae6503896ae86a168c6b79fa5 100644 (file)
@@ -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<float> estimated_error_per_cell(triangulation.n_active_cells());
 
-    FEValuesExtractors::Scalar pressure(dim);
+    const FEValuesExtractors::Scalar pressure(dim);
     KellyErrorEstimator<dim>::estimate(
       dof_handler,
       QGauss<dim - 1>(degree + 1),
index 036962ebd9b29e824931e4047936e80dbc3cf6bd..c50a6ad6654bb60dbc12db9fd730de12fbcca549 100644 (file)
@@ -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<float> estimated_error_per_cell(triangulation.n_active_cells());
 
-    FEValuesExtractors::Scalar pressure(dim);
+    const FEValuesExtractors::Scalar pressure(dim);
     KellyErrorEstimator<dim>::estimate(
       dof_handler,
       QGauss<dim - 1>(degree + 1),
index 3a7e0b941539d6c0572de12d91ea10235a901a6c..59d2f4b6e735aee5327ff1ea6890d28a1d338c61 100644 (file)
@@ -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<float> estimated_error_per_cell(triangulation.n_active_cells());
 
-    FEValuesExtractors::Scalar pressure(dim);
+    const FEValuesExtractors::Scalar pressure(dim);
     KellyErrorEstimator<dim>::estimate(
       dof_handler,
       QGauss<dim - 1>(degree + 1),
index 32e15a81c3c1b1b61c7cbb08e6728ac6b295419d..ba4e97f9639d427cc35a2fd054fcb768d40dee50 100644 (file)
@@ -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(
index ec5cbbf8023c4249c3ee6ac98ec5b7a489be02eb..26607e98ade61a0142ccae8e247fc597df0bd714 100644 (file)
@@ -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<Tensor<2, dim>> qp_Grad_u_t;
-  std::vector<Tensor<2, dim>> qp_Grad_u_t1;
-  std::vector<Tensor<2, dim>> qp_dot_Grad_u_t;
-  std::vector<Tensor<2, dim>> qp_dot_grad_u_t;
+  const FEValuesExtractors::Vector u_fe(0);
+  std::vector<Tensor<2, dim>>      qp_Grad_u_t;
+  std::vector<Tensor<2, dim>>      qp_Grad_u_t1;
+  std::vector<Tensor<2, dim>>      qp_dot_Grad_u_t;
+  std::vector<Tensor<2, dim>>      qp_dot_grad_u_t;
 
   FEValues<dim>         fe_values(fe, qf, update_gradients);
   MappingQEulerian<dim> q1_mapping(1, dof_handler, soln_t);
index df90da71314d4586c8db47b8e2be84b90c5b3e7c..9a512d2f6bad7c81a836514a1aa6382f083e3ff5 100644 (file)
@@ -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<types::global_dof_index, double> boundary_values;
     VectorTools::interpolate_boundary_values(dof_handler,
                                              0,
index 67cae3e4cddf5154c52200e96eec1a604bfe07f3..a3b1d9871e9851addb369439c801e2a12037ea37 100644 (file)
@@ -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<types::global_dof_index, double> boundary_values;
     VectorTools::interpolate_boundary_values(dof_handler,
                                              0,
index 1e854cbe2a174341b93b50a864996b645abc5e85..0b99edef37cd1722d19d6811d4b2753bec8e0b18 100644 (file)
@@ -638,7 +638,7 @@ namespace Step18
 #endif
 
 
-    FEValuesExtractors::Scalar                z_component(dim - 1);
+    const FEValuesExtractors::Scalar          z_component(dim - 1);
     std::map<types::global_dof_index, double> boundary_values;
     VectorTools::interpolate_boundary_values(mapping,
                                              dof_handler,
index 531eaa59eb142fd83fcf59dd418ef4eb82b2ee3f..c3dfdd47bc61b89ad731a7515329d8526a05c08e 100644 (file)
@@ -600,8 +600,8 @@ namespace Step31
     deallog << "   Rebuilding Stokes preconditioner..." << std::flush;
     assemble_stokes_preconditioner();
     Amg_preconditioner = std::make_shared<TrilinosWrappers::PreconditionAMG>();
-    std::vector<std::vector<bool>> constant_modes;
-    FEValuesExtractors::Vector     velocity_components(0);
+    std::vector<std::vector<bool>>   constant_modes;
+    const FEValuesExtractors::Vector velocity_components(0);
     DoFTools::extract_constant_modes(stokes_dof_handler,
                                      stokes_fe.component_mask(
                                        velocity_components),
index f873f3dae0f87d254d1b330746211646c1db5fe8..26cc9c9c325d9a4eeb5537c9caf025a4c5346d70 100644 (file)
@@ -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,

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.