]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Clear deletes matrix; initialize uses clear
authorGuido Kanschat <dr.guido.kanschat@gmail.com>
Mon, 19 Jun 2000 18:54:01 +0000 (18:54 +0000)
committerGuido Kanschat <dr.guido.kanschat@gmail.com>
Mon, 19 Jun 2000 18:54:01 +0000 (18:54 +0000)
git-svn-id: https://svn.dealii.org/trunk@3033 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/lac/include/lac/precondition_block.templates.h

index 04ac0056629d9bac29da768c6e3d8d2ee8c54d1b..69f75a238f8241086df2df501f28d16d38377fce 100644 (file)
@@ -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()));

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.