SparseMatrix::reinit(): accept other sparse matrix

diff --git a/include/deal.II/lac/sparse_matrix.h b/include/deal.II/lac/sparse_matrix.h
index d707f1f38ef5e7af95d969661914794a0216891c..2cd998e33cffb7ed4230bb44456e27c5e0d74200 100644 (file)
--- a/include/deal.II/lac/sparse_matrix.h
+++ b/include/deal.II/lac/sparse_matrix.h
@@ -696,6 +696,16 @@ public:
   virtual void
   reinit(const SparsityPattern &sparsity);
 
+  /**
+   * Reinitialize the sparse matrix with the sparsity pattern of the given
+   * @p sparse_matrix. See also comments of the function above.
+   *
+   * @note The elements of the matrix are set to zero by this function.
+   */
+  template <typename number2>
+  void
+  reinit(const SparseMatrix<number2> &sparse_matrix);
+
   /**
    * Release all memory and return to a state just like after having called
    * the default constructor. It also forgets the sparsity pattern it was
@@ -1787,6 +1797,16 @@ private:
 
 
 
+template <typename number>
+template <typename number2>
+void
+SparseMatrix<number>::reinit(const SparseMatrix<number2> &sparse_matrix)
+{
+  this->reinit(sparse_matrix.get_sparsity_pattern());
+}
+
+
+
 template <typename number>
 inline typename SparseMatrix<number>::size_type
 SparseMatrix<number>::m() const