Patch by Armin Ghajar Jazi: Also instantiate MappingQEulerian for Trilinos vectors.

git-svn-id: https://svn.dealii.org/trunk@30311 0785d39b-7218-0410-832d-ea1e28bc413d

diff --git a/deal.II/doc/news/changes.h b/deal.II/doc/news/changes.h
index f354c1a001d651f7a12df6a59496ee42bbb79313..75be54c8ea08acdcbcf158f5dbcffe26e9d8070a 100644 (file)
--- a/deal.II/doc/news/changes.h
+++ b/deal.II/doc/news/changes.h
@@ -56,6 +56,14 @@ inconvenience this causes.
 <h3>Specific improvements</h3>
 
 <ol>
+  <li>
+  New: MappingQEulerian is now also instantiated for vector elements
+  of type TrilinosWrappers::Vector as well as the MPI and block
+  variants.
+  <br>
+  (Armin Ghajar Jazi, 2013/08/14)
+  </li>
+
   <li>
   Fixed: The FiniteElement::face_to_cell_index() function had a bug
   that made it work incorrectly for elements that have more than one

diff --git a/deal.II/source/fe/mapping_q_eulerian.cc b/deal.II/source/fe/mapping_q_eulerian.cc
index 4e005c2f8f2deca4614ea6ae997db9981e03b75c..c8ba7374746cc74fef9a4a8b1b01953d25fd7886 100644 (file)
--- a/deal.II/source/fe/mapping_q_eulerian.cc
+++ b/deal.II/source/fe/mapping_q_eulerian.cc
@@ -18,6 +18,9 @@
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/lac/vector.h>
 #include <deal.II/lac/petsc_vector.h>
+#include <deal.II/lac/trilinos_vector.h>
+#include <deal.II/lac/trilinos_block_vector.h>
+#include <deal.II/lac/trilinos_parallel_block_vector.h>
 #include <deal.II/grid/tria_iterator.h>
 #include <deal.II/dofs/dof_handler.h>
 #include <deal.II/dofs/dof_accessor.h>

diff --git a/deal.II/source/fe/mapping_q_eulerian.inst.in b/deal.II/source/fe/mapping_q_eulerian.inst.in
index dd64c5497fdda71d773f7fa6d49d91992d2e8828..439ef92f5174050fc0b2e778f10408935be1e879 100644 (file)
--- a/deal.II/source/fe/mapping_q_eulerian.inst.in
+++ b/deal.II/source/fe/mapping_q_eulerian.inst.in
@@ -25,6 +25,17 @@ for (deal_II_dimension : DIMENSIONS; deal_II_space_dimension :  SPACE_DIMENSIONS
     template class MappingQEulerian<deal_II_dimension,
                                    PETScWrappers::MPI::Vector, deal_II_space_dimension>;
 #  endif   
+#  ifdef DEAL_II_WITH_TRILINOS
+    template class MappingQEulerian<deal_II_dimension,
+                                    TrilinosWrappers::Vector, deal_II_space_dimension>;
+    template class MappingQEulerian<deal_II_dimension,
+                                    TrilinosWrappers::MPI::Vector, deal_II_space_dimension>;
+    template class MappingQEulerian<deal_II_dimension,
+                                    TrilinosWrappers::BlockVector, deal_II_space_dimension>;
+    template class MappingQEulerian<deal_II_dimension,
+                                    TrilinosWrappers::MPI::BlockVector, deal_II_space_dimension>;
+#  endif
+
 #endif
   }