]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Fix Affineconstraints deprecations by initializing with `(locally_owned_dofs, locally...
authorMarc Fehling <mafehling.git@gmail.com>
Fri, 4 Jul 2025 14:33:10 +0000 (16:33 +0200)
committerMarc Fehling <mafehling.git@gmail.com>
Mon, 7 Jul 2025 18:09:12 +0000 (20:09 +0200)
197 files changed:
examples/step-90/step-90.cc
tests/arpack/parpack_advection_diffusion_petsc.cc
tests/arpack/parpack_advection_diffusion_trilinos.cc
tests/arpack/step-36_parpack.cc
tests/arpack/step-36_parpack_mf.cc
tests/arpack/step-36_parpack_mf_02.cc
tests/arpack/step-36_parpack_mf_03.cc
tests/arpack/step-36_parpack_trilinos.cc
tests/dofs/dof_tools_04a.cc
tests/dofs/dof_tools_extract_dofs_with_support_contained_within_02.cc
tests/dofs/dof_tools_extract_dofs_with_support_contained_within_04.cc
tests/fe/fe_enriched_color_07.cc
tests/fe/fe_enriched_step-36.cc
tests/fe/fe_enriched_step-36b.cc
tests/gla/block_mat_02.cc
tests/gla/block_mat_03.cc
tests/hp/laplace.h
tests/lac/diagonal_matrix_02.cc
tests/lac/linear_operator_12.cc
tests/lac/linear_operator_12a.cc
tests/lac/linear_operator_14.cc
tests/lac/step-40-linear_operator_01.cc
tests/lac/step-40-linear_operator_02.cc
tests/lac/step-40-linear_operator_03.cc
tests/lac/step-40-linear_operator_04.cc
tests/lac/step-40-linear_operator_05.cc
tests/lac/step-40-linear_operator_06.cc
tests/lac/utilities_02.cc
tests/mappings/mapping_q_eulerian_08.cc
tests/mappings/mapping_q_eulerian_09.cc
tests/mappings/mapping_q_eulerian_14.cc
tests/matrix_free/cell_categorization_02.cc
tests/matrix_free/compute_diagonal_02.cc
tests/matrix_free/dof_info_01.cc
tests/matrix_free/dof_info_02.cc
tests/matrix_free/interpolate_to_mg.cc
tests/matrix_free/laplace_operator_01.cc
tests/matrix_free/laplace_operator_02.cc
tests/matrix_free/laplace_operator_03.cc
tests/matrix_free/laplace_operator_05.cc
tests/matrix_free/mass_operator_01.cc
tests/matrix_free/mass_operator_02.cc
tests/matrix_free/mass_operator_03.cc
tests/matrix_free/mass_operator_04.cc
tests/matrix_free/matrix_vector_10.cc
tests/matrix_free/matrix_vector_11.cc
tests/matrix_free/matrix_vector_12.cc
tests/matrix_free/matrix_vector_13.cc
tests/matrix_free/matrix_vector_19.cc
tests/matrix_free/matrix_vector_22.cc
tests/matrix_free/matrix_vector_23.cc
tests/matrix_free/matrix_vector_26.cc
tests/matrix_free/matrix_vector_27.cc
tests/matrix_free/matrix_vector_blocks.cc
tests/matrix_free/matrix_vector_large_degree_02.cc
tests/matrix_free/parallel_multigrid_02.cc
tests/matrix_free/parallel_multigrid_adaptive_01.cc
tests/matrix_free/parallel_multigrid_adaptive_02.cc
tests/matrix_free/parallel_multigrid_adaptive_03.cc
tests/matrix_free/parallel_multigrid_adaptive_04.cc
tests/matrix_free/parallel_multigrid_adaptive_05.cc
tests/matrix_free/parallel_multigrid_adaptive_06.cc
tests/matrix_free/parallel_multigrid_adaptive_06ref.cc
tests/matrix_free/parallel_multigrid_adaptive_07.cc
tests/matrix_free/parallel_multigrid_adaptive_08.cc
tests/matrix_free/parallel_multigrid_mf_02.cc
tests/matrix_free/pre_and_post_loops_02.cc
tests/matrix_free/solver_cg_interleave.cc
tests/matrix_free/step-37-inhomogeneous-1.cc
tests/matrix_free/step-48.cc
tests/matrix_free/step-48c.cc
tests/matrix_free/stokes_computation.cc
tests/matrix_free_kokkos/matrix_free_device_initialize_vector.cc
tests/matrix_free_kokkos/matrix_free_device_matrix_vector_10.cc
tests/matrix_free_kokkos/matrix_free_device_matrix_vector_10a.cc
tests/matrix_free_kokkos/matrix_free_device_matrix_vector_19.cc
tests/matrix_free_kokkos/matrix_free_device_matrix_vector_25.cc
tests/matrix_free_kokkos/matrix_free_device_matrix_vector_26.cc
tests/matrix_free_kokkos/matrix_free_device_precondition_chebyshev.cc
tests/matrix_free_kokkos/matrix_free_device_reinit_01.cc
tests/meshworker/step-50-mesh_loop.cc
tests/mpi/constraint_matrix_condense_01.cc
tests/mpi/constraint_matrix_trilinos_bug.cc
tests/mpi/constraints_consistent_01.cc
tests/mpi/data_out_faces_01.cc
tests/mpi/hp_constraints_consistent_01.cc
tests/mpi/hp_constraints_consistent_02.cc
tests/mpi/hp_constraints_consistent_03.cc
tests/mpi/hp_hanging_node_constraints_01.cc
tests/mpi/hp_hanging_node_constraints_02.cc
tests/mpi/hp_step-4.cc
tests/mpi/hp_step-40.cc
tests/mpi/hp_step-40_variable_01.cc
tests/mpi/interpolate_to_different_mesh_02.cc
tests/mpi/make_zero_boundary_values.cc
tests/mpi/no_flux_constraints.cc
tests/mpi/no_flux_constraints_02.cc
tests/mpi/no_flux_constraints_03.cc
tests/mpi/normal_flux_constraints.cc
tests/mpi/p4est_2d_constraintmatrix_01.cc
tests/mpi/p4est_2d_constraintmatrix_02.cc
tests/mpi/p4est_2d_constraintmatrix_03.cc
tests/mpi/p4est_2d_constraintmatrix_04.cc
tests/mpi/p4est_2d_constraintmatrix_05.cc
tests/mpi/p4est_3d_constraintmatrix_01.cc
tests/mpi/p4est_3d_constraintmatrix_02.cc
tests/mpi/p4est_3d_constraintmatrix_03.cc
tests/mpi/p4est_3d_constraintmatrix_04.cc
tests/mpi/parallel_block_vector_05.cc
tests/mpi/parallel_block_vector_06.cc
tests/mpi/parallel_block_vector_07.cc
tests/mpi/parallel_block_vector_08.cc
tests/mpi/parallel_block_vector_09.cc
tests/mpi/parallel_block_vector_10.cc
tests/mpi/parallel_block_vector_11.cc
tests/mpi/parallel_block_vector_12.cc
tests/mpi/periodicity_01.cc
tests/mpi/periodicity_02.cc
tests/mpi/periodicity_03.cc
tests/mpi/periodicity_04.cc
tests/mpi/periodicity_06.cc
tests/mpi/periodicity_07.cc
tests/mpi/periodicity_08.cc
tests/mpi/petsc_distribute_01.cc
tests/mpi/petsc_distribute_01_inhomogeneous.cc
tests/mpi/petsc_step-27.cc
tests/mpi/project_bv_div_conf.cc
tests/mpi/step-37.cc
tests/mpi/step-40.cc
tests/mpi/step-40_cuthill_mckee.cc
tests/mpi/step-40_cuthill_mckee_MPI-subset.cc
tests/mpi/step-40_direct_solver.cc
tests/mpi/trilinos_distribute_01.cc
tests/mpi/trilinos_distribute_01_inhomogeneous.cc
tests/mpi/trilinos_distribute_04.cc
tests/mpi/trilinos_step-27.cc
tests/multigrid-global-coarsening/events_non_nested.cc
tests/multigrid-global-coarsening/interpolate_to_mg.cc
tests/multigrid-global-coarsening/interpolate_to_mg_01.cc
tests/multigrid-global-coarsening/multigrid_a_01.cc
tests/multigrid-global-coarsening/multigrid_a_02.cc
tests/multigrid-global-coarsening/multigrid_a_03.cc
tests/multigrid-global-coarsening/multigrid_p_01.cc
tests/multigrid-global-coarsening/multigrid_p_02.cc
tests/multigrid-global-coarsening/multigrid_p_03.cc
tests/multigrid-global-coarsening/multigrid_p_04.cc
tests/multigrid-global-coarsening/non_nested_multigrid_01.cc
tests/multigrid-global-coarsening/non_nested_multigrid_02.cc
tests/multigrid-global-coarsening/non_nested_multigrid_03.cc
tests/multigrid-global-coarsening/non_nested_multigrid_04.cc
tests/multigrid-global-coarsening/non_nested_multigrid_05.cc
tests/multigrid-global-coarsening/parallel_multigrid_02.cc
tests/multigrid-global-coarsening/parallel_multigrid_adaptive_01.cc
tests/multigrid-global-coarsening/parallel_multigrid_adaptive_02.cc
tests/multigrid-global-coarsening/parallel_multigrid_adaptive_03.cc
tests/multigrid-global-coarsening/parallel_multigrid_adaptive_04.cc
tests/multigrid-global-coarsening/parallel_multigrid_adaptive_05.cc
tests/multigrid-global-coarsening/parallel_multigrid_adaptive_06.cc
tests/multigrid-global-coarsening/parallel_multigrid_adaptive_06ref.cc
tests/multigrid-global-coarsening/parallel_multigrid_adaptive_07.cc
tests/multigrid-global-coarsening/parallel_multigrid_adaptive_08.cc
tests/multigrid-global-coarsening/parallel_multigrid_mf_02.cc
tests/multigrid-global-coarsening/step-37-inhomogeneous-1.cc
tests/multigrid-global-coarsening/stokes_computation.cc
tests/multigrid/constrained_dofs_01.cc
tests/multigrid/constrained_dofs_07.cc
tests/multigrid/events_01.cc
tests/multigrid/interface_matrix_entry_01.cc
tests/multigrid/interpolate_to_mg_01.cc
tests/multigrid/mg_coarse_01.cc
tests/multigrid/step-16-50-mpi-linear-operator.cc
tests/multigrid/step-16-50-mpi-smoother.cc
tests/multigrid/step-16-50-mpi.cc
tests/multigrid/step-50_01.cc
tests/multigrid/step-50_02.cc
tests/multigrid/transfer_04.cc
tests/multigrid/transfer_matrix_free_13.cc
tests/multigrid/transfer_prebuilt_05.cc
tests/non_matching/step-70.cc
tests/numerics/no_flux_18.cc
tests/numerics/nonlinear_solver_selector_03.cc
tests/numerics/project_parallel_common.h
tests/numerics/project_parallel_qp_common.h
tests/numerics/project_parallel_qpmf_common.h
tests/performance/timing_mg_glob_coarsen.cc
tests/performance/timing_step_37.cc
tests/performance/timing_step_37_gc.cc
tests/petsc/step-77-snes.cc
tests/petsc_complex/fe_get_function_values.cc
tests/petsc_complex/parallel_sparse_matrix_01.cc
tests/sharedtria/pbc_01.cc
tests/simplex/step-40.cc
tests/simplex/step-55.cc
tests/slepc/step-36_parallel.cc
tests/slepc/step-36_parallel_02.cc
tests/slepc/step-36_parallel_03.cc
tests/trilinos/solver_control_06.cc

index 4cc16e141e6ef9406d549472a165401a7eff3ada..8e156f93b77475514f2d1cb2c903ce39114056d4 100644 (file)
@@ -666,7 +666,8 @@ namespace Step90
     level_set_constraints.clear();
     const IndexSet level_set_locally_relevant_dofs =
       DoFTools::extract_locally_relevant_dofs(level_set_dof_handler);
-    level_set_constraints.reinit(level_set_locally_relevant_dofs);
+    level_set_constraints.reinit(level_set_dof_handler.locally_owned_dofs(),
+                                 level_set_locally_relevant_dofs);
     DoFTools::make_hanging_node_constraints(level_set_dof_handler,
                                             level_set_constraints);
     level_set_constraints.close();
index 2abfdb93cf24555ec04435300f70dda02547ee73..5720c1bb5130acea929bb0eca532cb165ffa8130 100644 (file)
@@ -190,7 +190,7 @@ test()
     dealii::DoFTools::extract_locally_relevant_dofs(dof_handler);
 
   constraints.clear();
-  constraints.reinit(locally_relevant_dofs);
+  constraints.reinit(locally_owned_dofs, locally_relevant_dofs);
   dealii::DoFTools::make_hanging_node_constraints(dof_handler, constraints);
   dealii::VectorTools::interpolate_boundary_values(
     dof_handler, 0, dealii::Functions::ZeroFunction<dim>(), constraints);
index ae4b9af7beee3649842514e8ccc72f674d5bf262..2ed14a4e8e59d2a1f3ad1e957e324a3278ab1567 100644 (file)
@@ -164,7 +164,7 @@ test()
   locally_relevant_dofs = DoFTools::extract_locally_relevant_dofs(dof_handler);
 
   constraints.clear();
-  constraints.reinit(locally_relevant_dofs);
+  constraints.reinit(locally_owned_dofs, locally_relevant_dofs);
   DoFTools::make_hanging_node_constraints(dof_handler, constraints);
   VectorTools::interpolate_boundary_values(dof_handler,
                                            0,
index cf8ca79f3901f7e3b68c905981d5ea7845fb7173..f642db32ee1071b5f52adb6e412a081d260ac34e 100644 (file)
@@ -192,7 +192,7 @@ test()
     dealii::DoFTools::extract_locally_relevant_dofs(dof_handler);
 
   constraints.clear();
-  constraints.reinit(locally_relevant_dofs);
+  constraints.reinit(locally_owned_dofs, locally_relevant_dofs);
   dealii::DoFTools::make_hanging_node_constraints(dof_handler, constraints);
   dealii::VectorTools::interpolate_boundary_values(
     dof_handler, 0, dealii::Functions::ZeroFunction<dim>(), constraints);
index ee2fda02532428c0638abbb6d5924fdcb5871c13..1dd7270dbc22bec2bec9bac21cfcb0c497958984 100644 (file)
@@ -85,7 +85,7 @@ test()
   const IndexSet locally_relevant_dofs =
     DoFTools::extract_locally_relevant_dofs(dof_handler);
   AffineConstraints<double> constraints;
-  constraints.reinit(locally_relevant_dofs);
+  constraints.reinit(dof_handler.locally_owned_dofs(), locally_relevant_dofs);
   DoFTools::make_hanging_node_constraints(dof_handler, constraints);
   VectorTools::interpolate_boundary_values(dof_handler,
                                            0,
index 6d82b81f974bf0fa0fca1fea1f9f7b9be05720b3..8ed61697839225b79f36175850104914a348d721 100644 (file)
@@ -85,7 +85,7 @@ test()
   const IndexSet locally_relevant_dofs =
     DoFTools::extract_locally_relevant_dofs(dof_handler);
   AffineConstraints<double> constraints;
-  constraints.reinit(locally_relevant_dofs);
+  constraints.reinit(dof_handler.locally_owned_dofs(), locally_relevant_dofs);
   DoFTools::make_hanging_node_constraints(dof_handler, constraints);
   VectorTools::interpolate_boundary_values(dof_handler,
                                            0,
index 86411478d2fa39c9b746140dfa55d6b6913f6552..d6a824b0d9f348c15ca8273a62c5d45d8f393019 100644 (file)
@@ -86,7 +86,7 @@ test()
   const IndexSet locally_relevant_dofs =
     DoFTools::extract_locally_relevant_dofs(dof_handler);
   AffineConstraints<double> constraints;
-  constraints.reinit(locally_relevant_dofs);
+  constraints.reinit(dof_handler.locally_owned_dofs(), locally_relevant_dofs);
   DoFTools::make_hanging_node_constraints(dof_handler, constraints);
   VectorTools::interpolate_boundary_values(dof_handler,
                                            0,
index 25d67a9545c0bdd1eb660ccdd1ad5ef05c81252f..8fd75373e91e121e243971da26f5eb19ceb66c67 100644 (file)
@@ -166,7 +166,7 @@ test()
   locally_relevant_dofs = DoFTools::extract_locally_relevant_dofs(dof_handler);
 
   constraints.clear();
-  constraints.reinit(locally_relevant_dofs);
+  constraints.reinit(locally_owned_dofs, locally_relevant_dofs);
   DoFTools::make_hanging_node_constraints(dof_handler, constraints);
   VectorTools::interpolate_boundary_values(dof_handler,
                                            0,
index f98b51ee9523ba9a9eefe08deb481c4e74a57157..8bb1d7305d7dd4154d08a8423efce8a0681fc19f 100644 (file)
@@ -55,7 +55,8 @@ check_this(const DoFHandler<dim> &dof_handler)
   const IndexSet is_hanging_node_constrained =
     DoFTools::extract_hanging_node_dofs(dof_handler);
 
-  AffineConstraints<double> constraints(locally_relevant_dofs);
+  AffineConstraints<double> constraints(dof_handler.locally_owned_dofs(),
+                                        locally_relevant_dofs);
   DoFTools::make_hanging_node_constraints(dof_handler, constraints);
   constraints.close();
 
index dce5f7e4bc650737f3452543f786f0bc4729ee30..f631c7aabe43ea07c612f1a78bde77c732ddfef8 100644 (file)
@@ -98,7 +98,7 @@ test(const unsigned int flag)
     DoFTools::extract_locally_relevant_dofs(dh);
 
   AffineConstraints<double> cm;
-  cm.reinit(locally_relevant_set);
+  cm.reinit(dh.locally_owned_dofs(), locally_relevant_set);
   DoFTools::make_hanging_node_constraints(dh, cm);
   cm.close();
 
index 68a825861337ab7db468272754b25f0f50ec0e23..d2431a41d56b74323aaadf4274d367542318b8f6 100644 (file)
@@ -101,7 +101,7 @@ test()
     DoFTools::extract_locally_relevant_dofs(dh);
 
   AffineConstraints<double> cm;
-  cm.reinit(locally_relevant_set);
+  cm.reinit(locally_owned_set, locally_relevant_set);
   DoFTools::make_hanging_node_constraints(dh, cm);
   cm.close();
 
index 6d2b36fa060a0f979136b0bd24cfb029dac17f61..fe4fa640ddd28724fdaac9ea5b17e01ca7850446 100644 (file)
@@ -1442,7 +1442,7 @@ LaplaceProblem<dim>::setup_system()
   locally_relevant_dofs = DoFTools::extract_locally_relevant_dofs(dof_handler);
 
   constraints.clear();
-  constraints.reinit(locally_relevant_dofs);
+  constraints.reinit(locally_owned_dofs, locally_relevant_dofs);
   DoFTools::make_hanging_node_constraints(dof_handler, constraints);
 
   SigmaFunction<dim> boundary_value_func;
index d3ee8395332dfb3b3135d87ad576ae91a9db05e0..156b29322218a18d3e5dc8c1a877a2e866de218a 100644 (file)
@@ -290,7 +290,7 @@ namespace Step36
       DoFTools::extract_locally_relevant_dofs(dof_handler);
 
     constraints.clear();
-    constraints.reinit(locally_relevant_dofs);
+    constraints.reinit(locally_owned_dofs, locally_relevant_dofs);
     DoFTools::make_hanging_node_constraints(dof_handler, constraints);
     VectorTools::interpolate_boundary_values(dof_handler,
                                              0,
index 5ccbb381a1acd1bbf7a27d4744fd8ec5be7f85ed..7ae746b6568d64bf9cfd1ac11e36b7218eb8fad3 100644 (file)
@@ -300,7 +300,7 @@ namespace Step36
       DoFTools::extract_locally_relevant_dofs(dof_handler);
 
     constraints.clear();
-    constraints.reinit(locally_relevant_dofs);
+    constraints.reinit(locally_owned_dofs, locally_relevant_dofs);
     DoFTools::make_hanging_node_constraints(dof_handler, constraints);
     VectorTools::interpolate_boundary_values(dof_handler,
                                              0,
index 61e94e486ec23cad1e4e62f07327f25fd71f83c5..e7e6f20a326ce27eb4a133b51996587bb4768854 100644 (file)
@@ -97,7 +97,7 @@ test()
 
   stokes_relevant_set.print(deallog.get_file_stream());
 
-  AffineConstraints<double> cm(stokes_relevant_set);
+  AffineConstraints<double> cm(stokes_index_set, stokes_relevant_set);
 
   DoFTools::make_hanging_node_constraints(stokes_dof_handler, cm);
 
@@ -190,7 +190,7 @@ test_LA_Trilinos()
     stokes_relevant_set.get_view(n_u, n_u + n_p));
   stokes_relevant_set.print(deallog.get_file_stream());
 
-  AffineConstraints<double> cm(stokes_relevant_set);
+  AffineConstraints<double> cm(stokes_index_set, stokes_relevant_set);
 
   DoFTools::make_hanging_node_constraints(stokes_dof_handler, cm);
 
index bc708ba1bf3ebee2e793d08777b78df78f99878a..74c9491e41870af7cf8309e5cad5d539da155c28 100644 (file)
@@ -103,7 +103,8 @@ test()
 
   using number = typename LA::MPI::BlockSparseMatrix::value_type;
 
-  AffineConstraints<number> constraints(locally_relevant_dofs);
+  AffineConstraints<number> constraints(locally_owned_dofs,
+                                        locally_relevant_dofs);
   constraints.close();
 
   BlockDynamicSparsityPattern bcsp(locally_relevant_partitioning);
@@ -207,7 +208,8 @@ test_alt()
   locally_relevant_dofs.print(deallog);
 
 
-  AffineConstraints<double> constraints(locally_relevant_dofs);
+  AffineConstraints<double> constraints(locally_owned_dofs,
+                                        locally_relevant_dofs);
   constraints.close();
 
   TrilinosWrappers::BlockSparsityPattern sp(locally_owned_partitioning,
index eaa770f399b3d50f7789977aa89ecd366bc1a6c5..f9f635af166ef8813e8434a49e824b970b133c24 100644 (file)
@@ -269,7 +269,7 @@ Laplace<dim>::setup_system()
   system_rhs = 0;
 
   constraints.clear();
-  // constraints.reinit(locally_relevant_dofs);
+  // constraints.reinit(locally_owned_dofs, locally_relevant_dofs);
   DoFTools::make_hanging_node_constraints(dof_handler, constraints);
   VectorTools::interpolate_boundary_values(dof_handler,
                                            0,
index 96e281c5d74e9967d1cb740f05d3683a888cccfa..097e2d418a54b8cc018787b366a30a22f4252f72 100644 (file)
@@ -69,7 +69,7 @@ test(const bool hanging_nodes = true)
   const IndexSet &owned_set    = dof.locally_owned_dofs();
   const IndexSet  relevant_set = DoFTools::extract_locally_relevant_dofs(dof);
 
-  AffineConstraints<double> constraints(relevant_set);
+  AffineConstraints<double> constraints(owned_set, relevant_set);
   DoFTools::make_hanging_node_constraints(dof, constraints);
   VectorTools::interpolate_boundary_values(dof,
                                            0,
index 675bfd291925ad82ab24376facb370b9b6b5f5ae..eeb210e56ea619f43e8aa43272cee742bce0fa46 100644 (file)
@@ -184,19 +184,18 @@ Step4<dim>::setup_system()
 {
   dof_handler.distribute_dofs(fe);
 
+  const IndexSet &locally_owned_dofs = dof_handler.locally_owned_dofs();
+  const IndexSet  locally_relevant_dofs =
+    DoFTools::extract_locally_relevant_dofs(dof_handler);
+
   constraints.clear();
-  constraints.reinit(dof_handler.locally_owned_dofs(),
-                     DoFTools::extract_locally_relevant_dofs(dof_handler));
+  constraints.reinit(locally_owned_dofs, locally_relevant_dofs);
   VectorTools::interpolate_boundary_values(dof_handler,
                                            0,
                                            BoundaryValues<dim>(),
                                            constraints);
   constraints.close();
 
-  const IndexSet &locally_owned_dofs = dof_handler.locally_owned_dofs();
-  const IndexSet  locally_relevant_dofs =
-    DoFTools::extract_locally_relevant_dofs(dof_handler);
-
   DynamicSparsityPattern dsp(locally_relevant_dofs);
   DoFTools::make_sparsity_pattern(dof_handler, dsp, constraints, false);
   SparsityTools::distribute_sparsity_pattern(dsp,
index 1e80b6298a686f3f94dedcbe7cde658e1ddd492a..578cfb9ca4d430292ad57547a944f37f8b8b906c 100644 (file)
@@ -185,19 +185,18 @@ Step4<dim>::setup_system()
 {
   dof_handler.distribute_dofs(fe);
 
+  const IndexSet &locally_owned_dofs = dof_handler.locally_owned_dofs();
+  const IndexSet  locally_relevant_dofs =
+    DoFTools::extract_locally_relevant_dofs(dof_handler);
+
   constraints.clear();
-  constraints.reinit(dof_handler.locally_owned_dofs(),
-                     DoFTools::extract_locally_relevant_dofs(dof_handler));
+  constraints.reinit(locally_owned_dofs, locally_relevant_dofs);
   VectorTools::interpolate_boundary_values(dof_handler,
                                            0,
                                            BoundaryValues<dim>(),
                                            constraints);
   constraints.close();
 
-  const IndexSet &locally_owned_dofs = dof_handler.locally_owned_dofs();
-  const IndexSet  locally_relevant_dofs =
-    DoFTools::extract_locally_relevant_dofs(dof_handler);
-
   DynamicSparsityPattern dsp(locally_relevant_dofs);
   DoFTools::make_sparsity_pattern(dof_handler, dsp, constraints, false);
   SparsityTools::distribute_sparsity_pattern(dsp,
index d7b3382c996c14f870635217e7abdacce9260d74..228d0f355d1b62f1d97eedd3e074c18a2c6cd0ff 100644 (file)
@@ -112,7 +112,7 @@ build_matrix_vector(TrilinosWrappers::BlockSparseMatrix &matrix,
                                    dofs_per_block[0] + dofs_per_block[1]));
 
   constraints.clear();
-  constraints.reinit(locally_relevant_dofs);
+  constraints.reinit(locally_owned_dofs, locally_relevant_dofs);
   constraints.close();
 
   // See
index 7ced2f48459b88e16978e300aa75255c7ceb0b67..c694f76e3b29c847f3d9372763008725761bc8ef 100644 (file)
@@ -132,7 +132,7 @@ namespace Step40
     system_rhs.reinit(locally_owned_dofs, mpi_communicator);
 
     constraints.clear();
-    constraints.reinit(locally_relevant_dofs);
+    constraints.reinit(locally_owned_dofs, locally_relevant_dofs);
     DoFTools::make_hanging_node_constraints(dof_handler, constraints);
     VectorTools::interpolate_boundary_values(dof_handler,
                                              0,
index 7d822280bc34905aed54d7eef3fafe5e447f79c2..a32d351deb3fc2701a33fd7be40de449b7499a3d 100644 (file)
@@ -132,7 +132,7 @@ namespace Step40
     system_rhs.reinit(locally_owned_dofs, mpi_communicator);
 
     constraints.clear();
-    constraints.reinit(locally_relevant_dofs);
+    constraints.reinit(locally_owned_dofs, locally_relevant_dofs);
     DoFTools::make_hanging_node_constraints(dof_handler, constraints);
     VectorTools::interpolate_boundary_values(dof_handler,
                                              0,
index 52da8df4882d2760bd427f0b02c85a5379d5afdf..6c8d731e87d7a5ed30ed16cfa6a040d7afbfac57 100644 (file)
@@ -133,7 +133,7 @@ namespace Step40
     system_rhs.reinit(locally_owned_dofs, mpi_communicator);
 
     constraints.clear();
-    constraints.reinit(locally_relevant_dofs);
+    constraints.reinit(locally_owned_dofs, locally_relevant_dofs);
     DoFTools::make_hanging_node_constraints(dof_handler, constraints);
     VectorTools::interpolate_boundary_values(dof_handler,
                                              0,
index 11edae8ed1cf0fc0cd680dad2d2aca4527279a7a..ff894c70db0b488fc5e9ef70ba1274ca9544ebaf 100644 (file)
@@ -133,7 +133,7 @@ namespace Step40
     system_rhs.reinit(locally_owned_dofs, mpi_communicator);
 
     constraints.clear();
-    constraints.reinit(locally_relevant_dofs);
+    constraints.reinit(locally_owned_dofs, locally_relevant_dofs);
     DoFTools::make_hanging_node_constraints(dof_handler, constraints);
     VectorTools::interpolate_boundary_values(dof_handler,
                                              0,
index c582dcb41c47e6a7521f18d98a0d30e886a49f13..c30fac176966cfc40f0aa6a12af9e298910bd5a8 100644 (file)
@@ -133,7 +133,7 @@ namespace Step40
     system_rhs.reinit(locally_owned_dofs, mpi_communicator);
 
     constraints.clear();
-    constraints.reinit(locally_relevant_dofs);
+    constraints.reinit(locally_owned_dofs, locally_relevant_dofs);
     DoFTools::make_hanging_node_constraints(dof_handler, constraints);
     VectorTools::interpolate_boundary_values(dof_handler,
                                              0,
index 59aaec64020543d4604ac4f8efdb17848b517837..68fc041efabd6ddd92b690e2ed2a4e810d68567c 100644 (file)
@@ -133,7 +133,7 @@ namespace Step40
     system_rhs.reinit(locally_owned_dofs, mpi_communicator);
 
     constraints.clear();
-    constraints.reinit(locally_relevant_dofs);
+    constraints.reinit(locally_owned_dofs, locally_relevant_dofs);
     DoFTools::make_hanging_node_constraints(dof_handler, constraints);
     VectorTools::interpolate_boundary_values(dof_handler,
                                              0,
index 86a8f14105d459fe41359c3d28b59e9e40ab306a..cd06a0df456abd3eecd6ba733335f27b15811fb1 100644 (file)
@@ -96,7 +96,7 @@ test()
   const IndexSet locally_relevant_dofs =
     DoFTools::extract_locally_relevant_dofs(dof_handler);
   AffineConstraints<double> constraints;
-  constraints.reinit(locally_relevant_dofs);
+  constraints.reinit(dof_handler.locally_owned_dofs(), locally_relevant_dofs);
   DoFTools::make_hanging_node_constraints(dof_handler, constraints);
   VectorTools::interpolate_boundary_values(dof_handler,
                                            0,
index 8a9a82f40b47bce99e08cde96d8e1bc41327d6dc..6e59c0d2803bc7b80fd332cb218581f35d48e3be 100644 (file)
@@ -166,12 +166,13 @@ test(const unsigned int n_ref = 0)
 
   // constraints:
   AffineConstraints<double> constraints_euler;
-  constraints_euler.reinit(locally_relevant_dofs_euler);
+  constraints_euler.reinit(locally_owned_dofs_euler,
+                           locally_relevant_dofs_euler);
   DoFTools::make_hanging_node_constraints(dof_handler_euler, constraints_euler);
   constraints_euler.close();
 
   AffineConstraints<double> constraints;
-  constraints.reinit(locally_relevant_dofs);
+  constraints.reinit(locally_owned_dofs, locally_relevant_dofs);
   DoFTools::make_hanging_node_constraints(dof_handler, constraints);
   constraints_euler.close();
 
@@ -292,7 +293,8 @@ test(const unsigned int n_ref = 0)
       AffineConstraints<double> level_constraints;
       const IndexSet            relevant_dofs =
         DoFTools::extract_locally_relevant_level_dofs(dof_handler, level);
-      level_constraints.reinit(relevant_dofs);
+      level_constraints.reinit(dof_handler.locally_owned_mg_dofs(level),
+                               relevant_dofs);
       level_constraints.add_lines(
         mg_constrained_dofs.get_boundary_indices(level));
       level_constraints.close();
index 67c67799cac7233c2847f7815c9405dd13371f45..30e362b601a8003a2e163512fd44c19914b87cfe 100644 (file)
@@ -110,12 +110,13 @@ test()
 
   // constraints:
   AffineConstraints<double> constraints_euler;
-  constraints_euler.reinit(locally_relevant_dofs_euler);
+  constraints_euler.reinit(locally_owned_dofs_euler,
+                           locally_relevant_dofs_euler);
   DoFTools::make_hanging_node_constraints(dof_handler_euler, constraints_euler);
   constraints_euler.close();
 
   AffineConstraints<double> constraints;
-  constraints.reinit(locally_relevant_dofs);
+  constraints.reinit(locally_owned_dofs, locally_relevant_dofs);
   DoFTools::make_hanging_node_constraints(dof_handler, constraints);
   constraints_euler.close();
 
index 47cb3eb1eaf1bc24065d2b287d3077c0343fc214..33405012e1f686050d47b7690a8c7dfdb192689b 100644 (file)
@@ -208,12 +208,13 @@ test(const unsigned int n_ref = 0)
 
   // constraints:
   AffineConstraints<double> constraints_euler;
-  constraints_euler.reinit(locally_relevant_dofs_euler);
+  constraints_euler.reinit(locally_owned_dofs_euler,
+                           locally_relevant_dofs_euler);
   DoFTools::make_hanging_node_constraints(dof_handler_euler, constraints_euler);
   constraints_euler.close();
 
   AffineConstraints<double> constraints;
-  constraints.reinit(locally_relevant_dofs);
+  constraints.reinit(locally_owned_dofs, locally_relevant_dofs);
   DoFTools::make_hanging_node_constraints(dof_handler, constraints);
   constraints.close();
 
@@ -333,7 +334,8 @@ test(const unsigned int n_ref = 0)
       AffineConstraints<double> level_constraints;
       const IndexSet            relevant_dofs =
         DoFTools::extract_locally_relevant_level_dofs(dof_handler, level);
-      level_constraints.reinit(relevant_dofs);
+      level_constraints.reinit(dof_handler.locally_owned_mg_dofs(level),
+                               relevant_dofs);
       level_constraints.add_lines(
         mg_constrained_dofs.get_boundary_indices(level));
       level_constraints.close();
index 62702d26b6ebb8a8f93e5fc4f643fd99dfa8c992..8715fa33915fbf81d24222eb6a746be8e542cc07 100644 (file)
@@ -147,7 +147,7 @@ test()
       AffineConstraints<double> level_constraints;
       const IndexSet            relevant_dofs =
         DoFTools::extract_locally_relevant_level_dofs(dof, level);
-      level_constraints.reinit(relevant_dofs);
+      level_constraints.reinit(dof.locally_owned_mg_dofs(level), relevant_dofs);
       level_constraints.add_lines(
         mg_constrained_dofs.get_boundary_indices(level));
       level_constraints.close();
index ce9e714a190b488e4c82736608fcc865a2829c95..53a9cecadf3e788d5b1336ea948f065b291ac441 100644 (file)
@@ -69,7 +69,7 @@ test()
   const IndexSet            locally_relevant_dofs =
     DoFTools::extract_locally_relevant_dofs(dof_handler);
 
-  constraints.reinit(locally_relevant_dofs);
+  constraints.reinit(dof_handler.locally_owned_dofs(), locally_relevant_dofs);
 
   DoFTools::make_hanging_node_constraints(dof_handler, constraints);
 
index e355c037d903b3ba160b92f11df17082a08e8964..cd44792177c7adbb93c376e1a5341e022955db9f 100644 (file)
@@ -75,7 +75,7 @@ test(const bool adaptive_ref = true)
   const IndexSet &owned_set    = dof.locally_owned_dofs();
   const IndexSet  relevant_set = DoFTools::extract_locally_relevant_dofs(dof);
 
-  AffineConstraints<double> constraints(relevant_set);
+  AffineConstraints<double> constraints(owned_set, relevant_set);
   DoFTools::make_hanging_node_constraints(dof, constraints);
   // constrain bottom part of the boundary (lower in y direction)
   VectorTools::interpolate_boundary_values(dof,
index 59d9cebc1286bf08bedda8d05ef8e142861e8fe8..302d5e6e2aca2698def5c94e4f0d1a6ac39bb653 100644 (file)
@@ -77,7 +77,7 @@ test(const bool adaptive_ref = true)
   const IndexSet &owned_set    = dof.locally_owned_dofs();
   const IndexSet  relevant_set = DoFTools::extract_locally_relevant_dofs(dof);
 
-  AffineConstraints<double> constraints(relevant_set);
+  AffineConstraints<double> constraints(owned_set, relevant_set);
   DoFTools::make_hanging_node_constraints(dof, constraints);
   // constrain bottom part of the boundary (lower in y direction)
   VectorTools::interpolate_boundary_values(dof,
index 1a738722e951ee3c57d7134c9f91bcbaf967da55..158174a86467afd36f058a7abb1ba78dd4ee409f 100644 (file)
@@ -133,7 +133,7 @@ test(const unsigned int n_glob_ref = 2, const unsigned int n_ref = 0)
 
   // constraints:
   AffineConstraints<double> constraints;
-  constraints.reinit(locally_relevant_dofs);
+  constraints.reinit(locally_owned_dofs, locally_relevant_dofs);
   DoFTools::make_hanging_node_constraints(dof_handler, constraints);
   constraints.close();
 
index 0615c840c3e6b9225c875f515309598b56d667c3..f34c0f9dd8c01e13938bf9cf0b7a42ceaa4948c6 100644 (file)
@@ -90,7 +90,7 @@ test()
   const IndexSet &owned_set    = dof.locally_owned_dofs();
   const IndexSet  relevant_set = DoFTools::extract_locally_relevant_dofs(dof);
 
-  AffineConstraints<double> constraints(relevant_set);
+  AffineConstraints<double> constraints(owned_set, relevant_set);
   DoFTools::make_hanging_node_constraints(dof, constraints);
   VectorTools::interpolate_boundary_values(dof,
                                            0,
index 4935a8798b6878eb2e3fb3b60fba9c9e55c86dc8..083c8807973fda0746e6c993ad8e276081f3ddc1 100644 (file)
@@ -131,7 +131,7 @@ test()
   const IndexSet &owned_set    = dof.locally_owned_dofs();
   const IndexSet  relevant_set = DoFTools::extract_locally_relevant_dofs(dof);
 
-  AffineConstraints<double> constraints(relevant_set);
+  AffineConstraints<double> constraints(owned_set, relevant_set);
   DoFTools::make_hanging_node_constraints(dof, constraints);
   VectorTools::interpolate_boundary_values(dof,
                                            0,
index a1fb3ca8b777c9703bd35943fdd86a02ae3e71e6..821a7b7d01f6890cc6ee0e5d101017dd043d7cba 100644 (file)
@@ -85,7 +85,7 @@ test()
   const IndexSet &owned_set    = dof.locally_owned_dofs();
   const IndexSet  relevant_set = DoFTools::extract_locally_relevant_dofs(dof);
 
-  AffineConstraints<double> constraints(relevant_set);
+  AffineConstraints<double> constraints(owned_set, relevant_set);
   DoFTools::make_hanging_node_constraints(dof, constraints);
   VectorTools::interpolate_boundary_values(dof,
                                            0,
index 1a3756e36b38a8f3dcabe9f7bdb27864ae8495f8..ebff63756641fd421f800cc53315b2e3fbaa59fb 100644 (file)
@@ -91,7 +91,7 @@ test()
   const IndexSet &owned_set    = dof.locally_owned_dofs();
   const IndexSet  relevant_set = DoFTools::extract_locally_relevant_dofs(dof);
 
-  AffineConstraints<double> constraints(relevant_set);
+  AffineConstraints<double> constraints(owned_set, relevant_set);
   DoFTools::make_hanging_node_constraints(dof, constraints);
   VectorTools::interpolate_boundary_values(dof,
                                            0,
index 3900c0f00e7d8f0ade095cdd9f93fc3f9c440a81..1598946c28a4ae6a6bef8db13e42477932ef1532 100644 (file)
@@ -97,7 +97,7 @@ test()
   const IndexSet &owned_set    = dof.locally_owned_dofs();
   const IndexSet  relevant_set = DoFTools::extract_locally_relevant_dofs(dof);
 
-  AffineConstraints<double> constraints(relevant_set);
+  AffineConstraints<double> constraints(owned_set, relevant_set);
   DoFTools::make_hanging_node_constraints(dof, constraints);
   VectorTools::interpolate_boundary_values(dof,
                                            0,
index 296c637ed77bcefa94f2673f39a613a0f3f6a1e3..30c24410a0448451209f6b07e473596c2c057297 100644 (file)
@@ -62,7 +62,7 @@ test()
   const IndexSet &owned_set    = dof.locally_owned_dofs();
   const IndexSet  relevant_set = DoFTools::extract_locally_relevant_dofs(dof);
 
-  AffineConstraints<double> constraints(relevant_set);
+  AffineConstraints<double> constraints(owned_set, relevant_set);
   DoFTools::make_hanging_node_constraints(dof, constraints);
   VectorTools::interpolate_boundary_values(dof,
                                            0,
index 7eb76845282b77ed36d94df441572010ad068b3b..8afe150a0e4f56283a90b04bd6ba654ef26cbaf5 100644 (file)
@@ -63,7 +63,7 @@ test()
   const IndexSet &owned_set    = dof.locally_owned_dofs();
   const IndexSet  relevant_set = DoFTools::extract_locally_relevant_dofs(dof);
 
-  AffineConstraints<double> constraints(relevant_set);
+  AffineConstraints<double> constraints(owned_set, relevant_set);
   DoFTools::make_hanging_node_constraints(dof, constraints);
   constraints.close();
 
index 1bf88214ba5571430b5bd91247cdf1a340e6aa6e..03c224db8ce3903a20d99bec28fff788386134cf 100644 (file)
@@ -58,8 +58,8 @@ test()
   const IndexSet &owned_set    = dof.locally_owned_dofs();
   const IndexSet  relevant_set = DoFTools::extract_locally_relevant_dofs(dof);
 
-  AffineConstraints<double> constraints_0(relevant_set),
-    constraints_1(relevant_set);
+  AffineConstraints<double> constraints_0(owned_set, relevant_set),
+    constraints_1(owned_set, relevant_set);
   VectorTools::interpolate_boundary_values(dof,
                                            0,
                                            Functions::ZeroFunction<dim>(),
index 07af8c270df38cdbe73d9f41610d2f01c01b08f4..273e544186121aae4f70fe914fb5f64c14e8416e 100644 (file)
@@ -97,7 +97,7 @@ test()
   IndexSet relevant_set;
   DoFTools::extract_locally_relevant_dofs(dof, relevant_set);
 
-  AffineConstraints<double> constraints(relevant_set);
+  AffineConstraints<double> constraints(owned_set, relevant_set);
   DoFTools::make_hanging_node_constraints(dof, constraints);
   VectorTools::interpolate_boundary_values(dof,
                                            0,
index f878620eb4023a1be6ec08efc50e4bdf31873325..9990cd7f16122f36b5086b9e2f7ce625a2480606 100644 (file)
@@ -89,7 +89,7 @@ test()
   const IndexSet &owned_set    = dof.locally_owned_dofs();
   const IndexSet  relevant_set = DoFTools::extract_locally_relevant_dofs(dof);
 
-  AffineConstraints<double> constraints(relevant_set);
+  AffineConstraints<double> constraints(owned_set, relevant_set);
   DoFTools::make_hanging_node_constraints(dof, constraints);
   VectorTools::interpolate_boundary_values(dof,
                                            0,
index 824a9e1451aac6d76a778c2a78c01825bce01d00..f8f5536743bb9a9d6badcde948eef1e82a888abd 100644 (file)
@@ -154,7 +154,7 @@ test()
   const IndexSet &owned_set    = dof.locally_owned_dofs();
   const IndexSet  relevant_set = DoFTools::extract_locally_relevant_dofs(dof);
 
-  AffineConstraints<double> constraints(relevant_set);
+  AffineConstraints<double> constraints(owned_set, relevant_set);
   DoFTools::make_hanging_node_constraints(dof, constraints);
   VectorTools::interpolate_boundary_values(dof,
                                            0,
index da1e600146ab8c1bd2c6c9c601924df2e16b112a..603785ba117cbdbcd9f274f927e7f64d0481f181 100644 (file)
@@ -152,7 +152,7 @@ test()
   const IndexSet &owned_set    = dof.locally_owned_dofs();
   const IndexSet  relevant_set = DoFTools::extract_locally_relevant_dofs(dof);
 
-  AffineConstraints<double> constraints(relevant_set);
+  AffineConstraints<double> constraints(owned_set, relevant_set);
   DoFTools::make_hanging_node_constraints(dof, constraints);
   VectorTools::interpolate_boundary_values(dof,
                                            0,
index 5dadfb924796345aec5f028d6d3819915fdc6540..6a884aa82251d658f6933fec7c08f609cfef6989 100644 (file)
@@ -90,7 +90,7 @@ test()
   const IndexSet &owned_set    = dof.locally_owned_dofs();
   const IndexSet  relevant_set = DoFTools::extract_locally_relevant_dofs(dof);
 
-  AffineConstraints<double> constraints(relevant_set);
+  AffineConstraints<double> constraints(owned_set, relevant_set);
   DoFTools::make_hanging_node_constraints(dof, constraints);
   VectorTools::interpolate_boundary_values(dof,
                                            0,
index 3cac8a93cd1621025355638464262d4a06471161..18c1fb3095ad6ed281564b3e4573270bec8536ec 100644 (file)
@@ -155,7 +155,7 @@ test()
   const IndexSet &owned_set    = dof.locally_owned_dofs();
   const IndexSet  relevant_set = DoFTools::extract_locally_relevant_dofs(dof);
 
-  AffineConstraints<double> constraints(relevant_set);
+  AffineConstraints<double> constraints(owned_set, relevant_set);
   DoFTools::make_hanging_node_constraints(dof, constraints);
   VectorTools::interpolate_boundary_values(dof,
                                            0,
index e71dd2ceb6b5ed92ff28538b9cecd7be136405ff..729eb5aa8a49256afb11643c5afbd40271f6c4cf 100644 (file)
@@ -167,7 +167,7 @@ test()
   const IndexSet &owned_set    = dof.locally_owned_dofs();
   const IndexSet  relevant_set = DoFTools::extract_locally_relevant_dofs(dof);
 
-  AffineConstraints<double> constraints(relevant_set);
+  AffineConstraints<double> constraints(owned_set, relevant_set);
   DoFTools::make_hanging_node_constraints(dof, constraints);
   VectorTools::interpolate_boundary_values(dof,
                                            0,
index 0505d6af532eb46e8c1377c838a28e225b423c6c..25a2a1c9274aee4520f4ca85ab9e2c9347bbe9e4 100644 (file)
@@ -91,7 +91,7 @@ test()
   IndexSet relevant_set;
   DoFTools::extract_locally_relevant_dofs(dof, relevant_set);
 
-  AffineConstraints<double> constraints(relevant_set);
+  AffineConstraints<double> constraints(owned_set, relevant_set);
   DoFTools::make_hanging_node_constraints(dof, constraints);
   VectorTools::interpolate_boundary_values(dof,
                                            0,
index 810b7d0c24bcef539bdea941752aab38dd05381a..6ceb8b7083d7e65a794f31c9b3a651a3f8173a05 100644 (file)
@@ -68,7 +68,7 @@ test()
   IndexSet relevant_set;
   DoFTools::extract_locally_relevant_dofs(dof, relevant_set);
 
-  AffineConstraints<double> constraints(relevant_set);
+  AffineConstraints<double> constraints(owned_set, relevant_set);
   DoFTools::make_hanging_node_constraints(dof, constraints);
   VectorTools::interpolate_boundary_values(dof,
                                            0,
index 161b9a6c1968e44b5cb25cb017372d084de7ac42..51c0f3f69accc67462b118a286ffb0dc3c91cf4d 100644 (file)
@@ -135,7 +135,7 @@ test()
   const IndexSet &owned_set    = dof.locally_owned_dofs();
   const IndexSet  relevant_set = DoFTools::extract_locally_relevant_dofs(dof);
 
-  AffineConstraints<double> constraints(relevant_set);
+  AffineConstraints<double> constraints(owned_set, relevant_set);
   DoFTools::make_hanging_node_constraints(dof, constraints);
   VectorTools::interpolate_boundary_values(dof,
                                            0,
index da47c9af52c9a4274ac9b3844310c0fc854d3de3..72c3f02c945c14e9cd5deca66c40d4a5c997b43f 100644 (file)
@@ -75,7 +75,7 @@ test()
   const IndexSet &owned_set    = dof.locally_owned_dofs();
   const IndexSet  relevant_set = DoFTools::extract_locally_relevant_dofs(dof);
 
-  AffineConstraints<double> constraints(relevant_set);
+  AffineConstraints<double> constraints(owned_set, relevant_set);
   DoFTools::make_hanging_node_constraints(dof, constraints);
   VectorTools::interpolate_boundary_values(dof,
                                            0,
index 047a0a8d86feb60b6e2eb7bf9a61e5cc5b091389..cc36986e3ae963d917591ae5b6422069d4e35fda 100644 (file)
@@ -107,7 +107,7 @@ do_test(const DoFHandler<dim> &dof)
 
   // fine-level constraints
   AffineConstraints<double> constraints;
-  constraints.reinit(locally_relevant_dofs);
+  constraints.reinit(dof.locally_owned_dofs(), locally_relevant_dofs);
   DoFTools::make_hanging_node_constraints(dof, constraints);
   VectorTools::interpolate_boundary_values(dof,
                                            dirichlet_boundary,
@@ -174,7 +174,7 @@ do_test(const DoFHandler<dim> &dof)
       AffineConstraints<double> level_constraints;
       const IndexSet            relevant_dofs =
         DoFTools::extract_locally_relevant_level_dofs(dof, level);
-      level_constraints.reinit(relevant_dofs);
+      level_constraints.reinit(dof.locally_owned_mg_dofs(level), relevant_dofs);
       level_constraints.add_lines(
         mg_constrained_dofs.get_boundary_indices(level));
       level_constraints.close();
index 74b0ca5b150f2ce4895a47128b8f0528c7bfd9bf..7c1a997d71371b475b9ebf9960fca4ac2ef00f84 100644 (file)
@@ -79,7 +79,7 @@ public:
       {
         const IndexSet relevant_dofs =
           DoFTools::extract_locally_relevant_dofs(dof_handler);
-        constraints.reinit(relevant_dofs);
+        constraints.reinit(dof_handler.locally_owned_dofs(), relevant_dofs);
         DoFTools::make_hanging_node_constraints(dof_handler, constraints);
         VectorTools::interpolate_boundary_values(dof_handler,
                                                  dirichlet_boundary,
@@ -87,9 +87,11 @@ public:
       }
     else
       {
+        const IndexSet &locally_owned =
+          dof_handler.locally_owned_mg_dofs(level);
         const IndexSet relevant_dofs =
           DoFTools::extract_locally_relevant_level_dofs(dof_handler, level);
-        constraints.reinit(relevant_dofs);
+        constraints.reinit(locally_owned, relevant_dofs);
         constraints.add_lines(mg_constrained_dofs.get_boundary_indices(level));
 
         const std::vector<types::global_dof_index> interface_indices =
@@ -98,8 +100,6 @@ public:
         edge_constrained_indices.clear();
         edge_constrained_indices.reserve(interface_indices.size());
         edge_constrained_values.resize(interface_indices.size());
-        const IndexSet &locally_owned =
-          dof_handler.locally_owned_mg_dofs(level);
         for (unsigned int i = 0; i < interface_indices.size(); ++i)
           if (locally_owned.is_element(interface_indices[i]))
             edge_constrained_indices.push_back(
@@ -478,7 +478,8 @@ do_test(const DoFHandler<dim> &dof)
   AffineConstraints<double> hanging_node_constraints;
   const IndexSet            locally_relevant_dofs =
     DoFTools::extract_locally_relevant_dofs(dof);
-  hanging_node_constraints.reinit(locally_relevant_dofs);
+  hanging_node_constraints.reinit(dof.locally_owned_dofs(),
+                                  locally_relevant_dofs);
   DoFTools::make_hanging_node_constraints(dof, hanging_node_constraints);
   hanging_node_constraints.close();
 
index c1dcbc6201e031c466c5d27428c6e551335ada70..e65ad7c017e4e4cf8a10bbaba6ab9403e435a99c 100644 (file)
@@ -108,7 +108,7 @@ do_test(const DoFHandler<dim> &dof)
 
   // fine-level constraints
   AffineConstraints<double> constraints;
-  constraints.reinit(locally_relevant_dofs);
+  constraints.reinit(dof.locally_owned_dofs(), locally_relevant_dofs);
   DoFTools::make_hanging_node_constraints(dof, constraints);
   VectorTools::interpolate_boundary_values(dof,
                                            dirichlet_boundary,
@@ -153,7 +153,8 @@ do_test(const DoFHandler<dim> &dof)
   {
     // this is to make it consistent with parallel_multigrid_adaptive.cc
     AffineConstraints<double> hanging_node_constraints;
-    hanging_node_constraints.reinit(locally_relevant_dofs);
+    hanging_node_constraints.reinit(dof.locally_owned_dofs(),
+                                    locally_relevant_dofs);
     DoFTools::make_hanging_node_constraints(dof, hanging_node_constraints);
     hanging_node_constraints.close();
 
@@ -188,7 +189,7 @@ do_test(const DoFHandler<dim> &dof)
       AffineConstraints<double> level_constraints;
       const IndexSet            relevant_dofs =
         DoFTools::extract_locally_relevant_level_dofs(dof, level);
-      level_constraints.reinit(relevant_dofs);
+      level_constraints.reinit(dof.locally_owned_mg_dofs(level), relevant_dofs);
       level_constraints.add_lines(
         mg_constrained_dofs.get_boundary_indices(level));
       level_constraints.close();
index 6d21864016012abb7066f93088ec30c4600216e3..f0fbf1b5cd9dbfbf90e22055a70a0ec8d027c67b 100644 (file)
@@ -86,7 +86,7 @@ public:
       {
         const IndexSet relevant_dofs =
           DoFTools::extract_locally_relevant_dofs(dof_handler);
-        constraints.reinit(relevant_dofs);
+        constraints.reinit(dof_handler.locally_owned_dofs(), relevant_dofs);
         DoFTools::make_hanging_node_constraints(dof_handler, constraints);
         VectorTools::interpolate_boundary_values(dof_handler,
                                                  dirichlet_boundary,
@@ -94,9 +94,11 @@ public:
       }
     else
       {
+        const IndexSet &locally_owned =
+          dof_handler.locally_owned_mg_dofs(level);
         const IndexSet relevant_dofs =
           DoFTools::extract_locally_relevant_level_dofs(dof_handler, level);
-        constraints.reinit(relevant_dofs);
+        constraints.reinit(locally_owned, relevant_dofs);
         constraints.add_lines(mg_constrained_dofs.get_boundary_indices(level));
 
         const std::vector<types::global_dof_index> interface_indices =
@@ -105,8 +107,6 @@ public:
         edge_constrained_indices.clear();
         edge_constrained_indices.reserve(interface_indices.size());
         edge_constrained_values.resize(interface_indices.size());
-        const IndexSet &locally_owned =
-          dof_handler.locally_owned_mg_dofs(level);
         for (unsigned int i = 0; i < interface_indices.size(); ++i)
           if (locally_owned.is_element(interface_indices[i]))
             edge_constrained_indices.push_back(
@@ -482,7 +482,8 @@ do_test(const DoFHandler<dim> &dof, const bool threaded)
   AffineConstraints<double> hanging_node_constraints;
   const IndexSet            locally_relevant_dofs =
     DoFTools::extract_locally_relevant_dofs(dof);
-  hanging_node_constraints.reinit(locally_relevant_dofs);
+  hanging_node_constraints.reinit(dof.locally_owned_dofs(),
+                                  locally_relevant_dofs);
   DoFTools::make_hanging_node_constraints(dof, hanging_node_constraints);
   hanging_node_constraints.close();
 
index 77d528149c18babf6ee9e10e36e871ce498097fe..606c59be453ec2177356a11c92391b297bf2c4bb 100644 (file)
@@ -108,7 +108,7 @@ do_test(const DoFHandler<dim> &dof)
 
   // fine-level constraints
   AffineConstraints<double> constraints;
-  constraints.reinit(locally_relevant_dofs);
+  constraints.reinit(dof.locally_owned_dofs(), locally_relevant_dofs);
   DoFTools::make_hanging_node_constraints(dof, constraints);
   VectorTools::interpolate_boundary_values(dof,
                                            dirichlet_boundary,
@@ -153,7 +153,8 @@ do_test(const DoFHandler<dim> &dof)
   {
     // this is to make it consistent with parallel_multigrid_adaptive.cc
     AffineConstraints<double> hanging_node_constraints;
-    hanging_node_constraints.reinit(locally_relevant_dofs);
+    hanging_node_constraints.reinit(dof.locally_owned_dofs(),
+                                    locally_relevant_dofs);
     DoFTools::make_hanging_node_constraints(dof, hanging_node_constraints);
     hanging_node_constraints.close();
 
@@ -188,7 +189,7 @@ do_test(const DoFHandler<dim> &dof)
       AffineConstraints<double> level_constraints;
       const IndexSet            relevant_dofs =
         DoFTools::extract_locally_relevant_level_dofs(dof, level);
-      level_constraints.reinit(relevant_dofs);
+      level_constraints.reinit(dof.locally_owned_mg_dofs(level), relevant_dofs);
       level_constraints.add_lines(
         mg_constrained_dofs.get_boundary_indices(level));
       level_constraints.close();
index ab2de8241aadadbf525a801616dfd00ac9f131f2..5e9cbb12b00176863cf2b634cb1fd21930b8a563 100644 (file)
@@ -83,7 +83,7 @@ public:
       {
         const IndexSet relevant_dofs =
           DoFTools::extract_locally_relevant_dofs(dof_handler);
-        constraints.reinit(relevant_dofs);
+        constraints.reinit(dof_handler.locally_owned_dofs(), relevant_dofs);
         DoFTools::make_hanging_node_constraints(dof_handler, constraints);
         VectorTools::interpolate_boundary_values(dof_handler,
                                                  dirichlet_boundary,
@@ -91,9 +91,11 @@ public:
       }
     else
       {
+        const IndexSet &locally_owned =
+          dof_handler.locally_owned_mg_dofs(level);
         const IndexSet relevant_dofs =
           DoFTools::extract_locally_relevant_level_dofs(dof_handler, level);
-        constraints.reinit(relevant_dofs);
+        constraints.reinit(locally_owned, relevant_dofs);
         constraints.add_lines(mg_constrained_dofs.get_boundary_indices(level));
 
         const std::vector<types::global_dof_index> interface_indices =
@@ -102,8 +104,6 @@ public:
         edge_constrained_indices.clear();
         edge_constrained_indices.reserve(interface_indices.size());
         edge_constrained_values.resize(interface_indices.size());
-        const IndexSet &locally_owned =
-          dof_handler.locally_owned_mg_dofs(level);
         for (unsigned int i = 0; i < interface_indices.size(); ++i)
           if (locally_owned.is_element(interface_indices[i]))
             edge_constrained_indices.push_back(
@@ -479,7 +479,8 @@ do_test(const DoFHandler<dim> &dof, const bool threaded)
   AffineConstraints<double> hanging_node_constraints;
   const IndexSet            locally_relevant_dofs =
     DoFTools::extract_locally_relevant_dofs(dof);
-  hanging_node_constraints.reinit(locally_relevant_dofs);
+  hanging_node_constraints.reinit(dof.locally_owned_dofs(),
+                                  locally_relevant_dofs);
   DoFTools::make_hanging_node_constraints(dof, hanging_node_constraints);
   hanging_node_constraints.close();
 
index 9083ac6fe709f8af0ba060079ac43232e244c2f6..44c1a6a40456fa2e67e42535615a94c12df15060 100644 (file)
@@ -230,7 +230,7 @@ do_test(const DoFHandler<dim> &dof, const unsigned int nb)
 
   // fine-level constraints
   AffineConstraints<double> constraints;
-  constraints.reinit(locally_relevant_dofs);
+  constraints.reinit(dof.locally_owned_dofs(), locally_relevant_dofs);
   DoFTools::make_hanging_node_constraints(dof, constraints);
   VectorTools::interpolate_boundary_values(dof,
                                            dirichlet_boundary,
@@ -279,7 +279,8 @@ do_test(const DoFHandler<dim> &dof, const unsigned int nb)
   {
     // this is to make it consistent with parallel_multigrid_adaptive.cc
     AffineConstraints<double> hanging_node_constraints;
-    hanging_node_constraints.reinit(locally_relevant_dofs);
+    hanging_node_constraints.reinit(dof.locally_owned_dofs(),
+                                    locally_relevant_dofs);
     DoFTools::make_hanging_node_constraints(dof, hanging_node_constraints);
     hanging_node_constraints.close();
 
@@ -316,7 +317,7 @@ do_test(const DoFHandler<dim> &dof, const unsigned int nb)
       AffineConstraints<double> level_constraints;
       const IndexSet            relevant_dofs =
         DoFTools::extract_locally_relevant_level_dofs(dof, level);
-      level_constraints.reinit(relevant_dofs);
+      level_constraints.reinit(dof.locally_owned_mg_dofs(level), relevant_dofs);
       level_constraints.add_lines(
         mg_constrained_dofs.get_boundary_indices(level));
       level_constraints.close();
index 77ddfadfbd123c5984cfeac8ccacd90b92767c34..925862c2bfeceb71ca8ec622c25bfc56195ff47f 100644 (file)
@@ -109,7 +109,7 @@ do_test(const DoFHandler<dim> &dof)
 
   // fine-level constraints
   AffineConstraints<double> constraints;
-  constraints.reinit(locally_relevant_dofs);
+  constraints.reinit(dof.locally_owned_dofs(), locally_relevant_dofs);
   DoFTools::make_hanging_node_constraints(dof, constraints);
   VectorTools::interpolate_boundary_values(dof,
                                            dirichlet_boundary,
@@ -154,7 +154,8 @@ do_test(const DoFHandler<dim> &dof)
   {
     // this is to make it consistent with parallel_multigrid_adaptive.cc
     AffineConstraints<double> hanging_node_constraints;
-    hanging_node_constraints.reinit(locally_relevant_dofs);
+    hanging_node_constraints.reinit(dof.locally_owned_dofs(),
+                                    locally_relevant_dofs);
     DoFTools::make_hanging_node_constraints(dof, hanging_node_constraints);
     hanging_node_constraints.close();
 
@@ -189,7 +190,7 @@ do_test(const DoFHandler<dim> &dof)
       AffineConstraints<double> level_constraints;
       const IndexSet            relevant_dofs =
         DoFTools::extract_locally_relevant_level_dofs(dof, level);
-      level_constraints.reinit(relevant_dofs);
+      level_constraints.reinit(dof.locally_owned_mg_dofs(level), relevant_dofs);
       level_constraints.add_lines(
         mg_constrained_dofs.get_boundary_indices(level));
       level_constraints.close();
index d3937c2c88b6e0f1b91d75848434c94cbb362afc..aa946558b099ca1d60468f77706b5e24abdd183c 100644 (file)
@@ -252,7 +252,8 @@ do_test(const std::vector<const DoFHandler<dim> *> &dof)
   std::vector<const AffineConstraints<double> *> constraints_ptrs(dof.size());
   for (unsigned int i = 0; i < dof.size(); ++i)
     {
-      constraints[i].reinit(locally_relevant_dofs[i]);
+      constraints[i].reinit(dof[i]->locally_owned_dofs(),
+                            locally_relevant_dofs[i]);
       DoFTools::make_hanging_node_constraints(*dof[i], constraints[i]);
       VectorTools::interpolate_boundary_values(*dof[i],
                                                dirichlet_boundary,
@@ -301,7 +302,8 @@ do_test(const std::vector<const DoFHandler<dim> *> &dof)
         // this is to make it consistent with parallel_multigrid_adaptive.cc
         AffineConstraints<double> hanging_node_constraints;
 
-        hanging_node_constraints.reinit(locally_relevant_dofs[b]);
+        hanging_node_constraints.reinit(dof[b]->locally_owned_dofs(),
+                                        locally_relevant_dofs[b]);
         DoFTools::make_hanging_node_constraints(*dof[b],
                                                 hanging_node_constraints);
         hanging_node_constraints.close();
@@ -350,7 +352,8 @@ do_test(const std::vector<const DoFHandler<dim> *> &dof)
         {
           const IndexSet relevant_dofs =
             DoFTools::extract_locally_relevant_level_dofs(*dof[i], level);
-          level_constraints[i].reinit(relevant_dofs);
+          level_constraints[i].reinit(dof[i]->locally_owned_mg_dofs(level),
+                                      relevant_dofs);
           level_constraints[i].add_lines(
             mg_constrained_dofs[i].get_boundary_indices(level));
           level_constraints[i].close();
index 4d0208ddf723d5699058b30737dfcf5c5f7e3675..168b2caf9fec0a0c5af164ab24667feec1ea4087 100644 (file)
@@ -226,7 +226,7 @@ do_test(const DoFHandler<dim> &dof)
 
   // fine-level constraints
   AffineConstraints<double> constraints;
-  constraints.reinit(locally_relevant_dofs);
+  constraints.reinit(dof.locally_owned_dofs(), locally_relevant_dofs);
   DoFTools::make_hanging_node_constraints(dof, constraints);
   VectorTools::interpolate_boundary_values(dof,
                                            dirichlet_boundary,
@@ -266,7 +266,8 @@ do_test(const DoFHandler<dim> &dof)
   {
     // this is to make it consistent with parallel_multigrid_adaptive.cc
     AffineConstraints<double> hanging_node_constraints;
-    hanging_node_constraints.reinit(locally_relevant_dofs);
+    hanging_node_constraints.reinit(dof.locally_owned_dofs(),
+                                    locally_relevant_dofs);
     DoFTools::make_hanging_node_constraints(dof, hanging_node_constraints);
     hanging_node_constraints.close();
 
@@ -296,7 +297,7 @@ do_test(const DoFHandler<dim> &dof)
       AffineConstraints<double> level_constraints;
       const IndexSet            relevant_dofs =
         DoFTools::extract_locally_relevant_level_dofs(dof, level);
-      level_constraints.reinit(relevant_dofs);
+      level_constraints.reinit(dof.locally_owned_mg_dofs(level), relevant_dofs);
       level_constraints.add_lines(
         mg_constrained_dofs.get_boundary_indices(level));
       level_constraints.close();
index 7bac0350a362ed4fc85880086dc6e972382d649e..223d0dfee49a2888b01667e92b6ab3ce375c7d02 100644 (file)
@@ -107,7 +107,7 @@ do_test(const DoFHandler<dim> &dof)
   const IndexSet locally_relevant_dofs =
     DoFTools::extract_locally_relevant_dofs(dof);
   AffineConstraints<double> constraints;
-  constraints.reinit(locally_relevant_dofs);
+  constraints.reinit(dof.locally_owned_dofs(), locally_relevant_dofs);
   VectorTools::interpolate_boundary_values(dof,
                                            dirichlet_boundary,
                                            constraints);
@@ -161,7 +161,7 @@ do_test(const DoFHandler<dim> &dof)
       AffineConstraints<double> level_constraints;
       const IndexSet            relevant_dofs =
         DoFTools::extract_locally_relevant_level_dofs(dof, level);
-      level_constraints.reinit(relevant_dofs);
+      level_constraints.reinit(dof.locally_owned_mg_dofs(level), relevant_dofs);
       level_constraints.add_lines(
         mg_constrained_dofs.get_boundary_indices(level));
       level_constraints.close();
index dd690aba27d1577b89c1f61bbd5882ecad2be409..def6fefd785c29df0cb4808181119cbf91fe902c 100644 (file)
@@ -141,7 +141,7 @@ test()
   dof.distribute_dofs(fe);
   AffineConstraints<double> constraints;
   const IndexSet relevant_dofs = DoFTools::extract_locally_relevant_dofs(dof);
-  constraints.reinit(relevant_dofs);
+  constraints.reinit(dof.locally_owned_dofs(), relevant_dofs);
   VectorTools::interpolate_boundary_values(dof,
                                            0,
                                            Functions::ZeroFunction<dim>(),
index 1539f7bc568b37695154851eb8656f3d3f7e86d6..bc1ce257613b83c5c53988cdcb28b0006f847ede 100644 (file)
@@ -192,23 +192,27 @@ test(const unsigned int fe_degree)
   DoFHandler<dim> dof(tria);
   dof.distribute_dofs(fe);
 
-  IndexSet                  owned_set = dof.locally_owned_dofs();
-  IndexSet                  relevant_set;
-  AffineConstraints<double> constraints(relevant_set);
   typename MatrixFree<dim, number>::AdditionalData addit_data;
   addit_data.tasks_parallel_scheme =
     MatrixFree<dim, number>::AdditionalData::none;
 
+  IndexSet                  relevant_set;
+  AffineConstraints<double> constraints;
   {
+    constraints.clear();
     relevant_set = DoFTools::extract_locally_relevant_dofs(dof);
+    constraints.reinit(dof.locally_owned_dofs(), relevant_set);
     constraints.close();
-
-    DoFRenumbering::matrix_free_data_locality(dof, constraints, addit_data);
   }
 
-  constraints.clear();
-  relevant_set = DoFTools::extract_locally_relevant_dofs(dof);
-  constraints.close();
+  DoFRenumbering::matrix_free_data_locality(dof, constraints, addit_data);
+
+  {
+    constraints.clear();
+    relevant_set = DoFTools::extract_locally_relevant_dofs(dof);
+    constraints.reinit(dof.locally_owned_dofs(), relevant_set);
+    constraints.close();
+  }
 
   const QGauss<1>         quadrature(dof.get_fe().degree + 1);
   MatrixFree<dim, number> mf_data;
index 65b44279a4c66a2016e6f8d56fb06b3e179f6be2..744016e34301f221415e1b366b110b89d30eff44 100644 (file)
@@ -402,7 +402,7 @@ namespace Step37
       DoFTools::extract_locally_relevant_dofs(dof_handler);
 
     constraints.clear();
-    constraints.reinit(locally_relevant_dofs);
+    constraints.reinit(dof_handler.locally_owned_dofs(), locally_relevant_dofs);
     DoFTools::make_hanging_node_constraints(dof_handler, constraints);
     VectorTools::interpolate_boundary_values(dof_handler,
                                              0,
@@ -444,7 +444,8 @@ namespace Step37
         const IndexSet relevant_dofs =
           DoFTools::extract_locally_relevant_level_dofs(dof_handler, level);
         AffineConstraints<double> level_constraints;
-        level_constraints.reinit(relevant_dofs);
+        level_constraints.reinit(dof_handler.locally_owned_mg_dofs(level),
+                                 relevant_dofs);
         level_constraints.add_lines(
           mg_constrained_dofs.get_boundary_indices(level));
         level_constraints.close();
index c7a737fb6cd2e08febef0ba128251d1b1ce6c1d9..90bc1605bd900f4443de11fa109aab1bfbeca452 100644 (file)
@@ -296,7 +296,7 @@ namespace Step48
     locally_relevant_dofs =
       DoFTools::extract_locally_relevant_dofs(dof_handler);
     constraints.clear();
-    constraints.reinit(locally_relevant_dofs);
+    constraints.reinit(dof_handler.locally_owned_dofs(), locally_relevant_dofs);
     DoFTools::make_hanging_node_constraints(dof_handler, constraints);
     constraints.close();
 
index 104255db82bf4db488347a961bab3db1e3e44f1f..ff860038c9cd9b116a0f93efde3d11fecce28bca 100644 (file)
@@ -266,7 +266,7 @@ namespace Step48
     locally_relevant_dofs =
       DoFTools::extract_locally_relevant_dofs(dof_handler);
     constraints.clear();
-    constraints.reinit(locally_relevant_dofs);
+    constraints.reinit(dof_handler.locally_owned_dofs(), locally_relevant_dofs);
     DoFTools::make_hanging_node_constraints(dof_handler, constraints);
     constraints.close();
 
index cf8f976ec9d6c27b4ba714b7b249548c9bc5a707..adfd4b9e79e596008f69e76f7e965c1c432ea8ac 100644 (file)
@@ -1040,7 +1040,8 @@ namespace StokesClass
 
     const IndexSet locally_relevant_dofs_u =
       DoFTools::extract_locally_relevant_dofs(dof_handler_u);
-    constraints_u.reinit(locally_relevant_dofs_u);
+    constraints_u.reinit(dof_handler_u.locally_owned_dofs(),
+                         locally_relevant_dofs_u);
     DoFTools::make_hanging_node_constraints(dof_handler_u, constraints_u);
 
     VectorTools::interpolate_boundary_values(
@@ -1049,7 +1050,8 @@ namespace StokesClass
 
     const IndexSet locally_relevant_dofs_p =
       DoFTools::extract_locally_relevant_dofs(dof_handler_p);
-    constraints_p.reinit(locally_relevant_dofs_p);
+    constraints_p.reinit(dof_handler_p.locally_owned_dofs(),
+                         locally_relevant_dofs_p);
     DoFTools::make_hanging_node_constraints(dof_handler_p, constraints_p);
     constraints_p.close();
 
@@ -1116,7 +1118,8 @@ namespace StokesClass
         const IndexSet relevant_dofs =
           DoFTools::extract_locally_relevant_level_dofs(dof_handler_u, level);
         AffineConstraints<double> level_constraints;
-        level_constraints.reinit(relevant_dofs);
+        level_constraints.reinit(dof_handler_u.locally_owned_mg_dofs(level),
+                                 relevant_dofs);
         level_constraints.add_lines(
           mg_constrained_dofs.get_boundary_indices(level));
         level_constraints.close();
@@ -1171,7 +1174,8 @@ namespace StokesClass
     AffineConstraints<double> constraints_u_no_dirchlet;
     const IndexSet            locally_relevant_dofs_u =
       DoFTools::extract_locally_relevant_dofs(dof_handler_u);
-    constraints_u_no_dirchlet.reinit(locally_relevant_dofs_u);
+    constraints_u_no_dirchlet.reinit(dof_handler_u.locally_owned_dofs(),
+                                     locally_relevant_dofs_u);
     DoFTools::make_hanging_node_constraints(dof_handler_u,
                                             constraints_u_no_dirchlet);
     constraints_u_no_dirchlet.close();
index 6eb2499267e9b40291065cf02bcea45410071254..5f764b548ed9b9e0341aed7ce1fb5ea5c0e009ba 100644 (file)
@@ -70,7 +70,7 @@ test()
   deallog << "locally relevant dofs :" << std::endl;
   relevant_set.print(deallog.get_file_stream());
 
-  AffineConstraints<double> constraints(relevant_set);
+  AffineConstraints<double> constraints(owned_set, relevant_set);
   constraints.close();
 
   MappingQ<dim>                     mapping(fe_degree);
index 65e2f78368f1c219f91bf6276e6680b13c1d68fd..080c28d35749f4f5e701a39245c91d1dc794e61f 100644 (file)
@@ -93,7 +93,7 @@ test()
   const IndexSet &owned_set    = dof.locally_owned_dofs();
   const IndexSet  relevant_set = DoFTools::extract_locally_relevant_dofs(dof);
 
-  AffineConstraints<double> constraints(relevant_set);
+  AffineConstraints<double> constraints(owned_set, relevant_set);
   DoFTools::make_hanging_node_constraints(dof, constraints);
   VectorTools::interpolate_boundary_values(dof,
                                            0,
index 418c8eb4f74ad9b42a6b6e77d995f13917a261be..b7b9f46bb703fed95c07041c0e0502da2964c1bf 100644 (file)
@@ -95,7 +95,7 @@ test()
   const IndexSet &owned_set    = dof.locally_owned_dofs();
   const IndexSet  relevant_set = DoFTools::extract_locally_relevant_dofs(dof);
 
-  AffineConstraints<double> constraints(relevant_set);
+  AffineConstraints<double> constraints(owned_set, relevant_set);
   DoFTools::make_hanging_node_constraints(dof, constraints);
   VectorTools::interpolate_boundary_values(dof,
                                            0,
index 30ca22b1b5f2de9d5448cfc6de9bbcb25a50fad3..ca73f861b44733c7a5a50586545a6885910cbacc 100644 (file)
@@ -89,7 +89,7 @@ test()
   const IndexSet &owned_set    = dof.locally_owned_dofs();
   const IndexSet  relevant_set = DoFTools::extract_locally_relevant_dofs(dof);
 
-  AffineConstraints<double> constraints(relevant_set);
+  AffineConstraints<double> constraints(owned_set, relevant_set);
   DoFTools::make_hanging_node_constraints(dof, constraints);
   VectorTools::interpolate_boundary_values(dof,
                                            0,
index 9ebdc1ec5f376df591d4298d5e2bc3046e801a35..e217f9bd8808fef401e8123c074fa8bcb51f60a8 100644 (file)
@@ -95,7 +95,7 @@ test()
   const IndexSet &owned_set    = dof.locally_owned_dofs();
   const IndexSet  relevant_set = DoFTools::extract_locally_relevant_dofs(dof);
 
-  AffineConstraints<double> constraints(relevant_set);
+  AffineConstraints<double> constraints(owned_set, relevant_set);
   DoFTools::make_hanging_node_constraints(dof, constraints);
   VectorTools::interpolate_boundary_values(dof,
                                            0,
index a68eb050c8133eec9659067741f11c0b39784b22..bbc3a1cc9e2f7fd8d956382e1920c81a079e9346 100644 (file)
@@ -101,7 +101,7 @@ test()
   const IndexSet &owned_set    = dof.locally_owned_dofs();
   const IndexSet  relevant_set = DoFTools::extract_locally_relevant_dofs(dof);
 
-  AffineConstraints<double> constraints(relevant_set);
+  AffineConstraints<double> constraints(owned_set, relevant_set);
   DoFTools::make_hanging_node_constraints(dof, constraints);
   VectorTools::interpolate_boundary_values(dof,
                                            0,
index 0e8391b8b7c28f367f6f2b7830ff4cf771feca66..050cf6755efe8728b854d257a8abfffb69ff4ba7 100644 (file)
@@ -94,7 +94,7 @@ test()
   const IndexSet &owned_set    = dof.locally_owned_dofs();
   const IndexSet  relevant_set = DoFTools::extract_locally_relevant_dofs(dof);
 
-  AffineConstraints<double> constraints(relevant_set);
+  AffineConstraints<double> constraints(owned_set, relevant_set);
   DoFTools::make_hanging_node_constraints(dof, constraints);
   VectorTools::interpolate_boundary_values(dof,
                                            0,
index fc33408e1dd25ddb0bfbdc38007d04c57997e919..895a6d4159a82ff4597c8df35a5486e8619cc7b9 100644 (file)
@@ -58,7 +58,7 @@ test()
   deallog << "locally relevant dofs :" << std::endl;
   owned_set.print(deallog.get_file_stream());
 
-  AffineConstraints<double> constraints(relevant_set);
+  AffineConstraints<double> constraints(owned_set, relevant_set);
   DoFTools::make_hanging_node_constraints(dof, constraints);
   VectorTools::interpolate_boundary_values(dof,
                                            0,
index 366dd2b6cb1947da63c3373927774a144f5890f5..738207ca2c6ed002909aa46b1f5f3321f5525011 100644 (file)
@@ -236,7 +236,8 @@ namespace Step50
     solution.reinit(mg_dof_handler.locally_owned_dofs(), MPI_COMM_WORLD);
     system_rhs.reinit(mg_dof_handler.locally_owned_dofs(), MPI_COMM_WORLD);
 
-    constraints.reinit(locally_relevant_set);
+    constraints.reinit(mg_dof_handler.locally_owned_dofs(),
+                       locally_relevant_set);
     DoFTools::make_hanging_node_constraints(mg_dof_handler, constraints);
 
     const std::set<types::boundary_id> dirichlet_boundary_ids = {0};
@@ -363,7 +364,8 @@ namespace Step50
       {
         const IndexSet dofset =
           DoFTools::extract_locally_relevant_level_dofs(mg_dof_handler, level);
-        boundary_constraints[level].reinit(dofset);
+        boundary_constraints[level].reinit(
+          mg_dof_handler.locally_owned_mg_dofs(level), dofset);
         boundary_constraints[level].add_lines(
           mg_constrained_dofs.get_refinement_edge_indices(level));
         boundary_constraints[level].add_lines(
index e87ea4e66257ad072e9e1827624554b1a238326d..c88ff86f053c9b76e160640d5742e72a8f9104ac 100644 (file)
@@ -54,7 +54,8 @@ test()
   force.reinit(locally_owned_dofs, mpi_communicator);
   Assert(!force.has_ghost_elements(), ExcInternalError());
 
-  AffineConstraints<PetscScalar> constraints(locally_relevant_dofs);
+  AffineConstraints<PetscScalar> constraints(locally_owned_dofs,
+                                             locally_relevant_dofs);
   constraints.clear();
   {
     const IndexSet boundary_dofs =
index bc9fd3acbb8619cb518b3e00d297df9d30c5d0b7..37b2e14e1f33fec7ae975134c095b3ec10fa0485 100644 (file)
@@ -88,7 +88,7 @@ test()
   TrilinosWrappers::MPI::Vector x_rel;
   x_rel.reinit(relevant_set, MPI_COMM_WORLD);
 
-  AffineConstraints<double> cm(relevant_set);
+  AffineConstraints<double> cm(owned_set, relevant_set);
   DoFTools::make_hanging_node_constraints(dofh, cm);
   ComponentMask velocity_mask(dim + 1, true);
 
index 3956d5acf50ded6848b8f59de74ad1d9daf7f1ea..42a8b3e0295532469c839a18a18b5d50731c2873 100644 (file)
@@ -58,7 +58,7 @@ check(parallel::distributed::Triangulation<dim> &tria)
   AffineConstraints<double> constraints;
 
   constraints.clear();
-  constraints.reinit(locally_relevant_dofs);
+  constraints.reinit(locally_owned_dofs, locally_relevant_dofs);
   DoFTools::make_hanging_node_constraints(dof_handler, constraints);
 
   for (unsigned int id = 0; id < 1; ++id)
index 46e856c8f203c396d8f1ca9e34453ace7d128a3e..5af3193f580e3dc53b1b7135e3a9e205b3ae66b4 100644 (file)
@@ -121,7 +121,7 @@ namespace pdd
 
     // Apply a constant function at the boundary
     constraints.clear();
-    constraints.reinit(locally_relevant_dofs);
+    constraints.reinit(locally_owned_dofs, locally_relevant_dofs);
     DoFTools::make_hanging_node_constraints(dof_handler, constraints);
     VectorTools::interpolate_boundary_values(
       dof_handler, 0, Functions::ConstantFunction<dim>(1.0), constraints);
index 6f5d98d1e8db4c1403d1c240629a2346b49e2fd7..9ea19099875c94759a313a53423115baef058a0b 100644 (file)
@@ -107,7 +107,7 @@ test(const unsigned int degree_center,
 
   AffineConstraints<double> constraints;
   constraints.clear();
-  constraints.reinit(locally_relevant_dofs);
+  constraints.reinit(dh.locally_owned_dofs(), locally_relevant_dofs);
 
   DoFTools::make_hanging_node_constraints(dh, constraints);
 
index f66b9a69818f4a233d2e602d03f9720d814c377f..69b54665ad153f782b58548d74cfd2dd23b5b8ff 100644 (file)
@@ -107,7 +107,7 @@ test(const unsigned int degree_center,
 
   AffineConstraints<double> constraints;
   constraints.clear();
-  constraints.reinit(locally_relevant_dofs);
+  constraints.reinit(dh.locally_owned_dofs(), locally_relevant_dofs);
 
   DoFTools::make_hanging_node_constraints(dh, constraints);
 
index 9f4d6703507396dbf890c92020c3ddc93b56aeef..fffb26091677dd63ad91d99b198fbd60f5fd603e 100644 (file)
@@ -122,7 +122,7 @@ test(const unsigned int degree_center,
 
   AffineConstraints<double> constraints;
   constraints.clear();
-  constraints.reinit(locally_relevant_dofs);
+  constraints.reinit(dh.locally_owned_dofs(), locally_relevant_dofs);
 
   DoFTools::make_hanging_node_constraints(dh, constraints);
 
index 2e3aabd9090e1adf63dd31ac9eef412c2ef6284a..3ae0ded75a59de2a0b2be0b9d5b47d1456752cbf 100644 (file)
@@ -77,7 +77,7 @@ test()
     DoFTools::extract_locally_relevant_dofs(dh);
 
   AffineConstraints<double> constraints;
-  constraints.reinit(locally_relevant_dofs);
+  constraints.reinit(dh.locally_owned_dofs(), locally_relevant_dofs);
   DoFTools::make_hanging_node_constraints(dh, constraints);
 
   deallog << "OK" << std::endl;
index 772e1b53b1ec0d7a39a88eb7b94827f09cbd8365..49ec94ae36ab14719f7fe83c675370384fc058d5 100644 (file)
@@ -97,7 +97,8 @@ test()
     DoFTools::extract_locally_relevant_dofs(dof_handler);
 
   AffineConstraints<double> hanging_node_constraints;
-  hanging_node_constraints.reinit(locally_relevant_dofs);
+  hanging_node_constraints.reinit(dof_handler.locally_owned_dofs(),
+                                  locally_relevant_dofs);
   DoFTools::make_hanging_node_constraints(dof_handler,
                                           hanging_node_constraints);
 
index 5725bb91f805b092e7839ae099ca5d7c23a5c279..74f667bebee381321c063fd787aa2ffc342b60b3 100644 (file)
@@ -191,7 +191,7 @@ namespace Step4
 
     constraints.clear();
 
-    constraints.reinit(locally_relevant_dofs);
+    constraints.reinit(locally_owned_dofs, locally_relevant_dofs);
 
     DoFTools::make_hanging_node_constraints(dof_handler, constraints);
 
index 2b0019cdfe54a523536c735a380fd207a98106a5..932d13c2eb6d386cc469db99c48d7280c4075c44 100644 (file)
@@ -146,7 +146,7 @@ namespace Step40
     system_rhs = PetscScalar();
 
     constraints.clear();
-    constraints.reinit(locally_relevant_dofs);
+    constraints.reinit(locally_owned_dofs, locally_relevant_dofs);
     DoFTools::make_hanging_node_constraints(dof_handler, constraints);
     VectorTools::interpolate_boundary_values(dof_handler,
                                              0,
index d145099bfbe9414ad7fecc2d4abf53b1dcd01150..1785714f382c02b4cde58a581859c771dcb02d58 100644 (file)
@@ -154,7 +154,7 @@ namespace Step40
     system_rhs = PetscScalar();
 
     constraints.clear();
-    constraints.reinit(locally_relevant_dofs);
+    constraints.reinit(locally_owned_dofs, locally_relevant_dofs);
     DoFTools::make_hanging_node_constraints(dof_handler, constraints);
     VectorTools::interpolate_boundary_values(
       dof_handler, 0, Functions::ZeroFunction<dim, PetscScalar>(), constraints);
index 59fb4e4c9c0522565391b81cbb11bf0273fa8f0f..391699fc710049cf3488e52c5b9d6103715df9d1 100644 (file)
@@ -140,7 +140,7 @@ SeventhProblem<dim>::setup_system()
   system_rhs.reinit(locally_owned_dofs, mpi_communicator);
   system_rhs = 0;
   constraints.clear();
-  constraints.reinit(locally_relevant_dofs);
+  constraints.reinit(locally_owned_dofs, locally_relevant_dofs);
   DoFTools::make_hanging_node_constraints(dof_handler, constraints);
   VectorTools::interpolate_boundary_values(dof_handler,
                                            0,
index 89fa94723623e4161c93263223733c05da181038..ec3f804e200c46e21fc8563f1f8714947b8abf08 100644 (file)
@@ -53,7 +53,8 @@ test()
   // locally relevant dofs, so do the test again with that
   {
     const IndexSet relevant_set = DoFTools::extract_locally_relevant_dofs(dofh);
-    AffineConstraints<double> boundary_values(relevant_set);
+    AffineConstraints<double> boundary_values(dofh.locally_owned_dofs(),
+                                              relevant_set);
     DoFTools::make_zero_boundary_constraints(dofh, boundary_values);
     if (Utilities::MPI::this_mpi_process(MPI_COMM_WORLD) == 0)
       boundary_values.print(deallog.get_file_stream());
index 4220b544c0d9bf628de9caefbe8e37df9e416a27..495264b003745caf7aea2c541f2f32e09e84d599 100644 (file)
@@ -106,7 +106,7 @@ test()
   const IndexSet relevant_set = DoFTools::extract_locally_relevant_dofs(dofh);
 
   AffineConstraints<double> constraints;
-  constraints.reinit(relevant_set);
+  constraints.reinit(dofh.locally_owned_dofs(), relevant_set);
   DoFTools::make_hanging_node_constraints(dofh, constraints);
   std::set<types::boundary_id> no_normal_flux_boundaries;
   no_normal_flux_boundaries.insert(1);
index db7060a00b2ddd86e3b10ccf3df4c91867e5b198..a89c7c7cf62fee9f86d30eae90e02be166accd52 100644 (file)
@@ -108,7 +108,7 @@ test()
   const IndexSet relevant_set = DoFTools::extract_locally_relevant_dofs(dofh);
 
   AffineConstraints<double> constraints;
-  constraints.reinit(relevant_set);
+  constraints.reinit(dofh.locally_owned_dofs(), relevant_set);
   DoFTools::make_hanging_node_constraints(dofh, constraints);
   std::set<types::boundary_id> no_normal_flux_boundaries;
   no_normal_flux_boundaries.insert(1);
index 374b6f75e2844fc8c1dc4329a5ea83cdb35a52ec..66a60d23fddd719ebca8acca0ead0a1f7f8ef89b 100644 (file)
@@ -84,7 +84,7 @@ test()
   const IndexSet relevant_set = DoFTools::extract_locally_relevant_dofs(dofh);
 
   AffineConstraints<double> constraints;
-  constraints.reinit(relevant_set);
+  constraints.reinit(dofh.locally_owned_dofs(), relevant_set);
   DoFTools::make_hanging_node_constraints(dofh, constraints);
   const std::set<types::boundary_id> no_normal_flux_boundaries = {0};
   const unsigned int                 degree                    = 1;
index 7054bc1113da79468a144d694112dcbaeab7596a..f6988f8c2dd63351a42be65d7c4ac837ab266367 100644 (file)
@@ -105,7 +105,7 @@ test()
   const IndexSet relevant_set = DoFTools::extract_locally_relevant_dofs(dofh);
 
   AffineConstraints<double> constraints;
-  constraints.reinit(relevant_set);
+  constraints.reinit(dofh.locally_owned_dofs(), relevant_set);
   DoFTools::make_hanging_node_constraints(dofh, constraints);
   std::set<types::boundary_id> no_normal_flux_boundaries;
   no_normal_flux_boundaries.insert(1);
index 07aedae22e857f0a37b80081d29f59088599ff2e..40a3d270c1d09b09880f928eac184b43ddb2ae1f 100644 (file)
@@ -62,7 +62,7 @@ test()
 
   AffineConstraints<double> cm;
   DoFTools::make_hanging_node_constraints(dofh, cm);
-  AffineConstraints<double> cm2(dof_set);
+  AffineConstraints<double> cm2(dofh.locally_owned_dofs(), dof_set);
   DoFTools::make_hanging_node_constraints(dofh, cm2);
 
   {
index 3a431fdbdbacb48ab57d85098d5bf2db2536d82a..d57dad20fb513d7585ece09adb01508116ab0d98 100644 (file)
@@ -64,7 +64,7 @@ test()
 
       AffineConstraints<double> cm;
       DoFTools::make_hanging_node_constraints(dofh, cm);
-      AffineConstraints<double> cm2(dof_set);
+      AffineConstraints<double> cm2(dofh.locally_owned_dofs(), dof_set);
       DoFTools::make_hanging_node_constraints(dofh, cm2);
 
       if (myid == 0)
index 4fb480f9862f58d7d3c9254cc9c79466c5ddaeed..ee298ef96766012602754e21bb3ee8fc6f2a33cc 100644 (file)
@@ -97,7 +97,7 @@ test(const bool use_manifold_for_normal, const Mapping<dim> &mapping)
   TrilinosWrappers::MPI::Vector x_rel;
   x_rel.reinit(relevant_set, MPI_COMM_WORLD);
 
-  AffineConstraints<double> cm(relevant_set);
+  AffineConstraints<double> cm(owned_set, relevant_set);
   DoFTools::make_hanging_node_constraints(dofh, cm);
   ComponentMask velocity_mask(dim + 1, true);
 
index 10cd178b85a116de75b43e4ae35d6e1f81a66f29..e0190d400875b63d0f6e87b9f653da72acb451d2 100644 (file)
@@ -204,7 +204,7 @@ test()
 
       x_rel.reinit(relevant_set, MPI_COMM_WORLD);
 
-      AffineConstraints<double> cm(relevant_set);
+      AffineConstraints<double> cm(owned_set, relevant_set);
       DoFTools::make_hanging_node_constraints(dofh, cm);
       /*  std::vector<bool> velocity_mask (dim+1, true);
 
index e135e2277090afb86a8d8c5d67e210d48d0a9b89..a6d9468c3c4789466fe23d1ab64b0408dffc461a 100644 (file)
@@ -68,7 +68,7 @@ test()
   // AffineConstraints<double>::distribute should not do
   // anything
   x1 = x_ref;
-  AffineConstraints<double> cm(relevant_set);
+  AffineConstraints<double> cm(owned_set, relevant_set);
   DoFTools::make_hanging_node_constraints(dofh, cm);
   cm.close();
   cm.distribute(x1);
index 17d544f779a1ac2d7d4c0df64a20b5b0f6958c86..0f2d193158e3329cb4e3a020192a4cd144fc6f8f 100644 (file)
@@ -72,7 +72,7 @@ test()
 
   AffineConstraints<double> cm;
   DoFTools::make_hanging_node_constraints(dofh, cm);
-  AffineConstraints<double> cm2(dof_set);
+  AffineConstraints<double> cm2(dofh.locally_owned_dofs(), dof_set);
   DoFTools::make_hanging_node_constraints(dofh, cm2);
 
   {
index 55776b50e1852caced76cd8ffacb66bbe9394ab6..a84622206b0cc957beb2dc35fbca54397665a012 100644 (file)
@@ -80,7 +80,7 @@ test()
   TrilinosWrappers::MPI::Vector x_rel;
   x_rel.reinit(relevant_set, MPI_COMM_WORLD);
 
-  AffineConstraints<double> cm(relevant_set);
+  AffineConstraints<double> cm(owned_set, relevant_set);
   DoFTools::make_hanging_node_constraints(dofh, cm);
   cm.close();
 
index 828623a90d19d5cf45c02b4451ef513cc412f3b0..08a2dfe71018406a23d5ffadf20416dcfdeaeb80 100644 (file)
@@ -204,7 +204,7 @@ test()
 
       x_rel.reinit(relevant_set, MPI_COMM_WORLD);
 
-      AffineConstraints<double> cm(relevant_set);
+      AffineConstraints<double> cm(owned_set, relevant_set);
       DoFTools::make_hanging_node_constraints(dofh, cm);
       /*  std::vector<bool> velocity_mask (dim+1, true);
 
index f9f653c0047594393cce576b06daedbd11332b9d..4edc7d8dc3f911e69fdd3114db23d1e15efbaa1c 100644 (file)
@@ -68,7 +68,8 @@ test()
   // build constraint matrix
   const IndexSet locally_relevant =
     DoFTools::extract_locally_relevant_dofs(dof);
-  AffineConstraints<double> constraints(locally_relevant);
+  AffineConstraints<double> constraints(dof.locally_owned_dofs(),
+                                        locally_relevant);
   DoFTools::make_hanging_node_constraints(dof, constraints);
   constraints.close();
 
index f61eeecc408c84ac21bf01bef3af17571a77f07b..12b0db0c91d29a2dec9771c4753cf1c7dff3eef2 100644 (file)
@@ -89,7 +89,7 @@ test(const unsigned int n_blocks = 5)
   const IndexSet &owned_set    = dof.locally_owned_dofs();
   const IndexSet  relevant_set = DoFTools::extract_locally_relevant_dofs(dof);
 
-  AffineConstraints<double> constraints(relevant_set);
+  AffineConstraints<double> constraints(owned_set, relevant_set);
   DoFTools::make_hanging_node_constraints(dof, constraints);
   VectorTools::interpolate_boundary_values(dof,
                                            0,
index 84fed2e3cc683d78fc3d41fbb1451fcf225841f8..9976cbab86d8ecef0e112fb44903379c20f4b307 100644 (file)
@@ -91,7 +91,7 @@ test(const unsigned int n_blocks = 5)
   const IndexSet &owned_set    = dof.locally_owned_dofs();
   const IndexSet  relevant_set = DoFTools::extract_locally_relevant_dofs(dof);
 
-  AffineConstraints<double> constraints(relevant_set);
+  AffineConstraints<double> constraints(owned_set, relevant_set);
   DoFTools::make_hanging_node_constraints(dof, constraints);
   VectorTools::interpolate_boundary_values(dof,
                                            0,
index e348579761210b037608610c6fcd2359962fe3e4..0b01b60c55519875210dab33981e89d2dd6f3a34 100644 (file)
@@ -89,7 +89,7 @@ test(const unsigned int n_blocks = 5)
   const IndexSet &owned_set    = dof.locally_owned_dofs();
   const IndexSet  relevant_set = DoFTools::extract_locally_relevant_dofs(dof);
 
-  AffineConstraints<double> constraints(relevant_set);
+  AffineConstraints<double> constraints(owned_set, relevant_set);
   DoFTools::make_hanging_node_constraints(dof, constraints);
   VectorTools::interpolate_boundary_values(dof,
                                            0,
index 1a9f0199904c9c4c88191f17b4c831afe151eb33..cb2f44df39911fec6fc3b99008b31cdac3771c0c 100644 (file)
@@ -89,7 +89,7 @@ test(const unsigned int n_blocks = 5)
   const IndexSet &owned_set    = dof.locally_owned_dofs();
   const IndexSet  relevant_set = DoFTools::extract_locally_relevant_dofs(dof);
 
-  AffineConstraints<double> constraints(relevant_set);
+  AffineConstraints<double> constraints(owned_set, relevant_set);
   DoFTools::make_hanging_node_constraints(dof, constraints);
   VectorTools::interpolate_boundary_values(dof,
                                            0,
index 28f04172a222c9d12babd3e5d3d20fe9f05ec29e..12104c5036aafa822e966069ff7ec60e60338f91 100644 (file)
@@ -92,7 +92,7 @@ test(const unsigned int n_blocks = 5)
   const IndexSet &owned_set    = dof.locally_owned_dofs();
   const IndexSet  relevant_set = DoFTools::extract_locally_relevant_dofs(dof);
 
-  AffineConstraints<double> constraints(relevant_set);
+  AffineConstraints<double> constraints(owned_set, relevant_set);
   DoFTools::make_hanging_node_constraints(dof, constraints);
   VectorTools::interpolate_boundary_values(dof,
                                            0,
index d0ab5186d19936e0e28fa2763fa43dca3fb04846..3f5421dbaedc6cfb2a268fcfc89855dcbcc84678 100644 (file)
@@ -88,7 +88,7 @@ test(const unsigned int n_blocks = 5)
   const IndexSet &owned_set    = dof.locally_owned_dofs();
   const IndexSet  relevant_set = DoFTools::extract_locally_relevant_dofs(dof);
 
-  AffineConstraints<double> constraints(relevant_set);
+  AffineConstraints<double> constraints(owned_set, relevant_set);
   DoFTools::make_hanging_node_constraints(dof, constraints);
   VectorTools::interpolate_boundary_values(dof,
                                            0,
index fbbe6966745dd2d3a97bb6159fc5ba2e07ffe638..5ed50a9113ec88a6265e5a152c6abef69c705e73 100644 (file)
@@ -90,7 +90,7 @@ test(const unsigned int n_blocks = 5)
   const IndexSet &owned_set    = dof.locally_owned_dofs();
   const IndexSet  relevant_set = DoFTools::extract_locally_relevant_dofs(dof);
 
-  AffineConstraints<double> constraints(relevant_set);
+  AffineConstraints<double> constraints(owned_set, relevant_set);
   DoFTools::make_hanging_node_constraints(dof, constraints);
   VectorTools::interpolate_boundary_values(dof,
                                            0,
index 16222bdd9b161be5f6b288f8c0d671db07797beb..9d2855facf1216a2b8a65c37f3a569b968363aff 100644 (file)
@@ -90,7 +90,7 @@ test(const unsigned int n = 5, const unsigned int m = 3)
   const IndexSet &owned_set    = dof.locally_owned_dofs();
   const IndexSet  relevant_set = DoFTools::extract_locally_relevant_dofs(dof);
 
-  AffineConstraints<double> constraints(relevant_set);
+  AffineConstraints<double> constraints(owned_set, relevant_set);
   DoFTools::make_hanging_node_constraints(dof, constraints);
   VectorTools::interpolate_boundary_values(dof,
                                            0,
index 8ef91fab051effd950677a5eb334e8ea2e86868d..dfcad39cd1846541c3cf0b5ec83705f8db0e71f4 100644 (file)
@@ -150,7 +150,7 @@ namespace Step40
 
     // Periodic Conditions
     constraints.clear();
-    constraints.reinit(locally_relevant_dofs);
+    constraints.reinit(locally_owned_dofs, locally_relevant_dofs);
     DoFTools::make_hanging_node_constraints(dof_handler, constraints);
     for (int i = 1; i < dim; ++i)
       DoFTools::make_periodicity_constraints(dof_handler,
index f2900e10601b8bbe64105c78d5bc9bf5a10e57ee..dfc3eebcce52aab5169fc9cd3a661f90a3c407a8 100644 (file)
@@ -340,7 +340,7 @@ namespace Step22
         locally_relevant_dofs.get_view(n_u, n_u + n_p));
 
       constraints.clear();
-      constraints.reinit(locally_relevant_dofs);
+      constraints.reinit(locally_owned_dofs, locally_relevant_dofs);
 
       const FEValuesExtractors::Vector velocities(0);
       const FEValuesExtractors::Scalar pressure(dim);
index 97145887fdbd1b5b8f562327a55452983675c4a2..cbcf7d1a6d483a4b0870a35bd792be63398b33cc 100644 (file)
@@ -267,7 +267,7 @@ namespace Step22
         locally_relevant_dofs.get_view(n_u, n_u + n_p));
 
       constraints.clear();
-      constraints.reinit(locally_relevant_dofs);
+      constraints.reinit(locally_owned_dofs, locally_relevant_dofs);
 
       const FEValuesExtractors::Vector velocities(0);
       const FEValuesExtractors::Scalar pressure(dim);
index 521844b0668820a598ef61ebc6139fe848efa63d..2115cfdc976ae122b133ce219a7a165d0e4b13ba 100644 (file)
@@ -188,7 +188,7 @@ check(const unsigned int orientation, bool reverse)
   const IndexSet  locally_relevant_dofs =
     DoFTools::extract_locally_relevant_dofs(dof_handler);
 
-  constraints.reinit(locally_relevant_dofs);
+  constraints.reinit(locally_owned_dofs, locally_relevant_dofs);
   {
     std::vector<
       GridTools::PeriodicFacePair<typename DoFHandler<dim>::cell_iterator>>
index 7c853c416b3f0cefe3242ffc0533d72ffa5802c7..2b1c9056a20d1f6f5c64e9e8f3e72fb3bf7f519c 100644 (file)
@@ -170,7 +170,7 @@ test(const unsigned numRefinementLevels = 2)
   /// creating combined hanging node and periodic constraint matrix
   AffineConstraints<double> constraints;
   constraints.clear();
-  constraints.reinit(locally_relevant_dofs);
+  constraints.reinit(dof_handler.locally_owned_dofs(), locally_relevant_dofs);
   DoFTools::make_hanging_node_constraints(dof_handler, constraints);
   std::vector<
     GridTools::PeriodicFacePair<typename DoFHandler<dim>::cell_iterator>>
index 58b9c0e1c9fc64fb1e6f2dd073cffea2f596f4c4..66e68dee4ef6122ac318f981c8451a81dfbebdf2 100644 (file)
@@ -148,7 +148,7 @@ test(const unsigned numRefinementLevels = 2)
   /// creating combined hanging node and periodic constraint matrix
   AffineConstraints<double> constraints;
   constraints.clear();
-  constraints.reinit(locally_relevant_dofs);
+  constraints.reinit(dof_handler.locally_owned_dofs(), locally_relevant_dofs);
 
   std::vector<
     GridTools::PeriodicFacePair<typename DoFHandler<dim>::cell_iterator>>
index 80e0acfc065ade8a9b25c5171a483d3ae9131f1a..58cf4fbe6d9e1854c6da4016e537b1594732efb2 100644 (file)
@@ -142,7 +142,7 @@ test()
   /// creating combined hanging node and periodic constraint matrix
   AffineConstraints<double> constraints;
   constraints.clear();
-  constraints.reinit(locally_relevant_dofs);
+  constraints.reinit(dof_handler.locally_owned_dofs(), locally_relevant_dofs);
 
   for (int d = 0; d < dim; ++d)
     DoFTools::make_periodicity_constraints(
index fffc5c50c529846002d2a7276b309787dce49d6e..da2e155748b3057d1c1da5b116b5a59ea6595d7c 100644 (file)
@@ -66,7 +66,8 @@ test()
     std::max<int>(100 * myid - 50, 0),
     std::min(static_cast<types::global_dof_index>(100 * myid + 150),
              vec.size()));
-  AffineConstraints<PetscScalar> cm(locally_relevant_range);
+  AffineConstraints<PetscScalar> cm(vec.locally_owned_elements(),
+                                    locally_relevant_range);
 
   // add constraints that constrain an element in the middle of the
   // local range of each processor against an element outside, both in
index 849762b943e5c00fd4e51d36a5be633f3523d122..231d1a5a6f0da9a0938d0d673acdb2f9a93098a0 100644 (file)
@@ -58,7 +58,8 @@ test()
     std::max<int>(100 * myid - 50, 0),
     std::min(static_cast<types::global_dof_index>(100 * myid + 150),
              vec.size()));
-  AffineConstraints<PetscScalar> cm(locally_relevant_range);
+  AffineConstraints<PetscScalar> cm(vec.locally_owned_elements(),
+                                    locally_relevant_range);
 
   // add constraints that constrain an element in the middle of the
   // local range of each processor against an element outside, both in
index 7fef9f96dc393a02d6672402e565f7464b85a4ed..c9571e51b462ee631c9800ec367f78fc4cc4c4a4 100644 (file)
@@ -216,7 +216,7 @@ namespace Step27
     system_rhs.reinit(locally_owned_dofs, mpi_communicator);
 
     constraints.clear();
-    constraints.reinit(locally_relevant_dofs);
+    constraints.reinit(locally_owned_dofs, locally_relevant_dofs);
     DoFTools::make_hanging_node_constraints(dof_handler, constraints);
     VectorTools::interpolate_boundary_values(dof_handler,
                                              0,
index 140173515019d10848bf15e7d9bb491f3e74b2b1..e2f9cff9ccc11024cad12c9c6bda701164c5155e 100644 (file)
@@ -197,7 +197,8 @@ namespace ResFlow
           << std::endl;
 
     {
-      constraints.reinit(locally_relevant_dofs);
+      constraints.reinit(dof_handler.locally_owned_dofs(),
+                         locally_relevant_dofs);
 
       const FEValuesExtractors::Vector velocities(0);
       DoFTools::make_hanging_node_constraints(dof_handler, constraints);
index cbe105a7c92ead5840e869393b27cdec925aa367..4bc0513ee0b28b93e73d6e09d40d5c8a52828aff 100644 (file)
@@ -216,7 +216,7 @@ namespace Step37
       DoFTools::extract_locally_relevant_dofs(dof_handler);
 
     constraints.clear();
-    constraints.reinit(locally_relevant_dofs);
+    constraints.reinit(dof_handler.locally_owned_dofs(), locally_relevant_dofs);
     DoFTools::make_hanging_node_constraints(dof_handler, constraints);
     VectorTools::interpolate_boundary_values(dof_handler,
                                              0,
@@ -258,7 +258,8 @@ namespace Step37
         const IndexSet relevant_dofs =
           DoFTools::extract_locally_relevant_level_dofs(dof_handler, level);
         AffineConstraints<double> level_constraints;
-        level_constraints.reinit(relevant_dofs);
+        level_constraints.reinit(dof_handler.locally_owned_mg_dofs(level),
+                                 relevant_dofs);
         level_constraints.add_lines(
           mg_constrained_dofs.get_boundary_indices(level));
         level_constraints.close();
@@ -322,8 +323,10 @@ namespace Step37
     non_homogeneous_constraints.clear();
     {
       AffineConstraints<double> hanging_nodes_laplace_constraints;
-      hanging_nodes_laplace_constraints.reinit(locally_relevant_dofs);
-      non_homogeneous_constraints.reinit(locally_relevant_dofs);
+      hanging_nodes_laplace_constraints.reinit(dof_handler.locally_owned_dofs(),
+                                               locally_relevant_dofs);
+      non_homogeneous_constraints.reinit(dof_handler.locally_owned_dofs(),
+                                         locally_relevant_dofs);
       DoFTools::make_hanging_node_constraints(
         dof_handler, hanging_nodes_laplace_constraints);
 
index 0ef74d70dc6c59de3a4bc7fc50d8a3fe8af67d9d..2c1db7ab070f074d7472ee5b079f6d6a4de2f1d5 100644 (file)
@@ -140,7 +140,7 @@ namespace Step40
     system_rhs = PetscScalar();
 
     constraints.clear();
-    constraints.reinit(locally_relevant_dofs);
+    constraints.reinit(locally_owned_dofs, locally_relevant_dofs);
     DoFTools::make_hanging_node_constraints(dof_handler, constraints);
     VectorTools::interpolate_boundary_values(dof_handler,
                                              0,
index 4d3c4540dbcc7a0de16364d9f411a35319b64bce..642b2cfaa55ee9e4ea752220beca391734823aa9 100644 (file)
@@ -202,7 +202,7 @@ namespace Step40
     system_rhs = PetscScalar();
 
     constraints.clear();
-    constraints.reinit(locally_relevant_dofs);
+    constraints.reinit(locally_owned_dofs, locally_relevant_dofs);
     DoFTools::make_hanging_node_constraints(dof_handler, constraints);
     VectorTools::interpolate_boundary_values(dof_handler,
                                              0,
index e77ecfbcf6a8b705e26783a0a64bed5c598e748c..46f9006da98928d7579dac1c4ae9a755c8798443 100644 (file)
@@ -203,7 +203,7 @@ namespace Step40
     system_rhs = PetscScalar();
 
     constraints.clear();
-    constraints.reinit(locally_relevant_dofs);
+    constraints.reinit(locally_owned_dofs, locally_relevant_dofs);
     DoFTools::make_hanging_node_constraints(dof_handler, constraints);
     VectorTools::interpolate_boundary_values(dof_handler,
                                              0,
index bef1f12ee1de0f572e4d99010c3b73e28f6f323b..3a692f897c8d508bcbd2bada4da425c26c7629ed 100644 (file)
@@ -140,7 +140,7 @@ namespace Step40
     system_rhs = PetscScalar();
 
     constraints.clear();
-    constraints.reinit(locally_relevant_dofs);
+    constraints.reinit(locally_owned_dofs, locally_relevant_dofs);
     DoFTools::make_hanging_node_constraints(dof_handler, constraints);
     VectorTools::interpolate_boundary_values(dof_handler,
                                              0,
index 4a689bb7d73930ee826df138322d19df74bb3619..90cd9d463b952da9528c432de82041c615011f5b 100644 (file)
@@ -71,7 +71,8 @@ test()
     std::max<int>(100 * myid - 50, 0),
     std::min(static_cast<types::global_dof_index>(100 * myid + 150),
              vec.size()));
-  AffineConstraints<double> cm(locally_relevant_range);
+  AffineConstraints<double> cm(vec.locally_owned_elements(),
+                               locally_relevant_range);
 
   // add constraints that constrain an element in the middle of the
   // local range of each processor against an element outside, both in
index 53afe910d0500eaa79a6d4e820d614a1931ad25d..3f7978cd69b62e10e62104693d46b78b6ba5c143 100644 (file)
@@ -63,7 +63,8 @@ test()
     std::max<int>(100 * myid - 50, 0),
     std::min(static_cast<types::global_dof_index>(100 * myid + 150),
              vec.size()));
-  AffineConstraints<double> cm(locally_relevant_range);
+  AffineConstraints<double> cm(vec.locally_owned_elements(),
+                               locally_relevant_range);
 
   // add constraints that constrain an element in the middle of the
   // local range of each processor against an element outside, both in
index 18e6264cbcdea6d3c2eab229c332574974aaeb21..926b5ab4ec552247f0c74100768e0ff817819198 100644 (file)
@@ -91,7 +91,7 @@ test()
     DoFTools::extract_locally_relevant_dofs(dofh);
 
   AffineConstraints<double> cm;
-  cm.reinit(locally_relevant_set);
+  cm.reinit(dofh.locally_owned_dofs(), locally_relevant_set);
   DoFTools::make_hanging_node_constraints(dofh, cm);
 
   std::map<types::boundary_id, const Function<dim> *> dirichlet_boundary;
index 70f54e8e2afec3b0a12f78fc882578d38e1489d7..87f9c3efab48f34a3a9558084da6102bae46ab7b 100644 (file)
@@ -220,7 +220,7 @@ namespace Step27
     system_rhs.reinit(locally_owned_dofs, mpi_communicator);
 
     constraints.clear();
-    constraints.reinit(locally_relevant_dofs);
+    constraints.reinit(locally_owned_dofs, locally_relevant_dofs);
     DoFTools::make_hanging_node_constraints(dof_handler, constraints);
     VectorTools::interpolate_boundary_values(dof_handler,
                                              0,
index ed47e797eed446062fb8ece185167ce1b355a916..d1515096d686087db81dfa9f7a554ac9fe4334cf 100644 (file)
@@ -89,7 +89,8 @@ test(const unsigned int n_refinements,
       // set up constraints
       const IndexSet locally_relevant_dofs =
         DoFTools::extract_locally_relevant_dofs(dof_handler);
-      constraint.reinit(locally_relevant_dofs);
+      constraint.reinit(dof_handler.locally_owned_dofs(),
+                        locally_relevant_dofs);
       VectorTools::interpolate_boundary_values(
         *mapping, dof_handler, 0, Functions::ZeroFunction<dim>(), constraint);
       constraint.close();
index febc6bcb413e51c95264fce2742cc0006d7a739a..de6bc7fe68a8735b81b1905464a987d156243775 100644 (file)
@@ -131,7 +131,7 @@ test(const unsigned int n_glob_ref = 2, const unsigned int n_ref = 0)
 
   // constraints:
   AffineConstraints<double> constraints;
-  constraints.reinit(locally_relevant_dofs);
+  constraints.reinit(locally_owned_dofs, locally_relevant_dofs);
   DoFTools::make_hanging_node_constraints(dof_handler, constraints);
   constraints.close();
 
index c13e8fae2d6085f7f6d329a80d3dff32a3a6a74c..542b2f6d90d15393dba246a074b0a333633daca8 100644 (file)
@@ -85,7 +85,7 @@ test()
     DoFTools::extract_locally_relevant_dofs(dof_handler);
 
   AffineConstraints<double> constraints;
-  constraints.reinit(locally_relevant_dofs);
+  constraints.reinit(dof_handler.locally_owned_dofs(), locally_relevant_dofs);
   DoFTools::make_hanging_node_constraints(dof_handler, constraints);
 
   std::vector<
index b8b8991feb5bd003a1d96702cc9a7fb375493c78..b69ac00a284dfc82ed670ffbe6a27475cb7c7ce9 100644 (file)
@@ -83,7 +83,8 @@ test(const unsigned int n_refinements,
       // set up constraints
       const IndexSet locally_relevant_dofs =
         DoFTools::extract_locally_relevant_dofs(dof_handler);
-      constraint.reinit(locally_relevant_dofs);
+      constraint.reinit(dof_handler.locally_owned_dofs(),
+                        locally_relevant_dofs);
       VectorTools::interpolate_boundary_values(
         *mapping, dof_handler, 0, Functions::ZeroFunction<dim>(), constraint);
       constraint.close();
index 695a486506f5afc115fa5d22ccc373272cce98eb..dac418b83b8089885b94c2d6604a44bc9d396f87 100644 (file)
@@ -68,7 +68,7 @@ test(const unsigned int n_refinements, const unsigned int fe_degree_fine)
       // set up constraints
       const IndexSet relevant_dofs =
         DoFTools::extract_locally_relevant_level_dofs(dof_handler, l);
-      constraint.reinit(relevant_dofs);
+      constraint.reinit(dof_handler.locally_owned_mg_dofs(l), relevant_dofs);
       constraint.add_lines(mg_constrained_dofs.get_boundary_indices(l));
       constraint.close();
 
index 4d749aa1a81c0232bde7877b5e9124cc677473a7..90daf33d4774724947e66bf70c3eab5f4fe60a83 100644 (file)
@@ -66,7 +66,7 @@ test(const unsigned int n_refinements, const unsigned int fe_degree_fine)
       // set up constraints
       const IndexSet relevant_dofs =
         DoFTools::extract_locally_relevant_level_dofs(dof_handler, l);
-      constraint.reinit(relevant_dofs);
+      constraint.reinit(dof_handler.locally_owned_mg_dofs(l), relevant_dofs);
       constraint.add_lines(mg_constrained_dofs.get_boundary_indices(l));
       constraint.close();
 
index 6a67d9115833fa0cb49fc51e5f8e8f28b8ac1814..22cea00ec0edb2478014ae4d35a8cf09792da1e4 100644 (file)
@@ -110,7 +110,8 @@ test(const unsigned int n_refinements,
       // set up constraints
       const IndexSet locally_relevant_dofs =
         DoFTools::extract_locally_relevant_dofs(dof_handler);
-      constraint.reinit(locally_relevant_dofs);
+      constraint.reinit(dof_handler.locally_owned_dofs(),
+                        locally_relevant_dofs);
       VectorTools::interpolate_boundary_values(
         *mapping, dof_handler, 0, Functions::ZeroFunction<dim>(), constraint);
       DoFTools::make_hanging_node_constraints(dof_handler, constraint);
index 3ee24f6575d588253c3a8e95af2d726cf68f334e..b0408f9ea1bd463a982d173dae3ad82529a6d8e7 100644 (file)
@@ -80,7 +80,7 @@ test(const unsigned int n_refinements, const unsigned int fe_degree_fine)
 
       const IndexSet relevant_dofs =
         DoFTools::extract_locally_relevant_level_dofs(dof_handler, 0 /*level*/);
-      constraint.reinit(relevant_dofs);
+      constraint.reinit(dof_handler.locally_owned_mg_dofs(0), relevant_dofs);
       constraint.add_lines(
         mg_constrained_dofs.get_boundary_indices(0 /*level*/));
       constraint.close();
index daa1a3260c5ee25facdbfbaef13709a19a4cfdfb..e377c8b4c9c2d223ace33fc31105cc11a0b7e214 100644 (file)
@@ -83,7 +83,8 @@ test(const unsigned int n_refinements,
 
       const IndexSet relevant_dofs =
         DoFTools::extract_locally_relevant_level_dofs(dof_handler, level);
-      constraint.reinit(relevant_dofs);
+      constraint.reinit(dof_handler.locally_owned_mg_dofs(level),
+                        relevant_dofs);
       constraint.add_lines(mg_constrained_dofs.get_boundary_indices(level));
       constraint.close();
 
index 47400165f7821f9587ac9f2acb451c7b5e31d15b..f4422a04664363142d5a6f888d6d4c899eb17d8d 100644 (file)
@@ -87,7 +87,8 @@ test(const unsigned int n_refinements,
   AffineConstraints<Number> fine_constraints;
   const IndexSet            locally_relevant_dofs =
     DoFTools::extract_locally_relevant_dofs(fine_dof_handler);
-  fine_constraints.reinit(locally_relevant_dofs);
+  fine_constraints.reinit(fine_dof_handler.locally_owned_dofs(),
+                          locally_relevant_dofs);
   VectorTools::interpolate_boundary_values(*mapping,
                                            fine_dof_handler,
                                            0,
@@ -152,7 +153,8 @@ test(const unsigned int n_refinements,
       // set up constraints
       const IndexSet locally_relevant_dofs =
         DoFTools::extract_locally_relevant_dofs(dof_handler);
-      constraint.reinit(locally_relevant_dofs);
+      constraint.reinit(dof_handler.locally_owned_dofs(),
+                        locally_relevant_dofs);
       VectorTools::interpolate_boundary_values(
         *mapping, dof_handler, 0, Functions::ZeroFunction<dim>(), constraint);
       DoFTools::make_hanging_node_constraints(dof_handler, constraint);
index 15b19a3edb3e2293a46251685cf218c834690478..dce4eb94b4424c63ce76722227a45a891109b5fc 100644 (file)
@@ -83,7 +83,8 @@ test(const unsigned int n_refinements,
       // set up constraints
       const IndexSet locally_relevant_dofs =
         DoFTools::extract_locally_relevant_dofs(dof_handler);
-      constraint.reinit(locally_relevant_dofs);
+      constraint.reinit(dof_handler.locally_owned_dofs(),
+                        locally_relevant_dofs);
       VectorTools::interpolate_boundary_values(
         *mapping, dof_handler, 0, Functions::ZeroFunction<dim>(), constraint);
       constraint.close();
index 4d64b9f27623cb5cc8ddc69a9201c0b40efdc4b5..11bfdcfa35c91a70817cf12f26eea01b392fbe58 100644 (file)
@@ -74,7 +74,8 @@ test(const unsigned int n_refinements,
       // set up constraints
       const IndexSet locally_relevant_dofs =
         DoFTools::extract_locally_relevant_dofs(dof_handler);
-      constraint.reinit(locally_relevant_dofs);
+      constraint.reinit(dof_handler.locally_owned_dofs(),
+                        locally_relevant_dofs);
       VectorTools::interpolate_boundary_values(
         *_mapping, dof_handler, 0, Functions::ZeroFunction<dim>(), constraint);
       constraint.close();
index 89d39e081a9bed5ed0ef6cf48d5f478ab70baf0d..6cd9ce42db8f59714424ae20162d6129218f0494 100644 (file)
@@ -65,7 +65,8 @@ test(const unsigned int n_refinements, const unsigned int fe_degree_fine)
       // set up constraints
       const IndexSet locally_relevant_dofs =
         DoFTools::extract_locally_relevant_dofs(dof_handler);
-      constraint.reinit(locally_relevant_dofs);
+      constraint.reinit(dof_handler.locally_owned_dofs(),
+                        locally_relevant_dofs);
       VectorTools::interpolate_boundary_values(
         mapping, dof_handler, 0, Functions::ZeroFunction<dim>(), constraint);
       constraint.close();
index bf28c8f8ee860d05e4d52f0f24789b3a2a22be02..99199e696dd3973a25b2477a7431711f5985647d 100644 (file)
@@ -140,7 +140,8 @@ test(const unsigned int n_refinements,
       // set up constraints
       const IndexSet locally_relevant_dofs =
         DoFTools::extract_locally_relevant_dofs(dof_handler);
-      constraint.reinit(locally_relevant_dofs);
+      constraint.reinit(dof_handler.locally_owned_dofs(),
+                        locally_relevant_dofs);
       VectorTools::interpolate_boundary_values(
         *mapping, dof_handler, 0, Functions::ZeroFunction<dim>(), constraint);
       DoFTools::make_hanging_node_constraints(dof_handler, constraint);
index 45660a5963ef3610b6254d17c7930a16ccb554f5..3616c68bd16bd5b8b46a023acbd19771b784319c 100644 (file)
@@ -79,10 +79,10 @@ test(const unsigned int n_refinements,
       dof_handler.distribute_dofs(*fe);
 
       // set up constraints
-      IndexSet locally_relevant_dofs;
-      DoFTools::extract_locally_relevant_dofs(dof_handler,
-                                              locally_relevant_dofs);
-      constraint.reinit(locally_relevant_dofs);
+      const IndexSet locally_relevant_dofs =
+        DoFTools::extract_locally_relevant_dofs(dof_handler);
+      constraint.reinit(dof_handler.locally_owned_dofs(),
+                        locally_relevant_dofs);
       VectorTools::interpolate_boundary_values(*_mapping,
                                                dof_handler,
                                                0,
index a87e4ba29f80ef0c3c94175358d0c45c05f45bf4..eb82a745a7ad1350d7df0c9d2eb363be22883feb 100644 (file)
@@ -107,7 +107,7 @@ do_test(const DoFHandler<dim> &dof)
 
   // fine-level constraints
   AffineConstraints<double> constraints;
-  constraints.reinit(locally_relevant_dofs);
+  constraints.reinit(dof.locally_owned_dofs(), locally_relevant_dofs);
   DoFTools::make_hanging_node_constraints(dof, constraints);
   VectorTools::interpolate_boundary_values(dof,
                                            dirichlet_boundary,
@@ -174,7 +174,7 @@ do_test(const DoFHandler<dim> &dof)
       AffineConstraints<double> level_constraints;
       const IndexSet            relevant_dofs =
         DoFTools::extract_locally_relevant_level_dofs(dof, level);
-      level_constraints.reinit(relevant_dofs);
+      level_constraints.reinit(dof.locally_owned_mg_dofs(level), relevant_dofs);
       level_constraints.add_lines(
         mg_constrained_dofs.get_boundary_indices(level));
       level_constraints.close();
index 10e6ba46c154f00afb64ece9d5b8272d11e98166..84a5be6cd0dea25b23c98648a7ef61cd334ae5cf 100644 (file)
@@ -79,7 +79,7 @@ public:
       {
         const IndexSet relevant_dofs =
           DoFTools::extract_locally_relevant_dofs(dof_handler);
-        constraints.reinit(relevant_dofs);
+        constraints.reinit(dof_handler.locally_owned_dofs(), relevant_dofs);
         DoFTools::make_hanging_node_constraints(dof_handler, constraints);
         VectorTools::interpolate_boundary_values(dof_handler,
                                                  dirichlet_boundary,
@@ -87,9 +87,11 @@ public:
       }
     else
       {
+        const IndexSet &locally_owned =
+          dof_handler.locally_owned_mg_dofs(level);
         const IndexSet relevant_dofs =
           DoFTools::extract_locally_relevant_level_dofs(dof_handler, level);
-        constraints.reinit(relevant_dofs);
+        constraints.reinit(locally_owned, relevant_dofs);
         constraints.add_lines(mg_constrained_dofs.get_boundary_indices(level));
 
         const std::vector<types::global_dof_index> interface_indices =
@@ -98,8 +100,6 @@ public:
         edge_constrained_indices.clear();
         edge_constrained_indices.reserve(interface_indices.size());
         edge_constrained_values.resize(interface_indices.size());
-        const IndexSet &locally_owned =
-          dof_handler.locally_owned_mg_dofs(level);
         for (unsigned int i = 0; i < interface_indices.size(); ++i)
           if (locally_owned.is_element(interface_indices[i]))
             edge_constrained_indices.push_back(
@@ -478,7 +478,8 @@ do_test(const DoFHandler<dim> &dof)
   AffineConstraints<double> hanging_node_constraints;
   const IndexSet            locally_relevant_dofs =
     DoFTools::extract_locally_relevant_dofs(dof);
-  hanging_node_constraints.reinit(locally_relevant_dofs);
+  hanging_node_constraints.reinit(dof.locally_owned_dofs(),
+                                  locally_relevant_dofs);
   DoFTools::make_hanging_node_constraints(dof, hanging_node_constraints);
   hanging_node_constraints.close();
 
index fa2f36cf9af2d53089e633655e8ff6ed95f46114..e27c3e23c9ed2dd1590b675108f267f7d5e79cf5 100644 (file)
@@ -108,7 +108,7 @@ do_test(const DoFHandler<dim> &dof)
 
   // fine-level constraints
   AffineConstraints<double> constraints;
-  constraints.reinit(locally_relevant_dofs);
+  constraints.reinit(dof.locally_owned_dofs(), locally_relevant_dofs);
   DoFTools::make_hanging_node_constraints(dof, constraints);
   VectorTools::interpolate_boundary_values(dof,
                                            dirichlet_boundary,
@@ -153,7 +153,8 @@ do_test(const DoFHandler<dim> &dof)
   {
     // this is to make it consistent with parallel_multigrid_adaptive.cc
     AffineConstraints<double> hanging_node_constraints;
-    hanging_node_constraints.reinit(locally_relevant_dofs);
+    hanging_node_constraints.reinit(dof.locally_owned_dofs(),
+                                    locally_relevant_dofs);
     DoFTools::make_hanging_node_constraints(dof, hanging_node_constraints);
     hanging_node_constraints.close();
 
@@ -188,7 +189,7 @@ do_test(const DoFHandler<dim> &dof)
       AffineConstraints<double> level_constraints;
       const IndexSet            relevant_dofs =
         DoFTools::extract_locally_relevant_level_dofs(dof, level);
-      level_constraints.reinit(relevant_dofs);
+      level_constraints.reinit(dof.locally_owned_mg_dofs(level), relevant_dofs);
       level_constraints.add_lines(
         mg_constrained_dofs.get_boundary_indices(level));
       level_constraints.close();
index 9fbff8713a8b13f50618a7573cbc58241f92f3e2..dba6722a676ad3548f1128ebf13e20a48c2c5b42 100644 (file)
@@ -86,7 +86,7 @@ public:
       {
         const IndexSet relevant_dofs =
           DoFTools::extract_locally_relevant_dofs(dof_handler);
-        constraints.reinit(relevant_dofs);
+        constraints.reinit(dof_handler.locally_owned_dofs(), relevant_dofs);
         DoFTools::make_hanging_node_constraints(dof_handler, constraints);
         VectorTools::interpolate_boundary_values(dof_handler,
                                                  dirichlet_boundary,
@@ -94,9 +94,11 @@ public:
       }
     else
       {
+        const IndexSet &locally_owned =
+          dof_handler.locally_owned_mg_dofs(level);
         const IndexSet relevant_dofs =
           DoFTools::extract_locally_relevant_level_dofs(dof_handler, level);
-        constraints.reinit(relevant_dofs);
+        constraints.reinit(locally_owned, relevant_dofs);
         constraints.add_lines(mg_constrained_dofs.get_boundary_indices(level));
 
         const std::vector<types::global_dof_index> interface_indices =
@@ -105,8 +107,6 @@ public:
         edge_constrained_indices.clear();
         edge_constrained_indices.reserve(interface_indices.size());
         edge_constrained_values.resize(interface_indices.size());
-        const IndexSet &locally_owned =
-          dof_handler.locally_owned_mg_dofs(level);
         for (unsigned int i = 0; i < interface_indices.size(); ++i)
           if (locally_owned.is_element(interface_indices[i]))
             edge_constrained_indices.push_back(
@@ -482,7 +482,8 @@ do_test(const DoFHandler<dim> &dof, const bool threaded)
   AffineConstraints<double> hanging_node_constraints;
   const IndexSet            locally_relevant_dofs =
     DoFTools::extract_locally_relevant_dofs(dof);
-  hanging_node_constraints.reinit(locally_relevant_dofs);
+  hanging_node_constraints.reinit(dof.locally_owned_dofs(),
+                                  locally_relevant_dofs);
   DoFTools::make_hanging_node_constraints(dof, hanging_node_constraints);
   hanging_node_constraints.close();
 
index 4d3d0e684f83df0a7d51f37e3a6a158715e8f0bd..28995351bd151415b870de78096e6b22f65edb08 100644 (file)
@@ -108,7 +108,7 @@ do_test(const DoFHandler<dim> &dof)
 
   // fine-level constraints
   AffineConstraints<double> constraints;
-  constraints.reinit(locally_relevant_dofs);
+  constraints.reinit(dof.locally_owned_dofs(), locally_relevant_dofs);
   DoFTools::make_hanging_node_constraints(dof, constraints);
   VectorTools::interpolate_boundary_values(dof,
                                            dirichlet_boundary,
@@ -153,7 +153,8 @@ do_test(const DoFHandler<dim> &dof)
   {
     // this is to make it consistent with parallel_multigrid_adaptive.cc
     AffineConstraints<double> hanging_node_constraints;
-    hanging_node_constraints.reinit(locally_relevant_dofs);
+    hanging_node_constraints.reinit(dof.locally_owned_dofs(),
+                                    locally_relevant_dofs);
     DoFTools::make_hanging_node_constraints(dof, hanging_node_constraints);
     hanging_node_constraints.close();
 
@@ -188,7 +189,7 @@ do_test(const DoFHandler<dim> &dof)
       AffineConstraints<double> level_constraints;
       const IndexSet            relevant_dofs =
         DoFTools::extract_locally_relevant_level_dofs(dof, level);
-      level_constraints.reinit(relevant_dofs);
+      level_constraints.reinit(dof.locally_owned_mg_dofs(level), relevant_dofs);
       level_constraints.add_lines(
         mg_constrained_dofs.get_boundary_indices(level));
       level_constraints.close();
index b3d63db9f019b4133cbb497a15ddb3803365a46b..7629f63fd84a33532ea17e51623ca468a5055716 100644 (file)
@@ -83,7 +83,7 @@ public:
       {
         const IndexSet relevant_dofs =
           DoFTools::extract_locally_relevant_dofs(dof_handler);
-        constraints.reinit(relevant_dofs);
+        constraints.reinit(dof_handler.locally_owned_dofs(), relevant_dofs);
         DoFTools::make_hanging_node_constraints(dof_handler, constraints);
         VectorTools::interpolate_boundary_values(dof_handler,
                                                  dirichlet_boundary,
@@ -91,9 +91,11 @@ public:
       }
     else
       {
+        const IndexSet &locally_owned =
+          dof_handler.locally_owned_mg_dofs(level);
         const IndexSet relevant_dofs =
           DoFTools::extract_locally_relevant_level_dofs(dof_handler, level);
-        constraints.reinit(relevant_dofs);
+        constraints.reinit(locally_owned, relevant_dofs);
         constraints.add_lines(mg_constrained_dofs.get_boundary_indices(level));
 
         const std::vector<types::global_dof_index> interface_indices =
@@ -102,8 +104,6 @@ public:
         edge_constrained_indices.clear();
         edge_constrained_indices.reserve(interface_indices.size());
         edge_constrained_values.resize(interface_indices.size());
-        const IndexSet &locally_owned =
-          dof_handler.locally_owned_mg_dofs(level);
         for (unsigned int i = 0; i < interface_indices.size(); ++i)
           if (locally_owned.is_element(interface_indices[i]))
             edge_constrained_indices.push_back(
@@ -479,7 +479,8 @@ do_test(const DoFHandler<dim> &dof, const bool threaded)
   AffineConstraints<double> hanging_node_constraints;
   const IndexSet            locally_relevant_dofs =
     DoFTools::extract_locally_relevant_dofs(dof);
-  hanging_node_constraints.reinit(locally_relevant_dofs);
+  hanging_node_constraints.reinit(dof.locally_owned_dofs(),
+                                  locally_relevant_dofs);
   DoFTools::make_hanging_node_constraints(dof, hanging_node_constraints);
   hanging_node_constraints.close();
 
index c5c587d52632822760a929a07e3c2c5dad327838..7be148f0507552ce8d9d2584113ec824c932cc15 100644 (file)
@@ -230,7 +230,7 @@ do_test(const DoFHandler<dim> &dof, const unsigned int nb)
 
   // fine-level constraints
   AffineConstraints<double> constraints;
-  constraints.reinit(locally_relevant_dofs);
+  constraints.reinit(dof.locally_owned_dofs(), locally_relevant_dofs);
   DoFTools::make_hanging_node_constraints(dof, constraints);
   VectorTools::interpolate_boundary_values(dof,
                                            dirichlet_boundary,
@@ -279,7 +279,8 @@ do_test(const DoFHandler<dim> &dof, const unsigned int nb)
   {
     // this is to make it consistent with parallel_multigrid_adaptive.cc
     AffineConstraints<double> hanging_node_constraints;
-    hanging_node_constraints.reinit(locally_relevant_dofs);
+    hanging_node_constraints.reinit(dof.locally_owned_dofs(),
+                                    locally_relevant_dofs);
     DoFTools::make_hanging_node_constraints(dof, hanging_node_constraints);
     hanging_node_constraints.close();
 
@@ -316,7 +317,7 @@ do_test(const DoFHandler<dim> &dof, const unsigned int nb)
       AffineConstraints<double> level_constraints;
       const IndexSet            relevant_dofs =
         DoFTools::extract_locally_relevant_level_dofs(dof, level);
-      level_constraints.reinit(relevant_dofs);
+      level_constraints.reinit(dof.locally_owned_mg_dofs(level), relevant_dofs);
       level_constraints.add_lines(
         mg_constrained_dofs.get_boundary_indices(level));
       level_constraints.close();
index 5ab225509195e84c4610c3af6e8b10c5be530280..c1726e71df5e7231c733c14132a1e3459f604da4 100644 (file)
@@ -109,7 +109,7 @@ do_test(const DoFHandler<dim> &dof)
 
   // fine-level constraints
   AffineConstraints<double> constraints;
-  constraints.reinit(locally_relevant_dofs);
+  constraints.reinit(dof.locally_owned_dofs(), locally_relevant_dofs);
   DoFTools::make_hanging_node_constraints(dof, constraints);
   VectorTools::interpolate_boundary_values(dof,
                                            dirichlet_boundary,
@@ -154,7 +154,8 @@ do_test(const DoFHandler<dim> &dof)
   {
     // this is to make it consistent with parallel_multigrid_adaptive.cc
     AffineConstraints<double> hanging_node_constraints;
-    hanging_node_constraints.reinit(locally_relevant_dofs);
+    hanging_node_constraints.reinit(dof.locally_owned_dofs(),
+                                    locally_relevant_dofs);
     DoFTools::make_hanging_node_constraints(dof, hanging_node_constraints);
     hanging_node_constraints.close();
 
@@ -189,7 +190,7 @@ do_test(const DoFHandler<dim> &dof)
       AffineConstraints<double> level_constraints;
       const IndexSet            relevant_dofs =
         DoFTools::extract_locally_relevant_level_dofs(dof, level);
-      level_constraints.reinit(relevant_dofs);
+      level_constraints.reinit(dof.locally_owned_mg_dofs(level), relevant_dofs);
       level_constraints.add_lines(
         mg_constrained_dofs.get_boundary_indices(level));
       level_constraints.close();
index 158fa98d311083b193a7839e913109b9b4b0bc39..23c02f6f5115816ed61f11420b6d0415a6e49e91 100644 (file)
@@ -252,7 +252,8 @@ do_test(const std::vector<const DoFHandler<dim> *> &dof)
   std::vector<const AffineConstraints<double> *> constraints_ptrs(dof.size());
   for (unsigned int i = 0; i < dof.size(); ++i)
     {
-      constraints[i].reinit(locally_relevant_dofs[i]);
+      constraints[i].reinit(dof[i]->locally_owned_dofs(),
+                            locally_relevant_dofs[i]);
       DoFTools::make_hanging_node_constraints(*dof[i], constraints[i]);
       VectorTools::interpolate_boundary_values(*dof[i],
                                                dirichlet_boundary,
@@ -301,7 +302,8 @@ do_test(const std::vector<const DoFHandler<dim> *> &dof)
         // this is to make it consistent with parallel_multigrid_adaptive.cc
         AffineConstraints<double> hanging_node_constraints;
 
-        hanging_node_constraints.reinit(locally_relevant_dofs[b]);
+        hanging_node_constraints.reinit(dof[b]->locally_owned_dofs(),
+                                        locally_relevant_dofs[b]);
         DoFTools::make_hanging_node_constraints(*dof[b],
                                                 hanging_node_constraints);
         hanging_node_constraints.close();
@@ -350,7 +352,8 @@ do_test(const std::vector<const DoFHandler<dim> *> &dof)
         {
           const IndexSet relevant_dofs =
             DoFTools::extract_locally_relevant_level_dofs(*dof[i], level);
-          level_constraints[i].reinit(relevant_dofs);
+          level_constraints[i].reinit(dof[i]->locally_owned_mg_dofs(level),
+                                      relevant_dofs);
           level_constraints[i].add_lines(
             mg_constrained_dofs[i].get_boundary_indices(level));
           level_constraints[i].close();
index 0f5a8cdccd24550f49be10696797cd4bd4824ac2..85392507e68070700e13dff6a0ab6152f99c0659 100644 (file)
@@ -226,7 +226,7 @@ do_test(const DoFHandler<dim> &dof)
 
   // fine-level constraints
   AffineConstraints<double> constraints;
-  constraints.reinit(locally_relevant_dofs);
+  constraints.reinit(dof.locally_owned_dofs(), locally_relevant_dofs);
   DoFTools::make_hanging_node_constraints(dof, constraints);
   VectorTools::interpolate_boundary_values(dof,
                                            dirichlet_boundary,
@@ -266,7 +266,8 @@ do_test(const DoFHandler<dim> &dof)
   {
     // this is to make it consistent with parallel_multigrid_adaptive.cc
     AffineConstraints<double> hanging_node_constraints;
-    hanging_node_constraints.reinit(locally_relevant_dofs);
+    hanging_node_constraints.reinit(dof.locally_owned_dofs(),
+                                    locally_relevant_dofs);
     DoFTools::make_hanging_node_constraints(dof, hanging_node_constraints);
     hanging_node_constraints.close();
 
@@ -296,7 +297,7 @@ do_test(const DoFHandler<dim> &dof)
       AffineConstraints<double> level_constraints;
       const IndexSet            relevant_dofs =
         DoFTools::extract_locally_relevant_level_dofs(dof, level);
-      level_constraints.reinit(relevant_dofs);
+      level_constraints.reinit(dof.locally_owned_mg_dofs(level), relevant_dofs);
       level_constraints.add_lines(
         mg_constrained_dofs.get_boundary_indices(level));
       level_constraints.close();
index 010c4cff2ce9a7c67e1149e3af0db7d33c0ea4d4..0496444b4e19e0246c5b44b25edd7f2e4079c4bf 100644 (file)
@@ -107,7 +107,7 @@ do_test(const DoFHandler<dim> &dof)
   const IndexSet locally_relevant_dofs =
     DoFTools::extract_locally_relevant_dofs(dof);
   AffineConstraints<double> constraints;
-  constraints.reinit(locally_relevant_dofs);
+  constraints.reinit(dof.locally_owned_dofs(), locally_relevant_dofs);
   VectorTools::interpolate_boundary_values(dof,
                                            dirichlet_boundary,
                                            constraints);
@@ -161,7 +161,7 @@ do_test(const DoFHandler<dim> &dof)
       AffineConstraints<double> level_constraints;
       const IndexSet            relevant_dofs =
         DoFTools::extract_locally_relevant_level_dofs(dof, level);
-      level_constraints.reinit(relevant_dofs);
+      level_constraints.reinit(dof.locally_owned_mg_dofs(level), relevant_dofs);
       level_constraints.add_lines(
         mg_constrained_dofs.get_boundary_indices(level));
       level_constraints.close();
index 12432f18f46eb043dfddba6193a88b45d8b3578d..7cdd947eaf6fdd90968eb97766a9b9ebc4628d24 100644 (file)
@@ -402,7 +402,7 @@ namespace Step37
       DoFTools::extract_locally_relevant_dofs(dof_handler);
 
     constraints.clear();
-    constraints.reinit(locally_relevant_dofs);
+    constraints.reinit(dof_handler.locally_owned_dofs(), locally_relevant_dofs);
     DoFTools::make_hanging_node_constraints(dof_handler, constraints);
     VectorTools::interpolate_boundary_values(dof_handler,
                                              0,
@@ -444,7 +444,8 @@ namespace Step37
         const IndexSet relevant_dofs =
           DoFTools::extract_locally_relevant_level_dofs(dof_handler, level);
         AffineConstraints<double> level_constraints;
-        level_constraints.reinit(relevant_dofs);
+        level_constraints.reinit(dof_handler.locally_owned_mg_dofs(level),
+                                 relevant_dofs);
         level_constraints.add_lines(
           mg_constrained_dofs.get_boundary_indices(level));
         level_constraints.close();
index ea9d2f50382cfaab832c8cc12512bc99d3f1eb89..e095b4a1d7727867cfa894b9b5be720e08af3557 100644 (file)
@@ -1040,7 +1040,8 @@ namespace StokesClass
 
     const IndexSet locally_relevant_dofs_u =
       DoFTools::extract_locally_relevant_dofs(dof_handler_u);
-    constraints_u.reinit(locally_relevant_dofs_u);
+    constraints_u.reinit(dof_handler_u.locally_owned_dofs(),
+                         locally_relevant_dofs_u);
     DoFTools::make_hanging_node_constraints(dof_handler_u, constraints_u);
 
     VectorTools::interpolate_boundary_values(
@@ -1049,7 +1050,8 @@ namespace StokesClass
 
     const IndexSet locally_relevant_dofs_p =
       DoFTools::extract_locally_relevant_dofs(dof_handler_p);
-    constraints_p.reinit(locally_relevant_dofs_p);
+    constraints_p.reinit(dof_handler_p.locally_owned_dofs(),
+                         locally_relevant_dofs_p);
     DoFTools::make_hanging_node_constraints(dof_handler_p, constraints_p);
     constraints_p.close();
 
@@ -1116,7 +1118,8 @@ namespace StokesClass
         const IndexSet relevant_dofs =
           DoFTools::extract_locally_relevant_level_dofs(dof_handler_u, level);
         AffineConstraints<double> level_constraints;
-        level_constraints.reinit(relevant_dofs);
+        level_constraints.reinit(dof_handler_u.locally_owned_mg_dofs(level),
+                                 relevant_dofs);
         level_constraints.add_lines(
           mg_constrained_dofs.get_boundary_indices(level));
         level_constraints.close();
@@ -1171,7 +1174,8 @@ namespace StokesClass
     AffineConstraints<double> constraints_u_no_dirchlet;
     const IndexSet            locally_relevant_dofs_u =
       DoFTools::extract_locally_relevant_dofs(dof_handler_u);
-    constraints_u_no_dirchlet.reinit(locally_relevant_dofs_u);
+    constraints_u_no_dirchlet.reinit(dof_handler_u.locally_owned_dofs(),
+                                     locally_relevant_dofs_u);
     DoFTools::make_hanging_node_constraints(dof_handler_u,
                                             constraints_u_no_dirchlet);
     constraints_u_no_dirchlet.close();
index dce5d4f1eaa2d04ebbe0614199d4438fc5766d4f..ae6431dd9fc2f3a6061cdb20e6474ce2f8399a76 100644 (file)
@@ -185,7 +185,8 @@ check_fe(FiniteElement<dim> &fe)
 
         const IndexSet relevant =
           DoFTools::extract_locally_relevant_level_dofs(dofh, level);
-        mg_boundary_constraints[level].reinit(relevant);
+        mg_boundary_constraints[level].reinit(dofh.locally_owned_mg_dofs(level),
+                                              relevant);
         mg_boundary_constraints[level].add_lines(boundary_indices[level]);
 
         mg_constrained_dofs_1.add_boundary_indices(dofh,
index d4a8012122203fcde051c47f3b878b979be53aa1..6b8c582e47045b4fc6e79a2480d90560d09d4596 100644 (file)
@@ -55,7 +55,8 @@ check()
       const IndexSet relevant_dofs =
         DoFTools::extract_locally_relevant_level_dofs(mgdof, level);
       AffineConstraints<double> level_constraints;
-      level_constraints.reinit(relevant_dofs);
+      level_constraints.reinit(mgdof.locally_owned_mg_dofs(level),
+                               relevant_dofs);
 
       typename DoFHandler<dim>::level_face_iterator face0 =
         mgdof.begin(level)->face(0);
index 9cf68327b06e30b51887cc452da8046082d10b8b..fa8f8b89e5067361d1db390025ea785a5229f0fa 100644 (file)
@@ -148,7 +148,8 @@ namespace Step50
 
     solution.reinit(mg_dof_handler.locally_owned_dofs(), MPI_COMM_WORLD);
     system_rhs.reinit(mg_dof_handler.locally_owned_dofs(), MPI_COMM_WORLD);
-    constraints.reinit(locally_relevant_set);
+    constraints.reinit(mg_dof_handler.locally_owned_dofs(),
+                       locally_relevant_set);
     DoFTools::make_hanging_node_constraints(mg_dof_handler, constraints);
 
     std::map<types::boundary_id, const Function<dim> *> dirichlet_boundary;
@@ -282,7 +283,8 @@ namespace Step50
       {
         const IndexSet dofset =
           DoFTools::extract_locally_relevant_level_dofs(mg_dof_handler, level);
-        boundary_constraints[level].reinit(dofset);
+        boundary_constraints[level].reinit(
+          mg_dof_handler.locally_owned_mg_dofs(level), dofset);
         boundary_constraints[level].add_lines(
           mg_constrained_dofs.get_refinement_edge_indices(level));
         boundary_constraints[level].add_lines(
index a563bcda8ce17db23c7600b0c853bd76799e86dc..c2076d47764a6c9e2743b607768e78d08ffd1df1 100644 (file)
@@ -91,7 +91,7 @@ test()
   const IndexSet locally_relevant_set =
     DoFTools::extract_locally_relevant_dofs(mg_dof_handler);
 
-  constraints.reinit(locally_relevant_set);
+  constraints.reinit(mg_dof_handler.locally_owned_dofs(), locally_relevant_set);
   DoFTools::make_hanging_node_constraints(mg_dof_handler, constraints);
 
   const std::set<types::boundary_id> dirichlet_boundary_ids = {0};
index 953382ced24347e9ab3bcc0648efaa23e0e01c4c..746c70ba33e0a0318b1115d67f55ecd579857a80 100644 (file)
@@ -85,7 +85,7 @@ test()
     DoFTools::extract_locally_relevant_dofs(dof_handler);
 
   AffineConstraints<double> constraints;
-  constraints.reinit(locally_relevant_dofs);
+  constraints.reinit(dof_handler.locally_owned_dofs(), locally_relevant_dofs);
   DoFTools::make_hanging_node_constraints(dof_handler, constraints);
 
   std::vector<
index 0224000798564043a5e80e67761aed72e7500b67..aedebd1e989befcbcbfb437026fde5ca51daea09 100644 (file)
@@ -205,7 +205,8 @@ namespace Step50
 
     solution.reinit(mg_dof_handler.locally_owned_dofs(), MPI_COMM_WORLD);
     system_rhs.reinit(mg_dof_handler.locally_owned_dofs(), MPI_COMM_WORLD);
-    constraints.reinit(locally_relevant_set);
+    constraints.reinit(mg_dof_handler.locally_owned_dofs(),
+                       locally_relevant_set);
     DoFTools::make_hanging_node_constraints(mg_dof_handler, constraints);
 
     std::map<types::boundary_id, const Function<dim> *> dirichlet_boundary;
@@ -349,7 +350,8 @@ namespace Step50
       {
         const IndexSet dofset =
           DoFTools::extract_locally_relevant_level_dofs(mg_dof_handler, level);
-        boundary_constraints[level].reinit(dofset);
+        boundary_constraints[level].reinit(
+          mg_dof_handler.locally_owned_mg_dofs(level), dofset);
         boundary_constraints[level].add_lines(
           mg_constrained_dofs.get_refinement_edge_indices(level));
         boundary_constraints[level].add_lines(
index 618e511cbe4cf390983f28cf49e736fc74891731..77223a3ec80adbae137c2a33f60cfcd7e02eef5a 100644 (file)
@@ -199,7 +199,8 @@ namespace Step50
     solution.reinit(mg_dof_handler.locally_owned_dofs(), MPI_COMM_WORLD);
     system_rhs.reinit(mg_dof_handler.locally_owned_dofs(), MPI_COMM_WORLD);
 
-    constraints.reinit(locally_relevant_set);
+    constraints.reinit(mg_dof_handler.locally_owned_dofs(),
+                       locally_relevant_set);
     DoFTools::make_hanging_node_constraints(mg_dof_handler, constraints);
 
     std::map<types::boundary_id, const Function<dim> *> dirichlet_boundary;
@@ -348,7 +349,8 @@ namespace Step50
       {
         const IndexSet dofset =
           DoFTools::extract_locally_relevant_level_dofs(mg_dof_handler, level);
-        boundary_constraints[level].reinit(dofset);
+        boundary_constraints[level].reinit(
+          mg_dof_handler.locally_owned_mg_dofs(level), dofset);
 
         for (const types::global_dof_index dof_index :
              mg_constrained_dofs.get_refinement_edge_indices(level))
index 6c3c688d600084b8103a8b59425ee0f1fa5c3ab1..1347c5019e140dd0bddc93d66bbf118928ee7b4c 100644 (file)
@@ -199,7 +199,8 @@ namespace Step50
     solution.reinit(mg_dof_handler.locally_owned_dofs(), MPI_COMM_WORLD);
     system_rhs.reinit(mg_dof_handler.locally_owned_dofs(), MPI_COMM_WORLD);
 
-    constraints.reinit(locally_relevant_set);
+    constraints.reinit(mg_dof_handler.locally_owned_dofs(),
+                       locally_relevant_set);
     DoFTools::make_hanging_node_constraints(mg_dof_handler, constraints);
 
     std::map<types::boundary_id, const Function<dim> *> dirichlet_boundary;
@@ -347,7 +348,8 @@ namespace Step50
       {
         const IndexSet dofset =
           DoFTools::extract_locally_relevant_level_dofs(mg_dof_handler, level);
-        boundary_constraints[level].reinit(dofset);
+        boundary_constraints[level].reinit(
+          mg_dof_handler.locally_owned_mg_dofs(level), dofset);
 
         for (const types::global_dof_index dof_index :
              mg_constrained_dofs.get_refinement_edge_indices(level))
index 61a2dc0213ab60b104b50eb511edbf8bcee86ab0..45e70401b4a70c5d3bd09873b23ed38064edae52 100644 (file)
@@ -199,7 +199,8 @@ namespace Step50
     solution.reinit(mg_dof_handler.locally_owned_dofs(), MPI_COMM_WORLD);
     system_rhs.reinit(mg_dof_handler.locally_owned_dofs(), MPI_COMM_WORLD);
 
-    constraints.reinit(locally_relevant_set);
+    constraints.reinit(mg_dof_handler.locally_owned_dofs(),
+                       locally_relevant_set);
     DoFTools::make_hanging_node_constraints(mg_dof_handler, constraints);
 
     std::map<types::boundary_id, const Function<dim> *> dirichlet_boundary;
@@ -348,7 +349,8 @@ namespace Step50
       {
         const IndexSet dofset =
           DoFTools::extract_locally_relevant_level_dofs(mg_dof_handler, level);
-        boundary_constraints[level].reinit(dofset);
+        boundary_constraints[level].reinit(
+          mg_dof_handler.locally_owned_mg_dofs(level), dofset);
         for (const types::global_dof_index dof_index :
              mg_constrained_dofs.get_refinement_edge_indices(level))
           boundary_constraints[level].constrain_dof_to_zero(dof_index);
index 56338239fa860df3bcc53e00878f268cba8d6f7c..3b45c23071cc622c6aeab539121c0cc110d9e7e7 100644 (file)
@@ -202,7 +202,8 @@ namespace Step50
 
     solution.reinit(mg_dof_handler.locally_owned_dofs(), MPI_COMM_WORLD);
     system_rhs.reinit(mg_dof_handler.locally_owned_dofs(), MPI_COMM_WORLD);
-    constraints.reinit(locally_relevant_set);
+    constraints.reinit(mg_dof_handler.locally_owned_dofs(),
+                       locally_relevant_set);
     DoFTools::make_hanging_node_constraints(mg_dof_handler, constraints);
 
     std::map<types::boundary_id, const Function<dim> *> dirichlet_boundary;
@@ -346,7 +347,8 @@ namespace Step50
       {
         const IndexSet dofset =
           DoFTools::extract_locally_relevant_level_dofs(mg_dof_handler, level);
-        boundary_constraints[level].reinit(dofset);
+        boundary_constraints[level].reinit(
+          mg_dof_handler.locally_owned_mg_dofs(level), dofset);
         boundary_constraints[level].add_lines(
           mg_constrained_dofs.get_refinement_edge_indices(level));
         boundary_constraints[level].add_lines(
index 710b3ad30a501e7c754b24f7339d301f04800236..47b97010226212aafcba1955a77335629a82496a 100644 (file)
@@ -203,7 +203,8 @@ namespace Step50
     solution.reinit(mg_dof_handler.locally_owned_dofs(), MPI_COMM_WORLD);
     system_rhs.reinit(mg_dof_handler.locally_owned_dofs(), MPI_COMM_WORLD);
 
-    constraints.reinit(locally_relevant_set);
+    constraints.reinit(mg_dof_handler.locally_owned_dofs(),
+                       locally_relevant_set);
     DoFTools::make_hanging_node_constraints(mg_dof_handler, constraints);
 
     const std::set<types::boundary_id> dirichlet_boundary_ids = {0};
@@ -354,7 +355,8 @@ namespace Step50
       {
         const IndexSet dofset =
           DoFTools::extract_locally_relevant_level_dofs(mg_dof_handler, level);
-        boundary_constraints[level].reinit(dofset);
+        boundary_constraints[level].reinit(
+          mg_dof_handler.locally_owned_mg_dofs(level), dofset);
         boundary_constraints[level].add_lines(
           mg_constrained_dofs.get_refinement_edge_indices(level));
         boundary_constraints[level].add_lines(
index e39defb95ed869bd51386f5bf6037e4aa88f199d..d131b47a7d765919376f86177297069c45631ba9 100644 (file)
@@ -130,7 +130,8 @@ check_fe(FiniteElement<dim> &fe)
   AffineConstraints<double> hanging_node_constraints;
   const IndexSet            locally_relevant_set =
     DoFTools::extract_locally_relevant_dofs(dofh);
-  hanging_node_constraints.reinit(locally_relevant_set);
+  hanging_node_constraints.reinit(dofh.locally_owned_dofs(),
+                                  locally_relevant_set);
   DoFTools::make_hanging_node_constraints(dofh, hanging_node_constraints);
   hanging_node_constraints.close();
 
index d6fd593072559a95afac4abc7aa520d42f22179c..85e561c98be0adb45ec9f49400784f6e9e256f41 100644 (file)
@@ -54,7 +54,7 @@ check()
   const IndexSet relevant_dofs =
     DoFTools::extract_locally_relevant_level_dofs(mgdof, 0);
   AffineConstraints<double> user_constraints;
-  user_constraints.reinit(relevant_dofs);
+  user_constraints.reinit(mgdof.locally_owned_mg_dofs(0), relevant_dofs);
 
   typename DoFHandler<dim>::level_face_iterator face0 = mgdof.begin(0)->face(0);
   std::vector<types::global_dof_index>          face_dofs(fe.dofs_per_face);
index dd79ecc498ab74893e88d335071a05e5cf81f6a8..2a38bca582410669e8c03a5745c3d24fc4193a85 100644 (file)
@@ -51,7 +51,7 @@ check()
   const IndexSet relevant_dofs =
     DoFTools::extract_locally_relevant_level_dofs(mgdof, 0);
   AffineConstraints<double> user_constraints;
-  user_constraints.reinit(relevant_dofs);
+  user_constraints.reinit(mgdof.locally_owned_mg_dofs(0), relevant_dofs);
 
   typename DoFHandler<dim>::level_face_iterator face0 = mgdof.begin(0)->face(0);
   std::vector<types::global_dof_index>          face_dofs(fe.dofs_per_face);
index de1d8bb5e72a3346d46fa5078eafd8ea1c567024..e39de8c2d1332521a37af97baeda92d14ba6e0d0 100644 (file)
@@ -1247,7 +1247,7 @@ namespace Step70
     fluid_relevant_dofs[1] = locally_relevant_dofs.get_view(n_u, n_u + n_p);
 
     {
-      constraints.reinit(locally_relevant_dofs);
+      constraints.reinit(fluid_dh.locally_owned_dofs(), locally_relevant_dofs);
 
       const FEValuesExtractors::Vector velocities(0);
       DoFTools::make_hanging_node_constraints(fluid_dh, constraints);
index bfad4b552ba2d5463fcbd204fe52c5ee09297f11..f62a4ec9ba6191620787da966566332a4bfba229 100644 (file)
@@ -1042,7 +1042,8 @@ namespace StokesClass
 
     const IndexSet locally_relevant_dofs_u =
       DoFTools::extract_locally_relevant_dofs(dof_handler_u);
-    constraints_u.reinit(locally_relevant_dofs_u);
+    constraints_u.reinit(dof_handler_u.locally_owned_dofs(),
+                         locally_relevant_dofs_u);
     DoFTools::make_hanging_node_constraints(dof_handler_u, constraints_u);
 
     // VectorTools::interpolate_boundary_values(
@@ -1058,7 +1059,8 @@ namespace StokesClass
 
     const IndexSet locally_relevant_dofs_p =
       DoFTools::extract_locally_relevant_dofs(dof_handler_p);
-    constraints_p.reinit(locally_relevant_dofs_p);
+    constraints_p.reinit(dof_handler_p.locally_owned_dofs(),
+                         locally_relevant_dofs_p);
     DoFTools::make_hanging_node_constraints(dof_handler_p, constraints_p);
     constraints_p.close();
 
@@ -1125,7 +1127,8 @@ namespace StokesClass
         const IndexSet relevant_dofs =
           DoFTools::extract_locally_relevant_level_dofs(dof_handler_u, level);
         AffineConstraints<double> level_constraints;
-        level_constraints.reinit(relevant_dofs);
+        level_constraints.reinit(dof_handler_u.locally_owned_mg_dofs(level),
+                                 relevant_dofs);
         // level_constraints.add_lines(
         //   mg_constrained_dofs.get_boundary_indices(level));
 
@@ -1198,7 +1201,8 @@ namespace StokesClass
     AffineConstraints<double> constraints_u_no_dirchlet;
     const IndexSet            locally_relevant_dofs_u =
       DoFTools::extract_locally_relevant_dofs(dof_handler_u);
-    constraints_u_no_dirchlet.reinit(locally_relevant_dofs_u);
+    constraints_u_no_dirchlet.reinit(dof_handler_u.locally_owned_dofs(),
+                                     locally_relevant_dofs_u);
     DoFTools::make_hanging_node_constraints(dof_handler_u,
                                             constraints_u_no_dirchlet);
     constraints_u_no_dirchlet.close();
index 940f980da95848a0aa48b03a7c9112831b6f221f..f6873a6e3eeb6fa1556e94e15889c12178669cf2 100644 (file)
@@ -170,7 +170,7 @@ namespace MPI_nonlinear_solver_selector_test
 
         {
           nonzero_constraints.clear();
-          nonzero_constraints.reinit(locally_relevant_dofs);
+          nonzero_constraints.reinit(locally_owned_dofs, locally_relevant_dofs);
           DoFTools::make_hanging_node_constraints(dof_handler,
                                                   nonzero_constraints);
 
@@ -192,7 +192,7 @@ namespace MPI_nonlinear_solver_selector_test
 
         {
           zero_constraints.clear();
-          zero_constraints.reinit(locally_relevant_dofs);
+          zero_constraints.reinit(locally_owned_dofs, locally_relevant_dofs);
           DoFTools::make_hanging_node_constraints(dof_handler,
                                                   zero_constraints);
           VectorTools::interpolate_boundary_values(
index 06dc39f2daab62e71bf41d435a6b9dc5cf88e576..417103564925440a9279fe80de18fa41a9b21f32 100644 (file)
@@ -118,7 +118,7 @@ do_project(const parallel::distributed::Triangulation<dim> &triangulation,
     DoFTools::extract_locally_relevant_dofs(dof_handler);
 
   AffineConstraints<double> constraints;
-  constraints.reinit(locally_relevant_dofs);
+  constraints.reinit(locally_owned_dofs, locally_relevant_dofs);
   DoFTools::make_hanging_node_constraints(dof_handler, constraints);
   constraints.close();
 
index 5d53c945eadc8038c03b99a6579f67a4cd28b302..e77ffff0ae22be51b518aa0663b088e0af7c9946 100644 (file)
@@ -114,7 +114,7 @@ do_project(const parallel::distributed::Triangulation<dim> &triangulation,
     DoFTools::extract_locally_relevant_dofs(dof_handler);
 
   AffineConstraints<double> constraints;
-  constraints.reinit(locally_relevant_dofs);
+  constraints.reinit(dof_handler.locally_owned_dofs(), locally_relevant_dofs);
   DoFTools::make_hanging_node_constraints(dof_handler, constraints);
   constraints.close();
 
index e9d650b934e77d2eb37af6d01805716ff5a58fec..11b52656aa34ad4753e41eef63dfb65e0c54196c 100644 (file)
@@ -126,7 +126,8 @@ do_project(const parallel::distributed::Triangulation<dim> &triangulation,
       locally_relevant_dofs[i] =
         DoFTools::extract_locally_relevant_dofs(*dof_handlers[i]);
 
-      constraints[i].reinit(locally_relevant_dofs[i]);
+      constraints[i].reinit(dof_handlers[i]->locally_owned_dofs(),
+                            locally_relevant_dofs[i]);
       DoFTools::make_hanging_node_constraints(*dof_handlers[i], constraints[i]);
       constraints[i].close();
 
index cc6313a7aa030ec6bc8a7a98ca5a86f97289b6cd..dd994254b03e6da464421aa70172e68d13445ae0 100644 (file)
@@ -504,7 +504,7 @@ LaplaceProblem<dim>::setup_dofs()
 
       IndexSet relevant_dofs = DoFTools::extract_locally_relevant_dofs(dof_h);
       AffineConstraints<float> &constraints = level_constraints[level];
-      constraints.reinit(relevant_dofs);
+      constraints.reinit(dof_h.locally_owned_dofs(), relevant_dofs);
       DoFTools::make_hanging_node_constraints(dof_h, constraints);
       constraints.close();
       typename MatrixFree<dim, float>::AdditionalData additional_data;
@@ -518,7 +518,7 @@ LaplaceProblem<dim>::setup_dofs()
       // now create the final constraints object
       relevant_dofs = DoFTools::extract_locally_relevant_dofs(dof_h);
       constraints.clear();
-      constraints.reinit(relevant_dofs);
+      constraints.reinit(dof_h.locally_owned_dofs(), relevant_dofs);
       DoFTools::make_hanging_node_constraints(dof_h, constraints);
       constraints.close();
     }
index d79b605f9fa74e244ca4b066864a47608eb49fee..d52651bf2c801c75dd288b5599bcce09387d32f8 100644 (file)
@@ -372,7 +372,7 @@ LaplaceProblem<dim>::setup_dofs()
     DoFTools::extract_locally_relevant_dofs(dof_handler);
 
   constraints.clear();
-  constraints.reinit(locally_relevant_dofs);
+  constraints.reinit(dof_handler.locally_owned_dofs(), locally_relevant_dofs);
   DoFTools::make_hanging_node_constraints(dof_handler, constraints);
   VectorTools::interpolate_boundary_values(
     mapping, dof_handler, 0, Functions::ZeroFunction<dim>(), constraints);
@@ -393,7 +393,7 @@ LaplaceProblem<dim>::setup_dofs()
       const IndexSet relevant_dofs =
         DoFTools::extract_locally_relevant_level_dofs(dof_handler, level);
       AffineConstraints<double> level_constraints;
-      level_constraints.reinit(relevant_dofs);
+      level_constraints.reinit(dof_handler.locally_owned_dofs(), relevant_dofs);
       level_constraints.add_lines(
         mg_constrained_dofs.get_boundary_indices(level));
       level_constraints.close();
@@ -442,7 +442,7 @@ LaplaceProblem<dim>::setup_matrix_free()
       const IndexSet relevant_dofs =
         DoFTools::extract_locally_relevant_level_dofs(dof_handler, level);
       AffineConstraints<double> level_constraints;
-      level_constraints.reinit(relevant_dofs);
+      level_constraints.reinit(dof_handler.locally_owned_dofs(), relevant_dofs);
       level_constraints.add_lines(
         mg_constrained_dofs.get_boundary_indices(level));
       level_constraints.close();
index f4dd71b5878d0814c8a291a651979b69e268dfaa..6f139ff2a2e3fe58bf6c82e578a9be9ebbdec790 100644 (file)
@@ -372,7 +372,7 @@ LaplaceProblem<dim>::setup_dofs()
     DoFTools::extract_locally_relevant_dofs(dof_handler);
 
   constraints.clear();
-  constraints.reinit(locally_relevant_dofs);
+  constraints.reinit(dof_handler.locally_owned_dofs(), locally_relevant_dofs);
   DoFTools::make_hanging_node_constraints(dof_handler, constraints);
   VectorTools::interpolate_boundary_values(
     mapping, dof_handler, 0, Functions::ZeroFunction<dim>(), constraints);
@@ -393,7 +393,7 @@ LaplaceProblem<dim>::setup_dofs()
       const IndexSet relevant_dofs =
         DoFTools::extract_locally_relevant_level_dofs(dof_handler, level);
       AffineConstraints<double> level_constraints;
-      level_constraints.reinit(relevant_dofs);
+      level_constraints.reinit(dof_handler.locally_owned_dofs(), relevant_dofs);
       level_constraints.add_lines(
         mg_constrained_dofs.get_boundary_indices(level));
       level_constraints.close();
@@ -442,7 +442,7 @@ LaplaceProblem<dim>::setup_matrix_free()
       const IndexSet relevant_dofs =
         DoFTools::extract_locally_relevant_level_dofs(dof_handler, level);
       AffineConstraints<double> level_constraints;
-      level_constraints.reinit(relevant_dofs);
+      level_constraints.reinit(dof_handler.locally_owned_dofs(), relevant_dofs);
       level_constraints.add_lines(
         mg_constrained_dofs.get_boundary_indices(level));
       level_constraints.close();
index c4f1bcd78525e462ecc378d56a8763a62e8b3d7d..9bc1317f353c72d20c2521359b3ad180982c6722 100644 (file)
@@ -182,7 +182,7 @@ namespace Step77
     // Specifically, we need two types of AffineConstraints.
     // One to handle homogeneous boundary conditions for the update step.
     zero_constraints.clear();
-    zero_constraints.reinit(locally_relevant_dofs);
+    zero_constraints.reinit(locally_owned_dofs, locally_relevant_dofs);
     DoFTools::make_hanging_node_constraints(dof_handler, zero_constraints);
     VectorTools::interpolate_boundary_values(dof_handler,
                                              0,
@@ -193,7 +193,7 @@ namespace Step77
     // And another one to handle non-homogeneous boundary conditions
     // when computing the residual function.
     bc_constraints.clear();
-    bc_constraints.reinit(locally_relevant_dofs);
+    bc_constraints.reinit(locally_owned_dofs, locally_relevant_dofs);
     DoFTools::make_hanging_node_constraints(dof_handler, bc_constraints);
     VectorTools::interpolate_boundary_values(dof_handler,
                                              0,
index 26d3134728c6ac58055904451e82adddee335acd..acc02a78cdfbdeeb7dc3cfc837919d67782a3c27 100644 (file)
@@ -104,7 +104,7 @@ test()
   const types::global_dof_index n_local_dofs = locally_owned_dofs.n_elements();
 
   AffineConstraints<double> constraints;
-  constraints.reinit(locally_relevant_dofs);
+  constraints.reinit(locally_owned_dofs, locally_relevant_dofs);
   DoFTools::make_hanging_node_constraints(dof_handler, constraints);
   constraints.close();
   // set vector:
index af9599bd543ff48cd357975fd069f4c1176125ef..cf736f7bf8ce7833d6b1bda9216d6b4f45bd99a1 100644 (file)
@@ -84,7 +84,7 @@ test(const unsigned int poly_degree = 1)
   vector.reinit(locally_owned_dofs, mpi_communicator);
 
   AffineConstraints<PetscScalar> constraints;
-  constraints.reinit(locally_relevant_dofs);
+  constraints.reinit(locally_owned_dofs, locally_relevant_dofs);
   DoFTools::make_hanging_node_constraints(dof_handler, constraints);
   constraints.close();
 
index a3d673b35b6df1116e0b5a2cc637b5e22b22fc1e..b628a13fac09a5395df4eafb4d4beec5f6211b9f 100644 (file)
@@ -52,7 +52,8 @@ test()
 
   AffineConstraints<double> constraints;
 
-  constraints.reinit(DoFTools::extract_locally_relevant_dofs(dof_handler));
+  constraints.reinit(dof_handler.locally_owned_dofs(),
+                     DoFTools::extract_locally_relevant_dofs(dof_handler));
   for (unsigned int d = 0; d < dim; ++d)
     DoFTools::make_periodicity_constraints(
       dof_handler, 2 * d, 2 * d + 1, d, constraints);
index 3fc0c210ac7554b5649005c3374f4a5ec1d34b7c..d4c96658b7ba0e09dc8c800c37b961b8021cc6fc 100644 (file)
@@ -172,7 +172,7 @@ namespace Step40
     system_rhs.reinit(locally_owned_dofs, mpi_communicator);
 
     constraints.clear();
-    constraints.reinit(locally_relevant_dofs);
+    constraints.reinit(locally_owned_dofs, locally_relevant_dofs);
     DoFTools::make_hanging_node_constraints(dof_handler, constraints);
     VectorTools::interpolate_boundary_values(
       mapping, dof_handler, 0, Functions::ZeroFunction<dim>(), constraints);
index e6dffbf5eccd31354495052eaa67bc267b5b4d27..0bf16b6f168c9251d68748ddcc97caa370ad92a6 100644 (file)
@@ -465,7 +465,8 @@ namespace Step55
     // to put this function call in, in case adaptive refinement gets
     // introduced later.
     {
-      constraints.reinit(locally_relevant_dofs);
+      constraints.reinit(dof_handler.locally_owned_dofs(),
+                         locally_relevant_dofs);
 
       const FEValuesExtractors::Vector velocities(0);
       DoFTools::make_hanging_node_constraints(dof_handler, constraints);
index bb3a3073d821d9bf140f41503a97acb449222914..03b25aa6fdaa45370ea2ccd5ad5fa3e663d09393 100644 (file)
@@ -107,7 +107,7 @@ test(std::string solver_name, std::string preconditioner_name)
     dealii::DoFTools::extract_locally_relevant_dofs(dof_handler);
 
   constraints.clear();
-  constraints.reinit(locally_relevant_dofs);
+  constraints.reinit(locally_owned_dofs, locally_relevant_dofs);
   dealii::DoFTools::make_hanging_node_constraints(dof_handler, constraints);
   dealii::VectorTools::interpolate_boundary_values(
     dof_handler, 0, dealii::Functions::ZeroFunction<dim>(), constraints);
index 3faba96cbb5d379b667552c2521b933fb00a7fa3..3c25f7bce4dfefa7fec51610e0fa294469101b08 100644 (file)
@@ -107,7 +107,7 @@ test(std::string solver_name, std::string preconditioner_name)
     dealii::DoFTools::extract_locally_relevant_dofs(dof_handler);
 
   constraints.clear();
-  constraints.reinit(locally_relevant_dofs);
+  constraints.reinit(locally_owned_dofs, locally_relevant_dofs);
   dealii::DoFTools::make_hanging_node_constraints(dof_handler, constraints);
   dealii::VectorTools::interpolate_boundary_values(
     dof_handler, 0, dealii::Functions::ZeroFunction<dim>(), constraints);
index 5a4e4f56548d038261b500e8252378359cf9bf8e..f916502696de45b2d9539b7cec7b06d37783305d 100644 (file)
@@ -107,7 +107,7 @@ test(std::string solver_name, std::string preconditioner_name)
     dealii::DoFTools::extract_locally_relevant_dofs(dof_handler);
 
   constraints.clear();
-  constraints.reinit(locally_relevant_dofs);
+  constraints.reinit(locally_owned_dofs, locally_relevant_dofs);
   dealii::DoFTools::make_hanging_node_constraints(dof_handler, constraints);
   dealii::VectorTools::interpolate_boundary_values(
     dof_handler, 0, dealii::Functions::ZeroFunction<dim>(), constraints);
index 4e8bbf84a0f309e24b35d022fa729080b04925ba..46733ab51de050255a6a05c14ecc53d1c68d55e6 100644 (file)
@@ -202,7 +202,7 @@ Test_Solver_Output::setup_system()
   system_rhs.reinit(locally_owned_dofs, mpi_comm);
 
   constraints.clear();
-  constraints.reinit(locally_relevant_dofs);
+  constraints.reinit(locally_owned_dofs, locally_relevant_dofs);
   DoFTools::make_hanging_node_constraints(dof_handler, constraints);
   VectorTools::interpolate_boundary_values(dof_handler,
                                            0,

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.