From e067a5006ba48ef0a0b78ee2f73cfb9dbf4a7184 Mon Sep 17 00:00:00 2001
From: guido <guido@0785d39b-7218-0410-832d-ea1e28bc413d>
Date: Thu, 22 Nov 2001 15:52:28 +0000
Subject: [PATCH] Bug with block vectors removed

git-svn-id: https://svn.dealii.org/trunk@5242 0785d39b-7218-0410-832d-ea1e28bc413d
---
 deal.II/lac/include/lac/filtered_matrix.templates.h | 12 ++++++++++--
 1 file changed, 10 insertions(+), 2 deletions(-)

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 ();
 };
 
-- 
2.39.5