From: guido Date: Thu, 22 Nov 2001 15:52:28 +0000 (+0000) Subject: Bug with block vectors removed X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=e067a5006ba48ef0a0b78ee2f73cfb9dbf4a7184;p=dealii-svn.git Bug with block vectors removed git-svn-id: https://svn.dealii.org/trunk@5242 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/deal.II/lac/include/lac/filtered_matrix.templates.h b/deal.II/lac/include/lac/filtered_matrix.templates.h index d26f3b59f3..15344f1305 100644 --- a/deal.II/lac/include/lac/filtered_matrix.templates.h +++ b/deal.II/lac/include/lac/filtered_matrix.templates.h @@ -404,8 +404,12 @@ void FilteredMatrix,BlockVector >:: allocate_tmp_vector () { + std::vector block_sizes (matrix->n_block_rows()); + for (unsigned int i=0; iblock(i,i).n(); + tmp_mutex.acquire (); - tmp_vector.reinit (matrix->n()); + tmp_vector.reinit (block_sizes); tmp_mutex.release (); }; @@ -416,8 +420,12 @@ void FilteredMatrix,BlockVector >:: allocate_tmp_vector () { + std::vector block_sizes (matrix->n_block_rows()); + for (unsigned int i=0; iblock(i,i).n(); + tmp_mutex.acquire (); - tmp_vector.reinit (matrix->n()); + tmp_vector.reinit (block_sizes); tmp_mutex.release (); };