From: Ralf Hartmann <Ralf.Hartmann@dlr.de>
Date: Wed, 19 May 2004 09:12:16 +0000 (+0000)
Subject: Revert last patch as it triggers strange PETSc error messages like
X-Git-Tag: v8.0.0~15165
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=833107f569042ff368389cdd660e21374219125c;p=dealii.git

Revert last patch as it triggers strange PETSc error messages like


git-svn-id: https://svn.dealii.org/trunk@9264 0785d39b-7218-0410-832d-ea1e28bc413d
---

diff --git a/deal.II/lac/source/petsc_parallel_sparse_matrix.cc b/deal.II/lac/source/petsc_parallel_sparse_matrix.cc
index 8b8d754573..a445c2ab24 100644
--- a/deal.II/lac/source/petsc_parallel_sparse_matrix.cc
+++ b/deal.II/lac/source/petsc_parallel_sparse_matrix.cc
@@ -115,9 +115,15 @@ namespace PETScWrappers
                                        // use the call sequence indicating only
                                        // a maximal number of elements per row
                                        // for all rows globally
-//TODO[WB]: We should do better by providing ways to tell PETSc how to partition the columns
+                                       //
+                                       // note that we partition the columns
+                                       // in the same way as we partition the
+                                       // rows, i.e. we assume that users
+                                       // specify the same number for
+                                       // local_rows as they do for local_size
+                                       // on parallel vectors
       const int ierr
-        = MatCreateMPIAIJ(communicator, local_rows, PETSC_DECIDE, m, n,
+        = MatCreateMPIAIJ(communicator, local_rows, local_rows, m, n,
                           n_nonzero_per_row, 0, 0, 0,
                           &matrix);
       AssertThrow (ierr == 0, ExcPETScError(ierr));
@@ -154,9 +160,14 @@ namespace PETScWrappers
       const std::vector<signed int> int_row_lengths (row_lengths.begin(),
                                                      row_lengths.end());
                                        
-//TODO[WB]: We should do better by providing ways to tell PETSc how to partition the columns
+                                       // note that we partition the columns
+                                       // in the same way as we partition the
+                                       // rows, i.e. we assume that users
+                                       // specify the same number for
+                                       // local_rows as they do for local_size
+                                       // on parallel vectors
       const int ierr
-        = MatCreateMPIAIJ(communicator, local_rows, PETSC_DECIDE, m, n,
+        = MatCreateMPIAIJ(communicator, local_rows, local_rows, m, n,
                           0, &int_row_lengths[0], 0, 0,
                           &matrix);
       AssertThrow (ierr == 0, ExcPETScError(ierr));