Flow123d
JS_before_hm-1755-g6249b9fc1
|
#include <mesh.h>
Classes | |
struct | ElementTmpData |
Public Types | |
enum | IntersectionSearch { BIHsearch = 1, BIHonly = 2, BBsearch = 3 } |
Types of search algorithm for finding intersection candidates. More... | |
enum | { x_coord =0, y_coord =1, z_coord =2 } |
typedef enum Mesh::IntersectionSearch | IntersectionSearch |
Types of search algorithm for finding intersection candidates. More... | |
Public Member Functions | |
TYPEDEF_ERR_INFO (EI_ElemLast, int) | |
TYPEDEF_ERR_INFO (EI_ElemNew, int) | |
TYPEDEF_ERR_INFO (EI_RegLast, std::string) | |
TYPEDEF_ERR_INFO (EI_RegNew, std::string) | |
TYPEDEF_ERR_INFO (EI_ElemId, int) | |
TYPEDEF_ERR_INFO (EI_ElemIdOther, int) | |
TYPEDEF_ERR_INFO (EI_Region, std::string) | |
TYPEDEF_ERR_INFO (EI_RegIdx, unsigned int) | |
TYPEDEF_ERR_INFO (EI_Dim, unsigned int) | |
TYPEDEF_ERR_INFO (EI_DimOther, unsigned int) | |
DECLARE_EXCEPTION (ExcDuplicateBoundary,<< "Duplicate boundary elements! \n"<< "Element id: "<< EI_ElemLast::val<< " on region name: "<< EI_RegLast::val<< "\n"<< "Element id: "<< EI_ElemNew::val<< " on region name: "<< EI_RegNew::val<< "\n") | |
DECLARE_EXCEPTION (ExcElmWrongOrder,<< "Element IDs in non-increasing order, ID: "<< EI_ElemId::val<< "\n") | |
DECLARE_EXCEPTION (ExcRegionElmDiffDim,<< "User defined region "<< EI_Region::qval<< " (id "<< EI_RegIdx::val<< ") by 'From_Elements' cannot have elements of different dimensions.\n"<< "Thrown due to: dim "<< EI_Dim::val<< " neq dim "<< EI_DimOther::val<< " (ele id "<< EI_ElemId::val<< ").\n"<< "Split elements by dim, create separate regions and then possibly use Union.\n") | |
DECLARE_EXCEPTION (ExcTooMatchingIds,<< "Mesh: Duplicate dim-join lower dim elements: "<< EI_ElemId::val<< ", "<< EI_ElemIdOther::val<< ".\n") | |
DECLARE_EXCEPTION (ExcBdrElemMatchRegular,<< "Boundary element (id: "<< EI_ElemId::val<< ") match a regular element (id: "<< EI_ElemIdOther::val<< ") of lower dimension.\n") | |
Mesh () | |
Mesh (Input::Record in_record, MPI_Comm com=MPI_COMM_WORLD) | |
virtual | ~Mesh () |
Destructor. More... | |
virtual unsigned int | n_nodes () const |
unsigned int | n_boundaries () const |
unsigned int | n_edges () const |
Edge | edge (uint edge_idx) const |
Boundary | boundary (uint edge_idx) const |
unsigned int | n_corners () |
const RegionDB & | region_db () const |
virtual Partitioning * | get_part () |
virtual const LongIdx * | get_local_part () |
Distribution * | get_el_ds () const |
LongIdx * | get_row_4_el () const |
LongIdx * | get_el_4_loc () const |
Distribution * | get_node_ds () const |
LongIdx * | get_node_4_loc () const |
unsigned int | n_local_nodes () const |
MPI_Comm | get_comm () const |
MixedMeshIntersections & | mixed_intersections () |
unsigned int | n_sides () const |
unsigned int | n_vb_neighbours () const |
unsigned int | max_edge_sides (unsigned int dim) const |
void | read_gmsh_from_stream (istream &in) |
void | init_from_input () |
void | setup_topology () |
void | elements_id_maps (vector< LongIdx > &bulk_elements_id, vector< LongIdx > &boundary_elements_id) const |
virtual std::shared_ptr< std::vector< LongIdx > > | check_compatible_mesh (Mesh &input_mesh) |
virtual std::shared_ptr< std::vector< LongIdx > > | check_compatible_discont_mesh (Mesh &input_mesh) |
virtual ElementAccessor< 3 > | element_accessor (unsigned int idx) const |
Create and return ElementAccessor to element of given idx. More... | |
NodeAccessor< 3 > | node (unsigned int idx) const |
Create and return NodeAccessor to node of given idx. More... | |
void | read_regions_from_input (Input::Array region_list) |
const vector< vector< unsigned int > > & | node_elements () |
void | check_and_finish () |
std::vector< BoundingBox > | get_element_boxes () |
Compute bounding boxes of elements contained in mesh. More... | |
const BIHTree & | get_bih_tree () |
Getter for BIH. Creates and compute BIH at first call. More... | |
void | intersect_element_lists (vector< unsigned int > const &nodes_list, vector< unsigned int > &intersection_element_list) |
void | add_node (unsigned int node_id, arma::vec3 coords) |
Add new node of given id and coordinates to mesh. More... | |
void | add_element (unsigned int elm_id, unsigned int dim, unsigned int region_id, unsigned int partition_id, std::vector< unsigned int > node_ids) |
Add new element of given id to mesh. More... | |
void | add_physical_name (unsigned int dim, unsigned int id, std::string name) |
Add new node of given id and coordinates to mesh. More... | |
FilePath | mesh_file () |
Return FilePath object representing "mesh_file" input key. More... | |
IntersectionSearch | get_intersection_search () |
Getter for input type selection for intersection search algorithm. More... | |
double | global_snap_radius () const |
Maximal distance of observe point from Mesh relative to its size. More... | |
void | init_element_vector (unsigned int size) |
Initialize element_vec_, set size and reset counters of boundary and bulk elements. More... | |
void | init_node_vector (unsigned int size) |
Initialize node_vec_, set size. More... | |
virtual Range< ElementAccessor< 3 > > | elements_range () const |
Returns range of bulk elements. More... | |
Range< NodeAccessor< 3 > > | node_range () const |
Returns range of nodes. More... | |
Range< Edge > | edge_range () const |
Returns range of edges. More... | |
virtual unsigned int | n_elements () const |
Returns count of boundary or bulk elements. More... | |
int | elem_index (int elem_id) const |
For element of given elem_id returns index in element_vec_ or (-1) if element doesn't exist. More... | |
int | find_elem_id (unsigned int pos) const |
Return element id (in GMSH file) of element of given position in element vector. More... | |
const std::vector< unsigned int > & | element_permutations () const |
Return permutation vector of elements. More... | |
int | node_index (int node_id) const |
For node of given node_id returns index in element_vec_ or (-1) if node doesn't exist. More... | |
int | find_node_id (unsigned int pos) const |
Return node id (in GMSH file) of node of given position in node vector. More... | |
const std::vector< unsigned int > & | node_permutations () const |
Return permutation vector of nodes. More... | |
void | check_element_size (unsigned int elem_idx) const |
Check if given index is in element_vec_. More... | |
unsigned int | create_boundary_elements () |
Create boundary elements from data of temporary structure, this method MUST be call after read mesh from file, return number of read boundary elements. More... | |
void | permute_tetrahedron (unsigned int elm_idx, std::vector< unsigned int > permutation_vec) |
Permute nodes of 3D elements of given elm_idx. More... | |
void | permute_triangle (unsigned int elm_idx, std::vector< unsigned int > permutation_vec) |
Permute nodes of 2D elements of given elm_idx. More... | |
BCMesh * | get_bc_mesh () |
Create boundary mesh if doesn't exist and return it. More... | |
Static Public Member Functions | |
static const Input::Type::Selection & | get_input_intersection_variant () |
The definition of input record for selection of variant of file format. More... | |
static const Input::Type::Record & | get_input_type () |
Public Attributes | |
vector< BoundaryData > | boundary_ |
std::shared_ptr< MixedMeshIntersections > | intersections |
vector< vector< unsigned int > > | master_elements |
DuplicateNodes * | tree |
vector< Neighbour > | vb_neighbours_ |
int | n_insides |
int | n_exsides |
int | n_sides_ |
int | n_lines |
int | n_triangles |
int | n_tetrahedras |
vector< vector< vector< unsigned int > > > | side_nodes |
vector< vector< unsigned int > > | node_elements_ |
For each node the vector contains a list of elements that use this node. More... | |
Static Public Attributes | |
static const unsigned int | undef_idx =-1 |
Protected Member Functions | |
void | init () |
void | make_neighbours_and_edges () |
void | make_edge_permutations () |
void | create_node_element_lists () |
bool | find_lower_dim_element (vector< unsigned int > &element_list, unsigned int dim, unsigned int &element_idx) |
bool | same_sides (const SideIter &si, vector< unsigned int > &side_nodes) |
void | element_to_neigh_vb () |
void | count_element_types () |
void | count_side_types () |
void | check_mesh_on_read () |
void | modify_element_ids (const RegionDB::MapElementIDToRegionID &map) |
Element * | add_element_to_vector (int id) |
Adds element to mesh data structures (element_vec_, element_ids_), returns pointer to this element. More... | |
void | init_element (Element *ele, unsigned int elm_id, unsigned int dim, RegionIdx region_idx, unsigned int partition_id, std::vector< unsigned int > node_ids) |
Initialize element. More... | |
void | output_internal_ngh_data () |
Output of neighboring data into raw output. More... | |
void | optimize () |
void | sort_permuted_nodes_elements (std::vector< int > new_node_ids, std::vector< int > new_elem_ids) |
Sort elements and nodes by order stored in permutation vectors. More... | |
Protected Attributes | |
unsigned int | n_bb_neigh |
unsigned int | n_vb_neigh |
unsigned int | max_edge_sides_ [3] |
Maximal number of sides per one edge in the actual mesh (set in make_neighbours_and_edges()). More... | |
bool | optimize_memory_locality |
RegionDB | region_db_ |
std::shared_ptr< Partitioning > | part_ |
std::shared_ptr< BIHTree > | bih_tree_ |
Input::Record | in_record_ |
MPI_Comm | comm_ |
vector< Element > | element_vec_ |
vector< ElementTmpData > | bc_element_tmp_ |
Hold data of boundary elements during reading mesh (allow to preserve correct order during reading of mix bulk-boundary element) More... | |
unsigned int | bulk_size_ |
Count of bulk elements. More... | |
unsigned int | boundary_loaded_size_ |
Count of boundary elements loaded from mesh file. More... | |
BidirectionalMap< int > | element_ids_ |
Maps element ids to indexes into vector element_vec_. More... | |
Armor::Array< double > | nodes_ |
BidirectionalMap< int > | node_ids_ |
Maps node ids to indexes into vector node_vec_. More... | |
std::vector< EdgeData > | edges |
Vector of MH edges, this should not be part of the geometrical mesh. More... | |
std::vector< unsigned int > | node_permutation_ |
Vector of node permutations of optimized mesh (see class MeshOptimizer) More... | |
std::vector< unsigned int > | elem_permutation_ |
Vector of element permutations of optimized mesh (see class MeshOptimizer) More... | |
Private Member Functions | |
void | distribute_nodes () |
Fill array node_4_loc_ and create object node_ds_ according to element distribution. More... | |
Private Attributes | |
LongIdx * | row_4_el |
Index set assigning to global element index the local index used in parallel vectors. More... | |
LongIdx * | el_4_loc |
Index set assigning to local element index its global index. More... | |
Distribution * | el_ds |
Parallel distribution of elements. More... | |
LongIdx * | node_4_loc_ |
Index set assigning to local node index its global index. More... | |
Distribution * | node_ds_ |
Parallel distribution of nodes. Depends on elements distribution. More... | |
unsigned int | n_local_nodes_ |
Hold number of local nodes (own + ghost), value is equal with size of node_4_loc array. More... | |
BCMesh * | bc_mesh_ |
Boundary mesh, object is created only if it's necessary. More... | |
Friends | |
class | Edge |
class | Side |
class | RegionSetBase |
class | Element |
class | BIHTree |
class | Boundary |
class | BCMesh |
template<int spacedim> | |
class | ElementAccessor |
template<int spacedim> | |
class | NodeAccessor |
typedef enum Mesh::IntersectionSearch Mesh::IntersectionSearch |
Types of search algorithm for finding intersection candidates.
anonymous enum |
Labels for coordinate indexes in arma::vec3 representing vectors and points.
Enumerator | |
---|---|
x_coord | |
y_coord | |
z_coord |
Mesh::Mesh | ( | ) |
Mesh::Mesh | ( | Input::Record | in_record, |
MPI_Comm | com = MPI_COMM_WORLD |
||
) |
void Mesh::add_element | ( | unsigned int | elm_id, |
unsigned int | dim, | ||
unsigned int | region_id, | ||
unsigned int | partition_id, | ||
std::vector< unsigned int > | node_ids | ||
) |
|
protected |
void Mesh::add_node | ( | unsigned int | node_id, |
arma::vec3 | coords | ||
) |
void Mesh::add_physical_name | ( | unsigned int | dim, |
unsigned int | id, | ||
std::string | name | ||
) |
void Mesh::check_and_finish | ( | ) |
|
virtual |
|
virtual |
|
inline |
|
protected |
|
protected |
|
protected |
unsigned int Mesh::create_boundary_elements | ( | ) |
|
protected |
Mesh::DECLARE_EXCEPTION | ( | ExcBdrElemMatchRegular | , |
<< "Boundary element (id: "<< EI_ElemId::val<< ") match a regular element (id: "<< EI_ElemIdOther::val<< ") of lower dimension.\n" | |||
) |
Mesh::DECLARE_EXCEPTION | ( | ExcDuplicateBoundary | , |
<< "Duplicate boundary elements! \n"<< "Element id: "<< EI_ElemLast::val<< " on region name: "<< EI_RegLast::val<< "\n"<< "Element id: "<< EI_ElemNew::val<< " on region name: "<< EI_RegNew::val<< "\n" | |||
) |
Mesh::DECLARE_EXCEPTION | ( | ExcElmWrongOrder | , |
<< "Element IDs in non-increasing | order, | ||
ID:"<< EI_ElemId::val<< "\n" | |||
) |
Mesh::DECLARE_EXCEPTION | ( | ExcRegionElmDiffDim | , |
<< "User defined region "<< EI_Region::qval<< " (id "<< EI_RegIdx::val<< ") by 'From_Elements' cannot have elements of different dimensions.\n"<< "Thrown due to: dim "<< EI_Dim::val<< " neq dim "<< EI_DimOther::val<< " (ele id "<< EI_ElemId::val<< ").\n"<< "Split elements by | dim, | ||
create separate regions and then possibly use Union.\n" | |||
) |
Mesh::DECLARE_EXCEPTION | ( | ExcTooMatchingIds | , |
<< "Mesh: Duplicate dim-join lower dim elements: "<< EI_ElemId::val<< " | , | ||
"<< EI_ElemIdOther::val<< ".\n" | |||
) |
|
private |
|
inline |
|
virtual |
Create and return ElementAccessor to element of given idx.
Reimplemented in BCMesh.
|
inline |
|
protected |
|
virtual |
|
inline |
|
protected |
Remove elements with dimension not equal to dim
from element_list
. Index of the first element of dimension dim-1
, is returned in element_idx
. If no such element is found the method returns false, if one such element is found the method returns true, if more elements are found we report an user input error.
Definition at line 513 of file mesh.cc.
|
inline |
BCMesh * Mesh::get_bc_mesh | ( | ) |
const BIHTree & Mesh::get_bih_tree | ( | ) |
|
inline |
|
inline |
|
inline |
std::vector< BoundingBox > Mesh::get_element_boxes | ( | ) |
|
static |
|
static |
Mesh::IntersectionSearch Mesh::get_intersection_search | ( | ) |
|
virtual |
|
inline |
|
inline |
|
virtual |
Returns pointer to partitioning object. Partitioning is created during setup_topology.
Reimplemented in BCMesh.
Definition at line 261 of file mesh.cc.
|
inline |
double Mesh::global_snap_radius | ( | ) | const |
|
protected |
Part of the constructor whichdoes not depedn on input record. Initializes node-side numbering according to RefElement.
Definition at line 144 of file mesh.cc.
|
protected |
void Mesh::init_element_vector | ( | unsigned int | size | ) |
void Mesh::init_from_input | ( | ) |
Reads input record, creates regions, read the mesh, setup topology. creates region sets.
void Mesh::init_node_vector | ( | unsigned int | size | ) |
void Mesh::intersect_element_lists | ( | vector< unsigned int > const & | nodes_list, |
vector< unsigned int > & | intersection_element_list | ||
) |
Find intersection of element lists given by Mesh::node_elements_ for elements givne by nodes_list
parameter. The result is placed into vector intersection_element_list
. If the node_list
is empty, and empty intersection is returned.
Definition at line 481 of file mesh.cc.
|
protected |
|
protected |
This replaces read_neighbours() in order to avoid using NGH preprocessor.
TODO:
Avoid maps:
5) need not to have temporary array for Edges, only postpone setting pointers in elements and set them after edges are found; we can temporary save Edge index instead of pointer in Neigbours and elements
6) Try replace Edge * by indexes in Neigbours and elements (anyway we have mesh pointer in elements so it is accessible also from Neigbours)
TODO:
Definition at line 549 of file mesh.cc.
|
inline |
|
inline |
MixedMeshIntersections & Mesh::mixed_intersections | ( | ) |
|
protected |
Possibly modify region id of elements sets by user in "regions" part of input file.
TODO: This method needs check in issue 'Review mesh setting'. Changes have been done during generalized region key and may be causing problems during the further development.
Definition at line 291 of file mesh.cc.
|
inline |
|
inlinevirtual |
|
inline |
|
inlinevirtual |
unsigned int Mesh::n_sides | ( | ) | const |
unsigned int Mesh::n_vb_neighbours | ( | ) | const |
NodeAccessor< 3 > Mesh::node | ( | unsigned int | idx | ) | const |
Create and return NodeAccessor to node of given idx.
Definition at line 919 of file mesh.cc.
|
inline |
|
inline |
Range< NodeAccessor< 3 > > Mesh::node_range | ( | ) | const |
|
protected |
Apply functionality of MeshOptimizer to sort nodes and elements.
Use Hilbert curve, need call sort_permuted_nodes_elements method.
Definition at line 417 of file mesh.cc.
|
protected |
void Mesh::permute_tetrahedron | ( | unsigned int | elm_idx, |
std::vector< unsigned int > | permutation_vec | ||
) |
void Mesh::permute_triangle | ( | unsigned int | elm_idx, |
std::vector< unsigned int > | permutation_vec | ||
) |
void Mesh::read_gmsh_from_stream | ( | istream & | in | ) |
Reads mesh from stream.
Method is especially used in unit tests.
void Mesh::read_regions_from_input | ( | Input::Array | region_list | ) |
Reads elements and their affiliation to regions and region sets defined by user in input file Format of input record is defined in method RegionSetBase::get_input_type()
region_list | Array input AbstractRecords which define regions, region sets and elements |
Definition at line 1171 of file mesh.cc.
void Mesh::setup_topology | ( | ) |
|
protected |
Mesh::TYPEDEF_ERR_INFO | ( | EI_Dim | , |
unsigned int | |||
) |
Mesh::TYPEDEF_ERR_INFO | ( | EI_DimOther | , |
unsigned int | |||
) |
Mesh::TYPEDEF_ERR_INFO | ( | EI_ElemId | , |
int | |||
) |
Mesh::TYPEDEF_ERR_INFO | ( | EI_ElemIdOther | , |
int | |||
) |
Mesh::TYPEDEF_ERR_INFO | ( | EI_ElemLast | , |
int | |||
) |
Mesh::TYPEDEF_ERR_INFO | ( | EI_ElemNew | , |
int | |||
) |
Mesh::TYPEDEF_ERR_INFO | ( | EI_RegIdx | , |
unsigned int | |||
) |
Mesh::TYPEDEF_ERR_INFO | ( | EI_Region | , |
std::string | |||
) |
Mesh::TYPEDEF_ERR_INFO | ( | EI_RegLast | , |
std::string | |||
) |
Mesh::TYPEDEF_ERR_INFO | ( | EI_RegNew | , |
std::string | |||
) |
|
friend |
|
friend |
|
friend |
|
protected |
|
private |
|
protected |
|
mutable |
|
protected |
|
protected |
|
protected |
|
protected |
|
private |
|
private |
|
protected |
Vector of element permutations of optimized mesh (see class MeshOptimizer)
|
protected |
|
protected |
std::shared_ptr<MixedMeshIntersections> Mesh::intersections |
|
protected |
Maximal number of sides per one edge in the actual mesh (set in make_neighbours_and_edges()).
|
private |
|
private |
|
private |
|
protected |
|
protected |
Vector of node permutations of optimized mesh (see class MeshOptimizer)
|
protected |
|
protected |
|
protected |
|
protected |
|
private |
DuplicateNodes* Mesh::tree |