#include <deal.II/distributed/shared_tria.h>

Inheritance diagram for parallel::shared::Triangulation< dim, spacedim >:

Public Types
enum	Settings { partition_auto = 0x0 , partition_metis = 0x1 , partition_zorder = 0x2 , partition_zoltan = 0x3 , partition_custom_signal = 0x4 , construct_multigrid_hierarchy = 0x8 }

using	active_cell_iterator = typename ::Triangulation< dim, spacedim >::active_cell_iterator

using	cell_iterator = typename ::Triangulation< dim, spacedim >::cell_iterator

enum	MeshSmoothing { none = 0x0 , limit_level_difference_at_vertices = 0x1 , eliminate_unrefined_islands = 0x2 , patch_level_1 = 0x4 , coarsest_level_1 = 0x8 , allow_anisotropic_smoothing = 0x10 , eliminate_refined_inner_islands = 0x100 , eliminate_refined_boundary_islands = 0x200 , do_not_produce_unrefined_islands = 0x400 , smoothing_on_refinement , smoothing_on_coarsening , maximum_smoothing = 0xffff ^ allow_anisotropic_smoothing }

using	level_cell_iterator = cell_iterator

using	face_iterator = TriaIterator< TriaAccessor< dim - 1, dim, spacedim > >

using	active_face_iterator = TriaActiveIterator< TriaAccessor< dim - 1, dim, spacedim > >

using	vertex_iterator = TriaIterator<::TriaAccessor< 0, dim, spacedim > >

using	active_vertex_iterator = TriaActiveIterator<::TriaAccessor< 0, dim, spacedim > >

using	line_iterator = typename IteratorSelector::line_iterator

using	active_line_iterator = typename IteratorSelector::active_line_iterator

using	quad_iterator = typename IteratorSelector::quad_iterator

using	active_quad_iterator = typename IteratorSelector::active_quad_iterator

using	hex_iterator = typename IteratorSelector::hex_iterator

using	active_hex_iterator = typename IteratorSelector::active_hex_iterator

Public Member Functions
	Triangulation (const MPI_Comm mpi_communicator, const typename ::Triangulation< dim, spacedim >::MeshSmoothing=(::Triangulation< dim, spacedim >::none), const bool allow_artificial_cells=false, const Settings settings=partition_auto)

virtual	~Triangulation () override=default

virtual bool	is_multilevel_hierarchy_constructed () const override

virtual void	execute_coarsening_and_refinement () override

virtual void	create_triangulation (const std::vector< Point< spacedim > > &vertices, const std::vector< CellData< dim > > &cells, const SubCellData &subcelldata) override

virtual void	create_triangulation (const TriangulationDescription::Description< dim, spacedim > &construction_data) override

virtual void	copy_triangulation (const ::Triangulation< dim, spacedim > &other_tria) override

template<class Archive >
void	load (Archive &ar, const unsigned int version)

const std::vector< types::subdomain_id > &	get_true_subdomain_ids_of_cells () const

const std::vector< types::subdomain_id > &	get_true_level_subdomain_ids_of_cells (const unsigned int level) const

bool	with_artificial_cells () const

virtual MPI_Comm	get_communicator () const override

virtual void	copy_triangulation (const Triangulation< dim, spacedim > &other_tria)

unsigned int	n_locally_owned_active_cells () const

virtual types::global_cell_index	n_global_active_cells () const override

virtual std::size_t	memory_consumption () const override

virtual unsigned int	n_global_levels () const override

types::subdomain_id	locally_owned_subdomain () const override

const std::set< types::subdomain_id > &	ghost_owners () const

const std::set< types::subdomain_id > &	level_ghost_owners () const

const std::weak_ptr< const Utilities::MPI::Partitioner >	global_active_cell_index_partitioner () const override

const std::weak_ptr< const Utilities::MPI::Partitioner >	global_level_cell_index_partitioner (const unsigned int level) const override

virtual std::vector< types::boundary_id >	get_boundary_ids () const override

virtual std::vector< types::manifold_id >	get_manifold_ids () const override

void	communicate_locally_moved_vertices (const std::vector< bool > &vertex_locally_moved)

virtual types::coarse_cell_id	n_global_coarse_cells () const override

virtual void	clear ()

virtual void	set_mesh_smoothing (const MeshSmoothing mesh_smoothing)

virtual const MeshSmoothing &	get_mesh_smoothing () const

void	set_manifold (const types::manifold_id number, const Manifold< dim, spacedim > &manifold_object)

void	reset_manifold (const types::manifold_id manifold_number)

void	reset_all_manifolds ()

void	set_all_manifold_ids (const types::manifold_id number)

void	set_all_manifold_ids_on_boundary (const types::manifold_id number)

void	set_all_manifold_ids_on_boundary (const types::boundary_id b_id, const types::manifold_id number)

const Manifold< dim, spacedim > &	get_manifold (const types::manifold_id number) const

virtual void	create_triangulation_compatibility (const std::vector< Point< spacedim > > &vertices, const std::vector< CellData< dim > > &cells, const SubCellData &subcelldata)

void	flip_all_direction_flags ()

bool	prepare_coarsening_and_refinement ()

bool	prepare_coarsening_and_refinement ()

bool	prepare_coarsening_and_refinement ()

unsigned int	n_quads () const

unsigned int	n_quads (const unsigned int) const

unsigned int	n_quads () const

unsigned int	n_quads (const unsigned int) const

unsigned int	n_quads () const

unsigned int	n_quads (const unsigned int) const

unsigned int	n_active_quads (const unsigned int) const

unsigned int	n_active_quads () const

unsigned int	n_active_quads (const unsigned int) const

unsigned int	n_active_quads () const

unsigned int	n_active_quads (const unsigned int) const

unsigned int	n_active_quads () const

unsigned int	n_hexs () const

unsigned int	n_hexs (const unsigned int level) const

unsigned int	n_active_hexs () const

unsigned int	n_active_hexs (const unsigned int level) const

unsigned int	max_adjacent_cells () const

unsigned int	max_adjacent_cells () const

unsigned int	max_adjacent_cells () const

unsigned int	n_raw_quads (const unsigned int) const

unsigned int	n_raw_quads (const unsigned int) const

unsigned int	n_raw_quads (const unsigned int) const

unsigned int	n_raw_quads (const unsigned int level) const

unsigned int	n_raw_quads (const unsigned int level) const

unsigned int	n_raw_quads (const unsigned int) const

unsigned int	n_raw_quads () const

unsigned int	n_raw_hexs (const unsigned int) const

unsigned int	n_raw_hexs (const unsigned int) const

unsigned int	n_raw_hexs (const unsigned int) const

unsigned int	n_raw_hexs (const unsigned int level) const

Mesh refinement
void	set_all_refine_flags ()

void	refine_global (const unsigned int times=1)

void	coarsen_global (const unsigned int times=1)

virtual bool	prepare_coarsening_and_refinement ()

History of a triangulation
void	save_refine_flags (std::ostream &out) const

void	save_refine_flags (std::vector< bool > &v) const

void	load_refine_flags (std::istream &in)

void	load_refine_flags (const std::vector< bool > &v)

void	save_coarsen_flags (std::ostream &out) const

void	save_coarsen_flags (std::vector< bool > &v) const

void	load_coarsen_flags (std::istream &out)

void	load_coarsen_flags (const std::vector< bool > &v)

bool	get_anisotropic_refinement_flag () const

User data
void	clear_user_flags ()

void	save_user_flags (std::ostream &out) const

void	save_user_flags (std::vector< bool > &v) const

void	load_user_flags (std::istream &in)

void	load_user_flags (const std::vector< bool > &v)

void	clear_user_flags_line ()

void	save_user_flags_line (std::ostream &out) const

void	save_user_flags_line (std::vector< bool > &v) const

void	load_user_flags_line (std::istream &in)

void	load_user_flags_line (const std::vector< bool > &v)

void	clear_user_flags_quad ()

void	save_user_flags_quad (std::ostream &out) const

void	save_user_flags_quad (std::vector< bool > &v) const

void	load_user_flags_quad (std::istream &in)

void	load_user_flags_quad (const std::vector< bool > &v)

void	clear_user_flags_hex ()

void	save_user_flags_hex (std::ostream &out) const

void	save_user_flags_hex (std::vector< bool > &v) const

void	load_user_flags_hex (std::istream &in)

void	load_user_flags_hex (const std::vector< bool > &v)

void	clear_user_data ()

void	save_user_indices (std::vector< unsigned int > &v) const

void	load_user_indices (const std::vector< unsigned int > &v)

void	save_user_pointers (std::vector< void * > &v) const

void	load_user_pointers (const std::vector< void * > &v)

void	save_user_indices_line (std::vector< unsigned int > &v) const

void	load_user_indices_line (const std::vector< unsigned int > &v)

void	save_user_indices_quad (std::vector< unsigned int > &v) const

void	load_user_indices_quad (const std::vector< unsigned int > &v)

void	save_user_indices_hex (std::vector< unsigned int > &v) const

void	load_user_indices_hex (const std::vector< unsigned int > &v)

void	save_user_pointers_line (std::vector< void * > &v) const

void	load_user_pointers_line (const std::vector< void * > &v)

void	save_user_pointers_quad (std::vector< void * > &v) const

void	load_user_pointers_quad (const std::vector< void * > &v)

void	save_user_pointers_hex (std::vector< void * > &v) const

void	load_user_pointers_hex (const std::vector< void * > &v)

Cell iterator functions
cell_iterator	begin (const unsigned int level=0) const

active_cell_iterator	begin_active (const unsigned int level=0) const

cell_iterator	end () const

cell_iterator	end (const unsigned int level) const

active_cell_iterator	end_active (const unsigned int level) const

cell_iterator	last () const

active_cell_iterator	last_active () const

cell_iterator	create_cell_iterator (const CellId &cell_id) const

Cell iterator functions returning ranges of iterators
IteratorRange< cell_iterator >	cell_iterators () const

IteratorRange< active_cell_iterator >	active_cell_iterators () const

IteratorRange< cell_iterator >	cell_iterators_on_level (const unsigned int level) const

IteratorRange< active_cell_iterator >	active_cell_iterators_on_level (const unsigned int level) const

Face iterator functions
face_iterator	begin_face () const

active_face_iterator	begin_active_face () const

face_iterator	end_face () const

IteratorRange< active_face_iterator >	active_face_iterators () const

Vertex iterator functions
vertex_iterator	begin_vertex () const

active_vertex_iterator	begin_active_vertex () const

vertex_iterator	end_vertex () const

Information about the triangulation
unsigned int	n_lines () const

unsigned int	n_lines (const unsigned int level) const

unsigned int	n_active_lines () const

unsigned int	n_active_lines (const unsigned int level) const

unsigned int	n_quads () const

unsigned int	n_quads (const unsigned int level) const

unsigned int	n_active_quads () const

unsigned int	n_active_quads (const unsigned int level) const

unsigned int	n_hexs () const

unsigned int	n_hexs (const unsigned int level) const

unsigned int	n_active_hexs () const

unsigned int	n_active_hexs (const unsigned int level) const

unsigned int	n_cells () const

unsigned int	n_cells (const unsigned int level) const

unsigned int	n_active_cells () const

unsigned int	n_active_cells (const unsigned int level) const

unsigned int	n_faces () const

unsigned int	n_active_faces () const

unsigned int	n_levels () const

virtual bool	has_hanging_nodes () const

unsigned int	n_vertices () const

const std::vector< Point< spacedim > > &	get_vertices () const

unsigned int	n_used_vertices () const

bool	vertex_used (const unsigned int index) const

const std::vector< bool > &	get_used_vertices () const

unsigned int	max_adjacent_cells () const

Triangulation< dim, spacedim > &	get_triangulation ()

const Triangulation< dim, spacedim > &	get_triangulation () const

Internal information about the number of objects
unsigned int	n_raw_lines () const

unsigned int	n_raw_lines (const unsigned int level) const

unsigned int	n_raw_quads () const

unsigned int	n_raw_quads (const unsigned int level) const

unsigned int	n_raw_hexs (const unsigned int level) const

unsigned int	n_raw_cells (const unsigned int level) const

unsigned int	n_raw_faces () const

template<class Archive >
void	save (Archive &ar, const unsigned int version) const

virtual void	add_periodicity (const std::vector< GridTools::PeriodicFacePair< cell_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

const std::vector< ReferenceCell > &	get_reference_cells () const

bool	all_reference_cells_are_hyper_cube () const

bool	all_reference_cells_are_simplex () const

bool	is_mixed_mesh () const

template<class Archive >
void	serialize (Archive &archive, const unsigned int version)

Subscriptor functionality
Classes derived from Subscriptor provide a facility to subscribe to this object. This is mostly used by the SmartPointer class.
void	subscribe (std::atomic< bool > *const validity, const std::string &identifier="") const

void	unsubscribe (std::atomic< bool > *const validity, const std::string &identifier="") const

unsigned int	n_subscriptions () const

template<typename StreamType >
void	list_subscribers (StreamType &stream) const

void	list_subscribers () const

Static Public Member Functions
static ::ExceptionBase &	ExcInUse (int arg1, std::string arg2, std::string arg3)

static ::ExceptionBase &	ExcNoSubscriber (std::string arg1, std::string arg2)

Static Public Attributes
static constexpr unsigned int	dimension = dim

static constexpr unsigned int	space_dimension = spacedim

Protected Member Functions
virtual void	update_number_cache ()

void	update_reference_cells () override

void	reset_global_cell_indices ()

Protected Attributes
const MPI_Comm	mpi_communicator

types::subdomain_id	my_subdomain

types::subdomain_id	n_subdomains

NumberCache	number_cache

Private Types
using	IteratorSelector = ::internal::TriangulationImplementation::Iterators< dim, spacedim >

using	map_value_type = decltype(counter_map)::value_type

using	map_iterator = decltype(counter_map)::iterator

Private Member Functions
void	partition ()

void	check_no_subscribers () const noexcept

Line iterator functions for internal use
raw_line_iterator	begin_raw_line (const unsigned int level=0) const

line_iterator	begin_line (const unsigned int level=0) const

active_line_iterator	begin_active_line (const unsigned int level=0) const

line_iterator	end_line () const

Quad iterator functions for internal use
raw_quad_iterator	begin_raw_quad (const unsigned int level=0) const

quad_iterator	begin_quad (const unsigned int level=0) const

active_quad_iterator	begin_active_quad (const unsigned int level=0) const

quad_iterator	end_quad () const

Private Attributes
const Settings	settings

const bool	allow_artificial_cells

std::vector< types::subdomain_id >	true_subdomain_ids_of_cells

std::vector< std::vector< types::subdomain_id > >	true_level_subdomain_ids_of_cells

std::atomic< unsigned int >	counter

std::map< std::string, unsigned int >	counter_map

std::vector< std::atomic< bool > * >	validity_pointers

const std::type_info *	object_info

Static Private Attributes
static std::mutex	mutex

Keeping up with what happens to a triangulation
enum	CellStatus { CELL_PERSIST , CELL_REFINE , CELL_COARSEN , CELL_INVALID }

Signals	signals

Exceptions
std::unique_ptr< ::internal::TriangulationImplementation::Policy< dim, spacedim > >	policy

std::vector< GridTools::PeriodicFacePair< cell_iterator > >	periodic_face_pairs_level_0

std::map< std::pair< cell_iterator, unsigned int >, std::pair< std::pair< cell_iterator, unsigned int >, std::bitset< 3 > > >	periodic_face_map

MeshSmoothing	smooth_grid

std::vector< ReferenceCell >	reference_cells

void	update_periodic_face_map ()

static ::ExceptionBase &	ExcInvalidLevel (int arg1, int arg2)

static ::ExceptionBase &	ExcTriangulationNotEmpty (int arg1, int arg2)

static ::ExceptionBase &	ExcGridReadError ()

static ::ExceptionBase &	ExcFacesHaveNoLevel ()

static ::ExceptionBase &	ExcEmptyLevel (int arg1)

static ::ExceptionBase &	ExcNonOrientableTriangulation ()

static ::ExceptionBase &	ExcBoundaryIdNotFound (types::boundary_id arg1)

static ::ExceptionBase &	ExcInconsistentCoarseningFlags ()

static void	write_bool_vector (const unsigned int magic_number1, const std::vector< bool > &v, const unsigned int magic_number2, std::ostream &out)

static void	read_bool_vector (const unsigned int magic_number1, std::vector< bool > &v, const unsigned int magic_number2, std::istream &in)

Hex iterator functions for internal use
std::vector< std::unique_ptr<::internal::TriangulationImplementation::TriaLevel > >	levels

std::unique_ptr<::internal::TriangulationImplementation::TriaFaces >	faces

std::vector< Point< spacedim > >	vertices

std::vector< bool >	vertices_used

std::map< types::manifold_id, std::unique_ptr< const Manifold< dim, spacedim > > >	manifolds

bool	anisotropic_refinement

const bool	check_for_distorted_cells

std::unique_ptr< std::map< unsigned int, types::boundary_id > >	vertex_to_boundary_id_map_1d

std::unique_ptr< std::map< unsigned int, types::manifold_id > >	vertex_to_manifold_id_map_1d

raw_hex_iterator	begin_raw_hex (const unsigned int level=0) const

hex_iterator	begin_hex (const unsigned int level=0) const

active_hex_iterator	begin_active_hex (const unsigned int level=0) const

hex_iterator	end_hex () const

void	clear_despite_subscriptions ()

void	reset_policy ()

void	reset_active_cell_indices ()

void	reset_cell_vertex_indices_cache ()

DistortedCellList	execute_refinement ()

void	execute_coarsening ()

void	fix_coarsen_flags ()

virtual unsigned int	coarse_cell_id_to_coarse_cell_index (const types::coarse_cell_id coarse_cell_id) const

virtual types::coarse_cell_id	coarse_cell_index_to_coarse_cell_id (const unsigned int coarse_cell_index) const

Cell iterator functions for internal use
raw_cell_iterator	begin_raw (const unsigned int level=0) const

raw_cell_iterator	end_raw (const unsigned int level) const

using	raw_cell_iterator = TriaRawIterator< CellAccessor< dim, spacedim > >

using	raw_face_iterator = TriaRawIterator< TriaAccessor< dim - 1, dim, spacedim > >

using	raw_vertex_iterator = TriaRawIterator<::TriaAccessor< 0, dim, spacedim > >

using	raw_line_iterator = typename IteratorSelector::raw_line_iterator

using	raw_quad_iterator = typename IteratorSelector::raw_quad_iterator

using	raw_hex_iterator = typename IteratorSelector::raw_hex_iterator

Detailed Description

template<int dim, int spacedim = dim>
class parallel::shared::Triangulation< dim, spacedim >

This class provides a parallel triangulation for which every processor knows about every cell of the global mesh (unlike for the parallel::distributed::Triangulation class) but in which cells are automatically partitioned when run with MPI so that each processor "owns" a subset of cells. The use of this class is demonstrated in step-18.

Different from the parallel::distributed::Triangulation and parallel::fullydistributed::Triangulation classes, this implies that the entire mesh is stored on each processor. While this is clearly a memory bottleneck that limits the use of this class to a few dozen or hundreds of MPI processes, the partitioning of the mesh can be used to partition work such as assembly or postprocessing between participating processors, and it can also be used to partition which processor stores which parts of matrices and vectors. As a consequence, using this class is often a gentler introduction to parallelizing a code than the more involved parallel::distributed::Triangulation class in which processors only know their own part of the mesh, but nothing about cells owned by other processors with the exception of a single layer of ghost cells around their own part of the domain.

Because every MPI process has a complete copy of the entire mesh as if it were stored on only this process, it needs to know for every active cell whether it is flagged for refinement or coarsening when doing mesh refinement via Triangulation::execute_coarsening_and_refinement(). In practice, each process only needs to set this information for its own "locally owned" cells; upon calling Triangulation::execute_coarsening_and_refinement(), the relevant information is then exchanged between processes internally, via MPI communication.

The class is also useful in cases where compute time and memory considerations dictate that the program needs to be run in parallel, but where algorithmic concerns require that every processor knows about the entire mesh. An example could be where an application has to have both volume and surface meshes that can then both be partitioned independently, but where it is difficult to ensure that the locally owned set of surface mesh cells is adjacent to the locally owned set of volume mesh cells and the other way around. In such cases, knowing the entirety of both meshes ensures that assembly of coupling terms can be implemented without also implementing overly complicated schemes to transfer information about adjacent cells from processor to processor.

The partitioning of cells between processors is done internally based on a number of different possibilities. By passing appropriate flags to the constructor of this class (see the parallel::shared::Triangulation::Settings enum), it is possible to select different ways of partitioning the mesh, including ways that are dictated by the application and not by the desire to minimize the length of the interface between subdomains owned by processors (as is done by the METIS and Zoltan packages, both of which are options for partitioning). The DoFHandler class knows how to enumerate degrees of freedom in ways appropriate for the partitioned mesh.

Note: This class, function, or variable is a template, and it can only be instantiated if the following condition is true:
(concepts::is_valid_dim_spacedim<dim, spacedim>)

concepts::is_valid_dim_spacedim
Definition template_constraints.h:720

If your compiler supports the C++20 standard, then this constraint will be enforced by a C++20 requires clause.

Definition at line 109 of file shared_tria.h.

Member Typedef Documentation

◆ active_cell_iterator

template<int dim, int spacedim = dim>

using parallel::shared::Triangulation< dim, spacedim >::active_cell_iterator = typename ::Triangulation<dim, spacedim>::active_cell_iterator

Definition at line 113 of file shared_tria.h.

◆ cell_iterator

template<int dim, int spacedim = dim>

using parallel::shared::Triangulation< dim, spacedim >::cell_iterator = typename ::Triangulation<dim, spacedim>::cell_iterator

Definition at line 115 of file shared_tria.h.

◆ IteratorSelector

template<int dim, int spacedim = dim>

using Triangulation< dim, spacedim >::IteratorSelector = ::internal::TriangulationImplementation::Iterators<dim, spacedim>

privateinherited

An internal alias to make the definition of the iterator classes simpler.

Definition at line 1143 of file tria.h.

◆ level_cell_iterator

template<int dim, int spacedim = dim>

using Triangulation< dim, spacedim >::level_cell_iterator = cell_iterator

inherited

The same as above to allow the usage of the "MeshType concept" also on the refinement levels.

Definition at line 1376 of file tria.h.

◆ raw_cell_iterator

template<int dim, int spacedim = dim>

using Triangulation< dim, spacedim >::raw_cell_iterator = TriaRawIterator<CellAccessor<dim, spacedim> >

privateinherited

Declare a number of iterator types for raw iterators, i.e., iterators that also iterate over holes in the list of cells left by cells that have been coarsened away in previous mesh refinement cycles.

Since users should never have to access these internal properties of how we store data, these iterator types are made private.

Definition at line 3821 of file tria.h.

◆ raw_face_iterator

template<int dim, int spacedim = dim>

using Triangulation< dim, spacedim >::raw_face_iterator = TriaRawIterator<TriaAccessor<dim - 1, dim, spacedim> >

privateinherited

Definition at line 3822 of file tria.h.

◆ raw_vertex_iterator

template<int dim, int spacedim = dim>

using Triangulation< dim, spacedim >::raw_vertex_iterator = TriaRawIterator<::TriaAccessor<0, dim, spacedim> >

privateinherited

Definition at line 3824 of file tria.h.

◆ raw_line_iterator

template<int dim, int spacedim = dim>

using Triangulation< dim, spacedim >::raw_line_iterator = typename IteratorSelector::raw_line_iterator

privateinherited

Definition at line 3826 of file tria.h.

◆ raw_quad_iterator

template<int dim, int spacedim = dim>

using Triangulation< dim, spacedim >::raw_quad_iterator = typename IteratorSelector::raw_quad_iterator

privateinherited

Definition at line 3827 of file tria.h.

◆ raw_hex_iterator

template<int dim, int spacedim = dim>

using Triangulation< dim, spacedim >::raw_hex_iterator = typename IteratorSelector::raw_hex_iterator

privateinherited

Definition at line 3828 of file tria.h.

◆ map_value_type

using Subscriptor::map_value_type = decltype(counter_map)::value_type

privateinherited

The data type used in counter_map.

Definition at line 230 of file subscriptor.h.

◆ map_iterator

using Subscriptor::map_iterator = decltype(counter_map)::iterator

privateinherited

The iterator type used in counter_map.

Definition at line 235 of file subscriptor.h.

Member Enumeration Documentation

◆ Settings

template<int dim, int spacedim = dim>

enum parallel::shared::Triangulation::Settings

Configuration flags for distributed Triangulations to be set in the constructor. Settings can be combined using bitwise OR.

The constructor requires that exactly one of partition_auto, partition_metis, partition_zorder, partition_zoltan and partition_custom_signal is set. If partition_auto is chosen, it will use partition_zoltan (if available), then partition_metis (if available) and finally partition_zorder.

Enumerator
partition_auto	Choose the partitioner depending on the enabled dependencies that were found when configuring deal.II. In particular, if the Trilinos package Zoltan was found, then use the `partition_zoltan` strategy. If Zoltan was not found but the METIS package was found, then use the partition_metis strategy. If neither of these were found, then use the partition_zorder partitioning strategy.
partition_metis	Use METIS partitioner to partition active cells.
partition_zorder	Partition active cells with the same scheme used in the p4est library. The term "Z-order" originates in the fact that cells are sorted using a space filling curve which in 2d connects the four children of a cell in the order bottom left, bottom right, top left, top right (i.e., with a curve that looks like a reverse "Z"), and does so recursively on all levels of a triangulation. This is also the order in which children are enumerated by the GeometryInfo class. The "Z-order" is also sometimes called "Morton ordering", see https://en.wikipedia.org/wiki/Z-order_curve . See also Z order glossary entry.
partition_zoltan	Use Zoltan to partition active cells.
partition_custom_signal	Partition cells using a custom, user defined function. This is accomplished by connecting the post_refinement signal to the triangulation whenever it is first created and passing the user defined function through the signal using `std::bind`. Here is an example: template <int dim> void mypartition(parallel::shared::Triangulation<dim> &tria) { // user defined partitioning scheme: assign subdomain_ids // round-robin in a mostly random way: std::vector<unsigned int> assignment = {0,0,1,2,0,0,2,1,0,2,2,1,2,2,0,0}; unsigned int index = 0; for (const auto &cell : tria.active_cell_iterators()) cell->set_subdomain_id(assignment[(index++)%16]); } int main () { parallel::shared::Triangulation<dim> tria( ..., parallel::shared::Triangulation<dim>::partition_custom_signal); tria.signals.post_refinement.connect(std::bind(&mypartition<dim>, std::ref(tria))); } Triangulation::signals Signals signals Definition tria.h:2479 parallel::shared::Triangulation Definition shared_tria.h:111 Triangulation::active_cell_iterators IteratorRange< active_cell_iterator > active_cell_iterators() const Triangulation::Signals::post_refinement boost::signals2::signal< void()> post_refinement Definition tria.h:2129 tria const ::Triangulation< dim, spacedim > & tria Definition tria_description.cc:704 An equivalent code using lambda functions would look like this: int main () { parallel::shared::Triangulation<dim> tria( ..., parallel::shared::Triangulation<dim>::partition_custom_signal); tria.signals.post_refinement.connect ( [&tria]() { // user defined partitioning scheme as above ... }); } Note If you plan to use a custom partition with geometric multigrid, you must manually partition the level cells in addition to the active cells.
construct_multigrid_hierarchy	This flag needs to be set to use the geometric multigrid functionality. This option requires additional computation and communication. Note: This flag should always be set alongside a flag for an active cell partitioning method.

Definition at line 131 of file shared_tria.h.

◆ MeshSmoothing

template<int dim, int spacedim = dim>

enum Triangulation::MeshSmoothing

inherited

Declare some symbolic names for mesh smoothing algorithms. The meaning of these flags is documented in the Triangulation class.

Enumerator
none	No mesh smoothing at all, except that meshes have to remain one- irregular.
limit_level_difference_at_vertices	It can be shown, that degradation of approximation occurs if the triangulation contains vertices which are member of cells with levels differing by more than one. One such example is the following: It would seem that in two space dimensions, the maximum jump in levels between cells sharing a common vertex is two (as in the example above). However, this is not true if more than four cells meet at a vertex. It is not uncommon that a coarse (initial) mesh contains vertices at which six or even eight cells meet, when small features of the domain have to be resolved even on the coarsest mesh. In that case, the maximum difference in levels is three or four, respectively. The problem gets even worse in three space dimensions. Looking at an interpolation of the second derivative of the finite element solution (assuming bilinear finite elements), one sees that the numerical solution is almost totally wrong, compared with the true second derivative. Indeed, on regular meshes, there exist sharp estimations that the H²-error is only of order one, so we should not be surprised; however, the numerical solution may show a value for the second derivative which may be a factor of ten away from the true value. These problems are located on the small cell adjacent to the center vertex, where cells of non-subsequent levels meet, as well as on the upper and right neighbor of this cell (but with a less degree of deviation from the true value). If the smoothing indicator given to the constructor contains the bit for limit_level_difference_at_vertices, situations as the above one are eliminated by also marking the upper right cell for refinement. In case of anisotropic refinement, the level of a cell is not linked to the refinement of a cell as directly as in case of isotropic refinement. Furthermore, a cell can be strongly refined in one direction and not or at least much less refined in another. Therefore, it is very difficult to decide, which cases should be excluded from the refinement process. As a consequence, when using anisotropic refinement, the limit_level_difference_at_vertices flag must not be set. On the other hand, the implementation of multigrid methods in deal.II requires that this bit be set.
eliminate_unrefined_islands	Single cells which are not refined and are surrounded by cells which are refined usually also lead to a sharp decline in approximation properties locally. The reason is that the nodes on the faces between unrefined and refined cells are not real degrees of freedom but carry constraints. The patch without additional degrees of freedom is thus significantly larger then the unrefined cell itself. If in the parameter passed to the constructor the bit for eliminate_unrefined_islands is set, all cells which are not flagged for refinement but which are surrounded by more refined cells than unrefined cells are flagged for refinement. Cells which are not yet refined but flagged for that are accounted for the number of refined neighbors. Cells on the boundary are not accounted for at all. An unrefined island is, by this definition also a cell which (in 2d) is surrounded by three refined cells and one unrefined one, or one surrounded by two refined cells, one unrefined one and is at the boundary on one side. It is thus not a true island, as the name of the flag may indicate. However, no better name came to mind to the author by now.
patch_level_1	A triangulation of patch level 1 consists of patches, i.e. of cells that are refined once. This flag ensures that a mesh of patch level 1 is still of patch level 1 after coarsening and refinement. It is, however, the user's responsibility to ensure that the mesh is of patch level 1 before calling Triangulation::execute_coarsening_and_refinement() the first time. The easiest way to achieve this is by calling global_refine(1) straight after creation of the triangulation. It follows that if at least one of the children of a cell is or will be refined than all children need to be refined. If the patch_level_1 flag is set, than the flags eliminate_unrefined_islands, eliminate_refined_inner_islands and eliminate_refined_boundary_islands will be ignored as they will be fulfilled automatically.
coarsest_level_1	Each coarse grid cell is refined at least once, i.e., the triangulation might have active cells on level 1 but not on level 0. This flag ensures that a mesh which has coarsest_level_1 has still coarsest_level_1 after coarsening and refinement. It is, however, the user's responsibility to ensure that the mesh has coarsest_level_1 before calling execute_coarsening_and_refinement the first time. The easiest way to achieve this is by calling global_refine(1) straight after creation of the triangulation. It follows that active cells on level 1 may not be coarsened. The main use of this flag is to ensure that each cell has at least one neighbor in each coordinate direction (i.e. each cell has at least a left or right, and at least an upper or lower neighbor in 2d). This is a necessary precondition for some algorithms that compute finite differences between cells. The DerivativeApproximation class is one of these algorithms that require that a triangulation is coarsest_level_1 unless all cells already have at least one neighbor in each coordinate direction on the coarsest level.
allow_anisotropic_smoothing	This flag is not included in `maximum_smoothing`. The flag is concerned with the following case: consider the case that an unrefined and a refined cell share a common face and that one of the children of the refined cell along the common face is flagged for further refinement. In that case, the resulting mesh would have more than one hanging node along one or more of the edges of the triangulation, a situation that is not allowed. Consequently, in order to perform the refinement, the coarser of the two original cells is also going to be refined. However, in many cases it is sufficient to refine the coarser of the two original cells in an anisotropic way to avoid the case of multiple hanging vertices on a single edge. Doing only the minimal anisotropic refinement can save cells and degrees of freedom. By specifying this flag, the library can produce these anisotropic refinements. The flag is not included by default since it may lead to anisotropically refined meshes even though no cell has ever been refined anisotropically explicitly by a user command. This surprising fact may lead to programs that do the wrong thing since they are not written for the additional cases that can happen with anisotropic meshes, see the discussion in the introduction to step-30.
eliminate_refined_inner_islands	This algorithm seeks for isolated cells which are refined or flagged for refinement. This definition is unlike that for eliminate_unrefined_islands, which would mean that an island is defined as a cell which is refined but more of its neighbors are not refined than are refined. For example, in 2d, a cell's refinement would be reverted if at most one of its neighbors is also refined (or refined but flagged for coarsening). The reason for the change in definition of an island is, that this option would be a bit dangerous, since if you consider a chain of refined cells (e.g. along a kink in the solution), the cells at the two ends would be coarsened, after which the next outermost cells would need to be coarsened. Therefore, only one loop of flagging cells like this could be done to avoid eating up the whole chain of refined cells (‘chain reaction’...). This algorithm also takes into account cells which are not actually refined but are flagged for refinement. If necessary, it takes away the refinement flag. Actually there are two versions of this flag, eliminate_refined_inner_islands and eliminate_refined_boundary_islands. The first eliminates islands defined by the definition above which are in the interior of the domain, while the second eliminates only those islands if the cell is at the boundary. The reason for this split of flags is that one often wants to eliminate such islands in the interior while those at the boundary may well be wanted, for example if one refines the mesh according to a criterion associated with a boundary integral or if one has rough boundary data.
eliminate_refined_boundary_islands	The result of this flag is very similar to eliminate_refined_inner_islands. See the documentation there.
do_not_produce_unrefined_islands	This flag prevents the occurrence of unrefined islands. In more detail: It prohibits the coarsening of a cell if 'most of the neighbors' will be refined after the step.
smoothing_on_refinement	This flag sums up all smoothing algorithms which may be performed upon refinement by flagging some more cells for refinement.
smoothing_on_coarsening	This flag sums up all smoothing algorithms which may be performed upon coarsening by flagging some more cells for coarsening.
maximum_smoothing	This flag includes all the above ones (therefore combines all smoothing algorithms implemented), with the exception of anisotropic smoothing.

Definition at line 1151 of file tria.h.

◆ CellStatus

template<int dim, int spacedim = dim>

enum Triangulation::CellStatus

inherited

Used to inform functions in derived classes how the cell with the given cell_iterator is going to change. Note that this may me different than the refine_flag() and coarsen_flag() in the cell_iterator in parallel calculations because of refinement constraints that this machine does not see.

Enumerator
CELL_PERSIST	The cell will not be refined or coarsened and might or might not move to a different processor.
CELL_REFINE	The cell will be or was refined.
CELL_COARSEN	The children of this cell will be or were coarsened into this cell.
CELL_INVALID	Invalid status. Will not occur for the user.

Definition at line 2056 of file tria.h.

Constructor & Destructor Documentation

◆ Triangulation()

template<int dim, int spacedim>

Triangulation< dim, spacedim >::Triangulation	(	const MPI_Comm	mpi_communicator,
		const typename ::Triangulation< dim, spacedim >::MeshSmoothing	smooth_grid = `(::Triangulation<dim, spacedim>::none)`,
		const bool	allow_artificial_cells = `false`,
		const Settings	settings = `partition_auto`
	)

Constructor.

The flag allow_artificial_cells can be used to enable artificial cells. If enabled, this class will behave similarly to parallel::distributed::Triangulation and parallel::fullydistributed::Triangulation in the sense that there will be locally owned cells, a single layer of ghost cells, and artificial cells. However, one should not forget that in contrast to those parallel triangulations all cells are duplicated on all processes, leading in most cases to significantly more artificial cells.

If artificial cells are disabled, all non-locally owned cells are considered ghost cells. This might lead to very expensive ghost-value update steps. While in the case of artificial cells, ghost-value updates lead to communication only with the direct process neighbors in a point-to-point fashion, these degenerate to an operation in which every process communicates with every other process (an "all-to-all" communication) if no artificial cells are available. If such ghost-value updates are the bottleneck in your code, you may want to consider enabling artificial cells.

Definition at line 43 of file shared_tria.cc.

◆ ~Triangulation()

template<int dim, int spacedim = dim>

virtual parallel::shared::Triangulation< dim, spacedim >::~Triangulation ( )

overridevirtualdefault

Destructor.

Reimplemented from Triangulation< dim, spacedim >.

Member Function Documentation

◆ is_multilevel_hierarchy_constructed()

template<int dim, int spacedim>

bool Triangulation< dim, spacedim >::is_multilevel_hierarchy_constructed

overridevirtual

Return if multilevel hierarchy is supported and has been constructed.

Implements parallel::TriangulationBase< dim, spacedim >.

Definition at line 78 of file shared_tria.cc.

◆ execute_coarsening_and_refinement()

template<int dim, int spacedim>

void Triangulation< dim, spacedim >::execute_coarsening_and_refinement

overridevirtual

Coarsen and refine the mesh according to refinement and coarsening flags set.

This step is equivalent to the Triangulation class with an addition of calling GridTools::partition_triangulation() at the end.

Reimplemented from Triangulation< dim, spacedim >.

Definition at line 359 of file shared_tria.cc.

◆ create_triangulation() [1/2]

template<int dim, int spacedim>

void Triangulation< dim, spacedim >::create_triangulation	(	const std::vector< Point< spacedim > > &	vertices,
		const std::vector< CellData< dim > > &	cells,
		const SubCellData &	subcelldata
	)

overridevirtual

Create a triangulation.

This function also partitions triangulation based on the MPI communicator provided to the constructor.

Reimplemented from Triangulation< dim, spacedim >.

Definition at line 430 of file shared_tria.cc.

◆ create_triangulation() [2/2]

template<int dim, int spacedim>

void Triangulation< dim, spacedim >::create_triangulation ( const TriangulationDescription::Description< dim, spacedim > & construction_data )

overridevirtual

Create a triangulation.

This function also partitions triangulation based on the MPI communicator provided to the constructor.

Note: Not implemented yet.

Reimplemented from Triangulation< dim, spacedim >.

Definition at line 456 of file shared_tria.cc.

◆ copy_triangulation() [1/2]

template<int dim, int spacedim>

void Triangulation< dim, spacedim >::copy_triangulation ( const ::Triangulation< dim, spacedim > & other_tria )

overridevirtual

Copy other_tria to this triangulation.

This function also partitions triangulation based on the MPI communicator provided to the constructor.

Note: This function can not be used with parallel::distributed::Triangulation, since it only stores those cells that it owns, one layer of ghost cells around the ones it locally owns, and a number of artificial cells.

Reimplemented from parallel::TriangulationBase< dim, spacedim >.

Definition at line 469 of file shared_tria.cc.

◆ load()

template<int dim, int spacedim>

template<class Archive >

void Triangulation< dim, spacedim >::load	(	Archive &	ar,
		const unsigned int	version
	)

Read the data of this object from a stream for the purpose of serialization. Throw away the previous content.

This function first does the same work as in Triangulation::load, then partitions the triangulation based on the MPI communicator provided to the constructor.

Definition at line 413 of file shared_tria.h.

◆ get_true_subdomain_ids_of_cells()

template<int dim, int spacedim>

const std::vector< types::subdomain_id > & Triangulation< dim, spacedim >::get_true_subdomain_ids_of_cells

Return a vector of length Triangulation::n_active_cells() where each element stores the subdomain id of the owner of this cell. The elements of the vector are obviously the same as the subdomain ids for locally owned and ghost cells, but are also correct for artificial cells that do not store who the owner of the cell is in their subdomain_id field.

Definition at line 334 of file shared_tria.cc.

◆ get_true_level_subdomain_ids_of_cells()

template<int dim, int spacedim>

const std::vector< types::subdomain_id > & Triangulation< dim, spacedim >::get_true_level_subdomain_ids_of_cells ( const unsigned int level ) const

Return a vector of length Triangulation::n_cells(level) where each element stores the level subdomain id of the owner of this cell. The elements of the vector are obviously the same as the level subdomain ids for locally owned and ghost cells, but are also correct for artificial cells that do not store who the owner of the cell is in their level_subdomain_id field.

Definition at line 344 of file shared_tria.cc.

◆ with_artificial_cells()

template<int dim, int spacedim>

bool Triangulation< dim, spacedim >::with_artificial_cells

Return allow_artificial_cells, namely true if artificial cells are allowed.

Definition at line 324 of file shared_tria.cc.

◆ partition()

template<int dim, int spacedim>

void Triangulation< dim, spacedim >::partition

private

This function calls GridTools::partition_triangulation () and if requested in the constructor of the class marks artificial cells.

Definition at line 88 of file shared_tria.cc.

◆ get_communicator()

template<int dim, int spacedim>

MPI_Comm parallel::TriangulationBase< dim, spacedim >::get_communicator

overridevirtualinherited

Return MPI communicator used by this triangulation.

Reimplemented from Triangulation< dim, spacedim >.

Definition at line 160 of file tria_base.cc.

◆ copy_triangulation() [2/2]

template<int dim, int spacedim = dim>

virtual void Triangulation< dim, spacedim >::copy_triangulation ( const Triangulation< dim, spacedim > & other_tria )

virtualinherited

Copy other_tria to this triangulation. This operation is not cheap, so you should be careful with using this. We do not implement this function as a copy constructor, since it makes it easier to maintain collections of triangulations if you can assign them values later on.

Keep in mind that this function also copies the pointer to the boundary descriptor previously set by the set_manifold function. You must therefore also guarantee that the Manifold objects describing the boundary have a lifetime at least as long as the copied triangulation.

This triangulation must be empty beforehand.

The function is made virtual since some derived classes might want to disable or extend the functionality of this function.

Note: Calling this function triggers the 'copy' signal on other_tria, i.e. the triangulation being copied from. It also triggers the 'create' signal of the current triangulation. See the section on signals in the general documentation for more information.; The list of connections to signals is not copied from the old to the new triangulation since these connections were established to monitor how the old triangulation changes, not how any triangulation it may be copied to changes.

Reimplemented in PersistentTriangulation< dim, spacedim >.

◆ n_locally_owned_active_cells()

template<int dim, int spacedim>

unsigned int parallel::TriangulationBase< dim, spacedim >::n_locally_owned_active_cells

inherited

Return the number of active cells in the triangulation that are locally owned, i.e. that have a subdomain_id equal to locally_owned_subdomain(). Note that there may be more active cells in the triangulation stored on the present processor, such as for example ghost cells, or cells further away from the locally owned block of cells but that are needed to ensure that the triangulation that stores this processor's set of active cells still remains balanced with respect to the 2:1 size ratio of adjacent cells.

As a consequence of the remark above, the result of this function is always smaller or equal to the result of the function with the same name in the Triangulation base class, which includes the active ghost and artificial cells (see also GlossArtificialCell and GlossGhostCell).

Definition at line 131 of file tria_base.cc.

◆ n_global_active_cells()

template<int dim, int spacedim>

types::global_cell_index parallel::TriangulationBase< dim, spacedim >::n_global_active_cells

overridevirtualinherited

Return the sum over all processors of the number of active cells owned by each processor. This equals the overall number of active cells in the triangulation.

Reimplemented from Triangulation< dim, spacedim >.

Definition at line 151 of file tria_base.cc.

◆ memory_consumption()

template<int dim, int spacedim>

std::size_t parallel::TriangulationBase< dim, spacedim >::memory_consumption

overridevirtualinherited

Return the local memory consumption in bytes.

Reimplemented from Triangulation< dim, spacedim >.

Reimplemented in parallel::fullydistributed::Triangulation< dim, spacedim >, and parallel::distributed::Triangulation< dim, spacedim >.

Definition at line 93 of file tria_base.cc.

◆ n_global_levels()

template<int dim, int spacedim>

unsigned int parallel::TriangulationBase< dim, spacedim >::n_global_levels

overridevirtualinherited

Return the global maximum level. This may be bigger than the number Triangulation::n_levels() (a function in this class's base class) returns if the current processor only stores cells in parts of the domain that are not very refined, but if other processors store cells in more deeply refined parts of the domain.

Reimplemented from Triangulation< dim, spacedim >.

Definition at line 141 of file tria_base.cc.

◆ locally_owned_subdomain()

template<int dim, int spacedim>

types::subdomain_id parallel::TriangulationBase< dim, spacedim >::locally_owned_subdomain

overridevirtualinherited

Return the subdomain id of those cells that are owned by the current processor. All cells in the triangulation that do not have this subdomain id are either owned by another processor or have children that only exist on other processors.

Reimplemented from Triangulation< dim, spacedim >.

Definition at line 345 of file tria_base.cc.

◆ ghost_owners()

template<int dim, int spacedim>

const std::set< types::subdomain_id > & parallel::TriangulationBase< dim, spacedim >::ghost_owners

inherited

Return a set of MPI ranks of the processors that have at least one ghost cell adjacent to the cells of the local processor. In other words, this is the set of subdomain_id() for all ghost cells.

The returned sets are symmetric, that is if i is contained in the list of processor j, then j will also be contained in the list of processor i.

Definition at line 355 of file tria_base.cc.

◆ level_ghost_owners()

template<int dim, int spacedim>

const std::set< types::subdomain_id > & parallel::TriangulationBase< dim, spacedim >::level_ghost_owners

inherited

Return a set of MPI ranks of the processors that have at least one level ghost cell adjacent to our cells used in geometric multigrid. In other words, this is the set of level_subdomain_id() for all level ghost cells.

The returned sets are symmetric, that is if i is contained in the list of processor j, then j will also be contained in the list of processor i.

Note: The level ghost owners can only be determined if the multigrid ownership has been assigned (by setting the construct_multigrid_hierarchy flag at construction time), otherwise the returned set will be empty.

Definition at line 365 of file tria_base.cc.

◆ global_active_cell_index_partitioner()

template<int dim, int spacedim>

const std::weak_ptr< const Utilities::MPI::Partitioner > parallel::TriangulationBase< dim, spacedim >::global_active_cell_index_partitioner ( ) const

overridevirtualinherited

Return the partitioner for the global indices of the cells on the active level of the triangulation, which is returned by the function CellAccessor::global_active_cell_index().

Reimplemented from Triangulation< dim, spacedim >.

Definition at line 639 of file tria_base.cc.

◆ global_level_cell_index_partitioner()

template<int dim, int spacedim>

const std::weak_ptr< const Utilities::MPI::Partitioner > parallel::TriangulationBase< dim, spacedim >::global_level_cell_index_partitioner ( const unsigned int level ) const

overridevirtualinherited

Return the partitioner for the global indices of the cells on the given level of the triangulation, which is returned by the function CellAccessor::global_level_cell_index().

Reimplemented from Triangulation< dim, spacedim >.

Definition at line 650 of file tria_base.cc.

◆ get_boundary_ids()

template<int dim, int spacedim>

std::vector< types::boundary_id > parallel::TriangulationBase< dim, spacedim >::get_boundary_ids

overridevirtualinherited

Note: This function involves a global communication gathering all current IDs from all processes.; This function involves a global communication gathering all current IDs from all processes.

Reimplemented from Triangulation< dim, spacedim >.

Definition at line 374 of file tria_base.cc.

◆ get_manifold_ids()

template<int dim, int spacedim>

std::vector< types::manifold_id > parallel::TriangulationBase< dim, spacedim >::get_manifold_ids

overridevirtualinherited

Note: This function involves a global communication gathering all current IDs from all processes.; This function involves a global communication gathering all current IDs from all processes.

Reimplemented from Triangulation< dim, spacedim >.

Definition at line 386 of file tria_base.cc.

◆ communicate_locally_moved_vertices()

template<int dim, int spacedim>

void parallel::TriangulationBase< dim, spacedim >::communicate_locally_moved_vertices ( const std::vector< bool > & vertex_locally_moved )

inherited

When vertices have been moved locally, for example using code like

cell->vertex(0) = new_location;

then this function can be used to update the location of vertices between MPI processes.

All the vertices that have been moved and might be in the ghost layer of a process have to be reported in the vertex_locally_moved argument. This ensures that that part of the information that has to be send between processes is actually sent. Additionally, it is quite important that vertices on the boundary between processes are reported on exactly one process (e.g. the one with the highest id). Otherwise we could expect undesirable results if multiple processes move a vertex differently. A typical strategy is to let processor \(i\) move those vertices that are adjacent to cells whose owners include processor \(i\) but no other processor \(j\) with \(j<i\); in other words, for vertices at the boundary of a subdomain, the processor with the lowest subdomain id "owns" a vertex.

Note: It only makes sense to move vertices that are either located on locally owned cells or on cells in the ghost layer. This is because you can be sure that these vertices indeed exist on the finest mesh aggregated over all processors, whereas vertices on artificial cells but not at least in the ghost layer may or may not exist on the globally finest mesh. Consequently, the vertex_locally_moved argument may not contain vertices that aren't at least on ghost cells.; This function moves vertices in such a way that on every processor, the vertices of every locally owned and ghost cell is consistent with the corresponding location of these cells on other processors. On the other hand, the locations of artificial cells will in general be wrong since artificial cells may or may not exist on other processors and consequently it is not possible to determine their location in any way. This is not usually a problem since one never does anything on artificial cells. However, it may lead to problems if the mesh with moved vertices is refined in a later step. If that's what you want to do, the right way to do it is to save the offset applied to every vertex, call this function, and before refining or coarsening the mesh apply the opposite offset and call this function again.

Parameters

vertex_locally_moved A bitmap indicating which vertices have been moved. The size of this array must be equal to Triangulation::n_vertices() and must be a subset of those vertices flagged by GridTools::get_locally_owned_vertices().

See also: This function is used, for example, in GridTools::distort_random().

Definition at line 582 of file tria_base.cc.

◆ n_global_coarse_cells()

template<int dim, int spacedim>

types::coarse_cell_id parallel::TriangulationBase< dim, spacedim >::n_global_coarse_cells ( ) const

overridevirtualinherited

Return the total number of coarse cells. If the coarse mesh is replicated on each process, this simply returns n_cells(0).

Reimplemented from Triangulation< dim, spacedim >.

Definition at line 664 of file tria_base.cc.

◆ update_number_cache()

template<int dim, int spacedim>

void parallel::TriangulationBase< dim, spacedim >::update_number_cache

protectedvirtualinherited

Update the number_cache variable after mesh creation or refinement.

Reimplemented in parallel::fullydistributed::Triangulation< dim, spacedim >.

Definition at line 170 of file tria_base.cc.

◆ update_reference_cells()

template<int dim, int spacedim>

void parallel::TriangulationBase< dim, spacedim >::update_reference_cells

overrideprotectedvirtualinherited

Reimplemented from Triangulation< dim, spacedim >.

Definition at line 315 of file tria_base.cc.

◆ reset_global_cell_indices()

template<int dim, int spacedim>

void parallel::TriangulationBase< dim, spacedim >::reset_global_cell_indices

protectedinherited

Reset global active cell indices and global level cell indices.

Definition at line 398 of file tria_base.cc.

◆ clear()

template<int dim, int spacedim = dim>

virtual void Triangulation< dim, spacedim >::clear ( )

virtualinherited

Reset this triangulation into a virgin state by deleting all data.

Note that this operation is only allowed if no subscriptions to this object exist any more, such as DoFHandler objects using it.

Reimplemented in parallel::distributed::Triangulation< dim, spacedim >, parallel::DistributedTriangulationBase< dim, spacedim >, parallel::DistributedTriangulationBase< 1, spacedim >, and parallel::DistributedTriangulationBase< dim, dim >.

◆ set_mesh_smoothing()

template<int dim, int spacedim = dim>

virtual void Triangulation< dim, spacedim >::set_mesh_smoothing ( const MeshSmoothing mesh_smoothing )

virtualinherited

Set the mesh smoothing to mesh_smoothing. This overrides the MeshSmoothing given to the constructor. It is allowed to call this function only if the triangulation is empty.

◆ get_mesh_smoothing()

template<int dim, int spacedim = dim>

virtual const MeshSmoothing & Triangulation< dim, spacedim >::get_mesh_smoothing ( ) const

virtualinherited

Return the mesh smoothing requirements that are obeyed.

◆ create_triangulation_compatibility()

template<int dim, int spacedim = dim>

virtual void Triangulation< dim, spacedim >::create_triangulation_compatibility	(	const std::vector< Point< spacedim > > &	vertices,
		const std::vector< CellData< dim > > &	cells,
		const SubCellData &	subcelldata
	)

virtualinherited

For backward compatibility, only. This function takes the cell data in the ordering as requested by deal.II versions up to 5.2, converts it to the new (lexicographic) ordering and calls create_triangulation().

Note: This function internally calls create_triangulation and therefore can throw the same exception as the other function.

Reimplemented in PersistentTriangulation< dim, spacedim >.

◆ flip_all_direction_flags()

template<int dim, int spacedim = dim>

void Triangulation< dim, spacedim >::flip_all_direction_flags ( )

inherited

Revert or flip the direction_flags of a dim<spacedim triangulation, see GlossDirectionFlag.

This function throws an exception if dim equals spacedim.

◆ set_all_refine_flags()

template<int dim, int spacedim = dim>

void Triangulation< dim, spacedim >::set_all_refine_flags ( )

inherited

Flag all active cells for refinement. This will refine all cells of all levels which are not already refined (i.e. only cells are refined which do not yet have children). The cells are only flagged, not refined, thus you have the chance to save the refinement flags.

◆ refine_global()

template<int dim, int spacedim = dim>

void Triangulation< dim, spacedim >::refine_global ( const unsigned int times = 1 )

inherited

Refine all cells times times. In other words, in each one of the times iterations, loop over all cells and refine each cell uniformly into \(2^\text{dim}\) children. In practice, this function repeats the following operations times times: call set_all_refine_flags() followed by execute_coarsening_and_refinement(). The end result is that the number of cells increases by a factor of \((2^\text{dim})^\text{times}=2^{\text{dim} \times \text{times}}\).

The execute_coarsening_and_refinement() function called in this loop may throw an exception if it creates cells that are distorted (see its documentation for an explanation). This exception will be propagated through this function if that happens, and you may not get the actual number of refinement steps in that case.

Note: This function triggers the pre- and post-refinement signals before and after doing each individual refinement cycle (i.e. more than once if times > 1) . See the section on signals in the general documentation of this class.

◆ coarsen_global()

template<int dim, int spacedim = dim>

void Triangulation< dim, spacedim >::coarsen_global ( const unsigned int times = 1 )

inherited

Coarsen all cells the given number of times.

In each of one of the times iterations, all cells will be marked for coarsening. If an active cell is already on the coarsest level, it will be ignored.

Note: This function triggers the pre- and post-refinement signals before and after doing each individual coarsening cycle (i.e. more than once if times > 1) . See the section on signals in the general documentation of this class.

◆ prepare_coarsening_and_refinement() [1/4]

template<int dim, int spacedim = dim>

virtual bool Triangulation< dim, spacedim >::prepare_coarsening_and_refinement ( )

virtualinherited

Do both preparation for refinement and coarsening as well as mesh smoothing.

Regarding the refinement process it fixes the closure of the refinement in dim>=2 (make sure that no two cells are adjacent with a refinement level differing with more than one), etc. It performs some mesh smoothing if the according flag was given to the constructor of this class. The function returns whether additional cells have been flagged for refinement.

See the general doc of this class for more information on smoothing upon refinement.

Regarding the coarsening part, flagging and deflagging cells in preparation of the actual coarsening step are done. This includes deleting coarsen flags from cells which may not be deleted (e.g. because one neighbor is more refined than the cell), doing some smoothing, etc.

The effect is that only those cells are flagged for coarsening which will actually be coarsened. This includes the fact that all flagged cells belong to parent cells of which all children are flagged.

The function returns whether some cells' flagging has been changed in the process.

This function uses the user flags, so store them if you still need them afterwards.

Reimplemented in parallel::fullydistributed::Triangulation< dim, spacedim >, and parallel::distributed::Triangulation< dim, spacedim >.

◆ prepare_coarsening_and_refinement() [2/4]

bool Triangulation< 1, 1 >::prepare_coarsening_and_refinement ( )

inherited

Definition at line 15168 of file tria.cc.

◆ prepare_coarsening_and_refinement() [3/4]

bool Triangulation< 1, 2 >::prepare_coarsening_and_refinement ( )

inherited

Definition at line 15186 of file tria.cc.

◆ prepare_coarsening_and_refinement() [4/4]

bool Triangulation< 1, 3 >::prepare_coarsening_and_refinement ( )

inherited

Definition at line 15204 of file tria.cc.

◆ save_refine_flags() [1/2]

template<int dim, int spacedim = dim>

void Triangulation< dim, spacedim >::save_refine_flags ( std::ostream & out ) const

inherited

Save the addresses of the cells which are flagged for refinement to out. For usage, read the general documentation for this class.

◆ save_refine_flags() [2/2]

template<int dim, int spacedim = dim>

void Triangulation< dim, spacedim >::save_refine_flags ( std::vector< bool > & v ) const

inherited

Same as above, but store the flags to a bitvector rather than to a file.

◆ load_refine_flags() [1/2]

template<int dim, int spacedim = dim>

void Triangulation< dim, spacedim >::load_refine_flags ( std::istream & in )

inherited

Read the information stored by save_refine_flags.

◆ load_refine_flags() [2/2]

template<int dim, int spacedim = dim>

void Triangulation< dim, spacedim >::load_refine_flags ( const std::vector< bool > & v )

inherited

Read the information stored by save_refine_flags.

◆ save_coarsen_flags() [1/2]

template<int dim, int spacedim = dim>

void Triangulation< dim, spacedim >::save_coarsen_flags ( std::ostream & out ) const

inherited

Analogue to save_refine_flags.

◆ save_coarsen_flags() [2/2]

template<int dim, int spacedim = dim>

void Triangulation< dim, spacedim >::save_coarsen_flags ( std::vector< bool > & v ) const

inherited

Same as above, but store the flags to a bitvector rather than to a file.

◆ load_coarsen_flags() [1/2]

template<int dim, int spacedim = dim>

void Triangulation< dim, spacedim >::load_coarsen_flags ( std::istream & out )

inherited

Analogue to load_refine_flags.

◆ load_coarsen_flags() [2/2]

template<int dim, int spacedim = dim>

void Triangulation< dim, spacedim >::load_coarsen_flags ( const std::vector< bool > & v )

inherited

Analogue to load_refine_flags.

◆ get_anisotropic_refinement_flag()

template<int dim, int spacedim = dim>

bool Triangulation< dim, spacedim >::get_anisotropic_refinement_flag ( ) const

inherited

Return whether this triangulation has ever undergone anisotropic (as opposed to only isotropic) refinement.

◆ clear_user_flags()

template<int dim, int spacedim = dim>

void Triangulation< dim, spacedim >::clear_user_flags ( )

inherited

Clear all user flags. See also GlossUserFlags.

◆ save_user_flags() [1/2]

template<int dim, int spacedim = dim>

void Triangulation< dim, spacedim >::save_user_flags ( std::ostream & out ) const

inherited

Save all user flags. See the general documentation for this class and the documentation for the save_refine_flags for more details. See also GlossUserFlags.

◆ save_user_flags() [2/2]

template<int dim, int spacedim = dim>

void Triangulation< dim, spacedim >::save_user_flags ( std::vector< bool > & v ) const

inherited

Same as above, but store the flags to a bitvector rather than to a file. The output vector is resized if necessary. See also GlossUserFlags.

◆ load_user_flags() [1/2]

template<int dim, int spacedim = dim>

void Triangulation< dim, spacedim >::load_user_flags ( std::istream & in )

inherited

Read the information stored by save_user_flags. See also GlossUserFlags.

◆ load_user_flags() [2/2]

template<int dim, int spacedim = dim>

void Triangulation< dim, spacedim >::load_user_flags ( const std::vector< bool > & v )

inherited

Read the information stored by save_user_flags. See also GlossUserFlags.

◆ clear_user_flags_line()

template<int dim, int spacedim = dim>

void Triangulation< dim, spacedim >::clear_user_flags_line ( )

inherited

Clear all user flags on lines. See also GlossUserFlags.

◆ save_user_flags_line() [1/2]

template<int dim, int spacedim = dim>

void Triangulation< dim, spacedim >::save_user_flags_line ( std::ostream & out ) const

inherited

Save the user flags on lines. See also GlossUserFlags.

◆ save_user_flags_line() [2/2]

template<int dim, int spacedim = dim>

void Triangulation< dim, spacedim >::save_user_flags_line ( std::vector< bool > & v ) const

inherited

Same as above, but store the flags to a bitvector rather than to a file. The output vector is resized if necessary. See also GlossUserFlags.

◆ load_user_flags_line() [1/2]

template<int dim, int spacedim = dim>

void Triangulation< dim, spacedim >::load_user_flags_line ( std::istream & in )

inherited

Load the user flags located on lines. See also GlossUserFlags.

◆ load_user_flags_line() [2/2]

template<int dim, int spacedim = dim>

void Triangulation< dim, spacedim >::load_user_flags_line ( const std::vector< bool > & v )

inherited

Load the user flags located on lines. See also GlossUserFlags.

◆ clear_user_flags_quad()

template<int dim, int spacedim = dim>

void Triangulation< dim, spacedim >::clear_user_flags_quad ( )

inherited

Clear all user flags on quads. See also GlossUserFlags.

◆ save_user_flags_quad() [1/2]

template<int dim, int spacedim = dim>

void Triangulation< dim, spacedim >::save_user_flags_quad ( std::ostream & out ) const

inherited

Save the user flags on quads. See also GlossUserFlags.

◆ save_user_flags_quad() [2/2]

template<int dim, int spacedim = dim>

void Triangulation< dim, spacedim >::save_user_flags_quad ( std::vector< bool > & v ) const

inherited

Same as above, but store the flags to a bitvector rather than to a file. The output vector is resized if necessary. See also GlossUserFlags.

◆ load_user_flags_quad() [1/2]

template<int dim, int spacedim = dim>

void Triangulation< dim, spacedim >::load_user_flags_quad ( std::istream & in )

inherited

Load the user flags located on quads. See also GlossUserFlags.

◆ load_user_flags_quad() [2/2]

template<int dim, int spacedim = dim>

void Triangulation< dim, spacedim >::load_user_flags_quad ( const std::vector< bool > & v )

inherited

Load the user flags located on quads. See also GlossUserFlags.

◆ clear_user_flags_hex()

template<int dim, int spacedim = dim>

void Triangulation< dim, spacedim >::clear_user_flags_hex ( )

inherited

Clear all user flags on quads. See also GlossUserFlags.

◆ save_user_flags_hex() [1/2]

template<int dim, int spacedim = dim>

void Triangulation< dim, spacedim >::save_user_flags_hex ( std::ostream & out ) const

inherited

Save the user flags on hexs. See also GlossUserFlags.

◆ save_user_flags_hex() [2/2]

template<int dim, int spacedim = dim>

void Triangulation< dim, spacedim >::save_user_flags_hex ( std::vector< bool > & v ) const

inherited

Same as above, but store the flags to a bitvector rather than to a file. The output vector is resized if necessary. See also GlossUserFlags.

◆ load_user_flags_hex() [1/2]

template<int dim, int spacedim = dim>

void Triangulation< dim, spacedim >::load_user_flags_hex ( std::istream & in )

inherited

Load the user flags located on hexs. See also GlossUserFlags.

◆ load_user_flags_hex() [2/2]

template<int dim, int spacedim = dim>

void Triangulation< dim, spacedim >::load_user_flags_hex ( const std::vector< bool > & v )

inherited

Load the user flags located on hexs. See also GlossUserFlags.

◆ clear_user_data()

template<int dim, int spacedim = dim>

void Triangulation< dim, spacedim >::clear_user_data ( )

inherited

Clear all user pointers and indices and allow the use of both for next access. See also GlossUserData.

◆ save_user_indices()

template<int dim, int spacedim = dim>

void Triangulation< dim, spacedim >::save_user_indices ( std::vector< unsigned int > & v ) const

inherited

Save all user indices. The output vector is resized if necessary. See also GlossUserData.

◆ load_user_indices()

template<int dim, int spacedim = dim>

void Triangulation< dim, spacedim >::load_user_indices ( const std::vector< unsigned int > & v )

inherited

Read the information stored by save_user_indices(). See also GlossUserData.

◆ save_user_pointers()

template<int dim, int spacedim = dim>

void Triangulation< dim, spacedim >::save_user_pointers ( std::vector< void * > & v ) const

inherited

Save all user pointers. The output vector is resized if necessary. See also GlossUserData.

◆ load_user_pointers()

template<int dim, int spacedim = dim>

void Triangulation< dim, spacedim >::load_user_pointers ( const std::vector< void * > & v )

inherited

Read the information stored by save_user_pointers(). See also GlossUserData.

◆ save_user_indices_line()

template<int dim, int spacedim = dim>

void Triangulation< dim, spacedim >::save_user_indices_line ( std::vector< unsigned int > & v ) const

inherited

Save the user indices on lines. The output vector is resized if necessary. See also GlossUserData.

◆ load_user_indices_line()

template<int dim, int spacedim = dim>

void Triangulation< dim, spacedim >::load_user_indices_line ( const std::vector< unsigned int > & v )

inherited

Load the user indices located on lines. See also GlossUserData.

◆ save_user_indices_quad()

template<int dim, int spacedim = dim>

void Triangulation< dim, spacedim >::save_user_indices_quad ( std::vector< unsigned int > & v ) const

inherited

Save the user indices on quads. The output vector is resized if necessary. See also GlossUserData.

◆ load_user_indices_quad()

template<int dim, int spacedim = dim>

void Triangulation< dim, spacedim >::load_user_indices_quad ( const std::vector< unsigned int > & v )

inherited

Load the user indices located on quads. See also GlossUserData.

◆ save_user_indices_hex()

template<int dim, int spacedim = dim>

void Triangulation< dim, spacedim >::save_user_indices_hex ( std::vector< unsigned int > & v ) const

inherited

Save the user indices on hexes. The output vector is resized if necessary. See also GlossUserData.

◆ load_user_indices_hex()

template<int dim, int spacedim = dim>

void Triangulation< dim, spacedim >::load_user_indices_hex ( const std::vector< unsigned int > & v )

inherited

Load the user indices located on hexs. See also GlossUserData.

◆ save_user_pointers_line()

template<int dim, int spacedim = dim>

void Triangulation< dim, spacedim >::save_user_pointers_line ( std::vector< void * > & v ) const

inherited

Save the user indices on lines. The output vector is resized if necessary. See also GlossUserData.

◆ load_user_pointers_line()

template<int dim, int spacedim = dim>

void Triangulation< dim, spacedim >::load_user_pointers_line ( const std::vector< void * > & v )

inherited

Load the user pointers located on lines. See also GlossUserData.

◆ save_user_pointers_quad()

template<int dim, int spacedim = dim>

void Triangulation< dim, spacedim >::save_user_pointers_quad ( std::vector< void * > & v ) const

inherited

Save the user pointers on quads. The output vector is resized if necessary. See also GlossUserData.

◆ load_user_pointers_quad()

template<int dim, int spacedim = dim>

void Triangulation< dim, spacedim >::load_user_pointers_quad ( const std::vector< void * > & v )

inherited

Load the user pointers located on quads. See also GlossUserData.

◆ save_user_pointers_hex()

template<int dim, int spacedim = dim>

void Triangulation< dim, spacedim >::save_user_pointers_hex ( std::vector< void * > & v ) const

inherited

Save the user pointers on hexes. The output vector is resized if necessary. See also GlossUserData.

◆ load_user_pointers_hex()

template<int dim, int spacedim = dim>

void Triangulation< dim, spacedim >::load_user_pointers_hex ( const std::vector< void * > & v )

inherited

Load the user pointers located on hexs. See also GlossUserData.

◆ begin()

template<int dim, int spacedim = dim>

cell_iterator Triangulation< dim, spacedim >::begin ( const unsigned int level = 0 ) const

inherited

Iterator to the first used cell on level level.

Note: The given level argument needs to correspond to a level of the triangulation, i.e., should be less than the value returned by n_levels(). On the other hand, for parallel computations using a parallel::distributed::Triangulation object, it is often convenient to write loops over the cells of all levels of the global mesh, even if the local portion of the triangulation does not actually have cells at one of the higher levels. In those cases, the level argument is accepted if it is less than what the n_global_levels() function returns. If the given level is between the values returned by n_levels() and n_global_levels(), then no cells exist in the local portion of the triangulation at this level, and the function simply returns what end() would return.

◆ begin_active()

template<int dim, int spacedim = dim>

active_cell_iterator Triangulation< dim, spacedim >::begin_active ( const unsigned int level = 0 ) const

inherited

Iterator to the first active cell on level level. If the given level does not contain any active cells (i.e., all cells on this level are further refined, then this function returns end_active(level) so that loops of the kind

for (const auto cell=tria.begin_active(level);
     cell!=tria.end_active(level);
     ++cell)
  {
    ...
  }

have zero iterations, as may be expected if there are no active cells on this level.

Note: The given level argument needs to correspond to a level of the triangulation, i.e., should be less than the value returned by n_levels(). On the other hand, for parallel computations using a parallel::distributed::Triangulation object, it is often convenient to write loops over the cells of all levels of the global mesh, even if the local portion of the triangulation does not actually have cells at one of the higher levels. In those cases, the level argument is accepted if it is less than what the n_global_levels() function returns. If the given level is between the values returned by n_levels() and n_global_levels(), then no cells exist in the local portion of the triangulation at this level, and the function simply returns what end() would return.

◆ end() [1/2]

template<int dim, int spacedim = dim>

cell_iterator Triangulation< dim, spacedim >::end ( ) const

inherited

Iterator past the end; this iterator serves for comparisons of iterators with past-the-end or before-the-beginning states.

◆ end() [2/2]

template<int dim, int spacedim = dim>

cell_iterator Triangulation< dim, spacedim >::end ( const unsigned int level ) const

inherited

Return an iterator which is the first iterator not on level. If level is the last level, then this returns end().

Note: The given level argument needs to correspond to a level of the triangulation, i.e., should be less than the value returned by n_levels(). On the other hand, for parallel computations using a parallel::distributed::Triangulation object, it is often convenient to write loops over the cells of all levels of the global mesh, even if the local portion of the triangulation does not actually have cells at one of the higher levels. In those cases, the level argument is accepted if it is less than what the n_global_levels() function returns. If the given level is between the values returned by n_levels() and n_global_levels(), then no cells exist in the local portion of the triangulation at this level, and the function simply returns what end() would return.

◆ end_active()

template<int dim, int spacedim = dim>

active_cell_iterator Triangulation< dim, spacedim >::end_active ( const unsigned int level ) const

inherited

Return an active iterator which is the first active iterator not on the given level. If level is the last level, then this returns end().

Note: The given level argument needs to correspond to a level of the triangulation, i.e., should be less than the value returned by n_levels(). On the other hand, for parallel computations using a parallel::distributed::Triangulation object, it is often convenient to write loops over the cells of all levels of the global mesh, even if the local portion of the triangulation does not actually have cells at one of the higher levels. In those cases, the level argument is accepted if it is less than what the n_global_levels() function returns. If the given level is between the values returned by n_levels() and n_global_levels(), then no cells exist in the local portion of the triangulation at this level, and the function simply returns what end() would return.

◆ last()

template<int dim, int spacedim = dim>

cell_iterator Triangulation< dim, spacedim >::last ( ) const

inherited

Return an iterator pointing to the last used cell.

◆ last_active()

template<int dim, int spacedim = dim>

active_cell_iterator Triangulation< dim, spacedim >::last_active ( ) const

inherited

Return an iterator pointing to the last active cell.

◆ create_cell_iterator()

template<int dim, int spacedim = dim>

cell_iterator Triangulation< dim, spacedim >::create_cell_iterator ( const CellId & cell_id ) const

inherited

Return an iterator to a cell of this Triangulation object constructed from an independent CellId object.

If the given argument corresponds to a valid cell in this triangulation, this operation will always succeed for sequential triangulations where the current processor stores all cells that are part of the triangulation. On the other hand, if this is a parallel triangulation, then the current processor may not actually know about this cell. In this case, this operation will succeed for locally relevant cells, but may not for artificial cells that are less refined on the current processor.

◆ begin_face()

template<int dim, int spacedim = dim>

face_iterator Triangulation< dim, spacedim >::begin_face ( ) const

inherited

Iterator to the first used face.

◆ begin_active_face()

template<int dim, int spacedim = dim>

active_face_iterator Triangulation< dim, spacedim >::begin_active_face ( ) const

inherited

Iterator to the first active face.

◆ end_face()

template<int dim, int spacedim = dim>

face_iterator Triangulation< dim, spacedim >::end_face ( ) const

inherited

Iterator past the end; this iterator serves for comparisons of iterators with past-the-end or before-the-beginning states.

◆ begin_vertex()

template<int dim, int spacedim = dim>

vertex_iterator Triangulation< dim, spacedim >::begin_vertex ( ) const

inherited

Iterator to the first used vertex. This function can only be used if dim is not one.

◆ begin_active_vertex()

template<int dim, int spacedim = dim>

active_vertex_iterator Triangulation< dim, spacedim >::begin_active_vertex ( ) const

inherited

Iterator to the first active vertex. Because all vertices are active, begin_vertex() and begin_active_vertex() return the same vertex. This function can only be used if dim is not one.

◆ end_vertex()

template<int dim, int spacedim = dim>

vertex_iterator Triangulation< dim, spacedim >::end_vertex ( ) const

inherited

Iterator past the end; this iterator serves for comparisons of iterators with past-the-end or before-the-beginning states. This function can only be used if dim is not one.

◆ n_lines() [1/2]

template<int dim, int spacedim = dim>

unsigned int Triangulation< dim, spacedim >::n_lines ( ) const

inherited

In the following, most functions are provided in two versions, with and without an argument describing the level. The versions with this argument are only applicable for objects describing the cells of the present triangulation. For example: in 2d n_lines(level) cannot be called, only n_lines(), as lines are faces in 2d and therefore have no level. Return the total number of used lines, active or not.

◆ n_lines() [2/2]

template<int dim, int spacedim = dim>

unsigned int Triangulation< dim, spacedim >::n_lines ( const unsigned int level ) const

inherited

Return the total number of used lines, active or not on level level.

◆ n_active_lines() [1/2]

template<int dim, int spacedim = dim>

unsigned int Triangulation< dim, spacedim >::n_active_lines ( ) const

inherited

Return the total number of active lines.

◆ n_active_lines() [2/2]

template<int dim, int spacedim = dim>

unsigned int Triangulation< dim, spacedim >::n_active_lines ( const unsigned int level ) const

inherited

Return the total number of active lines, on level level.

◆ n_quads() [1/8]

template<int dim, int spacedim = dim>

unsigned int Triangulation< dim, spacedim >::n_quads ( ) const

inherited

Return the total number of used quads, active or not.

◆ n_quads() [2/8]

template<int dim, int spacedim = dim>

unsigned int Triangulation< dim, spacedim >::n_quads ( const unsigned int level ) const

inherited

Return the total number of used quads, active or not on level level.

◆ n_quads() [3/8]

unsigned int Triangulation< 1, 1 >::n_quads ( ) const

inherited

Definition at line 14116 of file tria.cc.

◆ n_quads() [4/8]

unsigned int Triangulation< 1, 1 >::n_quads ( const unsigned int ) const

inherited

Definition at line 14124 of file tria.cc.

◆ n_quads() [5/8]

unsigned int Triangulation< 1, 2 >::n_quads ( ) const

inherited

Definition at line 14165 of file tria.cc.

◆ n_quads() [6/8]

unsigned int Triangulation< 1, 2 >::n_quads ( const unsigned int ) const

inherited

Definition at line 14173 of file tria.cc.

◆ n_quads() [7/8]

unsigned int Triangulation< 1, 3 >::n_quads ( ) const

inherited

Definition at line 14213 of file tria.cc.

◆ n_quads() [8/8]

unsigned int Triangulation< 1, 3 >::n_quads ( const unsigned int ) const

inherited

Definition at line 14221 of file tria.cc.

◆ n_active_quads() [1/8]

template<int dim, int spacedim = dim>

unsigned int Triangulation< dim, spacedim >::n_active_quads ( ) const

inherited

Return the total number of active quads, active or not.

◆ n_active_quads() [2/8]

template<int dim, int spacedim = dim>

unsigned int Triangulation< dim, spacedim >::n_active_quads ( const unsigned int level ) const

inherited

Return the total number of active quads, active or not on level level.

◆ n_active_quads() [3/8]

unsigned int Triangulation< 1, 1 >::n_active_quads ( const unsigned int ) const

inherited

Definition at line 14148 of file tria.cc.

◆ n_active_quads() [4/8]

unsigned int Triangulation< 1, 1 >::n_active_quads ( ) const

inherited

Definition at line 14156 of file tria.cc.

◆ n_active_quads() [5/8]

unsigned int Triangulation< 1, 2 >::n_active_quads ( const unsigned int ) const

inherited

Definition at line 14197 of file tria.cc.

◆ n_active_quads() [6/8]

unsigned int Triangulation< 1, 2 >::n_active_quads ( ) const

inherited

Definition at line 14205 of file tria.cc.

◆ n_active_quads() [7/8]

unsigned int Triangulation< 1, 3 >::n_active_quads ( const unsigned int ) const

inherited

Definition at line 14245 of file tria.cc.

◆ n_active_quads() [8/8]

unsigned int Triangulation< 1, 3 >::n_active_quads ( ) const

inherited

Definition at line 14253 of file tria.cc.

◆ n_hexs() [1/4]

template<int dim, int spacedim = dim>

unsigned int Triangulation< dim, spacedim >::n_hexs ( ) const

inherited

Return the total number of used hexahedra, active or not.

◆ n_hexs() [2/4]

template<int dim, int spacedim = dim>

unsigned int Triangulation< dim, spacedim >::n_hexs ( const unsigned int level ) const

inherited

Return the total number of used hexahedra, active or not on level level.

◆ n_hexs() [3/4]

unsigned int Triangulation< 3, 3 >::n_hexs ( ) const

inherited

Definition at line 14394 of file tria.cc.

◆ n_hexs() [4/4]

unsigned int Triangulation< 3, 3 >::n_hexs ( const unsigned int level ) const

inherited

Definition at line 14403 of file tria.cc.

◆ n_active_hexs() [1/4]

template<int dim, int spacedim = dim>

unsigned int Triangulation< dim, spacedim >::n_active_hexs ( ) const

inherited

Return the total number of active hexahedra, active or not.

◆ n_active_hexs() [2/4]

template<int dim, int spacedim = dim>

unsigned int Triangulation< dim, spacedim >::n_active_hexs ( const unsigned int level ) const

inherited

Return the total number of active hexahedra, active or not on level level.

◆ n_active_hexs() [3/4]

unsigned int Triangulation< 3, 3 >::n_active_hexs ( ) const

inherited

Definition at line 14423 of file tria.cc.

◆ n_active_hexs() [4/4]

unsigned int Triangulation< 3, 3 >::n_active_hexs ( const unsigned int level ) const

inherited

Definition at line 14432 of file tria.cc.

◆ n_cells() [1/2]

template<int dim, int spacedim = dim>

unsigned int Triangulation< dim, spacedim >::n_cells ( ) const

inherited

Return the total number of used cells, active or not. Maps to n_lines() in one space dimension and so on.

◆ n_cells() [2/2]

template<int dim, int spacedim = dim>

unsigned int Triangulation< dim, spacedim >::n_cells ( const unsigned int level ) const

inherited

Return the total number of used cells, active or not, on level level. Maps to n_lines(level) in one space dimension and so on.

◆ n_active_cells() [1/2]

template<int dim, int spacedim = dim>

unsigned int Triangulation< dim, spacedim >::n_active_cells ( ) const

inherited

Return the total number of active cells. Maps to n_active_lines() in one space dimension and so on.

◆ n_active_cells() [2/2]

template<int dim, int spacedim = dim>

unsigned int Triangulation< dim, spacedim >::n_active_cells ( const unsigned int level ) const

inherited

Return the total number of active cells on level level. Maps to n_active_lines(level) in one space dimension and so on.

◆ n_faces()

template<int dim, int spacedim = dim>

unsigned int Triangulation< dim, spacedim >::n_faces ( ) const

inherited

Return the total number of used faces, active or not. In 2d, the result equals n_lines(), in 3d it equals n_quads(), while in 1d it equals the number of used vertices.

◆ n_active_faces()

template<int dim, int spacedim = dim>

unsigned int Triangulation< dim, spacedim >::n_active_faces ( ) const

inherited

Return the total number of active faces. In 2d, the result equals n_active_lines(), in 3d it equals n_active_quads(), while in 1d it equals the number of used vertices.

◆ n_levels()

template<int dim, int spacedim = dim>

unsigned int Triangulation< dim, spacedim >::n_levels ( ) const

inherited

Return the number of levels in this triangulation.

Note: Internally, triangulations store data in levels, and there may be more levels in this data structure than one may think – for example, imagine a triangulation that we just got by coarsening the highest level so that it was completely depopulated. That level is not removed, since it will most likely be repopulated soon by the next refinement process. As a consequence, if you happened to run through raw cell iterators (which you can't do as a user of this class, but can internally), then the number of objects in the levels hierarchy is larger than the level of the most refined cell plus one. On the other hand, since this is rarely what a user of this class cares about, the function really just returns the level of the most refined active cell plus one. (The plus one is because in a coarse, unrefined mesh, all cells have level zero – making the number of levels equal to one.)

◆ has_hanging_nodes()

template<int dim, int spacedim = dim>

virtual bool Triangulation< dim, spacedim >::has_hanging_nodes ( ) const

virtualinherited

Return true if the triangulation has hanging nodes.

The function is made virtual since the result can be interpreted in different ways, depending on whether the triangulation lives only on a single processor, or may be distributed as done in the parallel::distributed::Triangulation class (see there for a description of what the function is supposed to do in the parallel context).

Reimplemented in parallel::DistributedTriangulationBase< dim, spacedim >, parallel::DistributedTriangulationBase< 1, spacedim >, and parallel::DistributedTriangulationBase< dim, dim >.

◆ n_vertices()

template<int dim, int spacedim = dim>

unsigned int Triangulation< dim, spacedim >::n_vertices ( ) const

inherited

Return the total number of vertices. Some of them may not be used, which usually happens upon coarsening of a triangulation when some vertices are discarded, but we do not want to renumber the remaining ones, leading to holes in the numbers of used vertices. You can get the number of used vertices using n_used_vertices function.

◆ get_vertices()

template<int dim, int spacedim = dim>

const std::vector< Point< spacedim > > & Triangulation< dim, spacedim >::get_vertices ( ) const

inherited

Return a constant reference to all the vertices present in this triangulation. Note that not necessarily all vertices in this array are actually used; for example, if you coarsen a mesh, then some vertices are deleted, but their positions in this array are unchanged as the indices of vertices are only allocated once. You can find out about which vertices are actually used by the function get_used_vertices().

◆ n_used_vertices()

template<int dim, int spacedim = dim>

unsigned int Triangulation< dim, spacedim >::n_used_vertices ( ) const

inherited

Return the number of vertices that are presently in use, i.e. belong to at least one used element.

◆ vertex_used()

template<int dim, int spacedim = dim>

bool Triangulation< dim, spacedim >::vertex_used ( const unsigned int index ) const

inherited

Return true if the vertex with this index is used.

◆ get_used_vertices()

template<int dim, int spacedim = dim>

const std::vector< bool > & Triangulation< dim, spacedim >::get_used_vertices ( ) const

inherited

Return a constant reference to the array of bools indicating whether an entry in the vertex array is used or not.

◆ max_adjacent_cells() [1/4]

template<int dim, int spacedim = dim>

unsigned int Triangulation< dim, spacedim >::max_adjacent_cells ( ) const

inherited

Return the maximum number of cells meeting at a common vertex. Since this number is an invariant under refinement, only the cells on the coarsest level are considered. The operation is thus reasonably fast. The invariance is only true for sufficiently many cells in the coarsest triangulation (e.g. for a single cell one would be returned), so a minimum of four is returned in two dimensions, 8 in three dimensions, etc, which is how many cells meet if the triangulation is refined.

In one space dimension, two is returned.

◆ max_adjacent_cells() [2/4]

unsigned int Triangulation< 1, 1 >::max_adjacent_cells ( ) const

inherited

Definition at line 14461 of file tria.cc.

◆ max_adjacent_cells() [3/4]

unsigned int Triangulation< 1, 2 >::max_adjacent_cells ( ) const

inherited

Definition at line 14470 of file tria.cc.

◆ max_adjacent_cells() [4/4]

unsigned int Triangulation< 1, 3 >::max_adjacent_cells ( ) const

inherited

Definition at line 14478 of file tria.cc.

◆ get_triangulation() [1/2]

template<int dim, int spacedim = dim>

Triangulation< dim, spacedim > & Triangulation< dim, spacedim >::get_triangulation ( )

inherited

Return a reference to the current object.

This doesn't seem to be very useful but allows to write code that can access the underlying triangulation for anything that satisfies the MeshType concept (which may not only be a triangulation, but also a DoFHandler, for example).

◆ get_triangulation() [2/2]

template<int dim, int spacedim = dim>

const Triangulation< dim, spacedim > & Triangulation< dim, spacedim >::get_triangulation ( ) const

inherited

Return a reference to the current object. This is the const-version of the previous function.

◆ n_raw_lines() [1/2]

template<int dim, int spacedim = dim>

unsigned int Triangulation< dim, spacedim >::n_raw_lines ( ) const

inherited

Total number of lines, used or unused.

Note: This function really exports internal information about the triangulation. It shouldn't be used in applications. The function is only part of the public interface of this class because it is used in some of the other classes that build very closely on it (in particular, the DoFHandler class).

◆ n_raw_lines() [2/2]

template<int dim, int spacedim = dim>

unsigned int Triangulation< dim, spacedim >::n_raw_lines ( const unsigned int level ) const

inherited

Number of lines, used or unused, on the given level.

Note: This function really exports internal information about the triangulation. It shouldn't be used in applications. The function is only part of the public interface of this class because it is used in some of the other classes that build very closely on it (in particular, the DoFHandler class).

◆ n_raw_quads() [1/9]

template<int dim, int spacedim = dim>

unsigned int Triangulation< dim, spacedim >::n_raw_quads ( ) const

inherited

Total number of quads, used or unused.

Note: This function really exports internal information about the triangulation. It shouldn't be used in applications. The function is only part of the public interface of this class because it is used in some of the other classes that build very closely on it (in particular, the DoFHandler class).

◆ n_raw_quads() [2/9]

template<int dim, int spacedim = dim>

unsigned int Triangulation< dim, spacedim >::n_raw_quads ( const unsigned int level ) const

inherited

Number of quads, used or unused, on the given level.

Note: This function really exports internal information about the triangulation. It shouldn't be used in applications. The function is only part of the public interface of this class because it is used in some of the other classes that build very closely on it (in particular, the DoFHandler class).

◆ n_raw_quads() [3/9]

unsigned int Triangulation< 1, 1 >::n_raw_quads ( const unsigned int ) const

inherited

Definition at line 14132 of file tria.cc.

◆ n_raw_quads() [4/9]

unsigned int Triangulation< 1, 2 >::n_raw_quads ( const unsigned int ) const

inherited

Definition at line 14181 of file tria.cc.

◆ n_raw_quads() [5/9]

unsigned int Triangulation< 1, 3 >::n_raw_quads ( const unsigned int ) const

inherited

Definition at line 14229 of file tria.cc.

◆ n_raw_quads() [6/9]

unsigned int Triangulation< 2, 2 >::n_raw_quads ( const unsigned int level ) const

inherited

Definition at line 14282 of file tria.cc.

◆ n_raw_quads() [7/9]

unsigned int Triangulation< 2, 3 >::n_raw_quads ( const unsigned int level ) const

inherited

Definition at line 14292 of file tria.cc.

◆ n_raw_quads() [8/9]

unsigned int Triangulation< 3, 3 >::n_raw_quads ( const unsigned int ) const

inherited

Definition at line 14301 of file tria.cc.

◆ n_raw_quads() [9/9]

unsigned int Triangulation< 3, 3 >::n_raw_quads ( ) const

inherited

Definition at line 14321 of file tria.cc.

◆ n_raw_hexs() [1/5]

template<int dim, int spacedim = dim>

unsigned int Triangulation< dim, spacedim >::n_raw_hexs ( const unsigned int level ) const

inherited

Number of hexs, used or unused, on the given level.

Note: This function really exports internal information about the triangulation. It shouldn't be used in applications. The function is only part of the public interface of this class because it is used in some of the other classes that build very closely on it (in particular, the DoFHandler class).

◆ n_raw_hexs() [2/5]

unsigned int Triangulation< 1, 1 >::n_raw_hexs ( const unsigned int ) const

inherited

Definition at line 14140 of file tria.cc.

◆ n_raw_hexs() [3/5]

unsigned int Triangulation< 1, 2 >::n_raw_hexs ( const unsigned int ) const

inherited

Definition at line 14189 of file tria.cc.

◆ n_raw_hexs() [4/5]

unsigned int Triangulation< 1, 3 >::n_raw_hexs ( const unsigned int ) const

inherited

Definition at line 14237 of file tria.cc.

◆ n_raw_hexs() [5/5]

unsigned int Triangulation< 3, 3 >::n_raw_hexs ( const unsigned int level ) const

inherited

Definition at line 14414 of file tria.cc.

◆ n_raw_cells()

template<int dim, int spacedim = dim>

unsigned int Triangulation< dim, spacedim >::n_raw_cells ( const unsigned int level ) const

inherited

Number of cells, used or unused, on the given level.

Note: This function really exports internal information about the triangulation. It shouldn't be used in applications. The function is only part of the public interface of this class because it is used in some of the other classes that build very closely on it (in particular, the DoFHandler class).

◆ n_raw_faces()

template<int dim, int spacedim = dim>

unsigned int Triangulation< dim, spacedim >::n_raw_faces ( ) const

inherited

Return the total number of faces, used or not. In 2d, the result equals n_raw_lines(), in 3d it equals n_raw_quads(), while in 1d it equals the number of vertices.

Note: This function really exports internal information about the triangulation. It shouldn't be used in applications. The function is only part of the public interface of this class because it is used in some of the other classes that build very closely on it (in particular, the DoFHandler class).

◆ save()

template<int dim, int spacedim = dim>

template<class Archive >

void Triangulation< dim, spacedim >::save	(	Archive &	ar,
		const unsigned int	version
	)		const

inherited

Write the data of this object to a stream for the purpose of serialization using the BOOST serialization library.

Note: This function does not save all member variables of the current triangulation. Rather, only certain kinds of information are stored. For more information see the general documentation of this class.

◆ add_periodicity()

template<int dim, int spacedim = dim>

virtual void Triangulation< dim, spacedim >::add_periodicity ( const std::vector< GridTools::PeriodicFacePair< cell_iterator > > & )

virtualinherited

Declare the (coarse) face pairs given in the argument of this function as periodic. This way it is possible to obtain neighbors across periodic boundaries.

The vector can be filled by the function GridTools::collect_periodic_faces.

For more information on periodic boundary conditions see GridTools::collect_periodic_faces, DoFTools::make_periodicity_constraints and step-45.

Note: Before this function can be used the Triangulation has to be initialized and must not be refined.

◆ get_periodic_face_map()

template<int dim, int spacedim = dim>

const std::map< std::pair< cell_iterator, unsigned int >, std::pair< std::pair< cell_iterator, unsigned int >, std::bitset< 3 > > > & Triangulation< dim, spacedim >::get_periodic_face_map ( ) const

inherited

Return the periodic_face_map.

◆ get_reference_cells()

template<int dim, int spacedim = dim>

const std::vector< ReferenceCell > & Triangulation< dim, spacedim >::get_reference_cells ( ) const

inherited

Return vector filled with the used reference-cell types of this triangulation.

◆ all_reference_cells_are_hyper_cube()

template<int dim, int spacedim = dim>

bool Triangulation< dim, spacedim >::all_reference_cells_are_hyper_cube ( ) const

inherited

Indicate if the triangulation only consists of hypercube-like cells, i.e., lines, quadrilaterals, or hexahedra.

◆ all_reference_cells_are_simplex()

template<int dim, int spacedim = dim>

bool Triangulation< dim, spacedim >::all_reference_cells_are_simplex ( ) const

inherited

Indicate if the triangulation only consists of simplex-like cells, i.e., lines, triangles, or tetrahedra.

◆ is_mixed_mesh()

template<int dim, int spacedim = dim>

bool Triangulation< dim, spacedim >::is_mixed_mesh ( ) const

inherited

Indicate if the triangulation consists of different cell types (mix of simplices, hypercubes, ...) or different face types, as in the case of pyramids or wedges..

◆ serialize()

template<int dim, int spacedim = dim>

template<class Archive >

void Triangulation< dim, spacedim >::serialize	(	Archive &	archive,
		const unsigned int	version
	)

inherited

Write and read the data of this object from a stream for the purpose of serialization. using the BOOST serialization library.

◆ write_bool_vector()

template<int dim, int spacedim = dim>

static void Triangulation< dim, spacedim >::write_bool_vector	(	const unsigned int	magic_number1,
		const std::vector< bool > &	v,
		const unsigned int	magic_number2,
		std::ostream &	out
	)

staticprotectedinherited

Write a bool vector to the given stream, writing a pre- and a postfix magic number. The vector is written in an almost binary format, i.e. the bool flags are packed but the data is written as ASCII text.

The flags are stored in a binary format: for each true, a 1 bit is stored, a 0 bit otherwise. The bits are stored as unsigned char, thus avoiding endianness. They are written to out in plain text, thus amounting to 3.6 bits in the output per bits in the input on the average. Other information (magic numbers and number of elements of the input vector) is stored as plain text as well. The format should therefore be interplatform compatible.

◆ read_bool_vector()

template<int dim, int spacedim = dim>

static void Triangulation< dim, spacedim >::read_bool_vector	(	const unsigned int	magic_number1,
		std::vector< bool > &	v,
		const unsigned int	magic_number2,
		std::istream &	in
	)

staticprotectedinherited

Re-read a vector of bools previously written by write_bool_vector and compare with the magic numbers.

◆ update_periodic_face_map()

template<int dim, int spacedim = dim>

void Triangulation< dim, spacedim >::update_periodic_face_map ( )

protectedinherited

Recreate information about periodic neighbors from periodic_face_pairs_level_0.

◆ begin_raw()

template<int dim, int spacedim = dim>

raw_cell_iterator Triangulation< dim, spacedim >::begin_raw ( const unsigned int level = 0 ) const

privateinherited

Iterator to the first cell, used or not, on level level. If a level has no cells, a past-the-end iterator is returned.

◆ end_raw()

template<int dim, int spacedim = dim>

raw_cell_iterator Triangulation< dim, spacedim >::end_raw ( const unsigned int level ) const

privateinherited

Return a raw iterator which is the first iterator not on level. If level is the last level, then this returns end().

◆ begin_raw_line()

template<int dim, int spacedim = dim>

raw_line_iterator Triangulation< dim, spacedim >::begin_raw_line ( const unsigned int level = 0 ) const

privateinherited

Iterator to the first line, used or not, on level level. If a level has no lines, a past-the-end iterator is returned. If lines are no cells, i.e. for dim>1 no level argument must be given. The same applies for all the other functions above, of course.

◆ begin_line()

template<int dim, int spacedim = dim>

line_iterator Triangulation< dim, spacedim >::begin_line ( const unsigned int level = 0 ) const

privateinherited

Iterator to the first used line on level level.

Note: The given level argument needs to correspond to a level of the triangulation, i.e., should be less than the value returned by n_levels(). On the other hand, for parallel computations using a parallel::distributed::Triangulation object, it is often convenient to write loops over the cells of all levels of the global mesh, even if the local portion of the triangulation does not actually have cells at one of the higher levels. In those cases, the level argument is accepted if it is less than what the n_global_levels() function returns. If the given level is between the values returned by n_levels() and n_global_levels(), then no cells exist in the local portion of the triangulation at this level, and the function simply returns what end() would return.

◆ begin_active_line()

template<int dim, int spacedim = dim>

active_line_iterator Triangulation< dim, spacedim >::begin_active_line ( const unsigned int level = 0 ) const

privateinherited

Iterator to the first active line on level level.

Note: The given level argument needs to correspond to a level of the triangulation, i.e., should be less than the value returned by n_levels(). On the other hand, for parallel computations using a parallel::distributed::Triangulation object, it is often convenient to write loops over the cells of all levels of the global mesh, even if the local portion of the triangulation does not actually have cells at one of the higher levels. In those cases, the level argument is accepted if it is less than what the n_global_levels() function returns. If the given level is between the values returned by n_levels() and n_global_levels(), then no cells exist in the local portion of the triangulation at this level, and the function simply returns what end() would return.

◆ end_line()

template<int dim, int spacedim = dim>

line_iterator Triangulation< dim, spacedim >::end_line ( ) const

privateinherited

Iterator past the end; this iterator serves for comparisons of iterators with past-the-end or before-the-beginning states.

◆ begin_raw_quad()

template<int dim, int spacedim = dim>

raw_quad_iterator Triangulation< dim, spacedim >::begin_raw_quad ( const unsigned int level = 0 ) const

privateinherited

Iterator to the first quad, used or not, on the given level. If a level has no quads, a past-the-end iterator is returned. If quads are no cells, i.e. for \(dim>2\) no level argument must be given.

Note: The given level argument needs to correspond to a level of the triangulation, i.e., should be less than the value returned by n_levels(). On the other hand, for parallel computations using a parallel::distributed::Triangulation object, it is often convenient to write loops over the cells of all levels of the global mesh, even if the local portion of the triangulation does not actually have cells at one of the higher levels. In those cases, the level argument is accepted if it is less than what the n_global_levels() function returns. If the given level is between the values returned by n_levels() and n_global_levels(), then no cells exist in the local portion of the triangulation at this level, and the function simply returns what end() would return.

◆ begin_quad()

template<int dim, int spacedim = dim>

quad_iterator Triangulation< dim, spacedim >::begin_quad ( const unsigned int level = 0 ) const

privateinherited

Iterator to the first used quad on level level.

Note: The given level argument needs to correspond to a level of the triangulation, i.e., should be less than the value returned by n_levels(). On the other hand, for parallel computations using a parallel::distributed::Triangulation object, it is often convenient to write loops over the cells of all levels of the global mesh, even if the local portion of the triangulation does not actually have cells at one of the higher levels. In those cases, the level argument is accepted if it is less than what the n_global_levels() function returns. If the given level is between the values returned by n_levels() and n_global_levels(), then no cells exist in the local portion of the triangulation at this level, and the function simply returns what end() would return.

◆ begin_active_quad()

template<int dim, int spacedim = dim>

active_quad_iterator Triangulation< dim, spacedim >::begin_active_quad ( const unsigned int level = 0 ) const

privateinherited

Iterator to the first active quad on level level.

Note: The given level argument needs to correspond to a level of the triangulation, i.e., should be less than the value returned by n_levels(). On the other hand, for parallel computations using a parallel::distributed::Triangulation object, it is often convenient to write loops over the cells of all levels of the global mesh, even if the local portion of the triangulation does not actually have cells at one of the higher levels. In those cases, the level argument is accepted if it is less than what the n_global_levels() function returns. If the given level is between the values returned by n_levels() and n_global_levels(), then no cells exist in the local portion of the triangulation at this level, and the function simply returns what end() would return.

◆ end_quad()

template<int dim, int spacedim = dim>

quad_iterator Triangulation< dim, spacedim >::end_quad ( ) const

privateinherited

Iterator past the end; this iterator serves for comparisons of iterators with past-the-end or before-the-beginning states.

◆ begin_raw_hex()

template<int dim, int spacedim = dim>

raw_hex_iterator Triangulation< dim, spacedim >::begin_raw_hex ( const unsigned int level = 0 ) const

privateinherited

Iterator to the first hex, used or not, on level level. If a level has no hexes, a past-the-end iterator is returned.

Note: The given level argument needs to correspond to a level of the triangulation, i.e., should be less than the value returned by n_levels(). On the other hand, for parallel computations using a parallel::distributed::Triangulation object, it is often convenient to write loops over the cells of all levels of the global mesh, even if the local portion of the triangulation does not actually have cells at one of the higher levels. In those cases, the level argument is accepted if it is less than what the n_global_levels() function returns. If the given level is between the values returned by n_levels() and n_global_levels(), then no cells exist in the local portion of the triangulation at this level, and the function simply returns what end() would return.

◆ begin_hex()

template<int dim, int spacedim = dim>

hex_iterator Triangulation< dim, spacedim >::begin_hex ( const unsigned int level = 0 ) const

privateinherited

Iterator to the first used hex on level level.

Note: The given level argument needs to correspond to a level of the triangulation, i.e., should be less than the value returned by n_levels(). On the other hand, for parallel computations using a parallel::distributed::Triangulation object, it is often convenient to write loops over the cells of all levels of the global mesh, even if the local portion of the triangulation does not actually have cells at one of the higher levels. In those cases, the level argument is accepted if it is less than what the n_global_levels() function returns. If the given level is between the values returned by n_levels() and n_global_levels(), then no cells exist in the local portion of the triangulation at this level, and the function simply returns what end() would return.

◆ begin_active_hex()

template<int dim, int spacedim = dim>

active_hex_iterator Triangulation< dim, spacedim >::begin_active_hex ( const unsigned int level = 0 ) const

privateinherited

Iterator to the first active hex on level level.

Note: The given level argument needs to correspond to a level of the triangulation, i.e., should be less than the value returned by n_levels(). On the other hand, for parallel computations using a parallel::distributed::Triangulation object, it is often convenient to write loops over the cells of all levels of the global mesh, even if the local portion of the triangulation does not actually have cells at one of the higher levels. In those cases, the level argument is accepted if it is less than what the n_global_levels() function returns. If the given level is between the values returned by n_levels() and n_global_levels(), then no cells exist in the local portion of the triangulation at this level, and the function simply returns what end() would return.

◆ end_hex()

template<int dim, int spacedim = dim>

hex_iterator Triangulation< dim, spacedim >::end_hex ( ) const

privateinherited

Iterator past the end; this iterator serves for comparisons of iterators with past-the-end or before-the-beginning states.

◆ clear_despite_subscriptions()

template<int dim, int spacedim = dim>

void Triangulation< dim, spacedim >::clear_despite_subscriptions ( )

privateinherited

The (public) function clear() will only work when the triangulation is not subscribed to by other users. The clear_despite_subscriptions() function now allows the triangulation being cleared even when there are subscriptions.

Make sure, you know what you do, when calling this function, as its use is reasonable in very rare cases, only. For example, when the subscriptions were for the initially empty Triangulation and the Triangulation object wants to release its memory before throwing an assertion due to input errors (e.g. in the create_triangulation() function).

◆ reset_policy()

template<int dim, int spacedim = dim>

void Triangulation< dim, spacedim >::reset_policy ( )

privateinherited

Reset triangulation policy.

◆ reset_active_cell_indices()

template<int dim, int spacedim = dim>

void Triangulation< dim, spacedim >::reset_active_cell_indices ( )

privateinherited

For all cells, set the active cell indices so that active cells know the how many-th active cell they are, and all other cells have an invalid value. This function is called after mesh creation, refinement, and serialization.

◆ reset_cell_vertex_indices_cache()

template<int dim, int spacedim = dim>

void Triangulation< dim, spacedim >::reset_cell_vertex_indices_cache ( )

privateinherited

Reset cache for the cells' vertex indices.

◆ execute_refinement()

template<int dim, int spacedim = dim>

DistortedCellList Triangulation< dim, spacedim >::execute_refinement ( )

privateinherited

Refine all cells on all levels which were previously flagged for refinement.

Note, that this function uses the line->user_flags for dim=2,3 and the quad->user_flags for dim=3.

The function returns a list of cells that have produced children that satisfy the criteria of distorted cells if the check_for_distorted_cells flag was specified upon creation of this object, at

◆ execute_coarsening()

template<int dim, int spacedim = dim>

void Triangulation< dim, spacedim >::execute_coarsening ( )

privateinherited

Coarsen all cells which were flagged for coarsening, or rather: delete all children of those cells of which all child cells are flagged for coarsening and several other constraints hold (see the general doc of this class).

◆ fix_coarsen_flags()

template<int dim, int spacedim = dim>

void Triangulation< dim, spacedim >::fix_coarsen_flags ( )

privateinherited

Make sure that either all or none of the children of a cell are tagged for coarsening.

◆ coarse_cell_id_to_coarse_cell_index()

template<int dim, int spacedim = dim>

virtual unsigned int Triangulation< dim, spacedim >::coarse_cell_id_to_coarse_cell_index ( const types::coarse_cell_id coarse_cell_id ) const

privatevirtualinherited

Translate the unique id of a coarse cell to its index. See the glossary entry on coarse cell IDs for more information.

Note: For serial and shared triangulation both id and index are the same. For distributed triangulations setting both might differ, since the id might correspond to a global id and the index to a local id.

Parameters

coarse_cell_id Unique id of the coarse cell.

Returns: Index of the coarse cell within the current triangulation.

Reimplemented in parallel::fullydistributed::Triangulation< dim, spacedim >, parallel::distributed::Triangulation< dim, spacedim >, and parallel::distributed::Triangulation< 1, spacedim >.

◆ coarse_cell_index_to_coarse_cell_id()

template<int dim, int spacedim = dim>

virtual types::coarse_cell_id Triangulation< dim, spacedim >::coarse_cell_index_to_coarse_cell_id ( const unsigned int coarse_cell_index ) const

privatevirtualinherited

Translate the index of coarse cell to its unique id. See the glossary entry on coarse cell IDs for more information.

Note: See the note of the method coarse_cell_id_to_coarse_cell_index().

Parameters

coarse_cell_index Index of the coarse cell.

Returns: Id of the coarse cell.

Reimplemented in parallel::fullydistributed::Triangulation< dim, spacedim >, parallel::distributed::Triangulation< dim, spacedim >, and parallel::distributed::Triangulation< 1, spacedim >.

◆ subscribe()

void Subscriptor::subscribe	(	std::atomic< bool > *const	validity,
		const std::string &	identifier = `""`
	)		const

inherited

Subscribes a user of the object by storing the pointer validity. The subscriber may be identified by text supplied as identifier.

Definition at line 136 of file subscriptor.cc.

◆ unsubscribe()

void Subscriptor::unsubscribe	(	std::atomic< bool > *const	validity,
		const std::string &	identifier = `""`
	)		const

inherited

Unsubscribes a user from the object.

Note: The identifier and the validity pointer must be the same as the one supplied to subscribe().

Definition at line 156 of file subscriptor.cc.

◆ n_subscriptions()

unsigned int Subscriptor::n_subscriptions ( ) const

inlineinherited

Return the present number of subscriptions to this object. This allows to use this class for reference counted lifetime determination where the last one to unsubscribe also deletes the object.

Definition at line 300 of file subscriptor.h.

◆ list_subscribers() [1/2]

template<typename StreamType >

void Subscriptor::list_subscribers ( StreamType & stream ) const

inlineinherited

List the subscribers to the input stream.

Definition at line 317 of file subscriptor.h.

◆ list_subscribers() [2/2]

void Subscriptor::list_subscribers ( ) const

inherited

List the subscribers to deallog.

Definition at line 204 of file subscriptor.cc.

◆ check_no_subscribers()

void Subscriptor::check_no_subscribers ( ) const

privatenoexceptinherited

Check that there are no objects subscribing to this object. If this check passes then it is safe to destroy the current object. It this check fails then this function will either abort or print an error message to deallog (by using the AssertNothrow mechanism), but will not throw an exception.

Note: Since this function is just a consistency check it does nothing in release mode.; If this function is called when there is an uncaught exception then, rather than aborting, this function prints an error message to the standard error stream and returns.

Definition at line 53 of file subscriptor.cc.

Member Data Documentation

◆ settings

template<int dim, int spacedim = dim>

const Settings parallel::shared::Triangulation< dim, spacedim >::settings

private

Settings

Definition at line 368 of file shared_tria.h.

◆ allow_artificial_cells

template<int dim, int spacedim = dim>

const bool parallel::shared::Triangulation< dim, spacedim >::allow_artificial_cells

private

A flag to decide whether or not artificial cells are allowed.

Definition at line 373 of file shared_tria.h.

◆ true_subdomain_ids_of_cells

template<int dim, int spacedim = dim>

std::vector<types::subdomain_id> parallel::shared::Triangulation< dim, spacedim >::true_subdomain_ids_of_cells

private

A vector containing subdomain IDs of cells obtained by partitioning using either zorder, METIS, or a user-defined partitioning scheme. In case allow_artificial_cells is false, this vector is consistent with IDs stored in cell->subdomain_id() of the triangulation class. When allow_artificial_cells is true, cells which are artificial will have cell->subdomain_id() == numbers::artificial;

The original partition information is stored to allow using sequential DoF distribution and partitioning functions with semi-artificial cells.

Definition at line 394 of file shared_tria.h.

◆ true_level_subdomain_ids_of_cells

template<int dim, int spacedim = dim>

std::vector<std::vector<types::subdomain_id> > parallel::shared::Triangulation< dim, spacedim >::true_level_subdomain_ids_of_cells

private

A vector containing level subdomain IDs of cells obtained by partitioning each level.

The original partition information is stored to allow using sequential DoF distribution and partitioning functions with semi-artificial cells.

Definition at line 405 of file shared_tria.h.

◆ mpi_communicator

template<int dim, int spacedim = dim>

const MPI_Comm parallel::TriangulationBase< dim, spacedim >::mpi_communicator

protectedinherited

MPI communicator to be used for the triangulation. We create a unique communicator for this class, which is a duplicate of the one passed to the constructor.

Definition at line 296 of file tria_base.h.

◆ my_subdomain

template<int dim, int spacedim = dim>

types::subdomain_id parallel::TriangulationBase< dim, spacedim >::my_subdomain

protectedinherited

The subdomain id to be used for the current processor. This is the MPI rank.

Definition at line 302 of file tria_base.h.

◆ n_subdomains

template<int dim, int spacedim = dim>

types::subdomain_id parallel::TriangulationBase< dim, spacedim >::n_subdomains

protectedinherited

The total number of subdomains (or the size of the MPI communicator).

Definition at line 307 of file tria_base.h.

◆ number_cache

template<int dim, int spacedim = dim>

NumberCache parallel::TriangulationBase< dim, spacedim >::number_cache

protectedinherited

Definition at line 365 of file tria_base.h.

◆ dimension

template<int dim, int spacedim = dim>

constexpr unsigned int Triangulation< dim, spacedim >::dimension = dim

staticconstexprinherited

Make the dimension available in function templates.

Definition at line 1557 of file tria.h.

◆ space_dimension

template<int dim, int spacedim = dim>

constexpr unsigned int Triangulation< dim, spacedim >::space_dimension = spacedim

staticconstexprinherited

Make the space-dimension available in function templates.

Definition at line 1562 of file tria.h.

◆ signals

template<int dim, int spacedim = dim>

Signals Triangulation< dim, spacedim >::signals

mutableinherited

Signals for the various actions that a triangulation can do to itself.

Definition at line 2479 of file tria.h.

◆ smooth_grid

template<int dim, int spacedim = dim>

MeshSmoothing Triangulation< dim, spacedim >::smooth_grid

protectedinherited

Do some smoothing in the process of refining the triangulation. See the general doc of this class for more information about this.

Definition at line 3731 of file tria.h.

◆ reference_cells

template<int dim, int spacedim = dim>

std::vector<ReferenceCell> Triangulation< dim, spacedim >::reference_cells

protectedinherited

Vector caching all reference-cell types of the given triangulation (also in the distributed case).

Definition at line 3737 of file tria.h.

◆ policy

template<int dim, int spacedim = dim>

std::unique_ptr< ::internal::TriangulationImplementation::Policy<dim, spacedim> > Triangulation< dim, spacedim >::policy

privateinherited

Policy with the Triangulation-specific tasks related to creation, refinement, and coarsening.

Definition at line 3789 of file tria.h.

◆ periodic_face_pairs_level_0

template<int dim, int spacedim = dim>

std::vector<GridTools::PeriodicFacePair<cell_iterator> > Triangulation< dim, spacedim >::periodic_face_pairs_level_0

privateinherited

If add_periodicity() is called, this variable stores the given periodic face pairs on level 0 for later access during the identification of ghost cells for the multigrid hierarchy and for setting up the periodic_face_map.

Definition at line 3798 of file tria.h.

◆ periodic_face_map

template<int dim, int spacedim = dim>

std::map<std::pair<cell_iterator, unsigned int>, std::pair<std::pair<cell_iterator, unsigned int>, std::bitset<3> > > Triangulation< dim, spacedim >::periodic_face_map

privateinherited

If add_periodicity() is called, this variable stores the active periodic face pairs.

Definition at line 3806 of file tria.h.

◆ levels

template<int dim, int spacedim = dim>

std::vector< std::unique_ptr<::internal::TriangulationImplementation::TriaLevel> > Triangulation< dim, spacedim >::levels

privateinherited

Array of pointers pointing to the objects storing the cell data on the different levels.

Definition at line 4184 of file tria.h.

◆ faces

template<int dim, int spacedim = dim>

std::unique_ptr<::internal::TriangulationImplementation::TriaFaces> Triangulation< dim, spacedim >::faces

privateinherited

Pointer to the faces of the triangulation. In 1d this contains nothing, in 2d it contains data concerning lines and in 3d quads and lines. All of these have no level and are therefore treated separately.

Definition at line 4192 of file tria.h.

◆ vertices

template<int dim, int spacedim = dim>

std::vector<Point<spacedim> > Triangulation< dim, spacedim >::vertices

privateinherited

Array of the vertices of this triangulation.

Definition at line 4198 of file tria.h.

◆ vertices_used

template<int dim, int spacedim = dim>

std::vector<bool> Triangulation< dim, spacedim >::vertices_used

privateinherited

Array storing a bit-pattern which vertices are used.

Definition at line 4203 of file tria.h.

◆ manifolds

template<int dim, int spacedim = dim>

std::map<types::manifold_id, std::unique_ptr<const Manifold<dim, spacedim> > > Triangulation< dim, spacedim >::manifolds

privateinherited

Collection of manifold objects. We store only objects, which are not of type FlatManifold.

Definition at line 4210 of file tria.h.

◆ anisotropic_refinement

template<int dim, int spacedim = dim>

bool Triangulation< dim, spacedim >::anisotropic_refinement

privateinherited

Flag indicating whether anisotropic refinement took place.

Definition at line 4215 of file tria.h.

◆ check_for_distorted_cells

template<int dim, int spacedim = dim>

const bool Triangulation< dim, spacedim >::check_for_distorted_cells

privateinherited

A flag that determines whether we are to check for distorted cells upon creation and refinement of a mesh.

Definition at line 4222 of file tria.h.

◆ vertex_to_boundary_id_map_1d

template<int dim, int spacedim = dim>

std::unique_ptr<std::map<unsigned int, types::boundary_id> > Triangulation< dim, spacedim >::vertex_to_boundary_id_map_1d

privateinherited

A map that relates the number of a boundary vertex to the boundary indicator. This field is only used in 1d. We have this field because we store boundary indicator information with faces in 2d and higher where we have space in the structures that store data for faces, but in 1d there is no such space for faces.

The field is declared as a pointer for a rather mundane reason: all other fields of this class that can be modified by the TriaAccessor hierarchy are pointers, and so these accessor classes store a const pointer to the triangulation. We could no longer do so for TriaAccessor<0,1,spacedim> if this field (that can be modified by TriaAccessor::set_boundary_id) were not a pointer.

Definition at line 4250 of file tria.h.

◆ vertex_to_manifold_id_map_1d

template<int dim, int spacedim = dim>

std::unique_ptr<std::map<unsigned int, types::manifold_id> > Triangulation< dim, spacedim >::vertex_to_manifold_id_map_1d

privateinherited

A map that relates the number of a boundary vertex to the manifold indicator. This field is only used in 1d. We have this field because we store manifold indicator information with faces in 2d and higher where we have space in the structures that store data for faces, but in 1d there is no such space for faces.

Note: Manifold objects are pretty useless for points since they are neither refined nor are their interiors mapped. We nevertheless allow storing manifold ids for points to be consistent in dimension-independent programs.

The field is declared as a pointer for a rather mundane reason: all other fields of this class that can be modified by the TriaAccessor hierarchy are pointers, and so these accessor classes store a const pointer to the triangulation. We could no longer do so for TriaAccessor<0,1,spacedim> if this field (that can be modified by TriaAccessor::set_manifold_id) were not a pointer.

Definition at line 4273 of file tria.h.

◆ counter

std::atomic<unsigned int> Subscriptor::counter

mutableprivateinherited

Store the number of objects which subscribed to this object. Initially, this number is zero, and upon destruction it shall be zero again (i.e. all objects which subscribed should have unsubscribed again).

The creator (and owner) of an object is counted in the map below if HE manages to supply identification.

We use the mutable keyword in order to allow subscription to constant objects also.

This counter may be read from and written to concurrently in multithreaded code: hence we use the std::atomic class template.

Definition at line 219 of file subscriptor.h.

◆ counter_map

std::map<std::string, unsigned int> Subscriptor::counter_map

mutableprivateinherited

In this map, we count subscriptions for each different identification string supplied to subscribe().

Definition at line 225 of file subscriptor.h.

◆ validity_pointers

std::vector<std::atomic<bool> *> Subscriptor::validity_pointers

mutableprivateinherited

In this vector, we store pointers to the validity bool in the SmartPointer objects that subscribe to this class.

Definition at line 241 of file subscriptor.h.

◆ object_info

const std::type_info* Subscriptor::object_info

mutableprivateinherited

Pointer to the typeinfo object of this object, from which we can later deduce the class name. Since this information on the derived class is neither available in the destructor, nor in the constructor, we obtain it in between and store it here.

Definition at line 249 of file subscriptor.h.

◆ mutex

std::mutex Subscriptor::mutex

staticprivateinherited

A mutex used to ensure data consistency when printing out the list of subscribers.

Definition at line 271 of file subscriptor.h.

The documentation for this class was generated from the following files:

include/deal.II/base/template_constraints.h
include/deal.II/distributed/shared_tria.h
source/distributed/shared_tria.cc

Public Types

Public Member Functions

Static Public Member Functions

Static Public Attributes

Protected Member Functions

Protected Attributes

Private Types

Private Member Functions

Private Attributes

Static Private Attributes

Keeping up with what happens to a triangulation

Exceptions

Hex iterator functions for internal use

Cell iterator functions for internal use

Detailed Description

Member Typedef Documentation

◆ active_cell_iterator

◆ cell_iterator

◆ IteratorSelector

◆ level_cell_iterator

◆ raw_cell_iterator

◆ raw_face_iterator

◆ raw_vertex_iterator

◆ raw_line_iterator

◆ raw_quad_iterator

◆ raw_hex_iterator

◆ map_value_type

◆ map_iterator

Member Enumeration Documentation

◆ Settings

◆ MeshSmoothing

◆ CellStatus

Constructor & Destructor Documentation

◆ Triangulation()

◆ ~Triangulation()

Member Function Documentation

◆ is_multilevel_hierarchy_constructed()

◆ execute_coarsening_and_refinement()

◆ create_triangulation() [1/2]

◆ create_triangulation() [2/2]

◆ copy_triangulation() [1/2]

◆ load()

◆ get_true_subdomain_ids_of_cells()

◆ get_true_level_subdomain_ids_of_cells()

◆ with_artificial_cells()

◆ partition()

◆ get_communicator()

◆ copy_triangulation() [2/2]

◆ n_locally_owned_active_cells()

◆ n_global_active_cells()

◆ memory_consumption()

◆ n_global_levels()

◆ locally_owned_subdomain()

◆ ghost_owners()

◆ level_ghost_owners()

◆ global_active_cell_index_partitioner()

◆ global_level_cell_index_partitioner()

◆ get_boundary_ids()

◆ get_manifold_ids()

◆ communicate_locally_moved_vertices()

◆ n_global_coarse_cells()

◆ update_number_cache()

◆ update_reference_cells()

◆ reset_global_cell_indices()

◆ clear()

◆ set_mesh_smoothing()

◆ get_mesh_smoothing()

◆ create_triangulation_compatibility()

◆ flip_all_direction_flags()

◆ set_all_refine_flags()

◆ refine_global()

◆ coarsen_global()

◆ prepare_coarsening_and_refinement() [1/4]

◆ prepare_coarsening_and_refinement() [2/4]

◆ prepare_coarsening_and_refinement() [3/4]

◆ prepare_coarsening_and_refinement() [4/4]

◆ save_refine_flags() [1/2]

◆ save_refine_flags() [2/2]

◆ load_refine_flags() [1/2]

◆ load_refine_flags() [2/2]