MFEM
v4.3.0
Finite element discretization library
|
A parallel extension of the NCMesh class. More...
#include <pncmesh.hpp>
Classes | |
class | ElementSet |
class | ElementValueMessage |
class | NeighborDerefinementMessage |
class | NeighborElementRankMessage |
class | NeighborRefinementMessage |
class | RebalanceDofMessage |
class | RebalanceMessage |
Public Types | |
typedef short | GroupId |
typedef std::vector< int > | CommGroup |
Public Types inherited from mfem::NCMesh | |
typedef std::int64_t | RefCoord |
Public Member Functions | |
ParNCMesh (MPI_Comm comm, const NCMesh &ncmesh, int *part=NULL) | |
Construct by partitioning a serial NCMesh. More... | |
ParNCMesh (MPI_Comm comm, std::istream &input, int version, int &curved, int &is_nc) | |
ParNCMesh (const ParNCMesh &other) | |
Deep copy of another instance. More... | |
virtual | ~ParNCMesh () |
virtual void | Refine (const Array< Refinement > &refinements) |
virtual void | LimitNCLevel (int max_nc_level) |
Parallel version of NCMesh::LimitNCLevel. More... | |
virtual void | CheckDerefinementNCLevel (const Table &deref_table, Array< int > &level_ok, int max_nc_level) |
virtual void | Derefine (const Array< int > &derefs) |
void | Rebalance (const Array< int > *custom_partition=NULL) |
int | GetNElements () const |
int | GetNGhostVertices () const |
int | GetNGhostEdges () const |
int | GetNGhostFaces () const |
int | GetNGhostElements () const |
const NCList & | GetSharedVertices () |
const NCList & | GetSharedEdges () |
const NCList & | GetSharedFaces () |
const NCList & | GetSharedList (int entity) |
Helper to get shared vertices/edges/faces ('entity' == 0/1/2 resp.). More... | |
int | GetFaceOrientation (int index) const |
Return (shared) face orientation relative to its owner element. More... | |
GroupId | GetEntityOwnerId (int entity, int index) |
Return vertex/edge/face ('entity' == 0/1/2, resp.) owner. More... | |
GroupId | GetEntityGroupId (int entity, int index) |
const CommGroup & | GetGroup (GroupId id) const |
Return a list of ranks contained in the group of the given ID. More... | |
bool | GroupContains (GroupId id, int rank) const |
Return true if group 'id' contains the given rank. More... | |
bool | IsGhost (int entity, int index) const |
Return true if the specified vertex/edge/face is a ghost. More... | |
int | ElementRank (int index) const |
int | GetMyRank () const |
void | SendRebalanceDofs (int old_ndofs, const Table &old_element_dofs, long old_global_offset, FiniteElementSpace *space) |
Use the communication pattern from last Rebalance() to send element DOFs. More... | |
void | RecvRebalanceDofs (Array< int > &elements, Array< long > &dofs) |
Receive element DOFs sent by SendRebalanceDofs(). More... | |
const Array< int > & | GetRebalanceOldIndex () const |
const Array< int > & | GetDerefineOldRanks () const |
template<typename Type > | |
void | SynchronizeDerefinementData (Array< Type > &elem_data, const Table &deref_table) |
virtual void | GetBoundaryClosure (const Array< int > &bdr_attr_is_ess, Array< int > &bdr_vertices, Array< int > &bdr_edges) |
virtual void | Trim () |
Save memory by releasing all non-essential and cached data. More... | |
long | MemoryUsage (bool with_base=true) const |
Return total number of bytes allocated. More... | |
int | PrintMemoryDetail (bool with_base=true) const |
void | GetDebugMesh (Mesh &debug_mesh) const |
Public Member Functions inherited from mfem::NCMesh | |
NCMesh (const Mesh *mesh) | |
NCMesh (std::istream &input, int version, int &curved, int &is_nc) | |
NCMesh (const NCMesh &other) | |
Deep copy of another instance. More... | |
virtual | ~NCMesh () |
int | Dimension () const |
int | SpaceDimension () const |
int | GetNVertices () const |
int | GetNEdges () const |
int | GetNFaces () const |
const Table & | GetDerefinementTable () |
const NCList & | GetFaceList () |
Return the current list of conforming and nonconforming faces. More... | |
const NCList & | GetEdgeList () |
Return the current list of conforming and nonconforming edges. More... | |
const NCList & | GetVertexList () |
const NCList & | GetNCList (int entity) |
Return vertex/edge/face list (entity = 0/1/2, respectively). More... | |
void | MarkCoarseLevel () |
const CoarseFineTransformations & | GetRefinementTransforms () |
const CoarseFineTransformations & | GetDerefinementTransforms () |
void | ClearTransforms () |
Free all internal data created by the above three functions. More... | |
void | GetEdgeVertices (const MeshId &edge_id, int vert_index[2], bool oriented=true) const |
Return Mesh vertex indices of an edge identified by 'edge_id'. More... | |
int | GetEdgeNCOrientation (const MeshId &edge_id) const |
int | GetFaceVerticesEdges (const MeshId &face_id, int vert_index[4], int edge_index[4], int edge_orientation[4]) const |
int | GetEdgeMaster (int v1, int v2) const |
Geometry::Type | GetElementGeometry (int index) const |
Return element geometry type. index is the Mesh element number. More... | |
Geometry::Type | GetFaceGeometry (int index) const |
Return face geometry type. index is the Mesh face number. More... | |
int | GetNumRootElements () |
Return the number of root elements. More... | |
int | GetElementDepth (int i) const |
Return the distance of leaf 'i' from the root. More... | |
int | GetElementSizeReduction (int i) const |
void | GetElementFacesAttributes (int i, Array< int > &faces, Array< int > &fattr) const |
Return the faces and face attributes of leaf element 'i'. More... | |
void | Print (std::ostream &out) const |
I/O: Print the mesh in "MFEM NC mesh v1.0" format. More... | |
bool | IsLegacyLoaded () const |
I/O: Return true if the mesh was loaded from the legacy v1.1 format. More... | |
void | LegacyToNewVertexOrdering (Array< int > &order) const |
I/O: Return a map from old (v1.1) vertex indices to new vertex indices. More... | |
long | MemoryUsage () const |
Return total number of bytes allocated. More... | |
int | PrintMemoryDetail () const |
void | DebugLeafOrder (std::ostream &out) const |
void | DebugDump (std::ostream &out) const |
Protected Types | |
typedef std::vector< CommGroup > | GroupList |
typedef std::map< CommGroup, GroupId > | GroupMap |
Protected Types inherited from mfem::NCMesh | |
typedef std::map< std::string, int > | RefPathMap |
Protected Member Functions | |
void | GetConformingSharedStructures (class ParMesh &pmesh) |
void | GetFaceNeighbors (class ParMesh &pmesh) |
virtual void | Update () |
int | Partition (long index, long total_elements) const |
Return the processor number for a global element number. More... | |
int | InitialPartition (int index) const |
Helper to get the partitioning when the serial mesh gets split initially. More... | |
long | PartitionFirstIndex (int rank, long total_elements) const |
Return the global index of the first element owned by processor 'rank'. More... | |
virtual void | BuildFaceList () |
virtual void | BuildEdgeList () |
virtual void | BuildVertexList () |
virtual void | ElementSharesFace (int elem, int local, int face) |
virtual void | ElementSharesEdge (int elem, int local, int enode) |
virtual void | ElementSharesVertex (int elem, int local, int vnode) |
GroupId | GetGroupId (const CommGroup &group) |
GroupId | GetSingletonGroup (int rank) |
void | InitOwners (int num, Array< GroupId > &entity_owner) |
void | MakeSharedList (const NCList &list, NCList &shared) |
void | AddConnections (int entity, int index, const Array< int > &ranks) |
void | CalculatePMatrixGroups () |
void | CreateGroups (int nentities, Array< Connection > &index_rank, Array< GroupId > &entity_group) |
void | CalcFaceOrientations () |
void | UpdateLayers () |
void | MakeSharedTable (int ngroups, int ent, Array< int > &shared_local, Table &group_shared, Array< char > *entity_geom=NULL, char geom=0) |
void | AdjustMeshIds (Array< MeshId > ids[], int rank) |
void | ChangeVertexMeshIdElement (NCMesh::MeshId &id, int elem) |
void | ChangeEdgeMeshIdElement (NCMesh::MeshId &id, int elem) |
void | ChangeRemainingMeshIds (Array< MeshId > &ids, int pos, const Array< Pair< int, int > > &find) |
void | EncodeMeshIds (std::ostream &os, Array< MeshId > ids[]) |
void | DecodeMeshIds (std::istream &is, Array< MeshId > ids[]) |
void | EncodeGroups (std::ostream &os, const Array< GroupId > &ids) |
void | DecodeGroups (std::istream &is, Array< GroupId > &ids) |
bool | CheckElementType (int elem, int type) |
void | ElementNeighborProcessors (int elem, Array< int > &ranks) |
void | NeighborProcessors (Array< int > &neighbors) |
void | Prune () |
bool | PruneTree (int elem) |
Internal. Recursive part of Prune(). More... | |
void | RedistributeElements (Array< int > &new_ranks, int target_elements, bool record_comm) |
void | ClearAuxPM () |
long | GroupsMemoryUsage () const |
Protected Member Functions inherited from mfem::NCMesh | |
void | GetMeshComponents (Mesh &mesh) const |
Fill Mesh::{vertices,elements,boundary} for the current finest level. More... | |
void | OnMeshUpdated (Mesh *mesh) |
void | MakeTopologyOnly () |
void | UpdateLeafElements () |
void | UpdateVertices () |
update Vertex::index and vertex_nodeId More... | |
void | CollectLeafElements (int elem, int state, Array< int > &ghosts, int &counter) |
void | InitRootState (int root_count) |
void | InitGeomFlags () |
bool | HavePrisms () const |
bool | HaveTets () const |
bool | IsGhost (const Element &el) const |
void | RefineElement (int elem, char ref_type) |
void | DerefineElement (int elem) |
int | AddElement (const Element &el) |
void | FreeElement (int id) |
int | NewHexahedron (int n0, int n1, int n2, int n3, int n4, int n5, int n6, int n7, int attr, int fattr0, int fattr1, int fattr2, int fattr3, int fattr4, int fattr5) |
int | NewWedge (int n0, int n1, int n2, int n3, int n4, int n5, int attr, int fattr0, int fattr1, int fattr2, int fattr3, int fattr4) |
int | NewTetrahedron (int n0, int n1, int n2, int n3, int attr, int fattr0, int fattr1, int fattr2, int fattr3) |
int | NewQuadrilateral (int n0, int n1, int n2, int n3, int attr, int eattr0, int eattr1, int eattr2, int eattr3) |
int | NewTriangle (int n0, int n1, int n2, int attr, int eattr0, int eattr1, int eattr2) |
int | NewSegment (int n0, int n1, int attr, int vattr1, int vattr2) |
mfem::Element * | NewMeshElement (int geom) const |
int | QuadFaceSplitType (int v1, int v2, int v3, int v4, int mid[5]=NULL) const |
bool | TriFaceSplit (int v1, int v2, int v3, int mid[3]=NULL) const |
void | ForceRefinement (int vn1, int vn2, int vn3, int vn4) |
void | FindEdgeElements (int vn1, int vn2, int vn3, int vn4, Array< MeshId > &prisms) const |
void | CheckAnisoPrism (int vn1, int vn2, int vn3, int vn4, const Refinement *refs, int nref) |
void | CheckAnisoFace (int vn1, int vn2, int vn3, int vn4, int mid12, int mid34, int level=0) |
void | CheckIsoFace (int vn1, int vn2, int vn3, int vn4, int en1, int en2, int en3, int en4, int midf) |
void | ReparentNode (int node, int new_p1, int new_p2) |
int | FindMidEdgeNode (int node1, int node2) const |
int | GetMidEdgeNode (int node1, int node2) |
int | GetMidFaceNode (int en1, int en2, int en3, int en4) |
void | ReferenceElement (int elem) |
void | UnreferenceElement (int elem, Array< int > &elemFaces) |
Face * | GetFace (Element &elem, int face_no) |
void | RegisterFaces (int elem, int *fattr=NULL) |
void | DeleteUnusedFaces (const Array< int > &elemFaces) |
void | CollectDerefinements (int elem, Array< Connection > &list) |
int | RetrieveNode (const Element &el, int index) |
Return el.node[index] correctly, even if the element is refined. More... | |
int | FindNodeExt (const Element &el, int node, bool abort=true) |
Extended version of find_node: works if 'el' is refined. More... | |
int | ReorderFacePointMat (int v0, int v1, int v2, int v3, int elem, const PointMatrix &pm, PointMatrix &reordered) const |
void | TraverseQuadFace (int vn0, int vn1, int vn2, int vn3, const PointMatrix &pm, int level, Face *eface[4], MatrixMap &matrix_map) |
bool | TraverseTriFace (int vn0, int vn1, int vn2, const PointMatrix &pm, int level, MatrixMap &matrix_map) |
void | TraverseTetEdge (int vn0, int vn1, const Point &p0, const Point &p1, MatrixMap &matrix_map) |
void | TraverseEdge (int vn0, int vn1, double t0, double t1, int flags, int level, MatrixMap &matrix_map) |
void | FindSetNeighbors (const Array< char > &elem_set, Array< int > *neighbors, Array< char > *neighbor_set=NULL) |
void | FindNeighbors (int elem, Array< int > &neighbors, const Array< int > *search_set=NULL) |
void | NeighborExpand (const Array< int > &elems, Array< int > &expanded, const Array< int > *search_set=NULL) |
void | CollectEdgeVertices (int v0, int v1, Array< int > &indices) |
void | CollectTriFaceVertices (int v0, int v1, int v2, Array< int > &indices) |
void | CollectQuadFaceVertices (int v0, int v1, int v2, int v3, Array< int > &indices) |
void | BuildElementToVertexTable () |
void | UpdateElementToVertexTable () |
int | GetVertexRootCoord (int elem, RefCoord coord[3]) const |
void | CollectIncidentElements (int elem, const RefCoord coord[3], Array< int > &list) const |
void | FindVertexCousins (int elem, int local, Array< int > &cousins) const |
void | GetPointMatrix (Geometry::Type geom, const char *ref_path, DenseMatrix &matrix) |
void | TraverseRefinements (int elem, int coarse_index, std::string &ref_path, RefPathMap &map) |
void | InitDerefTransforms () |
void | SetDerefMatrixCodes (int parent, Array< int > &fine_coarse) |
const double * | CalcVertexPos (int node) const |
int | GetEdgeMaster (int node) const |
void | FindFaceNodes (int face, int node[4]) |
int | EdgeSplitLevel (int vn1, int vn2) const |
int | TriFaceSplitLevel (int vn1, int vn2, int vn3) const |
void | QuadFaceSplitLevel (int vn1, int vn2, int vn3, int vn4, int &h_level, int &v_level) const |
void | CountSplits (int elem, int splits[3]) const |
void | GetLimitRefinements (Array< Refinement > &refinements, int max_level) |
int | PrintVertexParents (std::ostream *out) const |
Print the "vertex_parents" section of the mesh file. More... | |
void | LoadVertexParents (std::istream &input) |
Load the vertex parent hierarchy from a mesh file. More... | |
int | PrintBoundary (std::ostream *out) const |
void | LoadBoundary (std::istream &input) |
Load the "boundary" section of the mesh file. More... | |
void | PrintCoordinates (std::ostream &out) const |
Print the "coordinates" section of the mesh file. More... | |
void | LoadCoordinates (std::istream &input) |
Load the "coordinates" section of the mesh file. More... | |
void | InitRootElements () |
Count root elements and initialize root_state. More... | |
int | CountTopLevelNodes () const |
Return the index of the last top-level node plus one. More... | |
bool | ZeroRootStates () const |
Return true if all root_states are zero. More... | |
void | LoadCoarseElements (std::istream &input) |
Load the element refinement hierarchy from a legacy mesh file. More... | |
void | CopyElements (int elem, const BlockArray< Element > &tmp_elements) |
void | LoadLegacyFormat (std::istream &input, int &curved, int &is_nc) |
Load the deprecated MFEM mesh v1.1 format for backward compatibility. More... | |
Static Protected Member Functions | |
static int | get_face_orientation (Face &face, Element &e1, Element &e2, int local[2]=NULL) |
Static Protected Member Functions inherited from mfem::NCMesh | |
static int | find_node (const Element &el, int node) |
static int | find_element_edge (const Element &el, int vn0, int vn1, bool abort=true) |
static int | find_local_face (int geom, int a, int b, int c) |
static const PointMatrix & | GetGeomIdentity (Geometry::Type geom) |
Protected Attributes | |
MPI_Comm | MyComm |
int | NRanks |
GroupList | groups |
GroupMap | group_id |
Array< GroupId > | entity_owner [3] |
Array< GroupId > | entity_pmat_group [3] |
Array< GroupId > | entity_conf_group [3] |
Array< int > | entity_elem_local [3] |
NCList | shared_vertices |
NCList | shared_edges |
NCList | shared_faces |
Array< char > | face_orient |
Array< char > | element_type |
Array< int > | ghost_layer |
list of elements whose 'element_type' == 2. More... | |
Array< int > | boundary_layer |
list of type 3 elements More... | |
Array< int > | tmp_owner |
Array< char > | tmp_shared_flag |
Array< Connection > | entity_index_rank [3] |
Array< int > | tmp_neighbors |
RebalanceDofMessage::Map | send_rebalance_dofs |
RebalanceDofMessage::Map | recv_rebalance_dofs |
Array< int > | old_index_or_rank |
Array< DenseMatrix * > | aux_pm_store |
Stores modified point matrices created by GetFaceNeighbors. More... | |
Protected Attributes inherited from mfem::NCMesh | |
int | Dim |
int | spaceDim |
dimensions of the elements and the vertex coordinates More... | |
int | MyRank |
used in parallel, or when loading a parallel file in serial More... | |
bool | Iso |
true if the mesh only contains isotropic refinements More... | |
int | Geoms |
bit mask of element geometries present, see InitGeomFlags() More... | |
bool | Legacy |
true if the mesh was loaded from the legacy v1.1 format More... | |
HashTable< Node > | nodes |
HashTable< Face > | faces |
BlockArray< Element > | elements |
Array< int > | free_element_ids |
Array< int > | root_state |
Array< double > | coordinates |
int | NElements |
int | NVertices |
int | NEdges |
int | NFaces |
int | NGhostElements |
int | NGhostVertices |
int | NGhostEdges |
int | NGhostFaces |
Array< int > | leaf_elements |
finest elements, in Mesh ordering (+ ghosts) More... | |
Array< int > | leaf_sfc_index |
natural tree ordering of leaf elements More... | |
Array< int > | vertex_nodeId |
vertex-index to node-id map, see UpdateVertices More... | |
NCList | face_list |
lazy-initialized list of faces, see GetFaceList More... | |
NCList | edge_list |
lazy-initialized list of edges, see GetEdgeList More... | |
NCList | vertex_list |
lazy-initialized list of vertices, see GetVertexList More... | |
Array< int > | boundary_faces |
subset of all faces, set by BuildFaceList More... | |
Array< char > | face_geom |
face geometry by face index, set by OnMeshUpdated More... | |
Table | element_vertex |
leaf-element to vertex table, see FindSetNeighbors More... | |
Array< Refinement > | ref_stack |
stack of scheduled refinements (temporary) More... | |
HashTable< Node > | shadow |
temporary storage for reparented nodes More... | |
Array< Triple< int, int, int > > | reparents |
scheduled node reparents (tmp) More... | |
Table | derefinements |
possible derefinements, see GetDerefinementTable More... | |
CoarseFineTransformations | transforms |
storage for data returned by Get[De]RefinementTransforms() More... | |
Array< int > | coarse_elements |
state of leaf_elements before Refine(), set by MarkCoarseLevel() More... | |
TmpVertex * | tmp_vertex |
Friends | |
class | ParMesh |
class | NeighborRowMessage |
Additional Inherited Members | |
Static Public Member Functions inherited from mfem::NCMesh | |
static void | GridSfcOrdering2D (int width, int height, Array< int > &coords) |
static void | GridSfcOrdering3D (int width, int height, int depth, Array< int > &coords) |
Static Protected Attributes inherited from mfem::NCMesh | |
static PointMatrix | pm_seg_identity |
static PointMatrix | pm_tri_identity |
static PointMatrix | pm_quad_identity |
static PointMatrix | pm_tet_identity |
static PointMatrix | pm_prism_identity |
static PointMatrix | pm_hex_identity |
static GeomInfo | GI [Geometry::NumGeom] |
A parallel extension of the NCMesh class.
The basic idea (and assumption) is that all processors share the coarsest layer ("root elements"). This has the advantage that refinements can easily be exchanged between processors when rebalancing since individual elements can be uniquely identified by the index of the root element and a path in the refinement tree.
Each leaf element is owned by one of the processors (NCMesh::Element::rank). The underlying NCMesh stores not only elements for the current ('MyRank') processor, but also a minimal layer of adjacent "ghost" elements owned by other processors. The ghost layer is synchronized after refinement.
The ghost layer contains all vertex-, edge- and face-neighbors of the current processor's region. It is used to determine constraining relations and ownership of DOFs on the processor boundary. Ghost elements are never seen by the rest of MFEM as they are skipped when a Mesh is created from the NCMesh.
The processor that owns a vertex, edge or a face (and in turn its DOFs) is currently defined to be the one with the lowest rank in the group of processors that share the entity.
Vertices, edges and faces that are not owned by this ('MyRank') processor are ghosts, and are numbered after all real vertices/edges/faces, i.e., they have indices greater than NVertices, NEdges, NFaces, respectively.
A shared vertex/edge/face is identified in an interprocessor message by a pair of numbers. The first number specifies an element in an ElementSet (typically sent at the beginning of the message) that contains the v/e/f. The second number is the local index of the v/e/f in that element.
Definition at line 64 of file pncmesh.hpp.
typedef std::vector<int> mfem::ParNCMesh::CommGroup |
Definition at line 140 of file pncmesh.hpp.
typedef short mfem::ParNCMesh::GroupId |
Definition at line 139 of file pncmesh.hpp.
|
protected |
Definition at line 266 of file pncmesh.hpp.
|
protected |
Definition at line 267 of file pncmesh.hpp.
mfem::ParNCMesh::ParNCMesh | ( | MPI_Comm | comm, |
const NCMesh & | ncmesh, | ||
int * | part = NULL |
||
) |
Construct by partitioning a serial NCMesh.
SFC partitioning is used by default. A user-specified partition can be passed in 'part', where part[i] is the desired MPI rank for element i.
Definition at line 28 of file pncmesh.cpp.
mfem::ParNCMesh::ParNCMesh | ( | MPI_Comm | comm, |
std::istream & | input, | ||
int | version, | ||
int & | curved, | ||
int & | is_nc | ||
) |
Load from a stream, parallel version. See the serial NCMesh::NCMesh counterpart for a description of the parameters.
Definition at line 49 of file pncmesh.cpp.
mfem::ParNCMesh::ParNCMesh | ( | const ParNCMesh & | other | ) |
Deep copy of another instance.
Definition at line 76 of file pncmesh.cpp.
|
virtual |
Definition at line 85 of file pncmesh.cpp.
|
protected |
Definition at line 461 of file pncmesh.cpp.
Adjust some of the MeshIds before encoding for recipient 'rank', so that they only reference elements that exist in the recipient's ref. tree.
Definition at line 2242 of file pncmesh.cpp.
|
protectedvirtual |
Reimplemented from mfem::NCMesh.
Definition at line 208 of file pncmesh.cpp.
|
protectedvirtual |
Reimplemented from mfem::NCMesh.
Definition at line 141 of file pncmesh.cpp.
|
protectedvirtual |
Reimplemented from mfem::NCMesh.
Definition at line 269 of file pncmesh.cpp.
|
protected |
Definition at line 572 of file pncmesh.cpp.
|
protected |
Definition at line 469 of file pncmesh.cpp.
|
protected |
Definition at line 2343 of file pncmesh.cpp.
|
protected |
Definition at line 2369 of file pncmesh.cpp.
|
protected |
Definition at line 2325 of file pncmesh.cpp.
|
virtual |
Parallel version of NCMesh::CheckDerefinementNCLevel.
Reimplemented from mfem::NCMesh.
Definition at line 1594 of file pncmesh.cpp.
|
protected |
Definition at line 661 of file pncmesh.cpp.
|
protected |
Definition at line 1163 of file pncmesh.cpp.
|
protected |
Definition at line 427 of file pncmesh.cpp.
Definition at line 2523 of file pncmesh.cpp.
Definition at line 2424 of file pncmesh.cpp.
|
virtual |
Parallel reimplementation of NCMesh::Derefine, keeps ghost layers in sync. The interface is identical.
Reimplemented from mfem::NCMesh.
Definition at line 1342 of file pncmesh.cpp.
|
protected |
Return a list of processors that own elements in the immediate neighborhood of 'elem' (i.e., vertex, edge and face neighbors), and are not 'MyRank'.
Definition at line 678 of file pncmesh.cpp.
|
inline |
Returns owner processor for element 'index'. This is normally MyRank but for index >= NElements (i.e., for ghosts) it may be something else.
Definition at line 197 of file pncmesh.hpp.
|
protectedvirtual |
Reimplemented from mfem::NCMesh.
Definition at line 183 of file pncmesh.cpp.
|
protectedvirtual |
Reimplemented from mfem::NCMesh.
Definition at line 118 of file pncmesh.cpp.
|
protectedvirtual |
Reimplemented from mfem::NCMesh.
Definition at line 246 of file pncmesh.cpp.
Definition at line 2481 of file pncmesh.cpp.
Definition at line 2381 of file pncmesh.cpp.
|
staticprotected |
Definition at line 543 of file pncmesh.cpp.
|
virtual |
Extension of NCMesh::GetBoundaryClosure. Filters out ghost vertices and ghost edges from 'bdr_vertices' and 'bdr_edges'.
Reimplemented from mfem::NCMesh.
Definition at line 598 of file pncmesh.cpp.
|
protected |
For compatibility with conforming code in ParMesh and ParFESpace. Initializes shared structures in ParMesh: gtopo, shared_*, group_s*, s*_l*. The ParMesh then acts as a parallel mesh cut along the NC interfaces.
Definition at line 787 of file pncmesh.cpp.
void mfem::ParNCMesh::GetDebugMesh | ( | Mesh & | debug_mesh | ) | const |
Extract a debugging Mesh containing all leaf elements, including ghosts. The debug mesh will have element attributes set to element rank + 1.
Definition at line 2687 of file pncmesh.cpp.
|
inline |
Get previous (pre-Derefine) fine element ranks. This complements the CoarseFineTransformations::embeddings array in parallel.
Definition at line 220 of file pncmesh.hpp.
|
inline |
Return the P matrix communication group ID for a vertex/edge/face. The groups are calculated specifically to match the P matrix construction algorithm and its communication pattern.
Definition at line 157 of file pncmesh.hpp.
|
inline |
Return vertex/edge/face ('entity' == 0/1/2, resp.) owner.
Definition at line 143 of file pncmesh.hpp.
|
protected |
Populate face neighbor members of ParMesh from the ghost layer, without communication.
Definition at line 886 of file pncmesh.cpp.
|
inline |
Return (shared) face orientation relative to its owner element.
Definition at line 134 of file pncmesh.hpp.
Return a list of ranks contained in the group of the given ID.
Definition at line 169 of file pncmesh.hpp.
|
protected |
Definition at line 391 of file pncmesh.cpp.
|
inline |
Definition at line 205 of file pncmesh.hpp.
|
inline |
Definition at line 109 of file pncmesh.hpp.
|
inline |
Definition at line 112 of file pncmesh.hpp.
|
inline |
Definition at line 114 of file pncmesh.hpp.
|
inline |
Definition at line 113 of file pncmesh.hpp.
|
inline |
Definition at line 111 of file pncmesh.hpp.
|
inline |
Get previous indices (pre-Rebalance) of current elements. Index of -1 indicates that an element didn't exist in the mesh before.
Definition at line 216 of file pncmesh.hpp.
|
inline |
Definition at line 119 of file pncmesh.hpp.
|
inline |
Definition at line 120 of file pncmesh.hpp.
|
inline |
Helper to get shared vertices/edges/faces ('entity' == 0/1/2 resp.).
Definition at line 123 of file pncmesh.hpp.
|
inline |
Definition at line 118 of file pncmesh.hpp.
|
protected |
Definition at line 407 of file pncmesh.cpp.
bool mfem::ParNCMesh::GroupContains | ( | GroupId | id, |
int | rank | ||
) | const |
Return true if group 'id' contains the given rank.
Definition at line 416 of file pncmesh.cpp.
|
protected |
Definition at line 2745 of file pncmesh.cpp.
|
inlineprotected |
Helper to get the partitioning when the serial mesh gets split initially.
Definition at line 305 of file pncmesh.hpp.
Definition at line 301 of file pncmesh.cpp.
|
inline |
Return true if the specified vertex/edge/face is a ghost.
Definition at line 179 of file pncmesh.hpp.
|
virtual |
Parallel version of NCMesh::LimitNCLevel.
Reimplemented from mfem::NCMesh.
Definition at line 1324 of file pncmesh.cpp.
Definition at line 311 of file pncmesh.cpp.
|
protected |
Definition at line 727 of file pncmesh.cpp.
long mfem::ParNCMesh::MemoryUsage | ( | bool | with_base = true | ) | const |
Return total number of bytes allocated.
Definition at line 2768 of file pncmesh.cpp.
|
protected |
Get a list of ranks that own elements in the neighborhood of our region. NOTE: MyRank is not included.
Definition at line 710 of file pncmesh.cpp.
|
inlineprotected |
Return the processor number for a global element number.
Definition at line 301 of file pncmesh.hpp.
|
inlineprotected |
Return the global index of the first element owned by processor 'rank'.
Definition at line 309 of file pncmesh.hpp.
int mfem::ParNCMesh::PrintMemoryDetail | ( | bool | with_base = true | ) | const |
Definition at line 2802 of file pncmesh.cpp.
|
protected |
Traverse the (local) refinement tree and determine which subtrees are no longer needed, i.e., their leaves are not owned by us nor are they our ghosts. These subtrees are then derefined.
Definition at line 1212 of file pncmesh.cpp.
|
protected |
Internal. Recursive part of Prune().
Definition at line 1175 of file pncmesh.cpp.
void mfem::ParNCMesh::Rebalance | ( | const Array< int > * | custom_partition = NULL | ) |
Migrate leaf elements of the global refinement hierarchy (including ghost elements) so that each processor owns the same number of leaves (+-1). The default partitioning strategy is based on equal splitting of the space-filling sequence of leaf elements (custom_partition == NULL). Alternatively, a used-defined element-rank assignment array can be passed.
Definition at line 1652 of file pncmesh.cpp.
Receive element DOFs sent by SendRebalanceDofs().
Definition at line 2021 of file pncmesh.cpp.
|
protected |
Assign new Element::rank to leaf elements and send them to their new owners, keeping the ghost layer up to date. Used by Rebalance() and Derefine(). 'target_elements' is the number of elements this rank is supposed to own after the exchange. If this number is not known a priori, the parameter can be set to -1, but more expensive communication (synchronous sends and a barrier) will be used in that case.
Definition at line 1713 of file pncmesh.cpp.
|
virtual |
An override of NCMesh::Refine, which is called eventually, after making sure that refinements that occur on the processor boundary are sent to the neighbor processors so they can keep their ghost layers up to date.
Reimplemented from mfem::NCMesh.
Definition at line 1246 of file pncmesh.cpp.
void mfem::ParNCMesh::SendRebalanceDofs | ( | int | old_ndofs, |
const Table & | old_element_dofs, | ||
long | old_global_offset, | ||
FiniteElementSpace * | space | ||
) |
Use the communication pattern from last Rebalance() to send element DOFs.
Definition at line 1988 of file pncmesh.cpp.
template void mfem::ParNCMesh::SynchronizeDerefinementData< double > | ( | Array< Type > & | elem_data, |
const Table & | deref_table | ||
) |
Exchange element data for derefinements that straddle processor boundaries. 'elem_data' is enlarged and filled with ghost values.
Definition at line 1514 of file pncmesh.cpp.
|
virtual |
Save memory by releasing all non-essential and cached data.
Reimplemented from mfem::NCMesh.
Definition at line 2704 of file pncmesh.cpp.
|
protectedvirtual |
Apart from the primary data structure, which is the element/node/face hierarchy, there is secondary data that is derived from the primary data and needs to be updated when the primary data changes. Update() takes care of that and needs to be called after each refinement and derefinement.
Reimplemented from mfem::NCMesh.
Definition at line 90 of file pncmesh.cpp.
|
protected |
Definition at line 623 of file pncmesh.cpp.
|
friend |
Definition at line 541 of file pncmesh.hpp.
|
friend |
Definition at line 249 of file pncmesh.hpp.
|
protected |
Stores modified point matrices created by GetFaceNeighbors.
Definition at line 536 of file pncmesh.hpp.
|
protected |
list of type 3 elements
Definition at line 296 of file pncmesh.hpp.
|
protected |
Type of each leaf element: 1 - our element (rank == MyRank), 3 - our element, and neighbor to the ghost layer, 2 - ghost layer element (existing element, but rank != MyRank), 0 - element beyond the ghost layer, may not be a real element. Note: indexed by Element::index. See also UpdateLayers().
Definition at line 293 of file pncmesh.hpp.
Definition at line 278 of file pncmesh.hpp.
|
protected |
Definition at line 280 of file pncmesh.hpp.
|
protected |
Definition at line 325 of file pncmesh.hpp.
Definition at line 273 of file pncmesh.hpp.
Definition at line 275 of file pncmesh.hpp.
|
protected |
Definition at line 285 of file pncmesh.hpp.
|
protected |
list of elements whose 'element_type' == 2.
Definition at line 295 of file pncmesh.hpp.
|
protected |
Definition at line 270 of file pncmesh.hpp.
|
protected |
Definition at line 269 of file pncmesh.hpp.
|
protected |
Definition at line 263 of file pncmesh.hpp.
|
protected |
Definition at line 264 of file pncmesh.hpp.
|
protected |
After Rebalance, this array holds the old element indices, or -1 if an element didn't exist in the mesh previously. After Derefine, it holds the ranks of the old (potentially non-existent) fine elements.
Definition at line 533 of file pncmesh.hpp.
|
protected |
Definition at line 528 of file pncmesh.hpp.
|
protected |
Recorded communication pattern from last Rebalance. Used by Send/RecvRebalanceDofs to ship element DOFs.
Definition at line 527 of file pncmesh.hpp.
|
protected |
Definition at line 283 of file pncmesh.hpp.
|
protected |
Definition at line 283 of file pncmesh.hpp.
|
protected |
Definition at line 283 of file pncmesh.hpp.
|
protected |
Definition at line 404 of file pncmesh.hpp.
|
protected |
Definition at line 323 of file pncmesh.hpp.
|
protected |
Definition at line 324 of file pncmesh.hpp.