From d30b1bf412e9ada5ddae81ecb8ecd773898e77b6 Mon Sep 17 00:00:00 2001 From: "daniel.arndt" Date: Fri, 6 Sep 2013 09:05:25 +0000 Subject: [PATCH] identify_periodic_face_pairs creates periodicity information for two given boundary_ids on distributed meshes git-svn-id: https://svn.dealii.org/trunk@30625 0785d39b-7218-0410-832d-ea1e28bc413d --- deal.II/include/deal.II/grid/grid_tools.h | 26 ++++++++ deal.II/source/grid/grid_tools.cc | 81 +++++++++++++++++++++++ deal.II/source/grid/grid_tools.inst.in | 32 +++++++++ 3 files changed, 139 insertions(+) diff --git a/deal.II/include/deal.II/grid/grid_tools.h b/deal.II/include/deal.II/grid/grid_tools.h index eabd288760..83e40cf654 100644 --- a/deal.II/include/deal.II/grid/grid_tools.h +++ b/deal.II/include/deal.II/grid/grid_tools.h @@ -24,6 +24,7 @@ #include #include #include +#include #include #include @@ -1153,6 +1154,31 @@ namespace GridTools const int direction, const dealii::Tensor<1,DH::space_dimension> &offset); + /** + * Add periodicity information to the @p periodicity_vector for the + * boundaries with boundary id @p b_id1 and @p b_id2 in cartesian + * direction @p direction. + * + * This function tries to match all faces belonging to the first + * boundary with faces belonging to the second boundary + * by comparing the center of the cell faces. To find the correct + * corresponding faces, the direction argument indicates in which + * cartesian direction periodicity should be set. + * ((0,1,2) -> (x,y,z)-direction) + * + * The output of this function can be used in + * parallel::distributed::Triangulation::add_periodicity + */ + template + void + identify_periodic_face_pairs + (const DH &dof_handler, + const types::boundary_id b_id1, + const types::boundary_id b_id2, + const unsigned int direction, + std::vector > + &periodicity_vector); /** diff --git a/deal.II/source/grid/grid_tools.cc b/deal.II/source/grid/grid_tools.cc index a49ede55e2..4dfd9ec119 100644 --- a/deal.II/source/grid/grid_tools.cc +++ b/deal.II/source/grid/grid_tools.cc @@ -17,6 +17,7 @@ #include #include #include +#include #include #include #include @@ -2540,6 +2541,86 @@ next_cell: } + template + void + identify_periodic_face_pairs + (const DH &dof_handler, + const types::boundary_id b_id1, + const types::boundary_id b_id2, + const unsigned int direction, + std::vector > + &periodicity_vector) + { + static const unsigned int dim = DH::dimension; + static const unsigned int spacedim = DH::space_dimension; + + Assert (0<=direction && direction, + Point > face_locations; + + // Collect faces with boundary_indicator b_id1 + typename DH::cell_iterator cell = dof_handler.begin(); + typename DH::cell_iterator endc = dof_handler.end(); + for (; cell != endc; ++cell) + for (unsigned int f = 0; f < GeometryInfo::faces_per_cell; ++f) + if(cell->face(f)->boundary_indicator() == b_id1) + { + Point face_center (cell->face(f)->center()); + face_center(direction)=0.; + const std::pair + cell_face_pair = std::make_pair(cell, f); + face_locations[cell_face_pair]=face_center; + } + + // Match faces with boundary_indicator b_id2 to the ones in + //face_locations + cell = dof_handler.begin(); + for (; cell != endc; ++cell) + for(unsigned int f=0;f::faces_per_cell;++f) + if(cell->face(f)->boundary_indicator() == b_id2) + { + typename std::map, + Point >::const_iterator p + = face_locations.begin(); + + Point center2 (cell->face(f)->center()); + center2(direction)=0.; + + for (; p != face_locations.end(); ++p) + { + if (center2.distance(p->second) < 1e-4*cell->face(f)->diameter()) + { + const std_cxx1x::tuple + periodic_tuple + = std::make_tuple(p->first.first, + p->first.second, + cell, + f); + + periodicity_vector.push_back(periodic_tuple); + + face_locations.erase(p); + break; + } + Assert (p != face_locations.end(), + ExcMessage ("No corresponding face was found!")); + } + } + + Assert (face_locations.size() == 0, + ExcMessage ("There are unmatched faces!")); + } + } /* namespace GridTools */ diff --git a/deal.II/source/grid/grid_tools.inst.in b/deal.II/source/grid/grid_tools.inst.in index 4411897224..9b3bc3b7f6 100644 --- a/deal.II/source/grid/grid_tools.inst.in +++ b/deal.II/source/grid/grid_tools.inst.in @@ -271,9 +271,41 @@ for (X : TRIANGULATION_AND_DOFHANDLERS; deal_II_dimension : DIMENSIONS ; deal_II int, const Tensor<1,deal_II_space_dimension> &); + template + void + identify_periodic_face_pairs + (const X &, + const types::boundary_id, + const types::boundary_id, + const unsigned int, + std::vector > &); + #endif \} #endif } +for (deal_II_dimension : DIMENSIONS ; deal_II_space_dimension : SPACE_DIMENSIONS) +{ +#if deal_II_dimension <= deal_II_space_dimension + #if deal_II_dimension >= 2 + namespace GridTools \{ + template + void + identify_periodic_face_pairs + (const parallel::distributed::Triangulation &, + const types::boundary_id, + const types::boundary_id, + const unsigned int, + std::vector::cell_iterator, + unsigned int, + typename parallel::distributed::Triangulation::cell_iterator, + unsigned int> > &); + \} + #endif +#endif +} + + -- 2.39.5