]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
Revert r25212 and 25217 for the moment since it doesn't compile on the autotester.
authorbangerth <bangerth@0785d39b-7218-0410-832d-ea1e28bc413d>
Thu, 8 Mar 2012 19:04:03 +0000 (19:04 +0000)
committerbangerth <bangerth@0785d39b-7218-0410-832d-ea1e28bc413d>
Thu, 8 Mar 2012 19:04:03 +0000 (19:04 +0000)
git-svn-id: https://svn.dealii.org/trunk@25223 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/source/lac/petsc_parallel_sparse_matrix.cc

index 8681dc0438a0853b297ceb5490543761fffaaf28..8c5838729ce90a60ae87021e2cbc706b320c914a 100644 (file)
@@ -193,26 +193,12 @@ namespace PETScWrappers
                                        // use the call sequence indicating only
                                        // a maximal number of elements per row
                                        // for all rows globally
-      int ierr;
-
-#ifdef DEAL_II_PETSC_VERSION_DEV
-      ierr
-        = MatCreateAIJ (communicator,
-                       local_rows, local_columns,
-                       m, n,
-                       n_nonzero_per_row, 0, 0, 0,
-                       &matrix);
-      AssertThrow (ierr == 0, ExcPETScError(ierr));
-
-      ierr = MatSetOption (matrix, MAT_NEW_NONZERO_ALLOCATION_ERR, PETSC_FALSE);
-#else
-      ierr
-        = MatCreateMPIAIJ (communicator,
-                          local_rows, local_columns,
-                          m, n,
-                          n_nonzero_per_row, 0, 0, 0,
-                          &matrix);
-#endif
+      const int ierr
+        = MatCreateMPIAIJ(communicator,
+                         local_rows, local_columns,
+                         m, n,
+                          n_nonzero_per_row, 0, 0, 0,
+                          &matrix);
       AssertThrow (ierr == 0, ExcPETScError(ierr));
 
                                        // set symmetric flag, if so requested
@@ -248,7 +234,7 @@ namespace PETScWrappers
                                       // For the case that
                                       // local_columns is smaller
                                       // than one of the row lengths
-                                      // MatCreate(MPI)AIJ throws an
+                                      // MatCreateMPIAIJ throws an
                                       // error. In this case use a
                                       // PETScWrappers::SparseMatrix
       for (unsigned int i=0; i<row_lengths.size(); ++i)
@@ -273,27 +259,12 @@ namespace PETScWrappers
 
 //TODO: There must be a significantly better way to provide information about the off-diagonal blocks of the matrix. this way, petsc keeps allocating tiny chunks of memory, and gets completely hung up over this
 
-      int ierr;
-
-#ifdef DEAL_II_PETSC_VERSION_DEV
-      ierr
-        = MatCreateAIJ (communicator,
-                       local_rows, local_columns,
-                       m, n,
-                       0, &int_row_lengths[0], 0, 0,
-                       &matrix);
-      AssertThrow (ierr == 0, ExcPETScError(ierr));
-
-//TODO: Sometimes the actual number of nonzero entries allocated is greater than the number of nonzero entries, which petsc will complain about. There is probably a way to prevent a different number nonzero elements being allocated in the first place. (See also previous TODO).
-      ierr = MatSetOption (matrix, MAT_NEW_NONZERO_ALLOCATION_ERR, PETSC_FALSE);
-#else
-      ierr
-        = MatCreateMPIAIJ (communicator,
-                          local_rows, local_columns,
-                          m, n,
-                          0, &int_row_lengths[0], 0, 0,
-                          &matrix);
-#endif
+      const int ierr
+        = MatCreateMPIAIJ(communicator,
+                         local_rows, local_columns,
+                         m, n,
+                          0, &int_row_lengths[0], 0, 0,
+                          &matrix);
       AssertThrow (ierr == 0, ExcPETScError(ierr));
 
                                        // set symmetric flag, if so requested

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.