add has_begin detector

diff --git a/include/deal.II/matrix_free/fe_evaluation.h b/include/deal.II/matrix_free/fe_evaluation.h
index 8da043d70beeb7068552ff62a86c91844b460074..4b088f977d17bb6b022898310b8c9979d14ccb4c 100644 (file)
--- a/include/deal.II/matrix_free/fe_evaluation.h
+++ b/include/deal.II/matrix_free/fe_evaluation.h
@@ -3450,6 +3450,27 @@ namespace internal
   };
 
 
+
+  // same as above to check
+  // T::const_iterator T::begin() const
+  template <typename T>
+  struct has_begin
+  {
+  private:
+    static void
+    detect(...);
+
+    template <typename U>
+    static decltype(std::declval<U const>().begin())
+    detect(const U &);
+
+  public:
+    static constexpr bool value =
+      std::is_same<typename T::const_iterator,
+                   decltype(detect(std::declval<T>()))>::value;
+  };
+
+
   // access to generic vectors that have operator ().
   // FIXME: this is wrong for Trilinos/Petsc MPI vectors
   template <typename VectorType,

diff --git a/tests/matrix_free/fe_evaluation_type_traits.cc b/tests/matrix_free/fe_evaluation_type_traits.cc
index 12ab08476b5b76d6f496ed07ac3efebe27f3f15e..d3363649aa1864815118cb3c1f19143e658355a3 100644 (file)
--- a/tests/matrix_free/fe_evaluation_type_traits.cc
+++ b/tests/matrix_free/fe_evaluation_type_traits.cc
@@ -129,5 +129,14 @@ main()
                TrilinosWrappers::MPI::Vector>::value
           << std::endl;
 
+  // check has_begin():
+  deallog
+    << "has_begin:" << std::endl
+    << "LinearAlgebra::distributed::Vector = "
+    << internal::has_begin<LinearAlgebra::distributed::Vector<double>>::value
+    << std::endl
+    << "TrilinosWrappers::MPI::Vector = "
+    << internal::has_begin<TrilinosWrappers::MPI::Vector>::value << std::endl;
+
   deallog << "OK" << std::endl;
 }

diff --git a/tests/matrix_free/fe_evaluation_type_traits.with_trilinos=true.output b/tests/matrix_free/fe_evaluation_type_traits.with_trilinos=true.output
index 4e82d5220255e34b378ffc3ec424e0ded7f137c2..71bac87207dc316dc26a23a270da2cf1732cb6cb 100644 (file)
--- a/tests/matrix_free/fe_evaluation_type_traits.with_trilinos=true.output
+++ b/tests/matrix_free/fe_evaluation_type_traits.with_trilinos=true.output
@@ -10,4 +10,7 @@ DEAL::Dummy2::operator()
 DEAL::has_partitioners_are_compatible:
 DEAL::LinearAlgebra::distributed::Vector = 1
 DEAL::TrilinosWrappers::MPI::Vector = 0
+DEAL::has_begin:
+DEAL::LinearAlgebra::distributed::Vector = 1
+DEAL::TrilinosWrappers::MPI::Vector = 1
 DEAL::OK