]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Consistently format PETSc/SLEPc error codes.
authorDavid Wells <wellsd2@rpi.edu>
Thu, 17 Nov 2016 00:41:08 +0000 (19:41 -0500)
committerDavid Wells <wellsd2@rpi.edu>
Thu, 17 Nov 2016 15:35:04 +0000 (10:35 -0500)
12 files changed:
include/deal.II/lac/petsc_matrix_base.h
include/deal.II/lac/petsc_vector_base.h
source/lac/petsc_full_matrix.cc
source/lac/petsc_matrix_base.cc
source/lac/petsc_matrix_free.cc
source/lac/petsc_parallel_sparse_matrix.cc
source/lac/petsc_parallel_vector.cc
source/lac/petsc_solver.cc
source/lac/petsc_sparse_matrix.cc
source/lac/petsc_vector_base.cc
source/lac/slepc_solver.cc
source/lac/slepc_spectral_transformation.cc

index 0dd365bfff22683643c6f36ec39897a7e26ed3fc..3c3e7e1824d4768723dac7efdec122de23c457f8 100644 (file)
@@ -1256,9 +1256,9 @@ namespace PETScWrappers
         col_value_ptr = &column_values[0];
       }
 
-    const int ierr
-      = MatSetValues (matrix, 1, &petsc_i, n_columns, col_index_ptr,
-                      col_value_ptr, INSERT_VALUES);
+    const int ierr = MatSetValues (matrix, 1, &petsc_i, n_columns,
+                                   col_index_ptr,
+                                   col_value_ptr, INSERT_VALUES);
     AssertThrow (ierr == 0, ExcPETScError(ierr));
   }
 
@@ -1398,9 +1398,8 @@ namespace PETScWrappers
         col_value_ptr = &column_values[0];
       }
 
-    const int ierr
-      = MatSetValues (matrix, 1, &petsc_i, n_columns, col_index_ptr,
-                      col_value_ptr, ADD_VALUES);
+    const int ierr = MatSetValues (matrix, 1, &petsc_i, n_columns,
+                                   col_index_ptr, col_value_ptr, ADD_VALUES);
     AssertThrow (ierr == 0, ExcPETScError(ierr));
   }
 
index 62dbbc04db9b54721e012c55ed14fd12a7874d9f..066700bcc685f84cf992d61ee9a7364db2790f70 100644 (file)
@@ -878,8 +878,7 @@ namespace PETScWrappers
 
       const PetscInt petsc_i = index;
 
-      const int ierr
-        = VecSetValues (vector, 1, &petsc_i, &value, INSERT_VALUES);
+      const int ierr = VecSetValues (vector, 1, &petsc_i, &value, INSERT_VALUES);
       AssertThrow (ierr == 0, ExcPETScError(ierr));
 
       vector.last_action = VectorOperation::insert;
@@ -915,8 +914,7 @@ namespace PETScWrappers
 
       // use the PETSc function to add something
       const PetscInt petsc_i = index;
-      const int ierr
-        = VecSetValues (vector, 1, &petsc_i, &value, ADD_VALUES);
+      const int ierr = VecSetValues (vector, 1, &petsc_i, &value, ADD_VALUES);
       AssertThrow (ierr == 0, ExcPETScError(ierr));
 
 
@@ -953,8 +951,7 @@ namespace PETScWrappers
       // add something
       const PetscInt petsc_i = index;
       const PetscScalar subtractand = -value;
-      const int ierr
-        = VecSetValues (vector, 1, &petsc_i, &subtractand, ADD_VALUES);
+      const int ierr = VecSetValues (vector, 1, &petsc_i, &subtractand, ADD_VALUES);
       AssertThrow (ierr == 0, ExcPETScError(ierr));
 
       return *this;
@@ -990,8 +987,7 @@ namespace PETScWrappers
       const PetscScalar new_value
         = static_cast<PetscScalar>(*this) * value;
 
-      const int ierr
-        = VecSetValues (vector, 1, &petsc_i, &new_value, INSERT_VALUES);
+      const int ierr = VecSetValues (vector, 1, &petsc_i, &new_value, INSERT_VALUES);
       AssertThrow (ierr == 0, ExcPETScError(ierr));
 
       return *this;
@@ -1027,8 +1023,7 @@ namespace PETScWrappers
       const PetscScalar new_value
         = static_cast<PetscScalar>(*this) / value;
 
-      const int ierr
-        = VecSetValues (vector, 1, &petsc_i, &new_value, INSERT_VALUES);
+      const int ierr = VecSetValues (vector, 1, &petsc_i, &new_value, INSERT_VALUES);
       AssertThrow (ierr == 0, ExcPETScError(ierr));
 
       return *this;
index 9d2201c157a1cfbac945ed3d45f270c27d664cbd..334a34498c4ad30af5dce96e486444d261782350 100644 (file)
@@ -55,10 +55,8 @@ namespace PETScWrappers
   {
     // use the call sequence indicating only a maximal number of
     // elements per row for all rows globally
-    const int ierr
-      = MatCreateSeqDense (PETSC_COMM_SELF, m, n, PETSC_NULL,
-                           &matrix);
-
+    const int ierr = MatCreateSeqDense (PETSC_COMM_SELF, m, n, PETSC_NULL,
+                                        &matrix);
     AssertThrow (ierr == 0, ExcPETScError(ierr));
   }
 
index 6b1255cb0f1fe80c84d3cbb94f2fd0b05257a795..a94467b6c652d086f9292a917f61d5c08a8281c6 100644 (file)
@@ -158,11 +158,10 @@ namespace PETScWrappers
 #endif
 
 #if DEAL_II_PETSC_VERSION_LT(3,2,0)
-    const int ierr
-      = MatZeroRowsIS(matrix, index_set, new_diag_value);
+    const int ierr = MatZeroRowsIS(matrix, index_set, new_diag_value);
 #else
-    const int ierr
-      = MatZeroRowsIS(matrix, index_set, new_diag_value, PETSC_NULL, PETSC_NULL);
+    const int ierr = MatZeroRowsIS(matrix, index_set, new_diag_value, PETSC_NULL,
+                                   PETSC_NULL);
 #endif
     AssertThrow (ierr == 0, ExcPETScError(ierr));
 
@@ -183,9 +182,8 @@ namespace PETScWrappers
 
     PetscScalar value;
 
-    const int ierr
-      = MatGetValues (matrix, 1, &petsc_i, 1, &petsc_j,
-                      &value);
+    const int ierr = MatGetValues (matrix, 1, &petsc_i, 1, &petsc_j,
+                                   &value);
     AssertThrow (ierr == 0, ExcPETScError(ierr));
 
     return value;
@@ -247,7 +245,7 @@ namespace PETScWrappers
   {
     PetscInt n_rows, n_cols;
 
-    int ierr = MatGetSize (matrix, &n_rows, &n_cols);
+    const int ierr = MatGetSize (matrix, &n_rows, &n_cols);
     AssertThrow (ierr == 0, ExcPETScError(ierr));
 
     return n_rows;
@@ -260,7 +258,7 @@ namespace PETScWrappers
   {
     PetscInt n_rows, n_cols;
 
-    int ierr = MatGetSize (matrix, &n_rows, &n_cols);
+    const int ierr = MatGetSize (matrix, &n_rows, &n_cols);
     AssertThrow (ierr == 0, ExcPETScError(ierr));
 
     return n_cols;
@@ -273,7 +271,7 @@ namespace PETScWrappers
   {
     PetscInt n_rows, n_cols;
 
-    int ierr = MatGetLocalSize (matrix, &n_rows, &n_cols);
+    const int ierr = MatGetLocalSize (matrix, &n_rows, &n_cols);
     AssertThrow (ierr == 0, ExcPETScError(ierr));
 
     return n_rows;
@@ -299,8 +297,7 @@ namespace PETScWrappers
   MatrixBase::n_nonzero_elements () const
   {
     MatInfo mat_info;
-    const int ierr
-      = MatGetInfo (matrix, MAT_GLOBAL_SUM, &mat_info);
+    const int ierr = MatGetInfo (matrix, MAT_GLOBAL_SUM, &mat_info);
     AssertThrow (ierr == 0, ExcPETScError(ierr));
 
     return static_cast<size_type>(mat_info.nz_used);
@@ -328,8 +325,7 @@ namespace PETScWrappers
 
 //TODO: this is probably horribly inefficient; we should lobby for a way to
 //query this information from PETSc
-    int ierr;
-    ierr = MatGetRow(*this, row, &ncols, &colnums, &values);
+    int ierr = MatGetRow(*this, row, &ncols, &colnums, &values);
     AssertThrow (ierr == 0, ExcPETScError(ierr));
 
     // then restore the matrix and return the number of columns in this row as
@@ -350,8 +346,7 @@ namespace PETScWrappers
   {
     PetscReal result;
 
-    const int ierr
-      = MatNorm (matrix, NORM_1, &result);
+    const int ierr = MatNorm (matrix, NORM_1, &result);
     AssertThrow (ierr == 0, ExcPETScError(ierr));
 
     return result;
@@ -364,8 +359,7 @@ namespace PETScWrappers
   {
     PetscReal result;
 
-    const int ierr
-      = MatNorm (matrix, NORM_INFINITY, &result);
+    const int ierr = MatNorm (matrix, NORM_INFINITY, &result);
     AssertThrow (ierr == 0, ExcPETScError(ierr));
 
     return result;
@@ -378,8 +372,7 @@ namespace PETScWrappers
   {
     PetscReal result;
 
-    const int ierr
-      = MatNorm (matrix, NORM_FROBENIUS, &result);
+    const int ierr = MatNorm (matrix, NORM_FROBENIUS, &result);
     AssertThrow (ierr == 0, ExcPETScError(ierr));
 
     return result;
@@ -411,8 +404,7 @@ namespace PETScWrappers
   {
     PetscScalar result;
 
-    const int ierr
-      = MatGetTrace (matrix, &result);
+    const int ierr = MatGetTrace (matrix, &result);
     AssertThrow (ierr == 0, ExcPETScError(ierr));
 
     return result;
@@ -437,7 +429,6 @@ namespace PETScWrappers
   {
     const PetscScalar factor = 1./a;
     const int ierr = MatScale (matrix, factor);
-
     AssertThrow (ierr == 0, ExcPETScError(ierr));
 
     return *this;
@@ -450,9 +441,7 @@ namespace PETScWrappers
   {
     const int ierr = MatAXPY (matrix, factor,
                               other, DIFFERENT_NONZERO_PATTERN);
-    (void)ierr;
-
-    Assert (ierr == 0, ExcPETScError(ierr));
+    AssertThrow (ierr == 0, ExcPETScError(ierr));
 
     return *this;
   }
@@ -474,7 +463,6 @@ namespace PETScWrappers
     Assert (&src != &dst, ExcSourceEqualsDestination());
 
     const int ierr = MatMult (matrix, src, dst);
-    (void)ierr;
     AssertThrow (ierr == 0, ExcPETScError(ierr));
   }
 
@@ -487,7 +475,6 @@ namespace PETScWrappers
     Assert (&src != &dst, ExcSourceEqualsDestination());
 
     const int ierr = MatMultTranspose (matrix, src, dst);
-    (void)ierr;
     AssertThrow (ierr == 0, ExcPETScError(ierr));
   }
 
@@ -513,7 +500,6 @@ namespace PETScWrappers
     Assert (&src != &dst, ExcSourceEqualsDestination());
 
     const int ierr = MatMultTransposeAdd (matrix, src, dst, dst);
-    (void)ierr;
     AssertThrow (ierr == 0, ExcPETScError(ierr));
   }
 
@@ -543,8 +529,7 @@ namespace PETScWrappers
   void
   MatrixBase::transpose ()
   {
-    int ierr = MatTranspose(matrix, MAT_REUSE_MATRIX, &matrix);
-    (void)ierr;
+    const int ierr = MatTranspose(matrix, MAT_REUSE_MATRIX, &matrix);
     AssertThrow (ierr == 0, ExcPETScError(ierr));
   }
 
@@ -553,8 +538,7 @@ namespace PETScWrappers
   {
     PetscBooleanType truth;
     assert_is_compressed ();
-    int ierr = MatIsSymmetric (matrix, tolerance, &truth);
-    (void)ierr;
+    const int ierr = MatIsSymmetric (matrix, tolerance, &truth);
     AssertThrow (ierr == 0, ExcPETScError(ierr));
     return truth;
   }
@@ -565,8 +549,7 @@ namespace PETScWrappers
     PetscBooleanType truth;
 
     assert_is_compressed ();
-    int ierr = MatIsHermitian (matrix, tolerance, &truth);
-    (void)ierr;
+    const int ierr = MatIsHermitian (matrix, tolerance, &truth);
     AssertThrow (ierr == 0, ExcPETScError(ierr));
 
     return truth;
@@ -600,7 +583,6 @@ namespace PETScWrappers
     for (row = loc_range.first; row < loc_range.second; ++row)
       {
         int ierr = MatGetRow(*this, row, &ncols, &colnums, &values);
-        (void)ierr;
         AssertThrow (ierr == 0, ExcPETScError(ierr));
 
         for (PetscInt col = 0; col < ncols; ++col)
index 3a3e35d68fff10043ce18e91ba5e358dcd9ce827..da16fa776480bc4ffdf853d067c43c1ca7a02827 100644 (file)
@@ -239,7 +239,7 @@ namespace PETScWrappers
     // to the matrix-vector multiplication
     // of this MatrixFree object,
     void  *this_object;
-    int ierr = MatShellGetContext (A, &this_object);
+    const int ierr = MatShellGetContext (A, &this_object);
     AssertThrow (ierr == 0, ExcPETScError(ierr));
 
     // call vmult of this object:
index 2d0e365e119173a46ed59b7470b78bc7d516512b..9489b8b5d61a41767dc65b159f37a0cbaf9d2196 100644 (file)
@@ -132,7 +132,7 @@ namespace PETScWrappers
 
       this->communicator = other.communicator;
 
-      int ierr = MatCopy(other.matrix, matrix, SAME_NONZERO_PATTERN);
+      const int ierr = MatCopy(other.matrix, matrix, SAME_NONZERO_PATTERN);
       AssertThrow (ierr == 0, ExcPETScError(ierr));
     }
 
@@ -239,23 +239,19 @@ namespace PETScWrappers
       int ierr;
 
 #if DEAL_II_PETSC_VERSION_LT(3,3,0)
-      ierr
-        = MatCreateMPIAIJ (communicator,
+      ierr = MatCreateMPIAIJ (communicator,
+                              local_rows, local_columns,
+                              m, n,
+                              n_nonzero_per_row, 0,
+                              n_offdiag_nonzero_per_row, 0,
+                              &matrix);
+#else
+      ierr = MatCreateAIJ (communicator,
                            local_rows, local_columns,
                            m, n,
                            n_nonzero_per_row, 0,
                            n_offdiag_nonzero_per_row, 0,
                            &matrix);
-#else
-      ierr
-        = MatCreateAIJ (communicator,
-                        local_rows, local_columns,
-                        m, n,
-                        n_nonzero_per_row, 0,
-                        n_offdiag_nonzero_per_row, 0,
-                        &matrix);
-      AssertThrow (ierr == 0, ExcPETScError(ierr));
-
       set_matrix_option (matrix, MAT_NEW_NONZERO_ALLOCATION_ERR, PETSC_FALSE);
 #endif
       AssertThrow (ierr == 0, ExcPETScError(ierr));
@@ -311,22 +307,19 @@ namespace PETScWrappers
       int ierr;
 
 #if DEAL_II_PETSC_VERSION_LT(3,3,0)
-      ierr
-        = MatCreateMPIAIJ (communicator,
+      ierr = MatCreateMPIAIJ (communicator,
+                              local_rows, local_columns,
+                              m, n,
+                              0, &int_row_lengths[0],
+                              0, offdiag_row_lengths.size() ? &int_offdiag_row_lengths[0] : 0,
+                              &matrix);
+#else
+      ierr = MatCreateAIJ (communicator,
                            local_rows, local_columns,
                            m, n,
                            0, &int_row_lengths[0],
                            0, offdiag_row_lengths.size() ? &int_offdiag_row_lengths[0] : 0,
                            &matrix);
-#else
-      ierr
-        = MatCreateAIJ (communicator,
-                        local_rows, local_columns,
-                        m, n,
-                        0, &int_row_lengths[0],
-                        0, offdiag_row_lengths.size() ? &int_offdiag_row_lengths[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 unless explicitly disabled with MatSetOption. There is probably a way to prevent a different number nonzero elements being allocated in the first place. (See also previous TODO).
       set_matrix_option (matrix, MAT_NEW_NONZERO_ALLOCATION_ERR, PETSC_FALSE);
@@ -456,20 +449,20 @@ namespace PETScWrappers
           // then call the petsc function
           // that summarily allocates these
           // entries:
-          MatMPIAIJSetPreallocationCSR (matrix,
-                                        &rowstart_in_window[0],
-                                        &colnums_in_window[0],
-                                        0);
+          ierr = MatMPIAIJSetPreallocationCSR (matrix,
+                                               &rowstart_in_window[0],
+                                               &colnums_in_window[0],
+                                               0);
+          AssertThrow (ierr == 0, ExcPETScError(ierr));
         }
       else
         {
           PetscInt i=0;
-          MatMPIAIJSetPreallocationCSR (matrix,
-                                        &i,
-                                        &i,
-                                        0);
-
-
+          ierr = MatMPIAIJSetPreallocationCSR (matrix,
+                                               &i,
+                                               &i,
+                                               0);
+          AssertThrow (ierr == 0, ExcPETScError(ierr));
         }
       compress (dealii::VectorOperation::insert);
 
@@ -548,24 +541,21 @@ namespace PETScWrappers
       // the first _local_ row, i.e. it
       // doesn't index into an array for
       // _all_ rows.
-      const int ierr
-        = MatCreateMPIAIJ(communicator,
-                          local_rows_per_process[this_process],
-                          local_columns_per_process[this_process],
-                          sparsity_pattern.n_rows(),
-                          sparsity_pattern.n_cols(),
-                          0, &row_lengths_in_window[0],
-                          0, &row_lengths_out_of_window[0],
-                          &matrix);
+      const int ierr = MatCreateMPIAIJ(communicator,
+                                       local_rows_per_process[this_process],
+                                       local_columns_per_process[this_process],
+                                       sparsity_pattern.n_rows(),
+                                       sparsity_pattern.n_cols(),
+                                       0, &row_lengths_in_window[0],
+                                       0, &row_lengths_out_of_window[0],
+                                       &matrix);
       AssertThrow (ierr == 0, ExcPETScError(ierr));
 
 #else //PETSC_VERSION>=2.3.3
       // create the matrix. We
       // do not set row length but set the
       // correct SparsityPattern later.
-      int ierr;
-
-      ierr = MatCreate(communicator,&matrix);
+      int ierr = MatCreate(communicator,&matrix);
       AssertThrow (ierr == 0, ExcPETScError(ierr));
 
       ierr = MatSetSizes(matrix,
@@ -636,10 +626,11 @@ namespace PETScWrappers
           // then call the petsc function
           // that summarily allocates these
           // entries:
-          MatMPIAIJSetPreallocationCSR (matrix,
-                                        &rowstart_in_window[0],
-                                        &colnums_in_window[0],
-                                        0);
+          ierr = MatMPIAIJSetPreallocationCSR (matrix,
+                                               &rowstart_in_window[0],
+                                               &colnums_in_window[0],
+                                               0);
+          AssertThrow (ierr == 0, ExcPETScError(ierr));
 
 #if DEAL_II_PETSC_VERSION_LT(2,3,3)
           // this is only needed for old
index 50bc5f57f6c5b874b1468a48267781e61c5d019d..71f6fa376a51a41c9d42987b98cfa3724eba79da 100644 (file)
@@ -141,7 +141,6 @@ namespace PETScWrappers
 #else
           ierr = VecDestroy (&vector);
 #endif
-
           AssertThrow (ierr == 0, ExcPETScError(ierr));
 
           create_vector (n, local_sz);
@@ -164,7 +163,7 @@ namespace PETScWrappers
           reinit (v.locally_owned_elements(), v.ghost_indices, v.communicator);
           if (!omit_zeroing_entries)
             {
-              int ierr = VecSet(vector, 0.0);
+              const int ierr = VecSet(vector, 0.0);
               AssertThrow (ierr == 0, ExcPETScError(ierr));
             }
         }
@@ -272,9 +271,8 @@ namespace PETScWrappers
       Assert (local_size <= n, ExcIndexRange (local_size, 0, n));
       ghosted = false;
 
-      const int ierr
-        = VecCreateMPI (communicator, local_size, PETSC_DETERMINE,
-                        &vector);
+      const int ierr = VecCreateMPI (communicator, local_size, PETSC_DETERMINE,
+                                     &vector);
       AssertThrow (ierr == 0, ExcPETScError(ierr));
 
       Assert (size() == n,
@@ -303,14 +301,12 @@ namespace PETScWrappers
            :
            0);
 
-      int ierr
-        = VecCreateGhost(communicator,
-                         local_size,
-                         PETSC_DETERMINE,
-                         ghostindices.size(),
-                         ptr,
-                         &vector);
-
+      int ierr = VecCreateGhost(communicator,
+                                local_size,
+                                PETSC_DETERMINE,
+                                ghostindices.size(),
+                                ptr,
+                                &vector);
       AssertThrow (ierr == 0, ExcPETScError(ierr));
 
       Assert (size() == n,
@@ -387,15 +383,12 @@ namespace PETScWrappers
       PetscInt    nlocal, istart, iend;
 
       int ierr = VecGetArray (vector, &val);
-
       AssertThrow (ierr == 0, ExcPETScError(ierr));
 
       ierr = VecGetLocalSize (vector, &nlocal);
-
       AssertThrow (ierr == 0, ExcPETScError(ierr));
 
       ierr = VecGetOwnershipRange (vector, &istart, &iend);
-
       AssertThrow (ierr == 0, ExcPETScError(ierr));
 
       // save the state of out stream
index 0acf74c25b5dca7e220f6113dea76ea1292bab3d..a6107b6b84980d82faf3d83f22d4b905f3fa0a6c 100644 (file)
@@ -310,12 +310,10 @@ namespace PETScWrappers
     // completely pointless change in
     // spelling Chebyshev between PETSc 3.2
     // and 3.3...
-    int ierr;
-
 #if DEAL_II_PETSC_VERSION_LT(3,3,0)
-    ierr = KSPSetType (ksp, KSPCHEBYCHEV);
+    int ierr = KSPSetType (ksp, KSPCHEBYCHEV);
 #else
-    ierr = KSPSetType (ksp, KSPCHEBYSHEV);
+    int ierr = KSPSetType (ksp, KSPCHEBYSHEV);
 #endif
     AssertThrow (ierr == 0, ExcPETScError(ierr));
 
@@ -340,8 +338,7 @@ namespace PETScWrappers
   void
   SolverCG::set_solver_type (KSP &ksp) const
   {
-    int ierr;
-    ierr = KSPSetType (ksp, KSPCG);
+    int ierr = KSPSetType (ksp, KSPCG);
     AssertThrow (ierr == 0, ExcPETScError(ierr));
 
     // in the deal.II solvers, we always
@@ -365,8 +362,7 @@ namespace PETScWrappers
   void
   SolverBiCG::set_solver_type (KSP &ksp) const
   {
-    int ierr;
-    ierr = KSPSetType (ksp, KSPBICG);
+    int ierr = KSPSetType (ksp, KSPBICG);
     AssertThrow (ierr == 0, ExcPETScError(ierr));
 
     // in the deal.II solvers, we always
@@ -400,8 +396,7 @@ namespace PETScWrappers
   void
   SolverGMRES::set_solver_type (KSP &ksp) const
   {
-    int ierr;
-    ierr = KSPSetType (ksp, KSPGMRES);
+    int ierr = KSPSetType (ksp, KSPGMRES);
     AssertThrow (ierr == 0, ExcPETScError(ierr));
 
     // set the restart parameter from the
@@ -467,8 +462,7 @@ namespace PETScWrappers
   void
   SolverBicgstab::set_solver_type (KSP &ksp) const
   {
-    int ierr;
-    ierr = KSPSetType (ksp, KSPBCGS);
+    int ierr = KSPSetType (ksp, KSPBCGS);
     AssertThrow (ierr == 0, ExcPETScError(ierr));
 
     // in the deal.II solvers, we always
@@ -492,8 +486,7 @@ namespace PETScWrappers
   void
   SolverCGS::set_solver_type (KSP &ksp) const
   {
-    int ierr;
-    ierr = KSPSetType (ksp, KSPCGS);
+    int ierr = KSPSetType (ksp, KSPCGS);
     AssertThrow (ierr == 0, ExcPETScError(ierr));
 
     // in the deal.II solvers, we always
@@ -542,8 +535,7 @@ namespace PETScWrappers
   void
   SolverTCQMR::set_solver_type (KSP &ksp) const
   {
-    int ierr;
-    ierr = KSPSetType (ksp, KSPTCQMR);
+    int ierr = KSPSetType (ksp, KSPTCQMR);
     AssertThrow (ierr == 0, ExcPETScError(ierr));
 
     // in the deal.II solvers, we always
@@ -592,8 +584,7 @@ namespace PETScWrappers
   void
   SolverLSQR::set_solver_type (KSP &ksp) const
   {
-    int ierr;
-    ierr = KSPSetType (ksp, KSPLSQR);
+    int ierr = KSPSetType (ksp, KSPLSQR);
     AssertThrow (ierr == 0, ExcPETScError(ierr));
 
     // in the deal.II solvers, we always
@@ -664,8 +655,7 @@ namespace PETScWrappers
     * preconditioner.  Its use is due to SparseDirectMUMPS being a direct
     * (rather than iterative) solver
     */
-    int ierr;
-    ierr = KSPSetType (ksp, KSPPREONLY);
+    int ierr = KSPSetType (ksp, KSPPREONLY);
     AssertThrow (ierr == 0, ExcPETScError(ierr));
 
     /**
index 3efacaf7ea0bd36f0080cb64141e7da73b6af343..7acf5e0644ecd7dd9a796010ac60aa8a7045978a 100644 (file)
@@ -31,9 +31,8 @@ namespace PETScWrappers
   SparseMatrix::SparseMatrix ()
   {
     const int m=0, n=0, n_nonzero_per_row=0;
-    const int ierr
-      = MatCreateSeqAIJ(PETSC_COMM_SELF, m, n, n_nonzero_per_row,
-                        0, &matrix);
+    const int ierr = MatCreateSeqAIJ(PETSC_COMM_SELF, m, n, n_nonzero_per_row,
+                                     0, &matrix);
     AssertThrow (ierr == 0, ExcPETScError(ierr));
   }
 
@@ -145,9 +144,9 @@ namespace PETScWrappers
     // use the call sequence indicating only
     // a maximal number of elements per row
     // for all rows globally
-    const int ierr
-      = MatCreateSeqAIJ(PETSC_COMM_SELF, m, n, n_nonzero_per_row,
-                        0, &matrix);
+    const int ierr = MatCreateSeqAIJ(PETSC_COMM_SELF, m, n,
+                                     n_nonzero_per_row,
+                                     0, &matrix);
     AssertThrow (ierr == 0, ExcPETScError(ierr));
 
     // set symmetric flag, if so requested
@@ -179,9 +178,8 @@ namespace PETScWrappers
     const std::vector<PetscInt>
     int_row_lengths (row_lengths.begin(), row_lengths.end());
 
-    const int ierr
-      = MatCreateSeqAIJ(PETSC_COMM_SELF, m, n, 0,
-                        &int_row_lengths[0], &matrix);
+    const int ierr = MatCreateSeqAIJ(PETSC_COMM_SELF, m, n, 0,
+                                     &int_row_lengths[0], &matrix);
     AssertThrow (ierr == 0, ExcPETScError(ierr));
 
     // set symmetric flag, if so requested
@@ -246,7 +244,7 @@ namespace PETScWrappers
   SparseMatrix::m () const
   {
     PetscInt m,n;
-    PetscErrorCode ierr = MatGetSize(matrix, &m, &n);
+    const PetscErrorCode ierr = MatGetSize(matrix, &m, &n);
     AssertThrow (ierr == 0, ExcPETScError(ierr));
 
     return m;
@@ -256,7 +254,7 @@ namespace PETScWrappers
   SparseMatrix::n () const
   {
     PetscInt m,n;
-    PetscErrorCode ierr = MatGetSize(matrix, &m, &n);
+    const PetscErrorCode ierr = MatGetSize(matrix, &m, &n);
     AssertThrow (ierr == 0, ExcPETScError(ierr));
 
     return n;
index f329e78da12638153566c9cde608ac1f3fad41b8..71a2e7e6edf5a893b99248ecea84f0220b124adc 100644 (file)
@@ -1178,9 +1178,8 @@ namespace PETScWrappers
 #endif
 
         const InsertMode mode = (add_values ? ADD_VALUES : INSERT_VALUES);
-        const int ierr
-          = VecSetValues (vector, n_elements, petsc_indices, values,
-                          mode);
+        const int ierr = VecSetValues (vector, n_elements, petsc_indices,
+                                       values, mode);
         AssertThrow (ierr == 0, ExcPETScError(ierr));
       }
 
index 3178d5a30a77c804075e59bbdb6874784dd0f62d..0fdaf0254da830d7dbfc5ed1348ea23281f647e6 100644 (file)
@@ -63,9 +63,9 @@ namespace SLEPcWrappers
       {
         // Destroy the solver object.
 #if DEAL_II_PETSC_VERSION_LT(3,2,0)
-        int ierr = EPSDestroy (eps);
+        const int ierr = EPSDestroy (eps);
 #else
-        int ierr = EPSDestroy (&eps);
+        const int ierr = EPSDestroy (&eps);
 #endif
         AssertThrow (ierr == 0, ExcSLEPcError(ierr));
       }
@@ -75,7 +75,7 @@ namespace SLEPcWrappers
   SolverBase::set_matrices (const PETScWrappers::MatrixBase &A)
   {
     // standard eigenspectrum problem
-    int ierr = EPSSetOperators (eps, A, PETSC_NULL);
+    const int ierr = EPSSetOperators (eps, A, PETSC_NULL);
     AssertThrow (ierr == 0, ExcSLEPcError(ierr));
   }
 
@@ -84,7 +84,7 @@ namespace SLEPcWrappers
                             const PETScWrappers::MatrixBase &B)
   {
     // generalized eigenspectrum problem
-    int ierr = EPSSetOperators (eps, A, B);
+    const int ierr = EPSSetOperators (eps, A, B);
     AssertThrow (ierr == 0, ExcSLEPcError(ierr));
   }
 
@@ -93,7 +93,7 @@ namespace SLEPcWrappers
   {
     // set transformation type if any
     // STSetShift is called inside
-    int ierr = EPSSetST(eps,transformation.st);
+    const int ierr = EPSSetST(eps,transformation.st);
     AssertThrow (ierr == 0, SolverBase::ExcSLEPcError(ierr));
   }
 
@@ -103,12 +103,11 @@ namespace SLEPcWrappers
     Assert(this_initial_vector.l2_norm()>0.0,
            ExcMessage("Initial vector should be nonzero."));
 
-    int ierr;
     Vec vec = this_initial_vector;
 #if DEAL_II_PETSC_VERSION_LT(3,1,0)
-    ierr = EPSSetInitialVector (eps, &vec);
+    const int ierr = EPSSetInitialVector (eps, &vec);
 #else
-    ierr = EPSSetInitialSpace (eps, 1, &vec);
+    const int ierr = EPSSetInitialSpace (eps, 1, &vec);
 #endif
     AssertThrow (ierr == 0, ExcSLEPcError(ierr));
   }
@@ -119,7 +118,7 @@ namespace SLEPcWrappers
     // set target eigenvalues to solve for
     // in all transformation except STSHIFT there is a direct connection between
     // the target and the shift, read more on p41 of SLEPc manual.
-    int ierr = EPSSetTarget (eps, this_target );
+    const int ierr = EPSSetTarget (eps, this_target );
     AssertThrow (ierr == 0, ExcSLEPcError(ierr));
   }
 
@@ -127,14 +126,14 @@ namespace SLEPcWrappers
   SolverBase::set_which_eigenpairs (const EPSWhich eps_which)
   {
     // set which portion of the eigenspectrum to solve for
-    int ierr = EPSSetWhichEigenpairs (eps, eps_which);
+    const int ierr = EPSSetWhichEigenpairs (eps, eps_which);
     AssertThrow (ierr == 0, ExcSLEPcError(ierr));
   }
 
   void
   SolverBase::set_problem_type (const EPSProblemType eps_problem)
   {
-    int ierr = EPSSetProblemType (eps, eps_problem);
+    const int ierr = EPSSetProblemType (eps, eps_problem);
     AssertThrow (ierr == 0, ExcSLEPcError(ierr));
   }
 
@@ -142,11 +141,9 @@ namespace SLEPcWrappers
   SolverBase::solve (const unsigned int  n_eigenpairs,
                      unsigned int       *n_converged)
   {
-    int ierr;
-
     // set number of eigenvectors to compute
-    ierr = EPSSetDimensions (eps, n_eigenpairs,
-                             PETSC_DECIDE, PETSC_DECIDE);
+    int ierr = EPSSetDimensions (eps, n_eigenpairs,
+                                 PETSC_DECIDE, PETSC_DECIDE);
     AssertThrow (ierr == 0, ExcSLEPcError(ierr));
 
     // set the solve options to the eigenvalue problem solver context
@@ -233,9 +230,9 @@ namespace SLEPcWrappers
                              PETScWrappers::VectorBase &eigenvectors)
   {
     // get converged eigenpair
-    int ierr = EPSGetEigenpair (eps, index,
-                                &eigenvalues, PETSC_NULL,
-                                eigenvectors, PETSC_NULL);
+    const int ierr = EPSGetEigenpair (eps, index,
+                                      &eigenvalues, PETSC_NULL,
+                                      eigenvectors, PETSC_NULL);
     AssertThrow (ierr == 0, ExcSLEPcError(ierr));
   }
 
@@ -249,9 +246,9 @@ namespace SLEPcWrappers
   {
 #ifndef PETSC_USE_COMPLEX
     // get converged eigenpair
-    int ierr = EPSGetEigenpair (eps, index,
-                                &real_eigenvalues, &imag_eigenvalues,
-                                real_eigenvectors, imag_eigenvectors);
+    const int ierr = EPSGetEigenpair (eps, index,
+                                      &real_eigenvalues, &imag_eigenvalues,
+                                      real_eigenvectors, imag_eigenvectors);
     AssertThrow (ierr == 0, ExcSLEPcError(ierr));
 #else
     Assert ((false),
@@ -323,7 +320,7 @@ namespace SLEPcWrappers
     SolverBase (cn, mpi_communicator),
     additional_data (data)
   {
-    int ierr = EPSSetType (eps, const_cast<char *>(EPSKRYLOVSCHUR));
+    const int ierr = EPSSetType (eps, const_cast<char *>(EPSKRYLOVSCHUR));
     AssertThrow (ierr == 0, ExcSLEPcError(ierr));
   }
 
index 6e9805a0026df4a78b2ff3a03f7628ed383ebbb5..c10d6cdade55438202fb81554ce16f4eafe3d6b9 100644 (file)
@@ -33,7 +33,7 @@ namespace SLEPcWrappers
 {
   TransformationBase::TransformationBase (const MPI_Comm &mpi_communicator)
   {
-    int ierr = STCreate(mpi_communicator, &st);
+    const int ierr = STCreate(mpi_communicator, &st);
     AssertThrow (ierr == 0, SolverBase::ExcSLEPcError(ierr));
   }
 
@@ -41,14 +41,14 @@ namespace SLEPcWrappers
   {
     if (st!=NULL)
       {
-        int ierr = STDestroy(&st);
+        const int ierr = STDestroy(&st);
         AssertThrow (ierr == 0, SolverBase::ExcSLEPcError(ierr));
       }
   }
 
   void TransformationBase::set_matrix_mode(const STMatMode mode)
   {
-    int ierr = STSetMatMode(st,mode);
+    const int ierr = STSetMatMode(st,mode);
     AssertThrow (ierr == 0, SolverBase::ExcSLEPcError(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.