Flow123d
release_3.0.0-684-g928e266
|
#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) | |
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 () | |
Mesh (Input::Record in_record, MPI_Comm com=MPI_COMM_WORLD) | |
void | reinit (Input::Record in_record) |
virtual | ~Mesh () |
Destructor. More... | |
virtual unsigned int | n_nodes () const |
unsigned int | n_boundaries () const |
unsigned int | n_edges () 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 |
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 bool | check_compatible_mesh (Mesh &mesh, vector< LongIdx > &bulk_elements_id, vector< LongIdx > &boundary_elements_id) |
virtual ElementAccessor< 3 > | element_accessor (unsigned int idx) const |
Create and return ElementAccessor to element of given idx. More... | |
NodeAccessor< 3 > | node_accessor (unsigned int idx) const |
Create and return NodeAccessor to node of given idx. More... | |
void | read_regions_from_input (Input::Array region_list) |
vector< vector< unsigned int > > const & | 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... | |
virtual unsigned int | n_elements (bool boundary=false) 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... | |
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... | |
void | check_element_size (unsigned int elem_idx) const |
Check if given index is in element_vec_. More... | |
void | create_boundary_elements () |
Create boundary elements from data of temporary structure, this method MUST be call after read mesh from. 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< Boundary > | boundary_ |
std::vector< Edge > | edges |
Vector of MH edges, this should not be part of the geometrical mesh. More... | |
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 | 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 | modify_element_ids (const RegionDB::MapElementIDToRegionID &map) |
Element * | add_element_to_vector (int id, bool boundary=false) |
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... | |
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... | |
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... | |
vector< Node > | node_vec_ |
BidirectionalMap< int > | node_ids_ |
Maps node ids to indexes into vector node_vec_. 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... | |
BCMesh * | bc_mesh_ |
Boundary mesh, object is created only if it's necessary. More... | |
ofstream | raw_ngh_output_file |
Friends | |
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 | ( | ) |
|
inline |
|
protected |
|
protected |
void Mesh::create_boundary_elements | ( | ) |
|
protected |
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" | |||
) |
|
inline |
|
virtual |
Create and return ElementAccessor to element of given idx.
Reimplemented in BCMesh.
|
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 376 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 |
|
virtual |
Returns pointer to partitioning object. Partitioning is created during setup_topology.
Reimplemented in BCMesh.
Definition at line 245 of file mesh.cc.
|
inline |
double Mesh::global_snap_radius | ( | ) | const |
|
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 344 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 413 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 275 of file mesh.cc.
|
inline |
|
inlinevirtual |
|
inlinevirtual |
unsigned int Mesh::n_sides | ( | ) | const |
unsigned int Mesh::n_vb_neighbours | ( | ) | const |
NodeAccessor< 3 > Mesh::node_accessor | ( | unsigned int | idx | ) | const |
Create and return NodeAccessor to node of given idx.
Definition at line 725 of file mesh.cc.
|
inline |
Range< NodeAccessor< 3 > > Mesh::node_range | ( | ) | const |
|
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 884 of file mesh.cc.
void Mesh::reinit | ( | Input::Record | in_record | ) |
void Mesh::setup_topology | ( | ) |
Mesh::TYPEDEF_ERR_INFO | ( | EI_ElemLast | , |
int | |||
) |
Mesh::TYPEDEF_ERR_INFO | ( | EI_ElemNew | , |
int | |||
) |
Mesh::TYPEDEF_ERR_INFO | ( | EI_RegLast | , |
std::string | |||
) |
Mesh::TYPEDEF_ERR_INFO | ( | EI_RegNew | , |
std::string | |||
) |
|
friend |
|
friend |
|
friend |
|
protected |
|
private |
|
protected |
|
protected |
|
protected |
|
protected |
std::vector<Edge> Mesh::edges |
|
private |
|
private |
|
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()).
|
protected |
|
protected |
|
protected |
|
private |
DuplicateNodes* Mesh::tree |