FEDomainConnectivity

This data type handles the connectivity between two domains.

Structure

This class is designed for generating mapping between two domains or parts of two domains.

The structure of DomainConnectivity is given below.

TYPE :: DomainConnectivity_
  PRIVATE
  LOGICAL(LGT), PUBLIC :: isInitiated = .FALSE.
  LOGICAL(LGT), PUBLIC :: isFacetToCell = .FALSE.
  LOGICAL(LGT), PUBLIC :: isNodeToNode = .FALSE.
  LOGICAL(LGT), PUBLIC :: isCellToCell = .FALSE.
  INTEGER(I4B), ALLOCATABLE :: nodeToNode(:)
  INTEGER(I4B), ALLOCATABLE :: cellToCell(:)
  INTEGER(I4B), ALLOCATABLE :: cellToCellExtraData(:, :)
  TYPE(FacetConnectivity_), ALLOCATABLE :: facetToCell(:)
  TYPE(ElementConnectivity_), ALLOCATABLE :: elemToElem(:)
END TYPE

Two domains can have the following relation between them

p-type

pType has the following properties:

Same geometry
Same number of elements
CellToCell is an identity
Different number of nodes
Different order

h-type

hType has the following properties:

Same geometry
Different number of elements
Different number of nodes
order is same.
Same element topology for a given mesh
CellToCell will be a collection of elements in Domain 2.

r-type

rType has the following properties:

Same geometry
Same number of nodes
Same number of elements
CellToCell is identity matrix
NodeToNode is an identity matrix
Spatial position of nodes would be different

More details

isInitiated, is true if an instance of DomainConnectivity is initiated
isFacetToCell is true if facetToCell data is initiated
isNodeToNode is true if nodeToNode data is initiated
isCellToCell is true if cellToCell data is initiated
nodeToNode, contains the node-to-node mapping from one mesh/domain to another mesh/domain.
- We can create mapping between parts of domain (mesh) or domains. See InitiateNodeToNodeData
- Size of NodeToNode is equal to the largest node number in domain-1 or mesh-1 (depending upon how the data is initiated)
- NodeToNode(i) implies global node number in domain-2/mesh-2, corresponding to global node number i in domain-1/mesh-1.
cellToCell contains the cell-to-cell mapping from one mesh/domain to another mesh/domain.
- We can create mapping between parts of domain (mesh) or domains. See, InitiateCellToCellData
- Size of cellToCell is equal to the largest element number in domain-1 or mesh-1 (depending upon how the data is initiated)
- cellToCell(i) denotes global element number in domain-2/ mesh-2, corresponding to global element number i in domain-1/ mesh-1
cellToCellExtraData: Currently,
- cellToCellExtraData has two rows (a) the first row is dim, (b) the second row is entityNum.
- each column represents the element number, for example
  - iel1 in domain1, then cellToCell(iel1) gives iel2 in domain2
  - cellToCellExtraData(1, iel1) gives dimension of mesh which contains iel2
  - cellToCellExtraData(2, iel1) gives entityNum of mesh which contains iel2
  - In this way, domain2%getMeshPointer(dim, entityNum) can give us the pointer to the mesh which contains the iel2
facetToCell, facet to cell-connectivity data.
elemToElem, currently, we do not know how to implement it (TODO).

FacetConnectivity

TYPE :: FacetConnectivity_
  INTEGER(I4B) :: facetID = 0
    !! global element number of facet element in facet mesh
  INTEGER(I4B) :: GlobalCellData(4, 2) = 0
    !! 1,1 --> Global element number of master cell
    !! 2,1 --> master cell's local facet number connected to facet-elem
    !! 3,1 --> master mesh dimension
    !! 4,1 --> master mesh entity number
    !! 1,2 --> Global element number of slave cell
    !! 2,2 --> slave cell's local facet number connected to facet-elem
    !! 3,2 --> slave mesh dimension
    !! 4,2 --> slave mesh entity number
END TYPE FacetConnectivity_

ElementConnectivity

TYPE :: ElementConnectivity_
  INTEGER(I4B) :: masterGlobalElemNum = 0
  INTEGER(I4B) :: masterLocalFacetID = 0
  INTEGER(I4B) :: slaveGlobalElemNum = 0
  INTEGER(I4B) :: slaveLocalFacetID = 0
END TYPE ElementConnectivity_