DOLFINx 0.10.0.0
DOLFINx C++ interface
Loading...
Searching...
No Matches
Topology.h
1// Copyright (C) 2006-2024 Anders Logg and Garth N. Wells
2//
3// This file is part of DOLFINx (https://www.fenicsproject.org)
4//
5// SPDX-License-Identifier: LGPL-3.0-or-later
6
7#pragma once
8
9#include <array>
10#include <cstdint>
11#include <dolfinx/common/MPI.h>
12#include <map>
13#include <memory>
14#include <optional>
15#include <span>
16#include <tuple>
17#include <utility>
18#include <vector>
19
20namespace dolfinx::common
21{
22class IndexMap;
23}
24
25namespace dolfinx::graph
26{
27template <typename T>
28class AdjacencyList;
29}
30
31namespace dolfinx::mesh
32{
33enum class CellType;
34
46{
47public:
66 std::vector<CellType> cell_types,
67 std::shared_ptr<const common::IndexMap> vertex_map,
68 std::vector<std::shared_ptr<const common::IndexMap>> cell_maps,
69 std::vector<std::shared_ptr<graph::AdjacencyList<std::int32_t>>> cells,
70 const std::optional<std::vector<std::vector<std::int64_t>>>&
72 = std::nullopt);
73
75 Topology(const Topology& topology) = default;
76
78 Topology(Topology&& topology) = default;
79
81 ~Topology() = default;
82
84 Topology& operator=(const Topology& topology) = delete;
85
87 Topology& operator=(Topology&& topology) = default;
88
90 int dim() const noexcept;
91
95 const std::vector<CellType>& entity_types(int dim) const;
96
102 CellType cell_type() const;
103
109 std::vector<std::shared_ptr<const common::IndexMap>>
110 index_maps(int dim) const;
111
118 std::shared_ptr<const common::IndexMap> index_map(int dim) const;
119
133 std::shared_ptr<const graph::AdjacencyList<std::int32_t>>
134 connectivity(std::array<int, 2> d0, std::array<int, 2> d1) const;
135
145 std::shared_ptr<const graph::AdjacencyList<std::int32_t>>
146 connectivity(int d0, int d1) const;
147
149 const std::vector<std::uint32_t>& get_cell_permutation_info() const;
150
165 const std::vector<std::uint8_t>& get_facet_permutations() const;
166
169 std::vector<CellType> cell_types() const;
170
183 const std::vector<std::int32_t>& interprocess_facets(int index) const;
184
193 const std::vector<std::int32_t>& interprocess_facets() const;
194
199 bool create_entities(int dim);
200
205 void create_connectivity(int d0, int d1);
206
209
211 std::vector<std::vector<std::int64_t>> original_cell_index;
212
215 MPI_Comm comm() const;
216
217private:
218 // Cell types for entities in Topology, where _entity_types_new[d][i]
219 // is the ith entity type of dimension d
220 std::vector<std::vector<CellType>> _entity_types;
221
222 // Parallel layout of entities for each dimension and cell type
223 // flattened in the same layout as _entity_types above.
224 // std::vector<std::shared_ptr<const common::IndexMap>> _index_map;
225
226 // _index_maps[(d, i) is the index map for the ith entity type of
227 // dimension d
228 std::map<std::array<int, 2>, std::shared_ptr<const common::IndexMap>>
229 _index_maps;
230
231 // Connectivity between cell types _connectivity_new[(dim0, i0),
232 // (dim1, i1)] is the connection from (dim0, i0) -> (dim1, i1),
233 // where dim0 and dim1 are topological dimensions and i0 and i1
234 // are the indices of cell types (following the order in _entity_types).
235 std::map<std::pair<std::array<int, 2>, std::array<int, 2>>,
236 std::shared_ptr<graph::AdjacencyList<std::int32_t>>>
237 _connectivity;
238
239 // The facet permutations (local facet, cell))
240 // [cell0_0, cell0_1, ,cell0_2, cell1_0, cell1_1, ,cell1_2, ...,
241 // celln_0, celln_1, ,celln_2,]
242 std::vector<std::uint8_t> _facet_permutations;
243
244 // Cell permutation info. See the documentation for
245 // get_cell_permutation_info for documentation of how this is encoded.
246 std::vector<std::uint32_t> _cell_permutations;
247
248 // List of facets that are on the inter-process boundary for each
249 // facet type. _interprocess_facets[i] is the inter-process facets of
250 // facet type i.
251 std::vector<std::vector<std::int32_t>> _interprocess_facets;
252};
253
282create_topology(MPI_Comm comm, const std::vector<CellType>& cell_types,
283 std::vector<std::span<const std::int64_t>> cells,
284 std::vector<std::span<const std::int64_t>> original_cell_index,
285 std::vector<std::span<const int>> ghost_owners,
286 std::span<const std::int64_t> boundary_vertices);
287
311Topology create_topology(MPI_Comm comm, std::span<const std::int64_t> cells,
312 std::span<const std::int64_t> original_cell_index,
313 std::span<const int> ghost_owners, CellType cell_type,
314 std::span<const std::int64_t> boundary_vertices);
315
328std::tuple<Topology, std::vector<int32_t>, std::vector<int32_t>>
329create_subtopology(const Topology& topology, int dim,
330 std::span<const std::int32_t> entities);
331
341std::vector<std::int32_t>
342entities_to_index(const Topology& topology, int dim,
343 std::span<const std::int32_t> entities);
344
362std::vector<std::vector<std::int32_t>>
363compute_mixed_cell_pairs(const Topology& topology, mesh::CellType facet_type);
364
365} // namespace dolfinx::mesh
Definition IndexMap.h:94
Definition AdjacencyList.h:27
std::shared_ptr< const common::IndexMap > index_map(int dim) const
Get the IndexMap that described the parallel distribution of the mesh entities.
Definition Topology.cpp:823
void create_connectivity(int d0, int d1)
Create connectivity between given pair of dimensions, d0 / -> d1.
Definition Topology.cpp:950
std::vector< std::vector< std::int64_t > > original_cell_index
Original cell index for each cell type.
Definition Topology.h:211
void create_entity_permutations()
Compute entity permutations and reflections.
Definition Topology.cpp:990
const std::vector< std::uint32_t > & get_cell_permutation_info() const
Returns the permutation information.
Definition Topology.cpp:852
std::shared_ptr< const graph::AdjacencyList< std::int32_t > > connectivity(std::array< int, 2 > d0, std::array< int, 2 > d1) const
Get the connectivity from entities of topological dimension d0 to dimension d1.
Definition Topology.cpp:832
const std::vector< std::uint8_t > & get_facet_permutations() const
Get the numbers that encode the number of permutations to apply to facets.
Definition Topology.cpp:867
Topology(const Topology &topology)=default
Copy constructor.
~Topology()=default
Destructor.
Topology(std::vector< CellType > cell_types, std::shared_ptr< const common::IndexMap > vertex_map, std::vector< std::shared_ptr< const common::IndexMap > > cell_maps, std::vector< std::shared_ptr< graph::AdjacencyList< std::int32_t > > > cells, const std::optional< std::vector< std::vector< std::int64_t > > > &original_cell_index=std::nullopt)
Create a mesh topology.
Definition Topology.cpp:741
const std::vector< CellType > & entity_types(int dim) const
Entity types in the topology for a given dimension.
Definition Topology.cpp:790
int dim() const noexcept
Topological dimension of the mesh.
Definition Topology.cpp:785
const std::vector< std::int32_t > & interprocess_facets(int index) const
List of inter-process facets of a given type.
Definition Topology.cpp:881
Topology(Topology &&topology)=default
Move constructor.
std::vector< CellType > cell_types() const
Get the types of cells in the topology.
Definition Topology.cpp:805
Topology & operator=(const Topology &topology)=delete
Assignment.
MPI_Comm comm() const
Mesh MPI communicator.
Definition Topology.cpp:1010
CellType cell_type() const
Cell type.
Definition Topology.cpp:795
bool create_entities(int dim)
Create entities of given topological dimension.
Definition Topology.cpp:893
Topology & operator=(Topology &&topology)=default
Assignment.
std::vector< std::shared_ptr< const common::IndexMap > > index_maps(int dim) const
Get the index maps that described the parallel distribution of the mesh entities of a given topologic...
Definition Topology.cpp:811
Miscellaneous classes, functions and types.
Definition dolfinx_common.h:8
Graph data structures and algorithms.
Definition dofmapbuilder.h:26
Mesh data structures and algorithms on meshes.
Definition DofMap.h:32
Topology create_topology(MPI_Comm comm, const std::vector< CellType > &cell_types, std::vector< std::span< const std::int64_t > > cells, std::vector< std::span< const std::int64_t > > original_cell_index, std::vector< std::span< const int > > ghost_owners, std::span< const std::int64_t > boundary_vertices)
Create a mesh topology.
Definition Topology.cpp:1017
std::tuple< Topology, std::vector< int32_t >, std::vector< int32_t > > create_subtopology(const Topology &topology, int dim, std::span< const std::int32_t > entities)
Create a topology for a subset of entities of a given topological dimension.
Definition Topology.cpp:1303
CellType
Cell type identifier.
Definition cell_types.h:22
std::vector< std::vector< std::int32_t > > compute_mixed_cell_pairs(const Topology &topology, mesh::CellType facet_type)
Compute a list of cell-cell connections for each possible combination in the topology which have the ...
Definition Topology.cpp:1444
std::vector< std::int32_t > entities_to_index(const Topology &topology, int dim, std::span< const std::int32_t > entities)
Get entity indices for entities defined by their vertices.
Definition Topology.cpp:1386