From f19559daea231359ec54d274a5e6f99e40982665 Mon Sep 17 00:00:00 2001
From: Denis Davydov <davydden@gmail.com>
Date: Thu, 25 Oct 2018 09:10:30 +0200
Subject: [PATCH] add template parameter to
 DynamicSparsityPattern::compute_Tmmult_pattern()

---
 include/deal.II/lac/dynamic_sparsity_pattern.h |  6 +++---
 source/lac/dynamic_sparsity_pattern.cc         | 10 ++++++++--
 2 files changed, 11 insertions(+), 5 deletions(-)

diff --git a/include/deal.II/lac/dynamic_sparsity_pattern.h b/include/deal.II/lac/dynamic_sparsity_pattern.h
index c33fe8edcb..40ad75005e 100644
--- a/include/deal.II/lac/dynamic_sparsity_pattern.h
+++ b/include/deal.II/lac/dynamic_sparsity_pattern.h
@@ -32,7 +32,6 @@
 DEAL_II_NAMESPACE_OPEN
 
 class DynamicSparsityPattern;
-class SparsityPattern;
 
 /*! @addtogroup Sparsity
  *@{
@@ -450,9 +449,10 @@ public:
    * Construct and store in this object the sparsity pattern corresponding to
    * the product of transposed @p left and and non-transpose @p right sparsity pattern.
    */
+  template <typename SparsityPatternTypeLeft, typename SparsityPatternTypeRight>
   void
-  compute_Tmmult_pattern(const SparsityPattern &left,
-                         const SparsityPattern &right);
+  compute_Tmmult_pattern(const SparsityPatternTypeLeft & left,
+                         const SparsityPatternTypeRight &right);
 
   /**
    * Print the sparsity pattern. The output consists of one line per row of
diff --git a/source/lac/dynamic_sparsity_pattern.cc b/source/lac/dynamic_sparsity_pattern.cc
index 62ca57ff73..712e4706eb 100644
--- a/source/lac/dynamic_sparsity_pattern.cc
+++ b/source/lac/dynamic_sparsity_pattern.cc
@@ -405,9 +405,11 @@ DynamicSparsityPattern::symmetrize()
 
 
 
+template <typename SparsityPatternTypeLeft, typename SparsityPatternTypeRight>
 void
-DynamicSparsityPattern::compute_Tmmult_pattern(const SparsityPattern &sp_A,
-                                               const SparsityPattern &sp_B)
+DynamicSparsityPattern::compute_Tmmult_pattern(
+  const SparsityPatternTypeLeft & sp_A,
+  const SparsityPatternTypeRight &sp_B)
 {
   Assert(sp_A.n_rows() == sp_B.n_rows(),
          ExcDimensionMismatch(sp_A.n_rows(), sp_B.n_rows()));
@@ -604,4 +606,8 @@ template void
 DynamicSparsityPattern::compute_mmult_pattern(const SparsityPattern &,
                                               const SparsityPattern &);
 
+template void
+DynamicSparsityPattern::compute_Tmmult_pattern(const SparsityPattern &,
+                                               const SparsityPattern &);
+
 DEAL_II_NAMESPACE_CLOSE
-- 
2.39.5