Structure

The structure of FEDomain CLASS is given below.

TYPE, EXTENDS(AbstractDomain_) :: FEDomain_
  PRIVATE
  CLASS(AbstractMesh_), POINTER :: meshVolume => NULL()
    !! meshVolume list of meshes of volume entities
  CLASS(AbstractMesh_), POINTER :: meshSurface => NULL()
    !! meshSurface list of meshes of surface entities
  CLASS(AbstractMesh_), POINTER :: meshCurve => NULL()
    !! meshCurve list of meshes of curve entities
  CLASS(AbstractMesh_), POINTER :: meshPoint => NULL()
    !! meshPoint list of meshes of point entities
  CLASS(AbstractMesh_), POINTER :: mesh => NULL()
    !! mesh points to meshVolume for nsd = 3
    !! mesh points to meshSurface for nsd = 2
    !! mesh points to meshCurve for nsd = 1
    !! mesh points to meshPoint for nsd = 0
END TYPE FEDomain_

• isInitiated It is true if the domain is initiated.

• engine It is the name of the engine used to create the mesh

• majorVersion It is the major version of the mesh engine

• minorVersion It is the minor version of the mesh engine

• version this is the version of the mesh engine

• nsd this is the spatial dimension of the domain

• maxNptrs, the largest node number in the domain

• minNptrs, the smallest node number in the domain

• isElemNumberSparse True if the element numbers are sparse

• tEntitiesForNodes Total number of entities for the nodes

• tEntitiesForElements total number of entities for the mesh.

• tElements It is an integer vector of length 4. Lower bound is 0 and upper bound is 1.

  • tElements(0) denotes the total number of point elements
  • tElements(1) denotes the total number of curve elements
  • tElements(2) denotes the total number of surface elements
  • tElements(3) denotes the total number of volume elements.

• tEntities It is an integer vector of length 4. It denotes the total number of entities of a given co dimension.

  • tEntities(0) total number of entities (mesh) of point elements
  • tEntities(1) total number of entities (mesh) of curve elements
  • tEntities(2) total number of entities (mesh) of surface elements
  • tEntities(3) total number of entities (mesh) of volume elements.

• nodeCoord nodal coordinates. rows of nodeCoord represent the spatial dimension. columns of nodeCoord represents the spatial node number.

• local_nptrs It is the index mapping which converts global node number to local node number.

• global_nptrs, local to global node number mapping.

• meshList

  • meshList(0) list of meshes of point entities
  • meshList( 1 ) list of meshes of line entities
  • meshList( 2 ) list of meshes of surface entities
  • meshList( 3 ) list of meshes of volume entities

meshFacetData mesh facet data is given below.

TYPE MeshFacetData_
  INTEGER(I4B) :: masterMesh = 0
  INTEGER(I4B) :: slaveMesh = 0
  INTEGER(I4B), ALLOCATABLE :: masterCellNumber(:)
  INTEGER(I4B), ALLOCATABLE :: slaveCellNumber(:)
  INTEGER(I4B), ALLOCATABLE :: masterLocalFacetID(:)
  INTEGER(I4B), ALLOCATABLE :: slaveLocalFacetID(:)
  ! CLASS( Halo_ ), POINTER :: halo => NULL()
END TYPE MeshFacetData_

meshMap contains the connectivity of meshes of dimension nsd (that is, meshes of cell $\Omega$).

Let us say there are $n$ meshes of cell elements. That is,

• In 1D, cell means line elements
• In 2D, cell means surface elements
• In 3D, cell means volume elements

Then, we constructs a sparse matrix, $M$, of dimension $n \times n$. If $M(I,J), I,J=1,2,\cdots,n$ is 1 then mesh-I and mesh-J are connected with each other. Otherwise, they are not connected with each other. This information is stored in an instance of CSRSparsity in meshMap