From 853b48fb5e9f0ff47bee284356bab5495f2ae767 Mon Sep 17 00:00:00 2001
From: Wolfgang Bangerth <bangerth@math.tamu.edu>
Date: Tue, 3 Feb 2015 11:27:31 -0600
Subject: [PATCH] Remove deprecated function SparseMIC::decompose.

---
 doc/news/changes.h                         |  1 +
 include/deal.II/lac/sparse_mic.h           |  8 ---
 include/deal.II/lac/sparse_mic.templates.h | 60 +++++++++-------------
 source/lac/sparse_mic.cc                   | 10 +---
 4 files changed, 27 insertions(+), 52 deletions(-)

diff --git a/doc/news/changes.h b/doc/news/changes.h
index 6a3a881977..6f4719b35a 100644
--- a/doc/news/changes.h
+++ b/doc/news/changes.h
@@ -148,6 +148,7 @@ inconvenience this causes.
   - The deprecated constructor of SparseILU.
   - SparseILU::apply_decomposition.
   - The deprecated constructor of SparseMIC.
+  - SparseMIC::decompose.
   - The compress() functions without argument in the various vector
     classes. You should use the versions with a VectorOperation
     argument instead.
diff --git a/include/deal.II/lac/sparse_mic.h b/include/deal.II/lac/sparse_mic.h
index 80c276b94c..aefa5923c7 100644
--- a/include/deal.II/lac/sparse_mic.h
+++ b/include/deal.II/lac/sparse_mic.h
@@ -102,14 +102,6 @@ public:
   void initialize (const SparseMatrix<somenumber> &matrix,
                    const AdditionalData &parameters = AdditionalData());
 
-  /**
-   * @deprecated This method is deprecated, and left for backward compability.
-   * It will be removed in later versions. Use initialize() instead.
-   */
-  template <typename somenumber>
-  void decompose (const SparseMatrix<somenumber> &matrix,
-                  const double                   strengthen_diagonal=0.) DEAL_II_DEPRECATED;
-
   /**
    * Apply the incomplete decomposition, i.e. do one forward-backward step
    * $dst=(LU)^{-1}src$.
diff --git a/include/deal.II/lac/sparse_mic.templates.h b/include/deal.II/lac/sparse_mic.templates.h
index c5168be157..3ed279ba20 100644
--- a/include/deal.II/lac/sparse_mic.templates.h
+++ b/include/deal.II/lac/sparse_mic.templates.h
@@ -68,47 +68,15 @@ void SparseMIC<number>::initialize (const SparseMatrix<somenumber> &matrix,
 {
   SparseLUDecomposition<number>::initialize(matrix, data);
 
-  decompose(matrix, data.strengthen_diagonal);
-}
-
-
-
-template <typename number>
-void SparseMIC<number>::reinit (const SparsityPattern &sparsity)
-{
-  {
-    std::vector<number> tmp;
-    tmp.swap (diag);
-  }
-  {
-    std::vector<number> tmp;
-    tmp.swap (inv_diag);
-  }
-  {
-    std::vector<number> tmp;
-    tmp.swap (inner_sums);
-  }
-
-  SparseMatrix<number>::reinit(sparsity);
-  this->decomposed = false;
-}
-
-
-
-template <typename number>
-template <typename somenumber>
-void SparseMIC<number>::decompose (const SparseMatrix<somenumber> &matrix,
-                                   const double                    strengthen_diagonal)
-{
-  SparseLUDecomposition<number>::decompose(matrix, strengthen_diagonal);
+  SparseLUDecomposition<number>::decompose(matrix, data.strengthen_diagonal);
 
   Assert (matrix.m()==matrix.n(), ExcNotQuadratic ());
   Assert (this->m()==this->n(),   ExcNotQuadratic ());
   Assert (matrix.m()==this->m(),  ExcDimensionMismatch(matrix.m(), this->m()));
 
-  Assert (strengthen_diagonal>=0, ExcInvalidStrengthening (strengthen_diagonal));
+  Assert (data.strengthen_diagonal>=0, ExcInvalidStrengthening (data.strengthen_diagonal));
 
-  if (strengthen_diagonal > 0)
+  if (data.strengthen_diagonal > 0)
     this->strengthen_diagonal_impl ();
 
   // MIC implementation: (S. Margenov lectures)
@@ -150,6 +118,28 @@ void SparseMIC<number>::decompose (const SparseMatrix<somenumber> &matrix,
 
 
 
+template <typename number>
+void SparseMIC<number>::reinit (const SparsityPattern &sparsity)
+{
+  {
+    std::vector<number> tmp;
+    tmp.swap (diag);
+  }
+  {
+    std::vector<number> tmp;
+    tmp.swap (inv_diag);
+  }
+  {
+    std::vector<number> tmp;
+    tmp.swap (inner_sums);
+  }
+
+  SparseMatrix<number>::reinit(sparsity);
+  this->decomposed = false;
+}
+
+
+
 template <typename number>
 inline number
 SparseMIC<number>::get_rowsum (const size_type row) const
diff --git a/source/lac/sparse_mic.cc b/source/lac/sparse_mic.cc
index 89c1c92d96..eb3e2e6d13 100644
--- a/source/lac/sparse_mic.cc
+++ b/source/lac/sparse_mic.cc
@@ -1,6 +1,6 @@
 // ---------------------------------------------------------------------
 //
-// Copyright (C) 2002 - 2014 by the deal.II authors
+// Copyright (C) 2002 - 2015 by the deal.II authors
 //
 // This file is part of the deal.II library.
 //
@@ -23,28 +23,20 @@ DEAL_II_NAMESPACE_OPEN
 template class SparseMIC<double>;
 template void SparseMIC<double>::initialize<double> (const SparseMatrix<double> &,
                                                      const AdditionalData &data);
-template void SparseMIC<double>::decompose<double> (const SparseMatrix<double> &,
-                                                    const double);
 template void SparseMIC<double>::vmult<double> (Vector<double> &,
                                                 const Vector<double> &) const;
 template void SparseMIC<double>::initialize<float> (const SparseMatrix<float> &,
                                                     const AdditionalData &data);
-template void SparseMIC<double>::decompose<float> (const SparseMatrix<float> &,
-                                                   const double);
 template void SparseMIC<double>::vmult<float> (Vector<float> &,
                                                const Vector<float> &) const;
 
 template class SparseMIC<float>;
 template void SparseMIC<float>::initialize<double> (const SparseMatrix<double> &,
                                                     const AdditionalData &data);
-template void SparseMIC<float>::decompose<double> (const SparseMatrix<double> &,
-                                                   const double);
 template void SparseMIC<float>::vmult<double> (Vector<double> &,
                                                const Vector<double> &) const;
 template void SparseMIC<float>::initialize<float> (const SparseMatrix<float> &,
                                                    const AdditionalData &data);
-template void SparseMIC<float>::decompose<float> (const SparseMatrix<float> &,
-                                                  const double);
 template void SparseMIC<float>::vmult<float> (Vector<float> &,
                                               const Vector<float> &) const;
 
-- 
2.39.5