From 35a065efd1a67a36302732007ecce8ac35c9e9f7 Mon Sep 17 00:00:00 2001
From: Guido Kanschat <dr.guido.kanschat@gmail.com>
Date: Mon, 19 Jun 2000 18:54:01 +0000
Subject: [PATCH] Clear deletes matrix; initialize uses clear

git-svn-id: https://svn.dealii.org/trunk@3033 0785d39b-7218-0410-832d-ea1e28bc413d
---
 deal.II/lac/include/lac/precondition_block.templates.h | 8 +++++---
 1 file changed, 5 insertions(+), 3 deletions(-)

diff --git a/deal.II/lac/include/lac/precondition_block.templates.h b/deal.II/lac/include/lac/precondition_block.templates.h
index 04ac005662..69f75a238f 100644
--- a/deal.II/lac/include/lac/precondition_block.templates.h
+++ b/deal.II/lac/include/lac/precondition_block.templates.h
@@ -43,14 +43,16 @@ void PreconditionBlock<number,inverse_type>::clear ()
   if (_inverse.size()!=0)
     _inverse.erase(_inverse.begin(), _inverse.end());
   blocksize     = 0;
-  same_diagonal = false;  
+  same_diagonal = false;
+  A = 0;
 }
 
 
 template <typename number, typename inverse_type>
-void PreconditionBlock<number,inverse_type>
-::initialize (const SparseMatrix<number> &M, unsigned int bsize)
+void PreconditionBlock<number,inverse_type>::initialize (const SparseMatrix<number> &M,
+							 unsigned int bsize)
 {
+  clear();
   Assert (M.m() == M.n(), ExcMatrixNotSquare());
   A = &M;
   Assert (bsize>0, ExcIndexRange(bsize, 1, M.m()));
-- 
2.39.5