]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
step-32: use correct IndexSets
authorheister <heister@0785d39b-7218-0410-832d-ea1e28bc413d>
Sun, 9 Jun 2013 21:08:20 +0000 (21:08 +0000)
committerheister <heister@0785d39b-7218-0410-832d-ea1e28bc413d>
Sun, 9 Jun 2013 21:08:20 +0000 (21:08 +0000)
git-svn-id: https://svn.dealii.org/branches/branch_unify_linear_algebra@29801 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/examples/step-32/step-32.cc

index 0293847747d5ea232ea7d4d2015e69c493cf1c61..1841371d7ee461f642e2ecc878db9519f94c6e1d 100644 (file)
@@ -1846,7 +1846,7 @@ namespace Step32
     stokes_matrix.clear ();
 
 #ifdef USE_PETSC_LA
-    LA::MPI::CompressedBlockSparsityPattern sp(stokes_partitioning);
+    LA::MPI::CompressedBlockSparsityPattern sp(stokes_relevant_partitioning);
 #else
     LA::MPI::CompressedBlockSparsityPattern sp(stokes_partitioning, MPI_COMM_WORLD);    
 #endif
@@ -1898,7 +1898,7 @@ namespace Step32
     stokes_preconditioner_matrix.clear ();
 
 #ifdef USE_PETSC_LA
-    LA::MPI::CompressedBlockSparsityPattern sp (stokes_partitioning);
+    LA::MPI::CompressedBlockSparsityPattern sp (stokes_relevant_partitioning);
 #else
     TrilinosWrappers::BlockSparsityPattern sp(stokes_partitioning,
                                               MPI_COMM_WORLD);
@@ -1948,8 +1948,12 @@ namespace Step32
     temperature_matrix.clear ();
 
 #ifdef USE_PETSC_LA
+    IndexSet relevant;
+
+    DoFTools::extract_locally_relevant_dofs (temperature_dof_handler,
+                                               relevant);
     // TODO:    LA::MPI::CompressedSparsityPattern
-    CompressedSimpleSparsityPattern sp (temperature_partitioner);
+    CompressedSimpleSparsityPattern sp (relevant);
 #else
     TrilinosWrappers::SparsityPattern sp(temperature_partitioner,
                                          MPI_COMM_WORLD);
@@ -1959,10 +1963,6 @@ namespace Step32
                                      Utilities::MPI::
                                      this_mpi_process(MPI_COMM_WORLD));
 #ifdef USE_PETSC_LA
-    IndexSet relevant;
-    
-    DoFTools::extract_locally_relevant_dofs (temperature_dof_handler,
-                                               relevant);
 
     SparsityTools::distribute_sparsity_pattern(sp,
                                               temperature_dof_handler.n_locally_owned_dofs_per_processor(),
@@ -2059,7 +2059,7 @@ namespace Step32
     // various vectors we keep around in this program.
     stokes_partitioning.clear();
     stokes_relevant_partitioning.clear();   
-    IndexSet temperature_partitioning (n_T), temperature_relevant_partitioning (n_T);
+    IndexSet temperature_partitioning , temperature_relevant_partitioning;
     IndexSet stokes_relevant_set;
     {
       IndexSet stokes_index_set = stokes_dof_handler.locally_owned_dofs();

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.