Remove the SparseLUDecomposition::decomposed variable.

It is only ever written to, never read.

diff --git a/include/deal.II/lac/sparse_decomposition.h b/include/deal.II/lac/sparse_decomposition.h
index c08f1c8f4d424f22c63c8fd2186e27882063f48f..9ae6d10299072ddc991a0cd7a07e3ccd11966081 100644 (file)
--- a/include/deal.II/lac/sparse_decomposition.h
+++ b/include/deal.II/lac/sparse_decomposition.h
@@ -287,13 +287,6 @@ protected:
    */
   virtual number get_strengthen_diagonal(const number rowsum, const size_type row) const;
 
-  /**
-   * State flag. If not in decomposed state, it is illegal to apply the
-   * decomposition.  This flag is cleared when the underlaying SparseMatrix
-   * SparsityPattern is changed, and set by decompose().
-   */
-  bool decomposed;
-
   /**
    * The default strengthening value, returned by get_strengthen_diagonal().
    */

diff --git a/include/deal.II/lac/sparse_decomposition.templates.h b/include/deal.II/lac/sparse_decomposition.templates.h
index bcf772404bc163c402a56c23b18654d7a9442e6e..7baefcc11f8314359856056cf63089af035dd67e 100644 (file)
--- a/include/deal.II/lac/sparse_decomposition.templates.h
+++ b/include/deal.II/lac/sparse_decomposition.templates.h
@@ -31,7 +31,6 @@ template<typename number>
 SparseLUDecomposition<number>::SparseLUDecomposition()
   :
   SparseMatrix<number>(),
-  decomposed(false),
   own_sparsity(0)
 {}
 
@@ -47,8 +46,6 @@ SparseLUDecomposition<number>::~SparseLUDecomposition()
 template<typename number>
 void SparseLUDecomposition<number>::clear()
 {
-  decomposed = false;
-
   std::vector<const size_type *> tmp;
   tmp.swap (prebuilt_lower_bound);
 
@@ -122,7 +119,6 @@ void SparseLUDecomposition<number>::initialize (
   // now use this sparsity pattern
   Assert (sparsity_pattern_to_use->n_rows()==sparsity_pattern_to_use->n_cols(),
           typename SparsityPattern::ExcDiagonalNotOptimized());
-  decomposed = false;
   {
     std::vector<const size_type *> tmp;
     tmp.swap (prebuilt_lower_bound);

diff --git a/include/deal.II/lac/sparse_ilu.templates.h b/include/deal.II/lac/sparse_ilu.templates.h
index 4b6820acf437edca6d54824ef546066d04859286..bb6e47871b62b21ce109432dba151a8e7a8e3ee3 100644 (file)
--- a/include/deal.II/lac/sparse_ilu.templates.h
+++ b/include/deal.II/lac/sparse_ilu.templates.h
@@ -48,11 +48,9 @@ void SparseILU<number>::initialize (const SparseMatrix<somenumber> &matrix,
   Assert (data.strengthen_diagonal>=0,
           ExcInvalidStrengthening (data.strengthen_diagonal));
 
-  this->decomposed = false;
   this->strengthen_diagonal = data.strengthen_diagonal;
   this->prebuild_lower_bound ();
   this->copy_from (matrix);
-  this->decomposed = true;
 
   if (data.strengthen_diagonal>0)
     this->strengthen_diagonal_impl();

diff --git a/include/deal.II/lac/sparse_mic.templates.h b/include/deal.II/lac/sparse_mic.templates.h
index ae77a0e9588e1687077293d853b50f44c757429d..b097ab64e5155e3d57db4d6f70e75ee7805fecad 100644 (file)
--- a/include/deal.II/lac/sparse_mic.templates.h
+++ b/include/deal.II/lac/sparse_mic.templates.h
@@ -74,11 +74,9 @@ void SparseMIC<number>::initialize (const SparseMatrix<somenumber> &matrix,
 
   SparseLUDecomposition<number>::initialize(matrix, data);
 
-  this->decomposed = false;
   this->strengthen_diagonal = data.strengthen_diagonal;
   this->prebuild_lower_bound ();
   this->copy_from (matrix);
-  this->decomposed = true;
 
   if (data.strengthen_diagonal > 0)
     this->strengthen_diagonal_impl ();