From: Timo Heister <timo.heister@gmail.com>
Date: Thu, 30 Mar 2017 19:09:45 +0000 (-0400)
Subject: fix step-55
X-Git-Tag: v8.5.0-rc2~11^2
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=refs%2Fpull%2F4153%2Fhead;p=dealii.git

fix step-55

This causes crashes with PETSc and solver convergence issues with
Trilinos with more than one rank.
---

diff --git a/examples/step-55/step-55.cc b/examples/step-55/step-55.cc
index ab01da6bd1..652ae2f236 100644
--- a/examples/step-55/step-55.cc
+++ b/examples/step-55/step-55.cc
@@ -440,14 +440,14 @@ namespace Step55
 
       DoFTools::make_sparsity_pattern (dof_handler, coupling, dsp,
                                        constraints, false);
-      preconditioner_matrix.reinit (owned_partitioning,
-                                    //      owned_partitioning,
-                                    dsp,
-                                    mpi_communicator);
       SparsityTools::distribute_sparsity_pattern (dsp,
                                                   dof_handler.locally_owned_dofs_per_processor(),
                                                   mpi_communicator,
                                                   locally_relevant_dofs);
+      preconditioner_matrix.reinit (owned_partitioning,
+                                    //      owned_partitioning,
+                                    dsp,
+                                    mpi_communicator);
     }
 
     // Finally, we construct the block vectors with the right sizes. The