35 #include <boost/serialization/map.hpp>
36 #include <boost/serialization/split_member.hpp>
37 #include <boost/serialization/unique_ptr.hpp>
38 #include <boost/serialization/vector.hpp>
39 #include <boost/signals2.hpp>
54 "The name 'signals' is already defined. You are most likely using the QT library \
55 and using the 'signals' keyword. You can either #include the Qt headers (or any conflicting headers) \
56 *after* the deal.II headers or you can define the 'QT_NO_KEYWORDS' macro and use the 'Q_SIGNALS' macro."
61 template <
int dim,
int spacedim>
77 template <
typename CellIterator>
78 struct PeriodicFacePair;
81 template <
int,
int,
int>
83 template <
int spacedim>
85 template <
int,
int,
int>
90 namespace TriangulationImplementation
104 struct Implementation;
105 struct ImplementationMixedMesh;
108 namespace TriaAccessorImplementation
110 struct Implementation;
125 namespace TriangulationImplementation
184 std::shared_ptr<const Utilities::MPI::Partitioner>
190 std::vector<std::shared_ptr<const Utilities::MPI::Partitioner>>
210 template <
class Archive>
268 template <
class Archive>
327 template <
class Archive>
337 template <
int dim,
int spacedim = dim>
376 template <
int dim,
int spacedim = dim>
402 const std::vector<cell_relation_t> &cell_relations,
405 &pack_callbacks_fixed,
408 &pack_callbacks_variable,
409 const MPI_Comm &mpi_communicator);
419 unpack_cell_status(std::vector<cell_relation_t> &cell_relations)
const;
435 const std::vector<cell_relation_t> &cell_relations,
436 const unsigned int handle,
440 const boost::iterator_range<std::vector<char>::const_iterator> &)>
441 &unpack_callback)
const;
458 save(
const unsigned int global_first_cell,
459 const unsigned int global_num_cells,
460 const std::string &filename,
461 const MPI_Comm &mpi_communicator)
const;
482 load(
const unsigned int global_first_cell,
483 const unsigned int global_num_cells,
484 const unsigned int local_num_cells,
485 const std::string &filename,
486 const unsigned int n_attached_deserialize_fixed,
487 const unsigned int n_attached_deserialize_variable,
488 const MPI_Comm &mpi_communicator);
1335 template <
int dim,
int spacedim = dim>
1402 limit_level_difference_at_vertices = 0x1,
1423 eliminate_unrefined_islands = 0x2,
1439 patch_level_1 = 0x4,
1463 coarsest_level_1 = 0x8,
1488 allow_anisotropic_smoothing = 0x10,
1521 eliminate_refined_inner_islands = 0x100,
1526 eliminate_refined_boundary_islands = 0x200,
1532 do_not_produce_unrefined_islands = 0x400,
1538 smoothing_on_refinement =
1539 (limit_level_difference_at_vertices | eliminate_unrefined_islands),
1544 smoothing_on_coarsening =
1545 (eliminate_refined_inner_islands | eliminate_refined_boundary_islands |
1546 do_not_produce_unrefined_islands),
1553 maximum_smoothing = 0xffff ^ allow_anisotropic_smoothing
1758 static constexpr
unsigned int dimension = dim;
1763 static constexpr
unsigned int space_dimension = spacedim;
1780 const
bool check_for_distorted_cells = false);
1832 get_communicator() const;
1840 global_active_cell_index_partitioner() const;
1848 global_level_cell_index_partitioner(const
unsigned int level) const;
1861 get_mesh_smoothing() const;
1887 const
Manifold<dim, spacedim> &manifold_object);
1916 reset_all_manifolds();
1980 get_boundary_ids() const;
1995 get_manifold_ids() const;
2076 const std::vector<
CellData<dim>> &cells,
2094 &construction_data);
2103 flip_all_direction_flags();
2117 set_all_refine_flags();
2142 refine_global(const
unsigned int times = 1);
2157 coarsen_global(const
unsigned int times = 1);
2191 execute_coarsening_and_refinement();
2223 prepare_coarsening_and_refinement();
2275 template <typename
T>
2280 template <
typename InputIterator>
2284 return std::accumulate(
first, last,
T());
2348 boost::signals2::signal<void(
2358 boost::signals2::signal<void(
2368 boost::signals2::signal<void(
3367 virtual unsigned int
3400 const std::vector<Point<spacedim>> &
3420 const std::vector<bool> &
3581 template <
class Archive>
3583 save(Archive &ar,
const unsigned int version)
const;
3602 template <
class Archive>
3604 load(Archive &ar,
const unsigned int version);
3614 save(
const std::string &filename)
const;
3646 std::pair<cell_iterator, unsigned int>,
3647 std::pair<std::pair<cell_iterator, unsigned int>, std::bitset<3>>> &
3654 const std::vector<ReferenceCell> &
3685 template <
class Archive>
3691 BOOST_SERIALIZATION_SPLIT_MEMBER()
3808 const std::function<std::vector<char>(
const cell_iterator &,
3811 const bool returns_variable_size_data);
3864 const unsigned int handle,
3865 const std::function<
3868 const boost::iterator_range<std::vector<char>::const_iterator> &)>
3882 const unsigned int global_num_cells,
3883 const std::string &filename)
const;
3894 const unsigned int global_num_cells,
3895 const unsigned int local_num_cells,
3896 const std::string &filename,
3897 const unsigned int n_attached_deserialize_fixed,
3898 const unsigned int n_attached_deserialize_variable);
3919 std::vector<typename internal::CellAttachedDataSerializer<dim, spacedim>::
3942 <<
"You are requesting information from refinement level "
3944 <<
" of a triangulation, but this triangulation only has "
3945 << arg2 <<
" refinement levels. The given level " << arg1
3946 <<
" must be *less* than " << arg2 <<
'.');
3954 ExcTriangulationNotEmpty,
3957 <<
"You are trying to perform an operation on a triangulation "
3958 <<
"that is only allowed if the triangulation is currently empty. "
3959 <<
"However, it currently stores " << arg1 <<
" vertices and has "
3960 <<
"cells on " << arg2 <<
" levels.");
3979 <<
"You tried to do something on level " << arg1
3980 <<
", but this level is empty.");
3997 <<
"The given boundary_id " << arg1
3998 <<
" is not defined in this Triangulation!");
4006 ExcInconsistentCoarseningFlags,
4007 "A cell is flagged for coarsening, but either not all of its siblings "
4008 "are active or flagged for coarsening as well. Please clean up all "
4009 "coarsen flags on your triangulation via "
4010 "Triangulation::prepare_coarsening_and_refinement() beforehand!");
4042 const std::vector<bool> &v,
4043 const unsigned int magic_number2,
4052 std::vector<bool> &v,
4053 const unsigned int magic_number2,
4085 std::vector<GridTools::PeriodicFacePair<cell_iterator>>
4092 std::map<std::pair<cell_iterator, unsigned int>,
4093 std::pair<std::pair<cell_iterator, unsigned int>, std::bitset<3>>>
4439 virtual unsigned int
4458 const unsigned int coarse_cell_index)
const;
4465 std::unique_ptr<::internal::TriangulationImplementation::TriaLevel>>
4473 std::unique_ptr<::internal::TriangulationImplementation::TriaFaces>
4491 std::map<types::manifold_id, std::unique_ptr<const Manifold<dim, spacedim>>>
4531 std::unique_ptr<std::map<unsigned int, types::boundary_id>>
4554 std::unique_ptr<std::map<unsigned int, types::manifold_id>>
4558 template <
int,
int,
int>
4560 template <
int,
int,
int>
4566 friend struct ::internal::TriaAccessorImplementation::Implementation;
4568 friend struct ::internal::TriangulationImplementation::Implementation;
4569 friend struct ::internal::TriangulationImplementation::
4570 ImplementationMixedMesh;
4572 friend class ::internal::TriangulationImplementation::TriaObjects;
4576 #ifndef DEAL_II_MSVC
4577 static_assert(dim <= spacedim,
4578 "The dimension <dim> of a Triangulation must be less than or "
4579 "equal to the space dimension <spacedim> in which it lives.");
4590 namespace TriangulationImplementation
4592 template <
class Archive>
4597 ar &n_lines &n_lines_level;
4598 ar &n_active_lines &n_active_lines_level;
4602 template <
class Archive>
4608 ar &n_quads &n_quads_level;
4609 ar &n_active_quads &n_active_quads_level;
4613 template <
class Archive>
4619 ar &n_hexes &n_hexes_level;
4620 ar &n_active_hexes &n_active_hexes_level;
4627 template <
int dim,
int spacedim>
4630 const unsigned int index)
const
4638 template <
int dim,
int spacedim>
4645 template <
int dim,
int spacedim>
4653 template <
int dim,
int spacedim>
4662 template <
int dim,
int spacedim>
4664 inline const std::vector<Point<spacedim>>
4671 template <
int dim,
int spacedim>
4673 template <
class Archive>
4688 bool faces_is_nullptr = (
faces.get() ==
nullptr);
4689 ar &faces_is_nullptr;
4690 if (!faces_is_nullptr)
4710 template <
int dim,
int spacedim>
4712 template <
class Archive>
4725 for (
auto &level_ :
levels)
4727 std::unique_ptr<internal::TriangulationImplementation::TriaLevel>
level;
4729 level_ = std::move(
level);
4733 bool faces_is_nullptr =
true;
4734 ar &faces_is_nullptr;
4735 if (!faces_is_nullptr)
4750 level->active_cell_indices.resize(
level->refine_flags.size());
4751 level->global_active_cell_indices.resize(
level->refine_flags.size());
4752 level->global_level_cell_indices.resize(
level->refine_flags.size());
4761 bool my_check_for_distorted_cells;
4762 ar &my_check_for_distorted_cells;
4765 ExcMessage(
"The triangulation loaded into here must have the "
4766 "same setting with regard to reporting distorted "
4767 "cell as the one previously stored."));
4783 template <
int dim,
int spacedim>
4794 template <
int dim,
int spacedim>
4798 const unsigned int coarse_cell_index)
const
4800 return coarse_cell_index;
@ children_will_be_coarsened
unsigned int n_hexs(const unsigned int level) const
void save_user_flags_line(std::vector< bool > &v) const
virtual types::global_cell_index n_global_active_cells() const
const std::vector< bool > & get_used_vertices() const
quad_iterator begin_quad(const unsigned int level=0) const
typename IteratorSelector::raw_line_iterator raw_line_iterator
active_vertex_iterator begin_active_vertex() const
void load_user_indices_quad(const std::vector< unsigned int > &v)
unsigned int n_quads() const
void load_user_indices(const std::vector< unsigned int > &v)
std::vector< bool > vertices_used
bool anisotropic_refinement
active_quad_iterator begin_active_quad(const unsigned int level=0) const
bool get_anisotropic_refinement_flag() const
virtual types::coarse_cell_id n_global_coarse_cells() const
std::unique_ptr< std::map< unsigned int, types::manifold_id > > vertex_to_manifold_id_map_1d
void save_user_pointers_quad(std::vector< void * > &v) const
void save_user_flags_hex(std::ostream &out) const
void clear_user_flags_quad()
unsigned int n_faces() const
active_hex_iterator begin_active_hex(const unsigned int level=0) const
static void read_bool_vector(const unsigned int magic_number1, std::vector< bool > &v, const unsigned int magic_number2, std::istream &in)
unsigned int n_active_lines(const unsigned int level) const
bool all_reference_cells_are_hyper_cube() const
void load_user_flags_line(std::istream &in)
raw_hex_iterator begin_raw_hex(const unsigned int level=0) const
void save_user_flags_line(std::ostream &out) const
active_cell_iterator last_active() const
void save(Archive &ar, const unsigned int version) const
std::vector< typename internal::CellAttachedDataSerializer< dim, spacedim >::cell_relation_t > local_cell_relations
void reset_global_cell_indices()
face_iterator end_face() const
void reset_active_cell_indices()
cell_iterator create_cell_iterator(const CellId &cell_id) const
cell_iterator begin(const unsigned int level=0) const
std::map< std::pair< cell_iterator, unsigned int >, std::pair< std::pair< cell_iterator, unsigned int >, std::bitset< 3 > > > periodic_face_map
unsigned int n_active_cells(const unsigned int level) const
unsigned int n_cells(const unsigned int level) const
void save_user_pointers_line(std::vector< void * > &v) const
void save_user_flags_quad(std::vector< bool > &v) const
void load_refine_flags(std::istream &in)
void save_user_indices_line(std::vector< unsigned int > &v) const
raw_cell_iterator begin_raw(const unsigned int level=0) const
unsigned int n_lines() const
void load_coarsen_flags(const std::vector< bool > &v)
unsigned int n_raw_lines() const
virtual std::size_t memory_consumption() const
std::vector< Point< spacedim > > vertices
raw_quad_iterator begin_raw_quad(const unsigned int level=0) const
virtual types::subdomain_id locally_owned_subdomain() const
unsigned int n_raw_faces() const
unsigned int n_active_faces() const
const bool check_for_distorted_cells
raw_cell_iterator end_raw(const unsigned int level) const
line_iterator end_line() const
std::unique_ptr< std::map< unsigned int, types::boundary_id > > vertex_to_boundary_id_map_1d
void load_user_flags_quad(std::istream &in)
unsigned int n_active_cells() const
virtual void update_reference_cells()
std::vector< ReferenceCell > reference_cells
void update_periodic_face_map()
unsigned int n_raw_lines(const unsigned int level) const
void clear_despite_subscriptions()
void save_user_flags(std::ostream &out) const
void load_user_flags_hex(std::istream &in)
void load_user_pointers_quad(const std::vector< void * > &v)
std::unique_ptr<::internal::TriangulationImplementation::TriaFaces > faces
virtual void add_periodicity(const std::vector< GridTools::PeriodicFacePair< cell_iterator >> &)
unsigned int n_used_vertices() const
void reset_cell_vertex_indices_cache()
unsigned int n_active_lines() const
unsigned int n_levels() const
void load_user_flags_hex(const std::vector< bool > &v)
void load_user_indices_line(const std::vector< unsigned int > &v)
void clear_user_flags_hex()
void save_user_pointers_hex(std::vector< void * > &v) const
typename IteratorSelector::raw_quad_iterator raw_quad_iterator
void load_user_pointers(const std::vector< void * > &v)
unsigned int register_data_attach(const std::function< std::vector< char >(const cell_iterator &, const ::CellStatus)> &pack_callback, const bool returns_variable_size_data)
::internal::TriangulationImplementation::NumberCache< dim > number_cache
void load_user_flags_quad(const std::vector< bool > &v)
void load_user_flags_line(const std::vector< bool > &v)
void save_user_indices_hex(std::vector< unsigned int > &v) const
DistortedCellList execute_refinement()
active_line_iterator begin_active_line(const unsigned int level=0) const
void save_user_indices_quad(std::vector< unsigned int > &v) const
void load_user_pointers_hex(const std::vector< void * > &v)
cell_iterator end() const
virtual bool has_hanging_nodes() const
std::vector< GridTools::PeriodicFacePair< cell_iterator > > periodic_face_pairs_level_0
unsigned int n_raw_cells(const unsigned int level) const
cell_iterator end(const unsigned int level) const
bool contains_cell(const CellId &cell_id) const
virtual void save(const std::string &filename) const
virtual types::coarse_cell_id coarse_cell_index_to_coarse_cell_id(const unsigned int coarse_cell_index) const
unsigned int n_raw_quads(const unsigned int level) const
internal::CellAttachedData< dim, spacedim > cell_attached_data
void load_coarsen_flags(std::istream &out)
quad_iterator end_quad() const
line_iterator begin_line(const unsigned int level=0) const
unsigned int max_adjacent_cells() const
vertex_iterator begin_vertex() const
void save_attached_data(const unsigned int global_first_cell, const unsigned int global_num_cells, const std::string &filename) const
unsigned int n_hexs() const
vertex_iterator end_vertex() const
bool vertex_used(const unsigned int index) const
void load_user_pointers_line(const std::vector< void * > &v)
void save_user_flags(std::vector< bool > &v) const
hex_iterator end_hex() const
const Triangulation< dim, spacedim > & get_triangulation() const
virtual void load(const std::string &filename)
hex_iterator begin_hex(const unsigned int level=0) const
virtual unsigned int n_global_levels() const
active_cell_iterator end_active(const unsigned int level) const
bool is_mixed_mesh() const
virtual void update_cell_relations()
cell_iterator last() const
unsigned int n_active_quads() const
void save_coarsen_flags(std::vector< bool > &v) const
void load_user_indices_hex(const std::vector< unsigned int > &v)
unsigned int n_raw_quads() const
void save_user_pointers(std::vector< void * > &v) const
face_iterator begin_face() const
unsigned int n_cells() const
virtual bool prepare_coarsening_and_refinement()
unsigned int n_active_quads(const unsigned int level) const
void load_user_flags(const std::vector< bool > &v)
typename IteratorSelector::raw_hex_iterator raw_hex_iterator
const std::map< std::pair< cell_iterator, unsigned int >, std::pair< std::pair< cell_iterator, unsigned int >, std::bitset< 3 > > > & get_periodic_face_map() const
std::map< types::manifold_id, std::unique_ptr< const Manifold< dim, spacedim > > > manifolds
void serialize(Archive &archive, const unsigned int version)
void load_refine_flags(const std::vector< bool > &v)
MeshSmoothing smooth_grid
void save_refine_flags(std::ostream &out) const
std::unique_ptr< ::internal::TriangulationImplementation::Policy< dim, spacedim > > policy
Triangulation< dim, spacedim > & get_triangulation()
void save_user_flags_quad(std::ostream &out) const
internal::CellAttachedDataSerializer< dim, spacedim > data_serializer
unsigned int n_active_hexs(const unsigned int level) const
unsigned int n_vertices() const
void load(Archive &ar, const unsigned int version)
void save_user_indices(std::vector< unsigned int > &v) const
void notify_ready_to_unpack(const unsigned int handle, const std::function< void(const cell_iterator &, const ::CellStatus, const boost::iterator_range< std::vector< char >::const_iterator > &)> &unpack_callback)
void save_user_flags_hex(std::vector< bool > &v) const
bool all_reference_cells_are_simplex() const
std::vector< std::unique_ptr<::internal::TriangulationImplementation::TriaLevel > > levels
unsigned int n_raw_hexs(const unsigned int level) const
unsigned int n_lines(const unsigned int level) const
unsigned int n_active_hexs() const
virtual unsigned int coarse_cell_id_to_coarse_cell_index(const types::coarse_cell_id coarse_cell_id) const
unsigned int n_quads(const unsigned int level) const
void load_user_flags(std::istream &in)
const std::vector< ReferenceCell > & get_reference_cells() const
void save_coarsen_flags(std::ostream &out) const
active_face_iterator begin_active_face() const
void clear_user_flags_line()
raw_line_iterator begin_raw_line(const unsigned int level=0) const
static void write_bool_vector(const unsigned int magic_number1, const std::vector< bool > &v, const unsigned int magic_number2, std::ostream &out)
const std::vector< Point< spacedim > > & get_vertices() const
active_cell_iterator begin_active(const unsigned int level=0) const
void execute_coarsening()
void load_attached_data(const unsigned int global_first_cell, const unsigned int global_num_cells, const unsigned int local_num_cells, const std::string &filename, const unsigned int n_attached_deserialize_fixed, const unsigned int n_attached_deserialize_variable)
void save_refine_flags(std::vector< bool > &v) const
std::vector< char > dest_data_variable
std::vector< char > src_data_fixed
std::vector< char > dest_data_fixed
bool variable_size_data_stored
std::vector< int > src_sizes_variable
std::vector< unsigned int > sizes_fixed_cumulative
typename std::pair< cell_iterator, CellStatus > cell_relation_t
std::vector< int > dest_sizes_variable
std::vector< char > src_data_variable
#define DEAL_II_NAMESPACE_OPEN
#define DEAL_II_CXX20_REQUIRES(condition)
#define DEAL_II_NAMESPACE_CLOSE
IteratorRange< active_cell_iterator > active_cell_iterators() const
IteratorRange< active_face_iterator > active_face_iterators() const
IteratorRange< active_cell_iterator > active_cell_iterators_on_level(const unsigned int level) const
IteratorRange< cell_iterator > cell_iterators() const
IteratorRange< cell_iterator > cell_iterators_on_level(const unsigned int level) const
#define DeclException0(Exception0)
#define Assert(cond, exc)
#define DeclException2(Exception2, type1, type2, outsequence)
#define AssertIndexRange(index, range)
#define DeclExceptionMsg(Exception, defaulttext)
#define DeclException1(Exception1, type1, outsequence)
static ::ExceptionBase & ExcMessage(std::string arg1)
static ::ExceptionBase & ExcFacesHaveNoLevel()
typename IteratorSelector::hex_iterator hex_iterator
typename IteratorSelector::active_quad_iterator active_quad_iterator
typename IteratorSelector::active_hex_iterator active_hex_iterator
typename IteratorSelector::quad_iterator quad_iterator
typename IteratorSelector::line_iterator line_iterator
typename IteratorSelector::active_line_iterator active_line_iterator
void create_triangulation(Triangulation< dim, dim > &tria, const AdditionalData &additional_data=AdditionalData())
std::enable_if_t< std::is_fundamental_v< T >, std::size_t > memory_consumption(const T &t)
global_cell_index coarse_cell_id
unsigned int subdomain_id
unsigned int global_cell_index
T operator()(InputIterator first, InputIterator last) const
virtual ~DistortedCellList() noexcept override
boost::signals2::signal< void()> mesh_movement
boost::signals2::signal< void()> post_distributed_load
boost::signals2::signal< void()> pre_distributed_save
boost::signals2::signal< unsigned int(const cell_iterator &, const ::CellStatus), CellWeightSum< unsigned int > > weight
boost::signals2::signal< void(const Triangulation< dim, spacedim > &destination_tria)> copy
boost::signals2::signal< void()> pre_distributed_refinement
boost::signals2::signal< void()> post_distributed_refinement
boost::signals2::signal< void()> pre_refinement
boost::signals2::signal< void()> any_change
boost::signals2::signal< void()> create
boost::signals2::signal< void()> clear
boost::signals2::signal< void(const typename Triangulation< dim, spacedim >::cell_iterator &cell)> post_refinement_on_cell
boost::signals2::signal< void(const typename Triangulation< dim, spacedim >::cell_iterator &cell)> pre_coarsening_on_cell
boost::signals2::signal< void()> post_refinement
boost::signals2::signal< void()> pre_partition
boost::signals2::signal< void()> pre_distributed_repartition
boost::signals2::signal< void()> pre_distributed_load
boost::signals2::signal< void()> post_p4est_refinement
boost::signals2::signal< void()> post_distributed_repartition
boost::signals2::signal< void()> post_distributed_save
std::vector< pack_callback_t > pack_callbacks_variable
unsigned int n_attached_deserialize
std::vector< pack_callback_t > pack_callbacks_fixed
unsigned int n_attached_data_sets
std::function< std::vector< char >(cell_iterator, CellStatus)> pack_callback_t
std::vector< unsigned int > n_lines_level
unsigned int n_active_lines
std::shared_ptr< const Utilities::MPI::Partitioner > active_cell_index_partitioner
void serialize(Archive &ar, const unsigned int version)
std::vector< unsigned int > n_active_lines_level
std::vector< std::shared_ptr< const Utilities::MPI::Partitioner > > level_cell_index_partitioners
std::vector< unsigned int > n_active_quads_level
unsigned int n_active_quads
void serialize(Archive &ar, const unsigned int version)
std::vector< unsigned int > n_quads_level
unsigned int n_active_hexes
void serialize(Archive &ar, const unsigned int version)
std::vector< unsigned int > n_hexes_level
std::vector< unsigned int > n_active_hexes_level
const ::Triangulation< dim, spacedim > & tria