]> https://gitweb.dealii.org/ - dealii.git/commitdiff
fix reinit() of block matrix
authorTimo Heister <timo.heister@gmail.com>
Fri, 3 May 2013 01:26:47 +0000 (01:26 +0000)
committerTimo Heister <timo.heister@gmail.com>
Fri, 3 May 2013 01:26:47 +0000 (01:26 +0000)
git-svn-id: https://svn.dealii.org/branches/branch_unify_linear_algebra@29438 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/source/lac/petsc_parallel_block_sparse_matrix.cc

index 03c4d5628364b3252f97cbe6f93e0b2a8e33957f..69539924e0aea210a7deeac0dde25d06bff4ab20 100644 (file)
@@ -82,12 +82,12 @@ namespace PETScWrappers
           bcsp.n_block_cols());
 
       std::vector<unsigned int> row_sizes;
-      for (unsigned int r=0; r<this->n_block_rows(); ++r)
+      for (unsigned int r=0; r<bcsp.n_block_rows(); ++r)
               row_sizes.push_back( bcsp.block(r,0).n_rows() );
       this->row_block_indices.reinit (row_sizes);
 
       std::vector<unsigned int> col_sizes;
-      for (unsigned int c=0; c<this->n_block_cols(); ++c)
+      for (unsigned int c=0; c<bcsp.n_block_cols(); ++c)
               col_sizes.push_back( bcsp.block(0,c).n_cols() );
       this->column_block_indices.reinit (col_sizes);
 

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.