Correct instantiation of MappingQ1Euler.

git-svn-id: https://svn.dealii.org/branches/branch_petscscalar_complex@32701 0785d39b-7218-0410-832d-ea1e28bc413d

diff --git a/deal.II/source/fe/mapping_q1_eulerian.cc b/deal.II/source/fe/mapping_q1_eulerian.cc
index 49fd4d75177f8282fcc1bb33d62e9f23b9649541..8a9d8b3a431c57ae828a6b8d833802cd3906ea27 100644 (file)
--- a/deal.II/source/fe/mapping_q1_eulerian.cc
+++ b/deal.II/source/fe/mapping_q1_eulerian.cc
@@ -78,10 +78,7 @@ compute_mapping_support_points(const typename Triangulation<dim,spacedim>::cell_
   // now get the values of the shift
   // vectors at the vertices
   Vector<double> mapping_values (shiftmap_dof_handler->get_fe().dofs_per_cell);
-
-  // @whattodo Some of the problems in petsc_vector_base.h are instantiated here.
-  // dof_cell->get_dof_values (*euler_transform_vectors, mapping_values); <- complaint line 
-  Assert ((false), ExcMessage ("This function is corrupt: @whattodo"));
+  dof_cell->get_dof_values (*euler_transform_vectors, mapping_values); 
 
   for (unsigned int i=0; i<GeometryInfo<dim>::vertices_per_cell; ++i)
     {

diff --git a/deal.II/source/fe/mapping_q1_eulerian.inst.in b/deal.II/source/fe/mapping_q1_eulerian.inst.in
index 06d1ebc2906ab84c45aebd268682ce6ac49c54a0..c1b2e4a467f3b78f154452a12c47bae451e285a0 100644 (file)
--- a/deal.II/source/fe/mapping_q1_eulerian.inst.in
+++ b/deal.II/source/fe/mapping_q1_eulerian.inst.in
@@ -15,14 +15,16 @@
 // ---------------------------------------------------------------------
 
 
-
-for (deal_II_dimension : DIMENSIONS; deal_II_space_dimension :  SPACE_DIMENSIONS)
+// TODO[TY]: allow this to be instantiated by PETScWrappers::Vector and PETSC_WITH_COMPLEX
+for (deal_II_dimension : DIMENSIONS; deal_II_space_dimension : SPACE_DIMENSIONS)
   {
 #if deal_II_dimension <= deal_II_space_dimension       
     template class MappingQ1Eulerian<deal_II_dimension, Vector<double>,        deal_II_space_dimension>;
     template class MappingQ1Eulerian<deal_II_dimension, Vector<float>, deal_II_space_dimension>;
 #ifdef DEAL_II_WITH_PETSC
+#ifdef PETSC_WITH_COMPLEX
     template class MappingQ1Eulerian<deal_II_dimension, PETScWrappers::Vector, deal_II_space_dimension>;
 #endif
+#endif
 #endif
   }