From 947cf9d300570cdfbe785c7627d1ebd5cdb9787a Mon Sep 17 00:00:00 2001
From: young <young@0785d39b-7218-0410-832d-ea1e28bc413d>
Date: Thu, 27 Mar 2014 09:19:37 +0000
Subject: [PATCH] Clean up function interpolate and instantiate correctly.
 Leave a TODO.

git-svn-id: https://svn.dealii.org/branches/branch_petscscalar_complex@32699 0785d39b-7218-0410-832d-ea1e28bc413d
---
 .../deal.II/numerics/vector_tools.templates.h       | 13 +++++--------
 .../numerics/vector_tools_interpolate.inst.in       | 12 +++++++++++-
 2 files changed, 16 insertions(+), 9 deletions(-)

diff --git a/deal.II/include/deal.II/numerics/vector_tools.templates.h b/deal.II/include/deal.II/numerics/vector_tools.templates.h
index e41bd1be6e..7216a3c813 100644
--- a/deal.II/include/deal.II/numerics/vector_tools.templates.h
+++ b/deal.II/include/deal.II/numerics/vector_tools.templates.h
@@ -280,20 +280,17 @@ namespace VectorTools
 
 
 
-  // @whattodo TODO: To return functionality for complex numbers, this
-  // needs to be templated with Number and instantiated for
-  // REAL_SERIAL_VECTORS+double and
-  // COMPLEX_SERIAL_VECTORS+complex<double>
-  template <int dim, class InVector, class OutVector, int spacedim>
+  // @TODO[TY] also enable instantiation for COMPLEX_SERIAL_VECTORS and std::complex<double>
+  template <int dim, class InVector, class OutVector, int spacedim, typename number>
   void
   interpolate (const DoFHandler<dim,spacedim>  &dof_1,
                const DoFHandler<dim,spacedim>  &dof_2,
-               const FullMatrix<double>        &transfer, // <- FullMatrix<std::complex>
+               const FullMatrix<number>        &transfer,
                const InVector                  &data_1,
                OutVector                       &data_2)
   {
-    Vector<double> cell_data_1(dof_1.get_fe().dofs_per_cell); // <- Vector<std::complex>
-    Vector<double> cell_data_2(dof_2.get_fe().dofs_per_cell); // <- Vector<std::complex>
+    Vector<number> cell_data_1(dof_1.get_fe().dofs_per_cell);
+    Vector<number> cell_data_2(dof_2.get_fe().dofs_per_cell);
 
     std::vector<short unsigned int> touch_count (dof_2.n_dofs(), 0); //TODO: check on datatype... kinda strange (UK)
     std::vector<types::global_dof_index>       local_dof_indices (dof_2.get_fe().dofs_per_cell);
diff --git a/deal.II/source/numerics/vector_tools_interpolate.inst.in b/deal.II/source/numerics/vector_tools_interpolate.inst.in
index 9b46645411..00bb3370ee 100644
--- a/deal.II/source/numerics/vector_tools_interpolate.inst.in
+++ b/deal.II/source/numerics/vector_tools_interpolate.inst.in
@@ -33,11 +33,21 @@ for (VEC : REAL_SERIAL_VECTORS ; deal_II_dimension : DIMENSIONS; deal_II_space_d
          const Function<deal_II_space_dimension>&,
          VEC&);
 
+      \}
+#endif
+  }
+
+
+for (VEC : REAL_SERIAL_VECTORS ; deal_II_dimension : DIMENSIONS; deal_II_space_dimension :  SPACE_DIMENSIONS; NUMBER : double)
+  {
+#if deal_II_dimension <= deal_II_space_dimension
+    namespace VectorTools \{
+
       template
         void interpolate
         (const DoFHandler<deal_II_dimension,deal_II_space_dimension>&,
          const DoFHandler<deal_II_dimension,deal_II_space_dimension>&,
-         const FullMatrix<double>&,
+         const FullMatrix<NUMBER>&,
          const VEC&,
          VEC&);
 
-- 
2.39.5