]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
Fix a bug.
authorwolf <wolf@0785d39b-7218-0410-832d-ea1e28bc413d>
Mon, 10 May 2004 19:39:42 +0000 (19:39 +0000)
committerwolf <wolf@0785d39b-7218-0410-832d-ea1e28bc413d>
Mon, 10 May 2004 19:39:42 +0000 (19:39 +0000)
git-svn-id: https://svn.dealii.org/trunk@9229 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/lac/source/petsc_parallel_sparse_matrix.cc

index a445c2ab24a917a5eb69198bde351485d1b4290b..8b8d75457376fd3325e847b4a659fba369cc9a90 100644 (file)
@@ -115,15 +115,9 @@ namespace PETScWrappers
                                        // use the call sequence indicating only
                                        // a maximal number of elements per row
                                        // for all rows globally
-                                       //
-                                       // 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
+//TODO[WB]: We should do better by providing ways to tell PETSc how to partition the columns
       const int ierr
-        = MatCreateMPIAIJ(communicator, local_rows, local_rows, m, n,
+        = MatCreateMPIAIJ(communicator, local_rows, PETSC_DECIDE, m, n,
                           n_nonzero_per_row, 0, 0, 0,
                           &matrix);
       AssertThrow (ierr == 0, ExcPETScError(ierr));
@@ -160,14 +154,9 @@ namespace PETScWrappers
       const std::vector<signed int> int_row_lengths (row_lengths.begin(),
                                                      row_lengths.end());
                                        
-                                       // 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
+//TODO[WB]: We should do better by providing ways to tell PETSc how to partition the columns
       const int ierr
-        = MatCreateMPIAIJ(communicator, local_rows, local_rows, m, n,
+        = MatCreateMPIAIJ(communicator, local_rows, PETSC_DECIDE, m, n,
                           0, &int_row_lengths[0], 0, 0,
                           &matrix);
       AssertThrow (ierr == 0, ExcPETScError(ierr));

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.