33 #include <boost/serialization/map.hpp> 34 #include <boost/serialization/split_member.hpp> 35 #include <boost/serialization/unique_ptr.hpp> 36 #include <boost/serialization/vector.hpp> 37 #include <boost/signals2.hpp> 52 "The name 'signals' is already defined. You are most likely using the QT library \ 53 and using the 'signals' keyword. You can either #include the Qt headers (or any conflicting headers) \ 54 *after* the deal.II headers or you can define the 'QT_NO_KEYWORDS' macro and use the 'Q_SIGNALS' macro." 59 template <
int dim,
int spacedim>
75 template <
typename CellIterator>
76 struct PeriodicFacePair;
79 template <
int,
int,
int>
81 template <
int spacedim>
83 template <
int,
int,
int>
88 namespace TriangulationImplementation
102 struct Implementation;
103 struct ImplementationMixedMesh;
106 namespace TriaAccessorImplementation
108 struct Implementation;
114 template <
int dim,
int spacedim>
129 namespace TriangulationImplementation
202 template <
class Archive>
204 serialize(Archive &ar,
const unsigned int version);
260 template <
class Archive>
262 serialize(Archive &ar,
const unsigned int version);
319 template <
class Archive>
321 serialize(Archive &ar,
const unsigned int version);
1126 template <
int dim,
int spacedim = dim>
1192 limit_level_difference_at_vertices = 0x1,
1213 eliminate_unrefined_islands = 0x2,
1229 patch_level_1 = 0x4,
1253 coarsest_level_1 = 0x8,
1278 allow_anisotropic_smoothing = 0x10,
1311 eliminate_refined_inner_islands = 0x100,
1316 eliminate_refined_boundary_islands = 0x200,
1322 do_not_produce_unrefined_islands = 0x400,
1328 smoothing_on_refinement =
1329 (limit_level_difference_at_vertices | eliminate_unrefined_islands),
1334 smoothing_on_coarsening =
1335 (eliminate_refined_inner_islands | eliminate_refined_boundary_islands |
1336 do_not_produce_unrefined_islands),
1343 maximum_smoothing = 0xffff ^ allow_anisotropic_smoothing
1541 std::list<typename Triangulation<dim, spacedim>::cell_iterator>
1548 static const unsigned int dimension = dim;
1553 static const unsigned int space_dimension = spacedim;
1570 const bool check_for_distorted_cells =
false);
1629 get_mesh_smoothing()
const;
1684 reset_all_manifolds();
1747 virtual std::vector<types::boundary_id>
1748 get_boundary_ids()
const;
1762 virtual std::vector<types::manifold_id>
1763 get_manifold_ids()
const;
1861 &construction_data);
1872 create_triangulation_compatibility(
1884 flip_all_direction_flags();
1898 set_all_refine_flags();
1923 refine_global(
const unsigned int times = 1);
1938 coarsen_global(
const unsigned int times = 1);
1972 execute_coarsening_and_refinement();
2004 prepare_coarsening_and_refinement();
2049 template <
typename T>
2054 template <
typename InputIterator>
2058 return std::accumulate(first, last,
T());
2122 boost::signals2::signal<void(
2132 boost::signals2::signal<void(
2142 boost::signals2::signal<void(
2294 save_refine_flags(std::ostream &out)
const;
2300 save_refine_flags(std::vector<bool> &v)
const;
2306 load_refine_flags(std::istream &in);
2312 load_refine_flags(
const std::vector<bool> &v);
2318 save_coarsen_flags(std::ostream &out)
const;
2324 save_coarsen_flags(std::vector<bool> &v)
const;
2330 load_coarsen_flags(std::istream &out);
2336 load_coarsen_flags(
const std::vector<bool> &v);
2343 get_anisotropic_refinement_flag()
const;
2367 save_user_flags(std::ostream &out)
const;
2375 save_user_flags(std::vector<bool> &v)
const;
2382 load_user_flags(std::istream &in);
2389 load_user_flags(
const std::vector<bool> &v);
2396 clear_user_flags_line();
2403 save_user_flags_line(std::ostream &out)
const;
2411 save_user_flags_line(std::vector<bool> &v)
const;
2418 load_user_flags_line(std::istream &in);
2425 load_user_flags_line(
const std::vector<bool> &v);
2432 clear_user_flags_quad();
2439 save_user_flags_quad(std::ostream &out)
const;
2447 save_user_flags_quad(std::vector<bool> &v)
const;
2454 load_user_flags_quad(std::istream &in);
2461 load_user_flags_quad(
const std::vector<bool> &v);
2469 clear_user_flags_hex();
2476 save_user_flags_hex(std::ostream &out)
const;
2484 save_user_flags_hex(std::vector<bool> &v)
const;
2491 load_user_flags_hex(std::istream &in);
2498 load_user_flags_hex(
const std::vector<bool> &v);
2514 save_user_indices(std::vector<unsigned int> &v)
const;
2521 load_user_indices(
const std::vector<unsigned int> &v);
2529 save_user_pointers(std::vector<void *> &v)
const;
2536 load_user_pointers(
const std::vector<void *> &v);
2544 save_user_indices_line(std::vector<unsigned int> &v)
const;
2551 load_user_indices_line(
const std::vector<unsigned int> &v);
2559 save_user_indices_quad(std::vector<unsigned int> &v)
const;
2566 load_user_indices_quad(
const std::vector<unsigned int> &v);
2574 save_user_indices_hex(std::vector<unsigned int> &v)
const;
2581 load_user_indices_hex(
const std::vector<unsigned int> &v);
2588 save_user_pointers_line(std::vector<void *> &v)
const;
2595 load_user_pointers_line(
const std::vector<void *> &v);
2603 save_user_pointers_quad(std::vector<void *> &v)
const;
2610 load_user_pointers_quad(
const std::vector<void *> &v);
2618 save_user_pointers_hex(std::vector<void *> &v)
const;
2625 load_user_pointers_hex(
const std::vector<void *> &v);
2687 begin_active(
const unsigned int level = 0)
const;
2715 end(
const unsigned int level)
const;
2737 end_active(
const unsigned int level)
const;
2750 last_active()
const;
2765 create_cell_iterator(
const CellId &cell_id)
const;
2782 cell_iterators()
const;
2820 active_cell_iterators()
const;
2838 cell_iterators_on_level(
const unsigned int level)
const;
2856 active_cell_iterators_on_level(
const unsigned int level)
const;
2879 begin_active_face()
const;
2907 active_face_iterators()
const;
2925 begin_vertex()
const;
2933 begin_active_vertex()
const;
2971 n_lines(
const unsigned int level)
const;
2977 n_active_lines()
const;
2983 n_active_lines(
const unsigned int level)
const;
2995 n_quads(
const unsigned int level)
const;
3001 n_active_quads()
const;
3007 n_active_quads(
const unsigned int level)
const;
3020 n_hexs(
const unsigned int level)
const;
3026 n_active_hexs()
const;
3033 n_active_hexs(
const unsigned int level)
const;
3064 n_global_active_cells()
const;
3088 n_active_faces()
const;
3116 virtual unsigned int 3117 n_global_levels()
const;
3129 has_hanging_nodes()
const;
3149 const std::vector<Point<spacedim>> &
3150 get_vertices()
const;
3157 n_used_vertices()
const;
3163 vertex_used(
const unsigned int index)
const;
3169 const std::vector<bool> &
3170 get_used_vertices()
const;
3184 max_adjacent_cells()
const;
3193 locally_owned_subdomain()
const;
3205 get_triangulation();
3212 get_triangulation()
const;
3234 n_raw_lines()
const;
3246 n_raw_lines(
const unsigned int level)
const;
3258 n_raw_quads()
const;
3270 n_raw_quads(
const unsigned int level)
const;
3282 n_raw_hexs(
const unsigned int level)
const;
3294 n_raw_cells(
const unsigned int level)
const;
3308 n_raw_faces()
const;
3334 template <
class Archive>
3336 save(Archive &ar,
const unsigned int version)
const;
3355 template <
class Archive>
3357 load(Archive &ar,
const unsigned int version);
3383 std::pair<cell_iterator, unsigned int>,
3384 std::pair<std::pair<cell_iterator, unsigned int>, std::bitset<3>>> &
3385 get_periodic_face_map()
const;
3391 const std::vector<ReferenceCell> &
3392 get_reference_cells()
const;
3399 all_reference_cells_are_hyper_cube()
const;
3407 template <
class Archive>
3409 serialize(Archive &archive,
const unsigned int version);
3413 BOOST_SERIALIZATION_SPLIT_MEMBER()
3429 <<
"You are requesting information from refinement level " 3431 <<
" of a triangulation, but this triangulation only has " 3432 << arg2 <<
" refinement levels. The given level " << arg1
3433 <<
" must be *less* than " << arg2 <<
".");
3441 ExcTriangulationNotEmpty,
3444 <<
"You are trying to perform an operation on a triangulation " 3445 <<
"that is only allowed if the triangulation is currently empty. " 3446 <<
"However, it currently stores " << arg1 <<
" vertices and has " 3447 <<
"cells on " << arg2 <<
" levels.");
3466 <<
"You tried to do something on level " << arg1
3467 <<
", but this level is empty.");
3484 <<
"The given boundary_id " << arg1
3485 <<
" is not defined in this Triangulation!");
3493 ExcInconsistentCoarseningFlags,
3494 "A cell is flagged for coarsening, but either not all of its siblings " 3495 "are active or flagged for coarsening as well. Please clean up all " 3496 "coarsen flags on your triangulation via " 3497 "Triangulation::prepare_coarsening_and_refinement() beforehand!");
3530 write_bool_vector(
const unsigned int magic_number1,
3531 const std::vector<bool> &v,
3532 const unsigned int magic_number2,
3533 std::ostream & out);
3540 read_bool_vector(
const unsigned int magic_number1,
3541 std::vector<bool> &v,
3542 const unsigned int magic_number2,
3550 update_periodic_face_map();
3556 update_reference_cells();
3574 std::vector<GridTools::PeriodicFacePair<cell_iterator>>
3581 std::map<std::pair<cell_iterator, unsigned int>,
3582 std::pair<std::pair<cell_iterator, unsigned int>, std::bitset<3>>>
3612 begin_raw(
const unsigned int level = 0)
const;
3619 end_raw(
const unsigned int level)
const;
3637 begin_raw_line(
const unsigned int level = 0)
const;
3657 begin_line(
const unsigned int level = 0)
const;
3677 begin_active_line(
const unsigned int level = 0)
const;
3715 begin_raw_quad(
const unsigned int level = 0)
const;
3735 begin_quad(
const unsigned int level = 0)
const;
3755 begin_active_quad(
const unsigned int level = 0)
const;
3792 begin_raw_hex(
const unsigned int level = 0)
const;
3812 begin_hex(
const unsigned int level = 0)
const;
3832 begin_active_hex(
const unsigned int level = 0)
const;
3860 clear_despite_subscriptions();
3869 reset_active_cell_indices();
3875 reset_global_cell_indices();
3891 execute_refinement();
3900 execute_coarsening();
3907 fix_coarsen_flags();
3922 virtual unsigned int 3923 coarse_cell_id_to_coarse_cell_index(
3940 coarse_cell_index_to_coarse_cell_id(
3941 const unsigned int coarse_cell_index)
const;
3948 std::unique_ptr<::internal::TriangulationImplementation::TriaLevel>>
3956 std::unique_ptr<::internal::TriangulationImplementation::TriaFaces>
3974 std::map<types::manifold_id, std::unique_ptr<const Manifold<dim, spacedim>>>
4014 std::unique_ptr<std::map<unsigned int, types::boundary_id>>
4037 std::unique_ptr<std::map<unsigned int, types::manifold_id>>
4041 template <
int,
int,
int>
4043 template <
int,
int,
int>
4049 friend struct ::internal::TriaAccessorImplementation::Implementation;
4051 friend struct ::internal::TriangulationImplementation::Implementation;
4052 friend struct ::internal::TriangulationImplementation::
4053 ImplementationMixedMesh;
4055 friend class ::internal::TriangulationImplementation::TriaObjects;
4059 #ifndef DEAL_II_MSVC 4060 static_assert(dim <= spacedim,
4061 "The dimension <dim> of a Triangulation must be less than or " 4062 "equal to the space dimension <spacedim> in which it lives.");
4073 namespace TriangulationImplementation
4075 template <
class Archive>
4077 NumberCache<1>::serialize(Archive &ar,
const unsigned int)
4080 ar &n_lines &n_lines_level;
4081 ar &n_active_lines &n_active_lines_level;
4085 template <
class Archive>
4087 NumberCache<2>::serialize(Archive &ar,
const unsigned int version)
4089 this->NumberCache<1>::serialize(ar, version);
4091 ar &n_quads &n_quads_level;
4092 ar &n_active_quads &n_active_quads_level;
4096 template <
class Archive>
4098 NumberCache<3>::serialize(Archive &ar,
const unsigned int version)
4100 this->NumberCache<2>::serialize(ar, version);
4102 ar &n_hexes &n_hexes_level;
4103 ar &n_active_hexes &n_active_hexes_level;
4110 template <
int dim,
int spacedim>
4115 return vertices_used[index];
4120 template <
int dim,
int spacedim>
4124 return number_cache.n_levels;
4127 template <
int dim,
int spacedim>
4131 return number_cache.n_levels;
4135 template <
int dim,
int spacedim>
4144 template <
int dim,
int spacedim>
4145 inline const std::vector<Point<spacedim>> &
4152 template <
int dim,
int spacedim>
4153 template <
class Archive>
4161 unsigned int n_levels = levels.size();
4163 for (
const auto &
level : levels)
4169 bool faces_is_nullptr = (faces.get() ==
nullptr);
4170 ar & faces_is_nullptr;
4171 if (!faces_is_nullptr)
4177 ar &anisotropic_refinement;
4180 ar &check_for_distorted_cells;
4184 ar &vertex_to_boundary_id_map_1d;
4185 ar &vertex_to_manifold_id_map_1d;
4191 template <
int dim,
int spacedim>
4192 template <
class Archive>
4205 levels.resize(size);
4206 for (
auto &level_ : levels)
4208 std::unique_ptr<internal::TriangulationImplementation::TriaLevel>
level;
4210 level_ = std::move(level);
4214 bool faces_is_nullptr =
true;
4215 ar & faces_is_nullptr;
4216 if (!faces_is_nullptr)
4222 ar &anisotropic_refinement;
4229 for (
auto &level : levels)
4230 level->active_cell_indices.resize(level->refine_flags.size());
4231 reset_active_cell_indices();
4235 bool my_check_for_distorted_cells;
4236 ar & my_check_for_distorted_cells;
4238 Assert(my_check_for_distorted_cells == check_for_distorted_cells,
4239 ExcMessage(
"The triangulation loaded into here must have the " 4240 "same setting with regard to reporting distorted " 4241 "cell as the one previously stored."));
4245 ar &vertex_to_boundary_id_map_1d;
4246 ar &vertex_to_manifold_id_map_1d;
4257 template <
int dim,
int spacedim>
4267 template <
int dim,
int spacedim>
4270 const unsigned int coarse_cell_index)
const 4272 return coarse_cell_index;
boost::signals2::signal< void(const Triangulation< dim, spacedim > &destination_tria)> copy
::internal::TriangulationImplementation::NumberCache< dim > number_cache
virtual unsigned int coarse_cell_id_to_coarse_cell_index(const types::coarse_cell_id coarse_cell_id) const
boost::signals2::signal< void()> post_refinement
unsigned int n_vertices() const
std::vector< unsigned int > n_quads_level
boost::signals2::signal< void(const typename Triangulation< dim, spacedim >::cell_iterator &cell)> pre_coarsening_on_cell
#define DeclException2(Exception2, type1, type2, outsequence)
boost::signals2::signal< void()> post_distributed_repartition
typename IteratorSelector::line_iterator line_iterator
T operator()(InputIterator first, InputIterator last) const
unsigned int n_active_lines
MeshSmoothing smooth_grid
typename IteratorSelector::raw_line_iterator raw_line_iterator
std::vector< Point< spacedim > > vertices
std::map< types::manifold_id, std::unique_ptr< const Manifold< dim, spacedim > > > manifold
typename IteratorSelector::hex_iterator hex_iterator
typename IteratorSelector::active_line_iterator active_line_iterator
boost::signals2::signal< void()> pre_distributed_refinement
bool anisotropic_refinement
#define AssertIndexRange(index, range)
std::unique_ptr< ::internal::TriangulationImplementation::Policy< dim, spacedim > > policy
std::vector< GridTools::PeriodicFacePair< cell_iterator > > periodic_face_pairs_level_0
std::vector< ReferenceCell > reference_cells
static ::ExceptionBase & ExcFacesHaveNoLevel()
unsigned int n_levels() const
void create_triangulation(Triangulation< dim, dim > &tria, const AdditionalData &additional_data=AdditionalData())
unsigned int n_active_quads
static ::ExceptionBase & ExcMessage(std::string arg1)
boost::signals2::signal< void()> pre_refinement
const bool check_for_distorted_cells
#define DeclException1(Exception1, type1, outsequence)
unsigned int n_active_hexes
unsigned int n_raw_quads() const
#define Assert(cond, exc)
std::vector< std::unique_ptr<::internal::TriangulationImplementation::TriaLevel > > levels
boost::signals2::signal< void()> clear
unsigned int n_raw_hexs(const unsigned int level) const
unsigned int max_adjacent_cells() const
void load(Archive &ar, const unsigned int version)
std::list< typename Triangulation< dim, spacedim >::cell_iterator > distorted_cells
unsigned int n_quads() const
#define DeclExceptionMsg(Exception, defaulttext)
std::unique_ptr< std::map< unsigned int, types::manifold_id > > vertex_to_manifold_id_map_1d
std::vector< bool > vertices_used
#define DeclException0(Exception0)
virtual types::coarse_cell_id coarse_cell_index_to_coarse_cell_id(const unsigned int coarse_cell_index) const
boost::signals2::signal< void(const typename Triangulation< dim, spacedim >::cell_iterator &cell)> post_refinement_on_cell
typename IteratorSelector::active_quad_iterator active_quad_iterator
const std::vector< Point< spacedim > > & get_vertices() const
#define DEAL_II_NAMESPACE_CLOSE
boost::signals2::signal< void()> mesh_movement
typename IteratorSelector::raw_quad_iterator raw_quad_iterator
VectorType::value_type * end(VectorType &V)
std::unique_ptr<::internal::TriangulationImplementation::TriaFaces > faces
boost::signals2::signal< void()> pre_distributed_save
void save(Archive &ar, const unsigned int version) const
unsigned int n_active_cells(const internal::TriangulationImplementation::NumberCache< 1 > &c)
typename IteratorSelector::quad_iterator quad_iterator
TriangulationBase< dim, spacedim > Triangulation
boost::signals2::signal< void()> create
std::unique_ptr< std::map< unsigned int, types::boundary_id > > vertex_to_boundary_id_map_1d
unsigned int n_active_quads() const
unsigned int n_cells(const internal::TriangulationImplementation::NumberCache< 1 > &c)
std::vector< unsigned int > n_lines_level
boost::signals2::signal< void()> pre_partition
typename IteratorSelector::active_hex_iterator active_hex_iterator
boost::signals2::signal< void()> post_distributed_save
std::vector< unsigned int > n_active_hexes_level
std::map< std::pair< cell_iterator, unsigned int >, std::pair< std::pair< cell_iterator, unsigned int >, std::bitset< 3 > > > periodic_face_map
#define DEAL_II_NAMESPACE_OPEN
virtual unsigned int n_global_levels() const
VectorType::value_type * begin(VectorType &V)
typename IteratorSelector::raw_hex_iterator raw_hex_iterator
boost::signals2::signal< unsigned int(const cell_iterator &, const CellStatus), CellWeightSum< unsigned int > > cell_weight
boost::signals2::signal< void()> pre_distributed_load
std::vector< unsigned int > n_active_lines_level
std::vector< unsigned int > n_active_quads_level
boost::signals2::signal< void()> post_distributed_load
global_cell_index coarse_cell_id
std::vector< unsigned int > n_hexes_level
boost::signals2::signal< void()> pre_distributed_repartition
bool vertex_used(const unsigned int index) const
boost::signals2::signal< void()> any_change
std::enable_if< std::is_fundamental< T >::value, std::size_t >::type memory_consumption(const T &t)
boost::signals2::signal< void()> post_distributed_refinement