scale_mesh_base Module Reference

module FElib / Mesh / Base More...

Data Types
type	meshbase
	Base type to manage a computational mesh. More...
interface	meshbase_get_localmesh
type	meshdiminfo
	Derived type to manage an information of a mesh dimension. More...

Functions/Subroutines
subroutine, public	meshbase_init (this, ndimtype, refelem, nlocalmeshperprc, nsidetile, nprocs)
	Initialize an object to manage a computational mesh.
subroutine, public	meshbase_final (this)
	Finalize an object to manage a computational mesh.
subroutine, public	meshbase_setgeometricinfo (mesh, ndim)

Detailed Description

module FElib / Mesh / Base

Description

Base module to manage meshes for element-based methods

Author

Yuta Kawai, Team SCALE

Function/Subroutine Documentation

meshbase_init()

subroutine, public scale_mesh_base::meshbase_init	(	class(meshbase), intent(inout)	this,
		integer, intent(in)	ndimtype,
		class(elementbase), intent(in), target	refelem,
		integer, intent(in)	nlocalmeshperprc,
		integer, intent(in)	nsidetile,
		integer, intent(in), optional	nprocs )

Initialize an object to manage a computational mesh.

Parameters

[in]	ndimtype	Number of DIMTYPE
[in]	refelem	An object with a reference element
[in]	nlocalmeshperprc	Number of local meshes in each MPI process
[in]	nsidetile	Number of side tiles
[in]	nprocs	MPI processes (if not provided, it will be set to the value from PRC_nprocs)

Definition at line 94 of file scale_mesh_base.F90.

    
    use scale_prc, only: &
      prc_nprocs
    implicit none
 
    class(MeshBase), intent(inout) :: this
    integer, intent(in) :: ndimtype                      !< Number of DIMTYPE
    class(ElementBase), intent(in), target :: refElem    !< An object with a reference element
    integer, intent(in) :: NLocalMeshPerPrc              !< Number of local meshes in each MPI process
    integer, intent(in) :: NsideTile                     !< Number of side tiles
    integer, intent(in), optional :: nprocs              !< MPI processes (if not provided, it will be set to the value from PRC_nprocs)
 
    integer :: n
    !-----------------------------------------------------------------------------
    
    if ( present(nprocs) ) then
      this%PRC_NUM = nprocs
    else
      this%PRC_NUM = prc_nprocs
    end if
 
    this%LOCAL_MESH_NUM        = nlocalmeshperprc
    this%LOCAL_MESH_NUM_global = this%PRC_NUM * this%LOCAL_MESH_NUM
 
    this%refElem => refelem
        
    allocate( this%tileID_globalMap(nsidetile, this%LOCAL_MESH_NUM_global) )
    allocate( this%tileFaceID_globalMap(nsidetile, this%LOCAL_MESH_NUM_global) )
    allocate( this%tilePanelID_globalMap(nsidetile, this%LOCAL_MESH_NUM_global) )
    allocate( this%tileID_global2localMap(this%LOCAL_MESH_NUM_global) )
    allocate( this%PRCRank_globalMap(this%LOCAL_MESH_NUM_global) )
    allocate( this%dimInfo(ndimtype) )
    !$acc enter data create(this%tileID_globalMap, this%tileFaceID_globalMap, this%tilePanelID_globalMap, &
    !$acc   this%tileID_global2localMap, this%PRCRank_globalMap, this%dimInfo )
    
    this%isGenerated = .false.
 
    return

Referenced by scale_mesh_base1d::meshbase1d_init(), scale_mesh_base2d::meshbase2d_init(), scale_mesh_base3d::meshbase3d_init(), and scale_mesh_base::meshbase_get_localmesh::meshbase_get_localmesh().

meshbase_final()

subroutine, public scale_mesh_base::meshbase_final

(

class(meshbase), intent(inout)

this

)

Finalize an object to manage a computational mesh.

Definition at line 139 of file scale_mesh_base.F90.

    implicit none
    class(MeshBase), intent(inout) :: this
    !-----------------------------------------------------------------------------
  
    if ( allocated(this%tileID_globalMap) ) then
      !$acc exit data delete(this%tileID_globalMap, this%tileFaceID_globalMap, this%tilePanelID_globalMap, &
      !$acc   this%tileID_global2localMap, this%PRCRank_globalMap, this%dimInfo)
      deallocate( this%tileID_globalMap )
      deallocate( this%tileFaceID_globalMap )
      deallocate( this%tilePanelID_globalMap )
      deallocate( this%tileID_global2localMap )
      deallocate( this%PRCRank_globalMap )
      deallocate( this%dimInfo )
    end if
    return

Referenced by scale_mesh_base1d::meshbase1d_final(), scale_mesh_base2d::meshbase2d_final(), scale_mesh_base3d::meshbase3d_final(), and scale_mesh_base::meshbase_get_localmesh::meshbase_get_localmesh().

meshbase_setgeometricinfo()

subroutine, public scale_mesh_base::meshbase_setgeometricinfo	(	class(localmeshbase), intent(inout)	mesh,
		integer, intent(in)	ndim )

Definition at line 158 of file scale_mesh_base.F90.

    use scale_element_base, only: elementbase3d
    use scale_localmesh_3d, only: localmesh3d
    implicit none
    
    class(LocalMeshBase), intent(inout) :: mesh
    integer, intent(in) :: ndim
 
    class(ElementBase), pointer :: refElem
    !-----------------------------------------------------------------------------
 
    refelem => mesh%refElem
 
    allocate( mesh%pos_en(refelem%Np,mesh%Ne,ndim) )
    !allocate( mesh%fx(refElem%Nfaces*refElem%Nfp,mesh%Ne) )
    !allocate( mesh%fy(refElem%Nfaces*refElem%Nfp,mesh%Ne) )
    allocate( mesh%normal_fn(refelem%NfpTot,mesh%Ne,ndim) )
    allocate( mesh%sJ(refelem%NfpTot,mesh%Ne) )
    allocate( mesh%J(refelem%Np,mesh%Ne) )
    allocate( mesh%Fscale(refelem%NfpTot,mesh%Ne) )
    allocate( mesh%Escale(refelem%Np,mesh%Ne,ndim,ndim) )
    allocate( mesh%Gsqrt(refelem%Np,mesh%NeA) )
    !$acc enter data create(mesh%pos_en, mesh%normal_fn, mesh%sJ, mesh%J, &
    !$acc   mesh%Fscale, mesh%Escale, mesh%Gsqrt)
 
    return

Referenced by scale_mesh_base1d::meshbase1d_setgeometricinfo(), scale_mesh_base2d::meshbase2d_setgeometricinfo(), scale_mesh_base3d::meshbase3d_setgeometricinfo(), and scale_mesh_base::meshbase_get_localmesh::meshbase_get_localmesh().