From: Martin Kronbichler <kronbichler@lnm.mw.tum.de>
Date: Thu, 30 Apr 2020 05:39:29 +0000 (+0200)
Subject: Fix SFINAE detection of LinOp for matrix-free operators
X-Git-Tag: v9.2.0-rc1~146^2
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=3f2e10d8cbdb27dc25c0f5f90ebc2682d1a271ef;p=dealii.git

Fix SFINAE detection of LinOp for matrix-free operators
---

diff --git a/include/deal.II/lac/la_parallel_vector.h b/include/deal.II/lac/la_parallel_vector.h
index 7f4374099a..2841949469 100644
--- a/include/deal.II/lac/la_parallel_vector.h
+++ b/include/deal.II/lac/la_parallel_vector.h
@@ -1824,7 +1824,7 @@ namespace internal
 
         template <typename U>
         static decltype(std::declval<U>().initialize_dof_vector(
-          LinearAlgebra::distributed::Vector<double>()))
+          std::declval<LinearAlgebra::distributed::Vector<Number> &>()))
         detect(const U &);
 
       public:
@@ -1840,7 +1840,7 @@ namespace internal
                   MatrixType>::type * = nullptr>
       static void
       reinit_domain_vector(MatrixType &                                mat,
-                           LinearAlgebra::distributed::Vector<double> &vec,
+                           LinearAlgebra::distributed::Vector<Number> &vec,
                            bool /*omit_zeroing_entries*/)
       {
         vec.reinit(mat.locally_owned_domain_indices(),
@@ -1854,7 +1854,7 @@ namespace internal
                                 MatrixType>::type * = nullptr>
       static void
       reinit_domain_vector(MatrixType &                                mat,
-                           LinearAlgebra::distributed::Vector<double> &vec,
+                           LinearAlgebra::distributed::Vector<Number> &vec,
                            bool omit_zeroing_entries)
       {
         mat.initialize_dof_vector(vec);
@@ -1870,7 +1870,7 @@ namespace internal
                   MatrixType>::type * = nullptr>
       static void
       reinit_range_vector(MatrixType &                                mat,
-                          LinearAlgebra::distributed::Vector<double> &vec,
+                          LinearAlgebra::distributed::Vector<Number> &vec,
                           bool /*omit_zeroing_entries*/)
       {
         vec.reinit(mat.locally_owned_range_indices(),
@@ -1884,7 +1884,7 @@ namespace internal
                                 MatrixType>::type * = nullptr>
       static void
       reinit_range_vector(MatrixType &                                mat,
-                          LinearAlgebra::distributed::Vector<double> &vec,
+                          LinearAlgebra::distributed::Vector<Number> &vec,
                           bool omit_zeroing_entries)
       {
         mat.initialize_dof_vector(vec);