add template parameter to DynamicSparsityPattern::compute_Tmmult_pattern()

diff --git a/include/deal.II/lac/dynamic_sparsity_pattern.h b/include/deal.II/lac/dynamic_sparsity_pattern.h
index c33fe8edcbe34b33e76a1566cd727a8d86ef3543..40ad75005e9cd7305aaabb2435d1ccfd1aed7e65 100644 (file)
--- 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 62ca57ff7375616c47dc97e0d2af68c70848ba43..712e4706eb1e575c84bc15b1756011573f0cc33d 100644 (file)
--- 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