]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
Added instantiations for boundary extraction in codim one
authorheltai <heltai@0785d39b-7218-0410-832d-ea1e28bc413d>
Fri, 10 Jun 2011 10:13:26 +0000 (10:13 +0000)
committerheltai <heltai@0785d39b-7218-0410-832d-ea1e28bc413d>
Fri, 10 Jun 2011 10:13:26 +0000 (10:13 +0000)
git-svn-id: https://svn.dealii.org/trunk@23805 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/source/dofs/dof_tools.cc
deal.II/source/dofs/dof_tools.inst.in

index 3fe723ba722a8ee0bc4404121c4add616e8c2c3d..787d89041ead7e45cb030bd125de8995d6e16bc3 100644 (file)
@@ -750,7 +750,7 @@ namespace DoFTools
                                  const Table<2,Coupling> &int_mask,
                                  const Table<2,Coupling> &flux_mask)
       {
-       const FiniteElement<DH::dimension> &fe = dof.get_fe();
+       const FiniteElement<DH::dimension,DH::space_dimension> &fe = dof.get_fe();
 
        std::vector<unsigned int> dofs_on_this_cell(fe.dofs_per_cell);
        std::vector<unsigned int> dofs_on_other_cell(fe.dofs_per_cell);
@@ -1210,13 +1210,13 @@ namespace DoFTools
                                     // beginning of this function.
     std::vector<bool> user_flags;
     dof.get_tria().save_user_flags(user_flags);
-    const_cast<Triangulation<DH::dimension> &>(dof.get_tria()).clear_user_flags ();
+    const_cast<Triangulation<DH::dimension,DH::space_dimension> &>(dof.get_tria()).clear_user_flags ();
 
     internal::make_flux_sparsity_pattern (dof, sparsity,
                                          int_mask, flux_mask);
   
                                     // finally restore the user flags
-    const_cast<Triangulation<DH::dimension> &>(dof.get_tria()).load_user_flags(user_flags);
+    const_cast<Triangulation<DH::dimension,DH::space_dimension> &>(dof.get_tria()).load_user_flags(user_flags);
   }
 
 
@@ -3406,7 +3406,7 @@ namespace DoFTools
     Vector<double>       &dof_data,
     const unsigned int    component)
   {
-    const Triangulation<DH::dimension> &tria = dof_handler.get_tria();
+    const Triangulation<DH::dimension,DH::space_dimension> &tria = dof_handler.get_tria();
 
     Assert (cell_data.size()==tria.n_active_cells(),
            ExcWrongSize (cell_data.size(), tria.n_active_cells()));
@@ -3745,6 +3745,7 @@ namespace DoFTools
     Assert (boundary_indicators.find (255) == boundary_indicators.end(),
            ExcInvalidBoundaryIndicator());
     const unsigned int dim=DH::dimension;
+    const unsigned int spacedim=DH::space_dimension;
 
                                     // let's see whether we have to
                                     // check for certain boundary
@@ -3788,7 +3789,7 @@ namespace DoFTools
              (boundary_indicators.find (cell->face(face)->boundary_indicator())
               != boundary_indicators.end()))
            {
-             const FiniteElement<DH::dimension> &fe = cell->get_fe();
+             const FiniteElement<DH::dimension, DH::space_dimension> &fe = cell->get_fe();
 
              const unsigned int dofs_per_face = fe.dofs_per_face;
              face_dof_indices.resize (dofs_per_face);
@@ -3913,7 +3914,7 @@ namespace DoFTools
              (boundary_indicators.find (cell->face(face)->boundary_indicator())
               != boundary_indicators.end()))
            {
-             const FiniteElement<DH::dimension> &fe = cell->get_fe();
+             const FiniteElement<DH::dimension, DH::space_dimension> &fe = cell->get_fe();
 
              const unsigned int dofs_per_cell = fe.dofs_per_cell;
              cell_dof_indices.resize (dofs_per_cell);
index adbc04d862a1a512efefc2ef4cc6711d66acfeab..ddabf8f8c28a50d35cfbc2311e80732db0cfbfd6 100644 (file)
@@ -426,6 +426,13 @@ DoFTools::count_dofs_with_subdomain_association<hp::DoFHandler<deal_II_dimension
 
 #if deal_II_dimension < 3
 template
+void
+DoFTools::extract_boundary_dofs<DoFHandler<deal_II_dimension,deal_II_dimension+1> >
+  (const DoFHandler<deal_II_dimension,deal_II_dimension+1> &,
+   const std::vector<bool>                  &,
+   std::vector<bool>                        &,
+   const std::set<unsigned char> &);
+template
 unsigned int
 DoFTools::count_dofs_with_subdomain_association<DoFHandler<deal_II_dimension,deal_II_dimension+1> >
 (const DoFHandler<deal_II_dimension,deal_II_dimension+1> &,

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.