From: Guido Kanschat <dr.guido.kanschat@gmail.com>
Date: Thu, 22 Nov 2001 15:52:28 +0000 (+0000)
Subject: Bug with block vectors removed
X-Git-Tag: v8.0.0~18591
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=c0325e5145317ecc3af98a66e82e75caf661835a;p=dealii.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<BlockSparseMatrix<double>,BlockVector<double> >::
 allocate_tmp_vector () 
 {
+  std::vector<unsigned int> block_sizes (matrix->n_block_rows());
+  for (unsigned int i=0; i<block_sizes.size(); ++i)
+    block_sizes[i] = matrix->block(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<BlockSparseMatrix<float>,BlockVector<float> >::
 allocate_tmp_vector () 
 {
+  std::vector<unsigned int> block_sizes (matrix->n_block_rows());
+  for (unsigned int i=0; i<block_sizes.size(); ++i)
+    block_sizes[i] = matrix->block(i,i).n();
+  
   tmp_mutex.acquire ();
-  tmp_vector.reinit (matrix->n());
+  tmp_vector.reinit (block_sizes);
   tmp_mutex.release ();
 };