LAC: also allow ConstraintMatrix::distribute for std::complex

Fix a typetrait check in ConstraintMatrix::distribute such that said class
and method can be used with complex valued vector types.

diff --git a/include/deal.II/lac/constraint_matrix.templates.h b/include/deal.II/lac/constraint_matrix.templates.h
index 8c20a74d700c7dd651a17fb10d2fe1dc8af0440e..06cb05f08e4bea2ac067059d7453bdcfb2370f88 100644 (file)
--- a/include/deal.II/lac/constraint_matrix.templates.h
+++ b/include/deal.II/lac/constraint_matrix.templates.h
@@ -852,12 +852,17 @@ ConstraintMatrix::distribute (VectorType &vec) const
   // call compress() finally. the first case here is for the complicated case,
   // the last else is for the simple case (sequential vector)
   const IndexSet vec_owned_elements = vec.locally_owned_elements();
+
+  // FIXME: This has to be refactored into a typetrait
   if ((typeid(vec) != typeid(Vector<double>)) &&
       (typeid(vec) != typeid(Vector<float>)) &&
+      (typeid(vec) != typeid(Vector<std::complex<double> >)) &&
       (typeid(vec) != typeid(BlockVector<double>)) &&
       (typeid(vec) != typeid(BlockVector<float>)) &&
+      (typeid(vec) != typeid(BlockVector<std::complex<double> >)) &&
       (typeid(vec) != typeid(LinearAlgebra::Vector<double>)) &&
-      (typeid(vec) != typeid(LinearAlgebra::Vector<float>)))
+      (typeid(vec) != typeid(LinearAlgebra::Vector<float>)) &&
+      (typeid(vec) != typeid(LinearAlgebra::Vector<std::complex<double> >)))
     {
       // This processor owns only part of the vector. one may think that
       // every processor should be able to simply communicate those elements