34 #include <boost/serialization/map.hpp>
35 #include <boost/serialization/split_member.hpp>
36 #include <boost/serialization/unique_ptr.hpp>
37 #include <boost/serialization/vector.hpp>
38 #include <boost/signals2.hpp>
53 "The name 'signals' is already defined. You are most likely using the QT library \
54 and using the 'signals' keyword. You can either #include the Qt headers (or any conflicting headers) \
55 *after* the deal.II headers or you can define the 'QT_NO_KEYWORDS' macro and use the 'Q_SIGNALS' macro."
60 template <
int dim,
int spacedim>
76 template <
typename CellIterator>
77 struct PeriodicFacePair;
80 template <
int,
int,
int>
82 template <
int spacedim>
84 template <
int,
int,
int>
89 namespace TriangulationImplementation
103 struct Implementation;
104 struct ImplementationMixedMesh;
107 namespace TriaAccessorImplementation
109 struct Implementation;
124 namespace TriangulationImplementation
183 std::shared_ptr<const Utilities::MPI::Partitioner>
189 std::vector<std::shared_ptr<const Utilities::MPI::Partitioner>>
209 template <
class Archive>
267 template <
class Archive>
326 template <
class Archive>
1133 template <
int dim,
int spacedim = dim>
1200 limit_level_difference_at_vertices = 0x1,
1221 eliminate_unrefined_islands = 0x2,
1237 patch_level_1 = 0x4,
1261 coarsest_level_1 = 0x8,
1286 allow_anisotropic_smoothing = 0x10,
1319 eliminate_refined_inner_islands = 0x100,
1324 eliminate_refined_boundary_islands = 0x200,
1330 do_not_produce_unrefined_islands = 0x400,
1336 smoothing_on_refinement =
1337 (limit_level_difference_at_vertices | eliminate_unrefined_islands),
1342 smoothing_on_coarsening =
1343 (eliminate_refined_inner_islands | eliminate_refined_boundary_islands |
1344 do_not_produce_unrefined_islands),
1351 maximum_smoothing = 0xffff ^ allow_anisotropic_smoothing
1556 static constexpr
unsigned int dimension = dim;
1561 static constexpr
unsigned int space_dimension = spacedim;
1578 const
bool check_for_distorted_cells = false);
1630 get_communicator() const;
1638 global_active_cell_index_partitioner() const;
1646 global_level_cell_index_partitioner(const
unsigned int level) const;
1660 get_mesh_smoothing() const;
1686 const
Manifold<dim, spacedim> &manifold_object);
1715 reset_all_manifolds();
1779 get_boundary_ids() const;
1794 get_manifold_ids() const;
1874 const std::vector<
CellData<dim>> & cells,
1892 &construction_data);
1904 create_triangulation_compatibility(
1906 const std::vector<
CellData<dim>> & cells,
1916 flip_all_direction_flags();
1930 set_all_refine_flags();
1955 refine_global(const
unsigned int times = 1);
1970 coarsen_global(const
unsigned int times = 1);
2004 execute_coarsening_and_refinement();
2036 prepare_coarsening_and_refinement();
2081 template <
typename T>
2086 template <
typename InputIterator>
2090 return std::accumulate(
first, last,
T());
2154 boost::signals2::signal<void(
2164 boost::signals2::signal<void(
2174 boost::signals2::signal<void(
2250 : cell_weight(weight)
2265 boost::signals2::slot<signature_type, slot_function_type>;
2271 boost::signals2::signal<signature_type, combiner_type> &new_signal)
2272 : new_signal(new_signal)
2280 base_weight.disconnect();
2290 boost::signals2::connection
2293 boost::signals2::connect_position position = boost::signals2::at_back)
2295 if (base_weight.connected() ==
false)
2297 base_weight = new_signal.connect(
2301 Assert(base_weight.connected() && new_signal.num_slots() == 1,
2305 return new_signal.connect(slot, position);
2316 return new_signal.num_slots() -
2317 static_cast<std::size_t
>(base_weight.connected());
2327 if (num_slots() == 0)
2341 template <
typename S>
2345 new_signal.disconnect(connection);
2347 if (num_slots() == 0)
2349 Assert(base_weight.connected() && new_signal.num_slots() == 1,
2351 new_signal.disconnect(base_weight);
2362 return new_signal(iterator, status);
2374 boost::signals2::signal<signature_type, combiner_type> &
new_signal;
3323 virtual unsigned int
3356 const std::vector<Point<spacedim>> &
3376 const std::vector<bool> &
3541 template <
class Archive>
3543 save(Archive &ar,
const unsigned int version)
const;
3562 template <
class Archive>
3564 load(Archive &ar,
const unsigned int version);
3590 std::pair<cell_iterator, unsigned int>,
3591 std::pair<std::pair<cell_iterator, unsigned int>, std::bitset<3>>> &
3598 const std::vector<ReferenceCell> &
3629 template <
class Archive>
3635 BOOST_SERIALIZATION_SPLIT_MEMBER()
3651 <<
"You are requesting information from refinement level "
3653 <<
" of a triangulation, but this triangulation only has "
3654 << arg2 <<
" refinement levels. The given level " << arg1
3655 <<
" must be *less* than " << arg2 <<
'.');
3663 ExcTriangulationNotEmpty,
3666 <<
"You are trying to perform an operation on a triangulation "
3667 <<
"that is only allowed if the triangulation is currently empty. "
3668 <<
"However, it currently stores " << arg1 <<
" vertices and has "
3669 <<
"cells on " << arg2 <<
" levels.");
3688 <<
"You tried to do something on level " << arg1
3689 <<
", but this level is empty.");
3706 <<
"The given boundary_id " << arg1
3707 <<
" is not defined in this Triangulation!");
3715 ExcInconsistentCoarseningFlags,
3716 "A cell is flagged for coarsening, but either not all of its siblings "
3717 "are active or flagged for coarsening as well. Please clean up all "
3718 "coarsen flags on your triangulation via "
3719 "Triangulation::prepare_coarsening_and_refinement() beforehand!");
3753 const std::vector<bool> &v,
3754 const unsigned int magic_number2,
3755 std::ostream & out);
3763 std::vector<bool> &v,
3764 const unsigned int magic_number2,
3796 std::vector<GridTools::PeriodicFacePair<cell_iterator>>
3803 std::map<std::pair<cell_iterator, unsigned int>,
3804 std::pair<std::pair<cell_iterator, unsigned int>, std::bitset<3>>>
4156 virtual unsigned int
4175 const unsigned int coarse_cell_index)
const;
4182 std::unique_ptr<::internal::TriangulationImplementation::TriaLevel>>
4190 std::unique_ptr<::internal::TriangulationImplementation::TriaFaces>
4208 std::map<types::manifold_id, std::unique_ptr<const Manifold<dim, spacedim>>>
4248 std::unique_ptr<std::map<unsigned int, types::boundary_id>>
4271 std::unique_ptr<std::map<unsigned int, types::manifold_id>>
4275 template <
int,
int,
int>
4277 template <
int,
int,
int>
4283 friend struct ::internal::TriaAccessorImplementation::Implementation;
4285 friend struct ::internal::TriangulationImplementation::Implementation;
4286 friend struct ::internal::TriangulationImplementation::
4287 ImplementationMixedMesh;
4289 friend class ::internal::TriangulationImplementation::TriaObjects;
4293 #ifndef DEAL_II_MSVC
4294 static_assert(dim <= spacedim,
4295 "The dimension <dim> of a Triangulation must be less than or "
4296 "equal to the space dimension <spacedim> in which it lives.");
4307 namespace TriangulationImplementation
4309 template <
class Archive>
4314 ar &n_lines &n_lines_level;
4315 ar &n_active_lines &n_active_lines_level;
4319 template <
class Archive>
4325 ar &n_quads &n_quads_level;
4326 ar &n_active_quads &n_active_quads_level;
4330 template <
class Archive>
4336 ar &n_hexes &n_hexes_level;
4337 ar &n_active_hexes &n_active_hexes_level;
4344 template <
int dim,
int spacedim>
4347 const unsigned int index)
const
4355 template <
int dim,
int spacedim>
4362 template <
int dim,
int spacedim>
4370 template <
int dim,
int spacedim>
4379 template <
int dim,
int spacedim>
4381 inline const std::vector<Point<spacedim>>
4388 template <
int dim,
int spacedim>
4390 template <
class Archive>
4405 bool faces_is_nullptr = (
faces.get() ==
nullptr);
4406 ar & faces_is_nullptr;
4407 if (!faces_is_nullptr)
4427 template <
int dim,
int spacedim>
4429 template <
class Archive>
4442 for (
auto &level_ :
levels)
4444 std::unique_ptr<internal::TriangulationImplementation::TriaLevel>
level;
4446 level_ = std::move(
level);
4450 bool faces_is_nullptr =
true;
4451 ar & faces_is_nullptr;
4452 if (!faces_is_nullptr)
4467 level->active_cell_indices.resize(
level->refine_flags.size());
4468 level->global_active_cell_indices.resize(
level->refine_flags.size());
4469 level->global_level_cell_indices.resize(
level->refine_flags.size());
4478 bool my_check_for_distorted_cells;
4479 ar & my_check_for_distorted_cells;
4482 ExcMessage(
"The triangulation loaded into here must have the "
4483 "same setting with regard to reporting distorted "
4484 "cell as the one previously stored."));
4500 template <
int dim,
int spacedim>
4511 template <
int dim,
int spacedim>
4515 const unsigned int coarse_cell_index)
const
4517 return coarse_cell_index;
unsigned int operator()(const cell_iterator &iterator, const CellStatus status)
boost::signals2::slot< signature_type, slot_function_type > slot_type
std::size_t num_slots() const
boost::function< signature_type > slot_function_type
boost::signals2::connection connect(const slot_type &slot, boost::signals2::connect_position position=boost::signals2::at_back)
boost::signals2::signal< signature_type, combiner_type > & new_signal
LegacySignal(boost::signals2::signal< signature_type, combiner_type > &new_signal)
unsigned int(const cell_iterator &, const CellStatus) signature_type
boost::signals2::connection base_weight
void disconnect(const S &connection)
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
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)
::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
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
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
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
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
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
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 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 save_refine_flags(std::vector< bool > &v) const
#define DEAL_II_DEPRECATED
#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)
static ::ExceptionBase & ExcInternalError()
#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< T >::value, 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< unsigned int(const cell_iterator &, const CellStatus), CellWeightSum< unsigned int > > weight
boost::signals2::signal< void()> mesh_movement
boost::signals2::signal< void()> post_distributed_load
boost::signals2::signal< void()> pre_distributed_save
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< 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