From: Wolfgang Bangerth <bangerth@math.tamu.edu>
Date: Fri, 12 Jun 2015 01:00:03 +0000 (-0500)
Subject: Cleanup: only revert the global-to-local numbering in the parallel case.
X-Git-Tag: v8.3.0-rc1~120^2~2
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=918cddf9e7921cf779ac4d4b2de2085c82092c92;p=dealii.git

Cleanup: only revert the global-to-local numbering in the parallel case.
---

diff --git a/source/dofs/dof_renumbering.cc b/source/dofs/dof_renumbering.cc
index d197d5178a..d0cd8f6d4f 100644
--- a/source/dofs/dof_renumbering.cc
+++ b/source/dofs/dof_renumbering.cc
@@ -422,7 +422,7 @@ namespace DoFRenumbering
         for (unsigned int i=0; i<locally_owned.n_elements(); ++i)
           {
             const types::global_dof_index row = locally_owned.nth_index_in_set(i);
-            row_entries.resize(0);
+            row_entries.clear();
             for (DynamicSparsityPattern::iterator it =
                    dsp.begin(row); it != dsp.end(row); ++it)
               if (it->column() != row && locally_owned.is_element(it->column()))
@@ -431,23 +431,24 @@ namespace DoFRenumbering
                                  true);
           }
 
-        AssertDimension(new_indices.size(), sparsity.n_rows());
+        AssertDimension(new_indices.size(), locally_owned.n_elements());
         SparsityTools::reorder_Cuthill_McKee (sparsity, new_indices,
                                               starting_indices);
+        if (reversed_numbering)
+          new_indices = Utilities::reverse_permutation (new_indices);
+
+        // convert indices back to global index space
+        for (std::size_t i=0; i<new_indices.size(); ++i)
+          new_indices[i] = locally_owned.nth_index_in_set(new_indices[i]);
       }
     else
       {
         AssertDimension(new_indices.size(), dsp.n_rows());
         SparsityTools::reorder_Cuthill_McKee (dsp, new_indices,
                                               starting_indices);
+        if (reversed_numbering)
+          new_indices = Utilities::reverse_permutation (new_indices);
       }
-
-    if (reversed_numbering)
-      new_indices = Utilities::reverse_permutation (new_indices);
-
-    // convert indices back to global index space
-    for (std::size_t i=0; i<new_indices.size(); ++i)
-      new_indices[i] = locally_owned.nth_index_in_set(new_indices[i]);
   }