Skip to main content

SetMeshMap

This routine sets MeshMap.

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

Let us say there are nn 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, MM, of dimension n×nn \times n. If M(I,J),I,J=1,2,⋯ ,nM(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

In this way:

  • meshMap%IA(i) to meshMap%IA(i+1)-1 denotes the location of iith mesh in obj%meshFacetData
  • meshMap%JA() contains the jjth (column) indices.

This routine makes meshMap.

This routine needs information of facet element data stored inside each mesh of dimension obj%nsd. In this routine, if for any mesh, isFacetDataInitiated is false, then this routine make a call to InitiateFacetElements and intiates the required data.

INTERFACE
  MODULE SUBROUTINE SetMeshmap(obj)
    CLASS(Domain_), INTENT(INOUT) :: obj
  END SUBROUTINE SetMeshmap
END INTERFACE